Natural selection not only affects the actual loci under selection but also leaves “footprints” in patterns of genetic variation in linked genetic regions. This offers exciting opportunities for inferring selection and for understanding the processes shaping levels of genetic variation in natural populations. Here, we develop analytical approximations based on coalescent theory to characterize the genetic footprint of a complex, but potentially common type of natural selection: balancing selection with seasonally fluctuating allele frequencies. As we show analytically and confirm with stochastic simulations, seasonal allele frequency fluctuations can have important (and partly unexpected) consequences for the genetic footprint of balancing selection. Fluctuating balancing selection generally leads to an increase in genetic diversity close to the selected site, the effect of balancing selection, but reduces diversity further away from the selected site, which is a consequence of the allele-frequency fluctuations effectively producing recurrent bottlenecks of allelic backgrounds. This medium- and long-range reduction usually outweighs the short-range increase when averaging diversity levels across the entire chromosome. Strong fluctuating balancing selection even induces a loss of genetic variation in unlinked regions, e.g. on different chromosomes. If many loci in the genome are simultaneously under fluctuating balancing selection this can lead to substantial genome-wide reductions in genetic diversity, even when allele-frequency fluctuations are small and local footprints are difficult to detect. Thus, together with genetic drift, selective sweeps and background selection, fluctuating selection could be a major force shaping levels of genetic diversity in natural populations.

Here, we review the literature on sexual lability in dioecious angiosperm species with well-studied sex chromosomes. We distinguish three types of departures from strict dioecy, concerning either a minority of flowers in some individuals (leakiness) or the entire individual, which can constantly be bisexual or change sex. We found that for only four of the 22 species studied, reports of lability are lacking. The occurrence of lability is only weakly related to sex chromosome characteristics (number of sex-linked genes, age of the non-recombining region). These results contradict the naive idea that lability is an indication of the absence or the recent evolution of sex chromosomes, and thereby contribute to a growing consensus that sex chromosomes do not necessarily fix sex determination once and for all. We discuss some implications of these findings for the evolution of sex chromosomes, and suggest that more species with well-characterized lability should be studied with genomic data and tools. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Here, we review the literature on sexual lability in dioecious angiosperm species with well-studied sex chromosomes. We distinguish three types of departures from strict dioecy, concerning either a minority of flowers in some individuals (leakiness) or the entire individual, which can constantly be bisexual or change sex. We found that for only four of the 22 species studied, reports of lability are lacking. The occurrence of lability is only weakly related to sex chromosome characteristics (number of sex-linked genes, age of the non-recombining region). These results contradict the naive idea that lability is an indication of the absence or the recent evolution of sex chromosomes, and thereby contribute to a growing consensus that sex chromosomes do not necessarily fix sex determination once and for all. We discuss some implications of these findings for the evolution of sex chromosomes, and suggest that more species with well-characterized lability should be studied with genomic data and tools.This article is part of the theme issue ‘Sex determination, sex chromosome evolution, and the role of sexual differentiation in land plants’.

Sex chromosomes are typically comprised of a non-recombining region and a recombining pseudoautosomal region. Accurately quantifying the relative size of these regions is critical for sex-chromosome biology both from a functional and evolutionary perspective. The evolution of the pseudoautosomal boundary (PAB) is well documented in haplorrhines (apes and monkeys) but not in strepsirrhines (lemurs and lorises). Here, we studied the PAB of seven species representing the main strepsirrhine lineages by sequencing a male and a female genome in each species and using sex differences in coverage to identify the PAB. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events moved the PAB and shortened the pseudoautosomal region in haplorrhines. Strepsirrhines are well known to have much lower sexual dimorphism than haplorrhines. We suggest that mutations with antagonistic effects between males and females have driven recombination suppression and PAB evolution in haplorrhines

Recent research shows that introgression between closely-related species is an important source of adaptive alleles for a wide range of taxa. Typically, detection of adaptive introgression from genomic data relies on comparative analyses that require sequence data from both the recipient and the donor species. However, in many cases, the donor is unknown or the data is not currently available. Here, we introduce a genome-scan method-VolcanoFinder-to detect recent events of adaptive introgression using polymorphism data from the recipient species only. VolcanoFinder detects adaptive introgression sweeps from the pattern of excess intermediate-frequency polymorphism they produce in the flanking region of the genome, a pattern which appears as a volcano-shape in pairwise genetic diversity. Using coalescent theory, we derive analytical predictions for these patterns. Based on these results, we develop a composite-likelihood test to detect signatures of adaptive introgression relative to the genomic background. Simulation results show that VolcanoFinder has high statistical power to detect these signatures, even for older sweeps and for soft sweeps initiated by multiple migrant haplotypes. Finally, we implement VolcanoFinder to detect archaic introgression in European and sub-Saharan African human populations, and uncovered interesting candidates in both populations, such as TSHR in Europeans and TCHH-RPTN in Africans. We discuss their biological implications and provide guidelines for identifying and circumventing artifactual signals during empirical applications of VolcanoFinder.

