Are invasive marsh frogs (Pelophylax ridibundus) replacing the native P. lessonae /P. esculentus hybridogenetic complex in Western Europe? Genetic evidence from a field study.

The water-frog L-E system, widespread in Western Europe, comprises the pool frog Pelophylax lessonae and the hybridogenetic edible frog P. esculentus, which originated from hybridization between pool frogs and marsh frogs (P. ridibundus). In P. esculentus, the lessonae (L) genome is eliminated during meiosis and has to be gained anew each generation from a P. lessonae partner, while the ridibundus (R') genome is transmitted clonally. It therefore accumulates deleterious mutations, so that R'R' offspring from P. esculentus x P. esculentus crosses are normally unviable. This system is now threatened by invasive P. ridibundus (RR) imported from Eastern Europe and the Balkans. We investigated the genetic interactions between invasive marsh frogs and native water frogs in a Swiss wetland area, and used genetic data collected in the field to validate several components of a recently postulated mechanism of species replacement. We identified neo-ridibundus individuals derived from crosses between invasive ridibundus and native esculentus, as well as newly arisen hybridogenetic esculentus lineages stemming from crosses between invasive ridibundus (RR) and native lessonae (LL). As their ridibundus genomes are likely to carry less deleterious mutations, such lineages are expected to produce viable ridibundus offspring, contributing to species replacement. However, such crosses with invasive ridibundus only occurred at a limited scale; moreover, RR x LL crosses did not induce any introgression from the ridibundus to the lessonae genome. We did not find any ridibundus stemming from crosses between ancient esculentus lineages. Despite several decades of presence on the site, introduced ridibundus individuals only represent 15% of sampled frogs, and their spatial distribution seems shaped by specific ecological requirements rather than history of colonization. We therefore expect the three taxa to coexist stably in this area.

Joint evolution of dispersal and inbreeding load.

Inbreeding avoidance is often invoked to explain observed patterns of dispersal, and theoretical models indeed point to a possibly important role. However, while inbreeding load is usually assumed constant in these models, it is actually bound to vary dynamically under the combined influences of mutation, drift, and selection and thus to evolve jointly with dispersal. Here we report the results of individual-based stochastic simulations allowing such a joint evolution. We show that strongly deleterious mutations should play no significant role, owing to the low genomic mutation rate for such mutations. Mildly deleterious mutations, by contrast, may create enough heterosis to affect the evolution of dispersal as an inbreeding-avoidance mechanism, but only provided that they are also strongly recessive. If slightly recessive, they will spread among demes and accumulate at the metapopulation level, thus contributing to mutational load, but not to heterosis. The resulting loss of viability may then combine with demographic stochasticity to promote population fluctuations, which foster indirect incentives for dispersal. Our simulations suggest that, under biologically realistic parameter values, deleterious mutations have a limited impact on the evolution of dispersal, which on average exceeds by only one-third the values expected from kin-competition avoidance.

Distribution of the transposable elements bilbo and gypsy in original and colonizing populations of Drosophila subobscura

Background: Transposable elements (TEs) constitute a substantial amount of all eukaryotic genomes. They induce an important proportion of deleterious mutations by insertion into genes or gene regulatory regions. However, their mutational capabilities are not always adverse but can contribute to the genetic diversity and evolution of organisms. Knowledge of their distribution and activity in the genomes of populations under different environmental and demographic regimes, is important to understand their role in species evolution. In this work we study the chromosomaldistribution of two TEs, gypsy and bilbo, in original and colonizing populations of Drosophilasubobscura to reveal the putative effect of colonization on their insertion profile.Results: Chromosomal frequency distribution of two TEs in one original and three colonizingpopulations of D. subobscura, is different. Whereas the original population shows a low insertionfrequency in most TE sites, colonizing populations have a mixture of high (frequency ¿ 10%) andlow insertion sites for both TEs. Most highly occupied sites are coincident among colonizingpopulations and some of them are correlated to chromosomal arrangements. Comparisons of TEcopy number between the X chromosome and autosomes show that gypsy occupancy seems to becontrolled by negative selection, but bilbo one does not. Conclusion: These results are in accordance that TEs in Drosophila subobscura colonizing populations are submitted to a founder effect followed by genetic drift as a consequence of colonization. This would explain the high insertion frequencies of bilbo and gypsy in coincident sites of colonizing populations. High occupancy sites would represent insertion events prior to colonization. Sites of low frequency would be insertions that occurred after colonization and/orcopies from the original population whose frequency is decreasing in colonizing populations. Thiswork is a pioneer attempt to explain the chromosomal distribution of TEs in a colonizing specieswith high inversion polymorphism to reveal the putative effect of arrangements in TE insertionprofiles. In general no associations between arrangements and TE have been found, except in a fewcases where the association is very strong. Alternatively, founder drift effects, seem to play aleading role in TE genome distribution in colonizing populations.

Sex-chromosome turnovers: the hot-potato model.

Abstract Sex-determining systems often undergo high rates of turnover but for reasons that remain largely obscure. Two recent evolutionary models assign key roles, respectively, to sex-antagonistic (SA) mutations occurring on autosomes and to deleterious mutations accumulating on sex chromosomes. These two models capture essential but distinct key features of sex-chromosome evolution; accordingly, they make different predictions and present distinct limitations. Here we show that a combination of features from the two models has the potential to generate endless cycles of sex-chromosome transitions: SA alleles accruing on a chromosome after it has been co-opted for sex induce an arrest of recombination; the ensuing accumulation of deleterious mutations will soon make a new transition ineluctable. The dynamics generated by these interactions share several important features with empirical data, namely, (i) that patterns of heterogamety tend to be conserved during transitions and (ii) that autosomes are not recruited randomly, with some chromosome pairs more likely than others to be co-opted for sex.

A Y-like social chromosome causes alternative colony organization in fire ants.

Intraspecific variability in social organization is common, yet the underlying causes are rarely known. In the fire ant Solenopsis invicta, the existence of two divergent forms of social organization is under the control of a single Mendelian genomic element marked by two variants of an odorant-binding protein gene. Here we characterize the genomic region responsible for this important social polymorphism, and show that it is part of a pair of heteromorphic chromosomes that have many of the key properties of sex chromosomes. The two variants, hereafter referred to as the social B and social b (SB and Sb) chromosomes, are characterized by a large region of approximately 13 megabases (55% of the chromosome) in which recombination is completely suppressed between SB and Sb. Recombination seems to occur normally between the SB chromosomes but not between Sb chromosomes because Sb/Sb individuals are non-viable. Genomic comparisons revealed limited differentiation between SB and Sb, and the vast majority of the 616 genes identified in the non-recombining region are present in the two variants. The lack of recombination over more than half of the two heteromorphic social chromosomes can be explained by at least one large inversion of around 9 megabases, and this absence of recombination has led to the accumulation of deleterious mutations, including repetitive elements in the non-recombining region of Sb compared with the homologous region of SB. Importantly, most of the genes with demonstrated expression differences between individuals of the two social forms reside in the non-recombining region. These findings highlight how genomic rearrangements can maintain divergent adaptive social phenotypes involving many genes acting together by locally limiting recombination.

Sex reversal: a fountain of youth for sex chromosomes?

Nonrecombining Y chromosomes are expected to degenerate through the progressive accumulation of deleterious mutations. In lower vertebrates, however, most species display homomorphic sex chromosomes. To address this, paradox I propose a role for sex reversal, which occasionally occurs in ectotherms due to the general dependence of physiological processes on temperature. Because sex-specific recombination patterns depend on phenotypic, rather than genotypic sex, homomorphic X and Y chromosomes are expected to recombine in sex-reversed females. These rare events should generate bursts of new Y haplotypes, which will be quickly sorted out by natural or sexual selection. By counteracting Muller's ratchet, this regular purge should prevent the evolutionary decay of Y chromosomes. I review empirical data supporting this suggestion, and propose further investigations for testing it.

Evolutionary analysis of genetic variants involved in rare diseases

Next-generation sequencing techniques such as exome sequencing can successfully detect all genetic variants in a human exome and it has been useful together with the implementation of variant filters to identify causing-disease mutations. Two filters aremainly used for the mutations identification: low allele frequency and the computational annotation of the genetic variant. Bioinformatic tools to predict the effect of a givenvariant may have errors due to the existing bias in databases and sometimes show a limited coincidence among them. Advances in functional and comparative genomics are needed in order to properly annotate these variants.The goal of this study is to: first, functionally annotate Common Variable Immunodeficiency disease (CVID) variants with the available bioinformatic methods in order to assess the reliability of these strategies. Sencondly, as the development of new methods to reduce the number of candidate genetic variants is an active and necessary field of research, we are exploring the utility of gene function information at organism level as a filter for rare disease genes identification. Recently, it has been proposed that only 10-15% of human genes are essential and therefore we would expect that severe rare diseases are mostly caused by mutations on them. Our goal is to determine whether or not these rare and severe diseases are caused by deleterious mutations in these essential genes. If this hypothesis were true, taking into account essential genes as a filter would be an interesting parameter to identify causingdisease mutations.

A hyperactive quantitative trait locus allele of Arabidopsis BRX contributes to natural variation in root growth vigor.

Quantitative trait loci analysis of natural Arabidopsis thaliana accessions is increasingly exploited for gene isolation. However, to date this has mostly revealed deleterious mutations. Among them, a loss-of-function allele identified the root growth regulator BREVIS RADIX (BRX). Here we present evidence that BRX and the paralogous BRX-LIKE (BRXL) genes are under selective constraint in monocotyledons as well as dicotyledons. Unexpectedly, however, whereas none of the Arabidopsis orthologs except AtBRXL1 could complement brx null mutants when expressed constitutively, nearly all monocotyledon BRXLs tested could. Thus, BRXL proteins seem to be more diversified in dicotyledons than in monocotyledons. This functional diversification was correlated with accelerated rates of sequence divergence in the N-terminal regions. Population genetic analyses of 30 haplotypes are suggestive of an adaptive role of AtBRX and AtBRXL1. In two accessions, Lc-0 and Lov-5, seven amino acids are deleted in the variable region between the highly conserved C-terminal, so-called BRX domains. Genotyping of 42 additional accessions also found this deletion in Kz-1, Pu2-7, and Ws-0. In segregating recombinant inbred lines, the Lc-0 allele (AtBRX(Lc-0)) conferred significantly enhanced root growth. Moreover, when constitutively expressed in the same regulatory context, AtBRX(Lc-0) complemented brx mutants more efficiently than an allele without deletion. The same was observed for AtBRXL1, which compared with AtBRX carries a 13 amino acid deletion that encompasses the deletion found in AtBRX(Lc-0). Thus, the AtBRX(Lc-0) allele seems to contribute to natural variation in root growth vigor and provides a rare example of an experimentally confirmed, hyperactive allelic variant.

Genes, environment and sex : perspectives from ectothermic vertebrates

RESUMELes modèles classiques sur l'évolution des chromosomes sexuels supposent que des gènes sexe- antagonistes s'accumulent sur les chromosomes sexuels, entraînant ainsi l'apparition d'une région non- recombinante, qui se répand progressivement en favorisant l'accumulation de mutations délétères. En accord avec cette théorie, les chromosomes sexuels que l'on observe aujourd'hui chez les mammifères et les oiseaux sont considérablement différenciés. En revanche, chez la plupart des vertébrés ectothermes, les chromosomes sexuels sont indifférenciés et il existe une impressionnante diversité de mécanismes de détermination du sexe. Au cours de cette thèse, j'ai étudié l'évolution des chromosomes sexuels chez les vertébrés ectothermes, en outre pour mieux comprendre ce contraste avec les vertébrés endothermes. L'hypothèse « high-turnover » postule que les chromosomes sexuels sont remplacés régulièrement à partir d'autosomes afin d'éviter leur dégénérescence. L'hypothèse « fountain-of-youth » propose que la recombinaison entre le chromosome X et le chromosome Y au sein de femelles XY empêche la dégénérescence. Les résultats de ma thèse, basés sur des études théoriques et empiriques, suggèrent que les deux processus peuvent être entraînés par l'environnement et ainsi jouent un rôle important dans l'évolution des chromosomes sexuels chez les vertébrés ectothermes.SUMMARYClassical models of sex-chromosome evolution assume that sexually antagonistic genes accumulate on sex chromosomes leading to a non-recombining region, which progressively expands and favors the accumulation of deleterious mutations. Concordant with this theory, sex chromosomes in extant mammals and birds are considerably differentiated. In most ectothermic vertebrates, such as frogs, however, sex chromosomes are undifferentiated and a striking diversity of sex determination systems is observed. This thesis was aimed to investigate this apparent contrast of sex chromosome evolution between endothermic and ectothermic vertebrates. The "high-turnover" hypothesis holds that sex chromosomes arose regularly from autosomes preventing decay. The "fountain-of-youth" hypothesis posits that sex chromosomes undergo episodic X-Y recombination in sex-reversed XY females, thereby purging ("rejuvenating") the Y chromosome. We suggest that both processes likely played an important role in sex chromosome evolution of ectothermic vertebrates. The literature largely views sex determination as a dichotomous process: individual sex is assumed to be determined either by genetic (genotypic sex determination, GSD) or by environmental factors (environmental sex determination, ESD), most often temperature (temperature sex determination, TSD). We endorsed an alternative view, which sees GSD and TSD as the ends of a continuum. The conservatism of molecular processes among different systems of sex determination strongly supports the continuum view. We proposed to define sex as a threshold trait underlain by a liability factor, and reaction norms allowing modeling interactions between genotypic and temperature effects. We showed that temperature changes (due to e.g., climatic changes or range expansions) are expected to provoke turnovers in sex-determination mechanisms maintaining homomorphic sex chromosomes. The balanced lethal system of crested newts might be the result of such a sex determination turnover, originating from two variants of ancient Y-chromosomes. Observations from a group of tree frogs, on the other hand, supported the 'fountain of youth' hypothesis. We then showed that low rates of sex- reversals in species with GSD might actually be adaptive considering joint effects of deleterious mutation purging and sexually antagonistic selection. Ongoing climatic changes are expected to threaten species with TSD by biasing population sex ratios. In contrast, species with GSD are implicitly assumed immune against such changes, because genetic systems are thought to necessarily produce even sex ratios. We showed that this assumption may be wrong and that sex-ratio biases by climatic changes may represent a previously unrecognized extinction threat for some GSD species.

Demogenetics in metapopulations

SUMMARY : Human-induced habitat fragmentation constitutes a major threat to biodiversity. Small and isolated populations suffer from increased stochasticity and from limited rescue effects. These two factors may be sufficient to cause local extinctions but fragmentation induces some genetic consequences that can also contribute significantly to extinction risks. Increased genetic drift reduces the effectiveness of selection against deleterious mutations, leading to their progressive accumulation. Drift also decreases both the standing genetic variation and the rate of fixation of beneficial mutations, limiting the evolutionary potential of isolated populations. Demography and genetics further interact and feed back on each other, progressively driving fragmented populations into "extinction vortices". The aim of the thesis was to better understand the processes occurring in fragmented populations. For this, I combined simulation studies and empirical data from three species that live in structured habitats. Chapter 1 and 2 investigate the demography of two shrew species in fragmented habitats. I showed that connectivity and habitat quality strongly affect the demography of the greater white-tooted shrew, although demographic stochasticity was extremely high. I also demonstrated that habitat fragmentation is one of the leading factors allowing the local coexistence of two competing shrew species. Chapter 3 and 4 focus on measuring connectivity in fragmented populations based on genetic data. In particular, I showed that genetic data can be used to detect the landscape elements impeding dispersal. In Chapter 5 that deals with the accumulation of deleterious mutations in fragmented populations, I demonstrated that mutation accumulation, as well a time to extinction, can be predicted from simple demographic and genetic measures. In the last two chapters, I monitored individual reproductive success in an isolated tree frogs population. These data allowed quantifying the effective population size, a measure closely linked to population evolutionary potential. To conclude, this thesis brings some new insights into the processes occurring in fragmented populations, and I hope it will contribute to the improvement of the management and conservation of fragmented populations.

Evidence of recombination in putative ancient asexuals.

Ancient asexuals have been considered to be a contradiction of the basic tenets of evolutionary theory. Barred from rearranging genetic variation by recombination, their reduced number of gene arrangements is thought to hamper their response to changing environments. For the same reason, it should be difficult for them to avoid the build-up of deleterious mutations. Several groups of taxonomically diverse organisms are thought to be ancient asexuals, although clear evidence for or against the existence of recombination events is scarce. Several methods have recently been developed for predicting recombination events by analyzing aligned sequences of a given region of DNA that all originate from one species. The methods are based on phylogenetic, substitution, and compatibility analyses. Here we present the results of analyses of sequence data from different loci studied in several groups of evolutionarily distant species that are considered to be ancient asexuals, using seven different types of analysis. The groups of organisms were the arbuscular mycorrhizal fungi (Glomales), Darwinula stevensoni (Darwinuloidea crustacean ostracods) and the bdelloid rotifers (Bdelloidea), which are thought to have been asexual for the last 400, 25-100, and 35-40 Myr, respectively. The seven different analytical methods evaluated the evolutionary relationships among haplotypes, and these methods had previously been shown to be reliable for predicting the occurrence of recombination events. Despite the different degree of genetic variation among the different groups of organisms, at least some evidence for recombination was found in all species groups. In particular, predictions of recombination events in the arbuscular mycorrhizal fungi were frequent. Predictions of recombination were also found for sequence data that have previously been used to infer the absence of recombination in bdelloid rotifers. Although our results have to be taken with some caution because they could signal very ancient recombination events or possibly other genetic variation of nonrecombinant origin, they suggest that some cryptic recombination events may exist in these organisms.

Sex-Chromosome Homomorphy in Palearctic Tree Frogs Results from Both Turnovers and X-Y Recombination.

Contrasting with birds and mammals, poikilothermic vertebrates often have homomorphic sex chromosomes, possibly resulting from high rates of sex-chromosome turnovers and/or occasional X-Y recombination. Strong support for the latter mechanism was provided by four species of European tree frogs, which inherited from a common ancestor (∼5 Ma) the same pair of homomorphic sex chromosomes (linkage group 1, LG1), harboring the candidate sex-determining gene Dmrt1. Here, we test sex linkage of LG1 across six additional species of the Eurasian Hyla radiation with divergence times ranging from 6 to 40 Ma. LG1 turns out to be sex linked in six of nine resolved cases. Mapping the patterns of sex linkage to the Hyla phylogeny reveals several transitions in sex-determination systems within the last 10 My, including one switch in heterogamety. Phylogenetic trees of DNA sequences along LG1 are consistent with occasional X-Y recombination in all species where LG1 is sex linked. These patterns argue against one of the main potential causes for turnovers, namely the accumulation of deleterious mutations on nonrecombining chromosomes. Sibship analyses show that LG1 recombination is strongly reduced in males from most species investigated, including some in which it is autosomal. Intrinsically low male recombination might facilitate the evolution of male heterogamety, and the presence of important genes from the sex-determination cascade might predispose LG1 to become a sex chromosome.

FANCD2 Western blot as a diagnostic tool for Brazilian patients with Fanconi anemia

Fanconi anemia is a rare hereditary disease showing genetic heterogeneity due to a variety of mutations in genes involved in DNA repair pathways, which may lead to different clinical manifestations. Phenotypic variability makes diagnosis difficult based only on clinical manifestations, therefore laboratory tests are necessary. New advances in molecular pathogenesis of this disease led researchers to develop a diagnostic test based on Western blot for FANCD2. The objective of the present study was to determine the efficacy of this method for the diagnosis of 84 Brazilian patients with Fanconi anemia, all of whom tested positive for the diepoxybutane test, and 98 healthy controls. The FANCD2 monoubiquitinated isoform (FANCDS+/FANCD2L-) was not detected in 77 patients (91.7%). In 2 patients (2.4%), there was an absence of both the monoubiquitinated and the non-ubiquitinated proteins (FANCD2S-/FANCD2L-) and 5 patients (5.9%) had both isoforms (FANCD2S+/FANCD2L+). This last phenotype suggests downstream subtypes or mosaicism. All controls were diepoxybutane negative and were also negative on the FANCD2 Western blot. The Western blot for FANCD2 presented a sensitivity of 94% (79/84) and specificity of 100% (98/98). This method was confirmed as an efficient approach to screen Brazilian patients with deleterious mutations on FANCD2 (FANCD2S-/FANCD2L-) or other upstream genes of the FA/BRCA pathway (FANCDS+/FANCD2L-), to confirm the chromosome breakage test and to classify patients according to the level of FA/BRCA pathway defects. However, patients showing both FANCD2 isoforms (FANCD2S+/FANCD2L+) require additional studies to confirm mutations on downstream Fanconi anemia genes or the presence of mosaicism.

A hierarchical Bayesian model for predicting the functional consequences of amino-acid polymorphisms

Genetic polymorphisms in deoxyribonucleic acid coding regions may have a phenotypic effect on the carrier, e.g. by influencing susceptibility to disease. Detection of deleterious mutations via association studies is hampered by the large number of candidate sites; therefore methods are needed to narrow down the search to the most promising sites. For this, a possible approach is to use structural and sequence-based information of the encoded protein to predict whether a mutation at a particular site is likely to disrupt the functionality of the protein itself. We propose a hierarchical Bayesian multivariate adaptive regression spline (BMARS) model for supervised learning in this context and assess its predictive performance by using data from mutagenesis experiments on lac repressor and lysozyme proteins. In these experiments, about 12 amino-acid substitutions were performed at each native amino-acid position and the effect on protein functionality was assessed. The training data thus consist of repeated observations at each position, which the hierarchical framework is needed to account for. The model is trained on the lac repressor data and tested on the lysozyme mutations and vice versa. In particular, we show that the hierarchical BMARS model, by allowing for the clustered nature of the data, yields lower out-of-sample misclassification rates compared with both a BMARS and a frequen-tist MARS model, a support vector machine classifier and an optimally pruned classification tree.

Mutational meltdown in primary endosymbionts: selection limits Muller's ratchet

Background Primary bacterial endosymbionts of insects (p-endosymbionts) are thought to be undergoing the process of Muller's ratchet where they accrue slightly deleterious mutations due to genetic drift in small populations with negligible recombination rates. If this process were to go unchecked over time, theory predicts mutational meltdown and eventual extinction. Although genome degradation is common among p-endosymbionts, we do not observe widespread p-endosymbiont extinction, suggesting that Muller's ratchet may be slowed or even stopped over time. For example, selection may act to slow the effects of Muller's ratchet by removing slightly deleterious mutations before they go to fixation thereby causing a decrease in nucleotide substitutions rates in older p-endosymbiont lineages. Methodology/Principal Findings To determine whether selection is slowing the effects of Muller's ratchet, we determined the age of the Candidatus Riesia/sucking louse assemblage and analyzed the nucleotide substitution rates of several p-endosymbiont lineages that differ in the length of time that they have been associated with their insect hosts. We find that Riesia is the youngest p-endosymbiont known to date, and has been associated with its louse hosts for only 13–25 My. Further, it is the fastest evolving p-endosymbiont with substitution rates of 19–34% per 50 My. When comparing Riesia to other insect p-endosymbionts, we find that nucleotide substitution rates decrease dramatically as the age of endosymbiosis increases. Conclusions/Significance A decrease in nucleotide substitution rates over time suggests that selection may be limiting the effects of Muller's ratchet by removing individuals with the highest mutational loads and decreasing the rate at which new mutations become fixed. This countering effect of selection could slow the overall rate of endosymbiont extinction.