Sylvain Mousset

Power of neutrality tests to detect bottlenecks and hitchhiking.

The power of several neutrality tests to reject a simple bottleneck model is examined in a coalescent framework. Several tests are considered including some relying on the frequency spectrum of mutations and some reflecting the linkage disequilibrium structure of the data. We evaluate the effect of the age and of the strength of the bottleneck, and their interaction. We contrast two qualitatively different bottleneck effects depending on their strength. In genealogical terms, during severe bottlenecks, all lineages coalesce leading to a star-like gene genealogy of the sample. Some time after the bottleneck, once new mutations have arisen, they tend to show an excess of rare variants and a slight excess of haplotypes. On the contrary, more moderate bottlenecks allow several lineages to survive the demographic crash, leading to a balanced genealogy with long internal branches. Soon after the event, data tend to show an excess of intermediate frequency variants and a deficit of haplotypes. We show that for moderate sequencing efforts, severe bottlenecks can be detected only after an intermediate time period has allowed for mutations to occur, preferably by frequency spectrum statistics. Moderate bottlenecks can be more easily detected for more recent events, especially using haplotype statistics. Finally, for a single locus, the bottleneck results closely approximate those of a simple hitchhiking model. The main difference concerns the frequency distribution of mutations and haplotypes after moderate perturbations. Hitchhiking increases the number of rare ancestral mutations and leads to a more predominant major haplotype class. Thus, despite a number of common features between the two processes, hitchhiking cannot be strictly modeled by bottlenecks.

Rapid de novo evolution of X chromosome dosage compensation in Silene latifolia, a plant with young sex chromosomes.

Evidence for dosage compensation in Silene latifolia, a plant with 10-million-year-old sex chromosomes, reveals that dosage compensation can evolve rapidly in young XY systems and is not an animal-specific phenomenon.

Detecting positive selection within genomes: the problem of biased gene conversion

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 (dN/dS), particularly in shorter branches, could be due to gBGC. We demonstrate that in some cases, gBGC can lead to very high dN/dS (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.

Evidence for Widespread GC-biased Gene Conversion in Eukaryotes

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

Analysis of sequence variability in the macronuclear DNA of Paramecium tetraurelia: A somatic view of the germline

Ciliates are the only unicellular eukaryotes known to separate germinal and somatic functions. Diploid but silent micronuclei transmit the genetic information to the next sexual generation. Polyploid macronuclei express the genetic information from a streamlined version of the genome but are replaced at each sexual generation. The macronuclear genome of Paramecium tetraurelia was recently sequenced by a shotgun approach, providing access to the gene repertoire. The 72-Mb assembly represents a consensus sequence for the somatic DNA, which is produced after sexual events by reproducible rearrangements of the zygotic genome involving elimination of repeated sequences, precise excision of unique-copy internal eliminated sequences (IES), and amplification of the cellular genes to high copy number. We report use of the shotgun sequencing data (>10(6) reads representing 13 x coverage of a completely homozygous clone) to evaluate variability in the somatic DNA produced by these developmental genome rearrangements. Although DNA amplification appears uniform, both of the DNA elimination processes produce sequence heterogeneity. The variability that arises from IES excision allowed identification of hundreds of putative new IESs, compared to 42 that were previously known, and revealed cases of erroneous excision of segments of coding sequences. We demonstrate that IESs in coding regions are under selective pressure to introduce premature termination of translation in case of excision failure.

An efficient method to find potentially universal population genetic markers, applied to metazoans

BackgroundDespite 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.ResultsIn 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.ConclusionsOur 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.

An angiosperm-wide analysis of the gynodioecy–dioecy pathway

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.

Standard Sister Clade Comparison Fails when Testing Derived Character States

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

Selective sweep near the In(2L)t inversion breakpoint in an African population of Drosophila melanogaster.

Chromosomal inversions largely inhibit recombination and may be associated with selective forces, such as hitch-hiking effects: the effect of positive selection on linked loci. A West African population of Drosophila melanogaster showed a high frequency (0.61) of the In(2L)t inversion. Departure from neutrality statistically associated with the inversion polymorphism was previously recorded at Su(H), a locus distant from the proximal breakpoint of the inversion. These results were consistent with hitch-hiking effects with recombination. The present sequence polymorphism survey involves a 1 kb fragment of the Vha68-1 locus located closer to the proximal breakpoint of the inversion. It shows a significant deficit of polymorphism with respect to divergence when compared with other loci studied in the same population, thus suggesting selective effects. Only 11 polymorphic sites are present in a sample of 20 chromosomes and these sites present a significant excess of rare-frequency variants. The major haplotype shows an unexpectedly high frequency. Our estimate of the background selection effect is not sufficient to account for the observed reduction of polymorphism. Intraspecific variation is structured between inverted and standard chromosomes; there are no shared polymorphisms but also no fixed differences between them. This pattern, together with that found on other loci previously studied near this inversion breakpoint, suggests hitch-hiking effects enhanced by the inversion.

SEX-DETector: A Probabilistic Approach to Study Sex Chromosomes in Non-Model Organisms

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