L13 ‘Moved’ seems a bit unspecific in this sentence, I would replace with extended PAR. P3 L12 I’m not aware of any direct evidence that recombination between the primate X and Y stopped instantaneously as ‘at once’ suggests. I would remove this. P3 L14 I am not sure what process is referred to in this sentence. P7 L12 I believe Wright et al showed extension of recombination suppression not additional strata P7 L17

ABSTRACT The human sex chromosomes provide a striking example of how sex chromosomes can become modified in gene content and structure (i.e. become heteromorphic). Why sex chromosomes in humans and other lineages can reach such levels of heteromorphy is not clear. One hypothesis posits that mutations with antagonistic effects between males and females can drive recombination suppression and subsequent differentiation between X and Y chromosomes. Here, we tested this hypothesis by focusing on strepsirrhine primates (lemurs and lorises) that have much lower sexual dimorphism and opportunities for sexually antagonistic mutations than haplorrhines (apes and monkeys). We sequenced seven strepsirrhine species and identified the pseudoautosomal boundary (PAB), the border between recombining and non-recombining regions on the sex chromosomes. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events have occurred in haplorrhines, supporting the view that sexually antagonistic mutations have driven sex chromosomes heteromorphy in primates.

The potential of using the latitudinal distributions of fossil wood genera as a proxy for Jurassic terrestrial climate conditions in western Laurasia was tested using a database of 386 entries. The fossil wood genera Brachyoxylon, Protaxodioxylon and Simplicioxylon are southern taxa, which have only occasionally extended north of 45°N (all latitude references are to palaeolatitudes). Conversely, Protocedroxylon and Xenoxylon are mostly found at high latitudes. The standard deviation of the data on the palaeolatitudes of the northern genera is, however, much larger than that of the southern genera. This is interpreted as resulting from brief periodic southward extensions of the Protocedroxylon and Xenoxylon distributions, while the southern genera remained confined south of 40°N. The interval during which the northern genera replaced the southern genera at mid-latitudes can be hypothesized to be times of global terrestrial climate cooling or increased rainfall. Comparisons between the spatiotemporal distributions of these ‘northern’ taxa and independent climate-sensitive indicators (clay mineral assemblages, oxygen isotopes, and neritic carbonate distribution) suggest that these genera were adapted to more humid conditions characterizing the mid-latitudes of the Northern Hemisphere. Accordingly, southward migrations of this climatic zone occurred during the Early and Late Pliensbachian, Early Toarcian, and Early Middle Oxfordian time intervals. These reconstructions are limited by the low resolution and restricted spatial coverage of the wood record. Therefore, further sampling at the ammonite zone level is needed. Nevertheless, the present study suggests that wood taxa distribution represents an underexploited but potentially powerful record of terrestrial environmental changes during the Jurassic Period.

We propose a probabilistic framework to infer autosomal and sex-linked genes from RNA-seq data of a cross for any sex chromosome type (XY, ZW, and UV). Sex chromosomes (especially the non-recombining and repeat-dense Y, W, U, and V) are notoriously difficult to sequence. Strategies have been developed to obtain partially assembled sex chromosome sequences. Most of them remain difficult to apply to numerous non-model organisms, either because they require a reference genome, or because they are designed for evolutionarily old systems. Sequencing a cross (parents and progeny) by RNA-seq to study the segregation of alleles and infer sex-linked genes is a cost-efficient strategy, which also provides expression level estimates. However, the lack of a proper statistical frame work has limited a broader application of this approach. Tests on empirical Silene data show that our method identifies 20-35% more sex-linked genes than existing pipelines, while making reliable inferences for downstream analyses. Approximately 12 individuals are needed for optimal results based on simulations. For species with an unknown sex-determination system, the method can assess the presence and type (XY vs. ZW) of sex chromosomes through a model comparison strategy. The method is particularly well optimized for sex chromosomes of young or intermediate age, which are expected in thousands of yet unstudied lineages. Any organisms, including non-model ones for which nothing is known a priori, that can be bred in the lab, are suitable for our method. SEX-DETector and its implementation in a Galaxy workflow are made freely available.

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%.

Comparing species richness in sister clades that differ in a character state is one of the ways to study factors influencing diversification. While most of its applications have focussed on traits that increase diversification, some have been used to study the association of a trait with lower species richness, e.g., the occurrence of dioecy in flowering plants. We show here, using simulations and an analytical model, that the null expectation of equal species richness that is generally used in sister clade comparisons is wrong in the case of a derived trait occurring independently from speciation: one should expect fewer species in the clade with the derived character state when there is no difference in diversification rates. This is due to the waiting time for the derived state to appear, which causes it to occur more often on longer branches. This has the important implication that the probability for a clade to possess the derived state depends on the tree geometry, and thus on species richness: species-poorer clades are more likely to possess the derived state. We develop a statistical test for sister clade comparisons to study the effect of a derived character state. Applying it to a data set of dioecious clades, we find that we cannot confirm earlier work that concluded that dioecy decreases diversification; on the contrary, it seems to be associated to higher species richness than expected. [angiosperms; dioecy; diversification; sister clades; species richness.]

Large tandem repeat sequences have been poorly investigated as severe technical limitations and their frequent absence from the genome reference hinder their analysis. Extensive allelotyping of this class of variation has not been possible until now and their mutational dynamics are still poorly known. In order to estimate the mutation rate of a macrosatellite, we analysed in detail the RNU2 locus, which displays at least 50 different alleles containing 5-82 copies of a 6.1 kb repeat unit. Mining data from the 1000 Genomes Project allowed us to precisely estimate copy numbers of the RNU2 repeat unit using read depth of coverage. This further revealed significantly different mean values in various recent modern human populations, favoring a scenario of fast evolution of this locus. Its proximity to a disease gene with numerous founder mutations, BRCA1, within the same linkage disequilibrium block, offered the unique opportunity to trace RNU2 arrays over a large timescale. Analysis of the transmission of RNU2 arrays associated with one 'private' mutation in an extended kindred and four founder mutations in multiple kindreds gave an estimation by maximum likelihood of 5 × 10(-3) mutations per generation, which is close to that of microsatellites.

In angiosperms, dioecious clades tend to have fewer species than their nondioecious sister clades. This departure from the expected equal species richness in the standard sister clade test has been interpreted as implying that dioecious clades diversify less and has initiated a series of studies suggesting that dioecy might be an 'evolutionary dead end'. However, two of us recently showed that the 'equal species richness' null hypothesis is not valid in the case of derived char acters, such as dioecy, and proposed a new test for sister clade comparisons; preliminary results, using a data set available in the litterature, indicated that dioecious clades migth diversify more than expected. However, it is crucial for this new test to distinguish between ancestral and derived cases of dioecy, a criterion that was not taken into account in the available data set. Here, we present a new data set that was obtained by searching the phylogenetic literature on more than 600 completely dioecious angiosperm genera and identifying 115 sister clade pairs for which dioecy is likely to be derived (including > 50% of the dioecious species). Applying the new sister clade test to this new dataset, we confirm the preliminary result that dioecy is associated with an increased diversification rate, a result that does not support the idea that dioecy is an evolutionary dead end in angiosperms. The traits usually associated with dioecy, that is, an arborescent growth form, abiotic pollination, fleshy fruits or a tropical distribution, do not influence the diversification rate. Rather than a low diversification rate, the observed species richness patterns of dioecious clades seem to be better explained by a low transition rate to dioecy and frequent losses.

BACKGROUND AND AIMS:About 6 % of an estimated total of 240 000 species of angiosperms are dioecious. The main precursors of this sexual system are thought to be monoecy and gynodioecy. A previous angiosperm-wide study revealed that many dioecious species have evolved through the monoecy pathway; some case studies and a large body of theoretical research also provide evidence in support of the gynodioecy pathway. If plants have evolved through the gynodioecy pathway, gynodioecious and dioecious species should co-occur in the same genera. However, to date, no large-scale analysis has been conducted to determine the prevalence of the gynodioecy pathway in angiosperms. In this study, this gap in knowledge was addressed by performing an angiosperm-wide survey in order to test for co-occurrence as evidence of the gynodioecy pathway.METHODS:Data from different sources were compiled to obtain (to our knowledge) the largest dataset on gynodioecy available, with 275 genera that include at least one gynodioecious species. This dataset was combined with a dioecy dataset from the literature, and a study was made of how often dioecious and gynodioecious species could be found in the same genera using a contingency table framework.KEY RESULTS:It was found that, overall, angiosperm genera with both gynodioecious and dioecious species occur more frequently than expected, in agreement with the gynodioecy pathway. Importantly, this trend holds when studying different classes separately (or sub-classes, orders and families), suggesting that the gynodioecy pathway is not restricted to a few taxa but may instead be widespread in angiosperms.CONCLUSIONS:This work complements that previously carried out on the monoecy pathway and suggests that gynodioecy is also a common pathway in angiosperms. The results also identify angiosperm families where some (or all) dioecious species may have evolved from gynodioecious precursors. These families could be the targets of future small-scale studies on transitions to dioecy taking phylogeny explicitly into account.

Genetic diversity plays a major role in ensuring the adaptive capacity and long-term viability of populations. Recent studies have shown that evaluations of the ecological benefits of ecosystem restoration should incorporate the impact of restoration on population genetic variability to accurately guide restoration priorities. This study aimed to analyse the population genetic diversity of an aquatic plant, Berula erecta, in wetlands restored by sediment dredging using natural local wetlands as a reference to assess the impact of ecological restoration on genetic diversity. The analyses were performed using microsatellite markers. Restoration led to a considerable increase in the open water surface area at the site, and the restored areas were efficiently colonised by B. erecta. The populations exhibited a high level of genetic diversity, even though clonal propagation is the major reproductive mode of this species. Two years after restoration, population genetic diversity (AR allelic richness, Fis and PG) was very similar between the restored and natural wetlands. For one restored wetland, the individuals in the restored area of the wetland showed greater allelic richness than those in the unrestored areas of the same wetland. In addition, most of the multilocus genotypes observed in the restored wetlands only occurred as single individuals, suggesting that clonal propagation is not the only mode of recolonisation after restoration. The source of the new genotypes that appeared after restoration may have been a local relict population, seeds remaining in the seed bank after restoration (cryptic genetic diversity), or propagules (seed or vegetative fragments) that migrated into the wetland. These results suggest that the benefits of restoration may depend strongly on the genetic diversity of established relict and dormant populations. Population diversity, in turn, may partly depend on the ecological characteristics of the site (i.e., connectivity of the water body to the surrounding wetlands, successional stage and site history). (C) 2013 Elsevier B.V. All rights reserved.

Exon Primed Intron Crossing (EPIC) markers provide molecular tools that are susceptible to be variable within species while remaining amplifiable by PCR using potentially universal primers. In this study we tested the possibility of obtaining PCR products from 50 EPIC markers on 23 species belonging to seven different phyla (Porifera, Cnidaria, Arthropoda, Nematoda, Mollusca, Annelida, Echinodermata) using 70 new primer pairs. A previous study had identified and tested those loci in a dozen species, including another phylum, Urochordata (Chenuil et al., 2010). Results were contrasted among species. The best results were achieved with the oyster (Mollusca) where 28 loci provided amplicons susceptible to contain an intron according to their size. This was however not the case with the other mollusk Crepidula fornicata, which seems to have undergone a reduction in intron number or intron size. In the Porifera, 13 loci appeared susceptible to contain an intron, a surprisingly high number for this phylum considering its phylogenetic distance with genomic data used to design the primers. For two cnidarian species, numerous loci (24) were obtained. Ecdysozoan phyla (arthropods and nematodes) proved less successful than others as expected considering reports of their rapid rate of genome evolution and the worst results were obtained for several arthropods. Some general patterns among phyla arose, and we discuss how the results of this EPIC survey may give new insights into genome evolution of the study species. This work confirms that this set of EPIC loci provides an easy-to-use toolbox to identify genetic markers potentially useful for population genetics, phylogeography or phylogenetic studies for a large panel of metazoan species. We then argue that obtaining diploid sequence genotypes for these loci became simple and affordable owing to Next-Generation Sequencing development. Species surveyed in this study belong to several genera (Acanthaster, Alvinocaris, Aplysina, Aurelia, Crepidula, Eunicella, Hediste, Hemimysis, Litoditis, Lophelia, Mesopodopsis, Mya, Ophiocten, Ophioderma, Ostrea, Pelagia, Platynereis, Rhizostoma, Rimicaris), two of them, belonging to the family Vesicomydae and Eunicidae, could not be determined at the genus level.

GC-biased gene conversion (gBGC) is a process that tends to increase the GC content of recombining DNA over evolutionary time and is thought to explain the evolution of GC content in mammals and yeasts. Evidence for gBGC outside these two groups is growing but is still limited. Here, we analyzed 36 completely sequenced genomes representing four of the five major groups in eukaryotes (Unikonts, Excavates, Chromalveolates and Plantae). gBGC was investigated by directly comparing GC content and recombination rates in species where recombination data are available, that is, half of them. To study all species of our dataset, we used chromosome size as a proxy for recombination rate and compared it with GC content. Among the 17 species showing a significant relationship between GC content and chromosome size, 15 are consistent with the predictions of the gBGC model. Importantly, the species showing a pattern consistent with gBGC are found in all the four major groups of eukaryotes studied, which suggests that gBGC may be widespread in eukaryotes.

Silene latifolia is a dioecious plant with heteromorphic sex chromosomes that have originated only ∼10 million years ago and is a promising model organism to study sex chromosome evolution in plants. Previous work suggests that S. latifolia XY chromosomes have gradually stopped recombining and the Y chromosome is undergoing degeneration as in animal sex chromosomes. However, this work has been limited by the paucity of sex-linked genes available. Here, we used 35 Gb of RNA-seq data from multiple males (XY) and females (XX) of an S. latifolia inbred line to detect sex-linked SNPs and identified more than 1,700 sex-linked contigs (with X-linked and Y-linked alleles). Analyses using known sex-linked and autosomal genes, together with simulations indicate that these newly identified sex-linked contigs are reliable. Using read numbers, we then estimated expression levels of X-linked and Y-linked alleles in males and found an overall trend of reduced expression of Y-linked alleles, consistent with a widespread ongoing degeneration of the S. latifolia Y chromosome. By comparing expression intensities of X-linked alleles in males and females, we found that X-linked allele expression increases as Y-linked allele expression decreases in males, which makes expression of sex-linked contigs similar in both sexes. This phenomenon is known as dosage compensation and has so far only been observed in evolutionary old animal sex chromosome systems. Our results suggest that dosage compensation has evolved in plants and that it can quickly evolve de novo after the origin of sex chromosomes.

BACKGROUND:Despite the impressive growth of sequence databases, the limited availability of nuclear markers that are sufficiently polymorphic for population genetics and phylogeography and applicable across various phyla restricts many potential studies, particularly in non-model organisms. Numerous introns have invariant positions among kingdoms, providing a potential source for such markers. Unfortunately, most of the few known EPIC (Exon Primed Intron Crossing) loci are restricted to vertebrates or belong to multigenic families.RESULTS:In order to develop markers with broad applicability, we designed a bioinformatic approach aimed at avoiding multigenic families while identifying intron positions conserved across metazoan phyla. We developed a program facilitating the identification of EPIC loci which allowed slight variation in intron position. From the Homolens databases we selected 29 gene families which contained 52 promising introns for which we designed 93 primer pairs. PCR tests were performed on several ascidians, echinoderms, bivalves and cnidarians. On average, 24 different introns per genus were amplified in bilaterians. Remarkably, five of the introns successfully amplified in all of the metazoan genera tested (a dozen genera, including cnidarians). The influence of several factors on amplification success was investigated. Success rate was not related to the phylogenetic relatedness of a taxon to the groups that most influenced primer design, showing that these EPIC markers are extremely conserved in animals.CONCLUSIONS:Our new method now makes it possible to (i) rapidly isolate a set of EPIC markers for any phylum, even outside the animal kingdom, and thus, (ii) compare genetic diversity at potentially homologous polymorphic loci between divergent taxa.

The identification of loci influenced by positive selection is a major goal of evolutionary genetics. A popular approach is to perform scans of alignments on a genome-wide scale in order to find regions evolving at accelerated rates on a particular branch of a phylogenetic tree. However, positive selection is not the only process that can lead to accelerated evolution. Notably, GC-biased gene conversion (gBGC) is a recombination-associated process that results in the biased fixation of G and C nucleotides. This process can potentially generate bursts of nucleotide substitutions within hotspots of meiotic recombination. Here, we analyse the results of a scan for positive selection on genes on branches across the primate phylogeny. We show that genes identified as targets of positive selection have a significant tendency to exhibit the genomic signature of gBGC. Using a maximum-likelihood framework, we estimate that more than 20 per cent of cases of significantly elevated non-synonymous to synonymous substitution rates ratio ( d N / d S ), particularly in shorter branches, could be due to gBGC. We demonstrate that in some cases, gBGC can lead to very high d N / d S (more than 2). Our results indicate that gBGC significantly affects the evolution of coding sequences in primates, often leading to patterns of evolution that can be mistaken for positive selection.