Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria

In populations that are small and asexual, mutations with slight negative effects on fitness will drift to fixation more often than in large or sexual populations in which they will be eliminated by selection. If such mutations occur in substantial numbers, the combined effects of long-term asexuality and small population size may result in substantial accumulation of mildly deleterious substitutions. Prokaryotic endosymbionts of animals that are transmitted maternally for very long periods are effectively asexual and experience smaller effective population size than their free-living relatives. The contrast between such endosymbionts and related free-living bacteria allows us to test whether a population structure imposing frequent bottlenecks and asexuality does lead to an accumulation of slightly deleterious substitutions. Here we show that several independently derived insect endosymbionts, each with a long history of maternal transmission, have accumulated destabilizing base substitutions in the highly conserved 16S rRNA. Stabilities of Domain I of this subunit are 15–25% lower in endosymbionts than in closely related free-living bacteria. By mapping destabilizing substitutions onto a reconstructed phylogeny, we show that decreased ribosomal stability has evolved separately in each endosymbiont lineage. Our phylogenetic approach allows us to demonstrate statistical significance for this pattern: becoming endosymbiotic predictably results in decreased stability of rRNA secondary structure.

The evolution of XY recombination: sexually antagonistic selection versus deleterious mutation load.

Recombination arrest between X and Y chromosomes, driven by sexually antagonistic genes, is expected to induce their progressive differentiation. However, in contrast to birds and mammals (which display the predicted pattern), most cold-blooded vertebrates have homomorphic sex chromosomes. Two main hypotheses have been proposed to account for this, namely high turnover rates of sex-determining systems and occasional XY recombination. Using individual-based simulations, we formalize the evolution of XY recombination (here mediated by sex reversal; the "fountain-of-youth" model) under the contrasting forces of sexually antagonistic selection and deleterious mutations. The shift between the domains of elimination and accumulation occurs at much lower selection coefficients for the Y than for the X. In the absence of dosage compensation, mildly deleterious mutations accumulating on the Y depress male fitness, thereby providing incentives for XY recombination. Under our settings, this occurs via "demasculinization" of the Y, allowing recombination in XY (sex-reversed) females. As we also show, this generates a conflict with the X, which coevolves to oppose sex reversal. The resulting rare events of XY sex reversal are enough to purge the Y from its load of deleterious mutations. Our results support the "fountain of youth" as a plausible mechanism to account for the maintenance of sex-chromosome homomorphy.

Deleterious mutation accumulation in asexual Timema stick insects.

Sexual reproduction is extremely widespread in spite of its presumed costs relative to asexual reproduction, indicating that it must provide significant advantages. One postulated benefit of sex and recombination is that they facilitate the purging of mildly deleterious mutations, which would accumulate in asexual lineages and contribute to their short evolutionary life span. To test this prediction, we estimated the accumulation rate of coding (nonsynonymous) mutations, which are expected to be deleterious, in parts of one mitochondrial (COI) and two nuclear (Actin and Hsp70) genes in six independently derived asexual lineages and related sexual species of Timema stick insects. We found signatures of increased coding mutation accumulation in all six asexual Timema and for each of the three analyzed genes, with 3.6- to 13.4-fold higher rates in the asexuals as compared with the sexuals. In addition, because coding mutations in the asexuals often resulted in considerable hydrophobicity changes at the concerned amino acid positions, coding mutations in the asexuals are likely associated with more strongly deleterious effects than in the sexuals. Our results demonstrate that deleterious mutation accumulation can differentially affect sexual and asexual lineages and support the idea that deleterious mutation accumulation plays an important role in limiting the long-term persistence of all-female lineages.

Predicting the deleterious effects of mutation load in fragmented populations.

Human-induced habitat fragmentation constitutes a major threat to biodiversity. Both genetic and demographic factors combine to drive small and isolated populations into extinction vortices. Nevertheless, the deleterious effects of inbreeding and drift load may depend on population structure, migration patterns, and mating systems and are difficult to predict in the absence of crossing experiments. We performed stochastic individual-based simulations aimed at predicting the effects of deleterious mutations on population fitness (offspring viability and median time to extinction) under a variety of settings (landscape configurations, migration models, and mating systems) on the basis of easy-to-collect demographic and genetic information. Pooling all simulations, a large part (70%) of variance in offspring viability was explained by a combination of genetic structure (F(ST)) and within-deme heterozygosity (H(S)). A similar part of variance in median time to extinction was explained by a combination of local population size (N) and heterozygosity (H(S)). In both cases the predictive power increased above 80% when information on mating systems was available. These results provide robust predictive models to evaluate the viability prospects of fragmented populations.

Sex-chromosome turnovers induced by deleterious mutation load.

In sharp contrast with mammals and birds, many cold-blooded vertebrates present homomorphic sex chromosomes. Empirical evidence supports a role for frequent turnovers, which replace nonrecombining sex chromosomes before they have time to decay. Three main mechanisms have been proposed for such turnovers, relying either on neutral processes, sex-ratio selection, or intrinsic benefits of the new sex-determining genes (due, e.g., to linkage with sexually antagonistic mutations). Here, we suggest an additional mechanism, arising from the load of deleterious mutations that accumulate on nonrecombining sex chromosomes. In the absence of dosage compensation, this load should progressively lower survival rate in the heterogametic sex. Turnovers should occur when this cost outweighs the benefits gained from any sexually antagonistic genes carried by the nonrecombining sex chromosome. We use individual-based simulations of a Muller's ratchet process to test this prediction, and investigate how the relevant parameters (effective population size, strength and dominance of deleterious mutations, size of nonrecombining segment, and strength of sexually antagonistic selection) are expected to affect the rate of turnovers.

Expansion load: recessive mutations and the role of standing genetic variation

Expanding populations incur a mutation burden – the so-called expansion load. Previous studies of expansion load have focused on codominant mutations. An important consequence of this assumption is that expansion load stems exclusively from the accumulation of new mutations occurring in individuals living at the wave front. Using individual-based simulations, we study here the dynamics of standing genetic variation at the front of expansions, and its consequences on mean fitness if mutations are recessive. We find that deleterious genetic diversity is quickly lost at the front of the expansion, but the loss of deleterious mutations at some loci is compensated by an increase of their frequencies at other loci. The frequency of deleterious homozygotes therefore increases along the expansion axis, whereas the average number of deleterious mutations per individual remains nearly constant across the species range. This reveals two important differences to codominant models: (i) mean fitness at the front of the expansion drops much faster if mutations are recessive, and (ii) mutation load can increase during the expansion even if the total number of deleterious mutations per individual remains constant. We use our model to make predictions about the shape of the site frequency spectrum at the front of range expansion, and about correlations between heterozygosity and fitness in different parts of the species range. Importantly, these predictions provide opportunities to empirically validate our theoretical results. We discuss our findings in the light of recent results on the distribution of deleterious genetic variation across human populations and link them to empirical results on the correlation of heterozygosity and fitness found in many natural range expansions.

Phenotypic effect of mutations in evolving populations of RNA molecules

Background. The secondary structure of folded RNA sequences is a good model to map phenotype onto genotype, as represented by the RNA sequence. Computational studies of the evolution of ensembles of RNA molecules towards target secondary structures yield valuable clues to the mechanisms behind adaptation of complex populations. The relationship between the space of sequences and structures, the organization of RNA ensembles at mutation-selection equilibrium, the time of adaptation as a function of the population parameters, the presence of collective effects in quasispecies, or the optimal mutation rates to promote adaptation all are issues that can be explored within this framework. Results. We investigate the effect of microscopic mutations on the phenotype of RNA molecules during their in silico evolution and adaptation. We calculate the distribution of the effects of mutations on fitness, the relative fractions of beneficial and deleterious mutations and the corresponding selection coefficients for populations evolving under different mutation rates. Three different situations are explored: the mutation-selection equilibrium (optimized population) in three different fitness landscapes, the dynamics during adaptation towards a goal structure (adapting population), and the behavior under periodic population bottlenecks (perturbed population). Conclusions. The ratio between the number of beneficial and deleterious mutations experienced by a population of RNA sequences increases with the value of the mutation rate µ at which evolution proceeds. In contrast, the selective value of mutations remains almost constant, independent of µ, indicating that adaptation occurs through an increase in the amount of beneficial mutations, with little variations in the average effect they have on fitness. Statistical analyses of the distribution of fitness effects reveal that small effects, either beneficial or deleterious, are well described by a Pareto distribution. These results are robust under changes in the fitness landscape, remarkably when, in addition to selecting a target secondary structure, specific subsequences or low-energy folds are required. A population perturbed by bottlenecks behaves similarly to an adapting population, struggling to return to the optimized state. Whether it can survive in the long run or whether it goes extinct depends critically on the length of the time interval between bottlenecks. © 2010 Stich et al; licensee BioMed Central Ltd.

Structural aspects of the p.P222Q homozygous mutation of HSD3B2 gene in a patient with congenital adrenal hyperplasia

Type II 3β-hydroxysteroid dehydrogenase/Δ5-Δ4-isomerase (3β-HSD2), encoded by the HSD3B2 gene, is a key enzyme involved in the biosynthesis of all the classes of steroid hormones. Deleterious mutations in the HSD3B2 gene cause the classical deficiency of 3β-HSD2, which is a rare autosomal recessive disease that leads to congenital adrenal hyperplasia (CAH). CAH is the most frequent cause of ambiguous genitalia and adrenal insufficiency in newborn infants with variable degrees of salt losing. Here we report the molecular and structural analysis of the HSD3B2 gene in a 46,XY child, who was born from consanguineous parents, and presented with ambiguous genitalia and salt losing. The patient carries a homozygous nucleotide c.665C>A change in exon 4 that putatively substitutes the proline at codon 222 for glutamine. Molecular homology modeling of normal and mutant 3β-HSD2 enzymes emphasizes codon 222 as an important residue for the folding pattern of the enzyme and validates a suitable model for analysis of new mutations.

The constancy of global regulation across a species: the concentrations of ppGpp and RpoS are strain-specific in Escherichia coli

Background: Sigma factors and the alarmone ppGpp control the allocation of RNA polymerase to promoters under stressful conditions. Both ppGpp and the sigma factor sigma(S) (RpoS) are potentially subject to variability across the species Escherichia coli. To find out the extent of strain variation we measured the level of RpoS and ppGpp using 31 E. coli strains from the ECOR collection and one reference K-12 strain. Results: Nine ECORs had highly deleterious mutations in rpoS, 12 had RpoS protein up to 7-fold above that of the reference strain MG1655 and the remainder had comparable or lower levels. Strain variation was also evident in ppGpp accumulation under carbon starvation and spoT mutations were present in several low-ppGpp strains. Three relationships between RpoS and ppGpp levels were found: isolates with zero RpoS but various ppGpp levels, strains where RpoS levels were proportional to ppGpp and a third unexpected class in which RpoS was present but not proportional to ppGpp concentration. High-RpoS and high-ppGpp strains accumulated rpoS mutations under nutrient limitation, providing a source of polymorphisms. Conclusions: The ppGpp and sigma(S) variance means that the expression of genes involved in translation, stress and other traits affected by ppGpp and/or RpoS are likely to be strain-specific and suggest that influential components of regulatory networks are frequently reset by microevolution. Different strains of E. coli have different relationships between ppGpp and RpoS levels and only some exhibit a proportionality between increasing ppGpp and RpoS levels as demonstrated for E. coli K-12.

Positive Epistasis Drives the Acquisition of Multidrug Resistance

The evolution of multiple antibiotic resistance is an increasing global problem. Resistance mutations are known to impair fitness, and the evolution of resistance to multiple drugs depends both on their costs individually and on how they interact-epistasis. Information on the level of epistasis between antibiotic resistance mutations is of key importance to understanding epistasis amongst deleterious alleles, a key theoretical question, and to improving public health measures. Here we show that in an antibiotic-free environment the cost of multiple resistance is smaller than expected, a signature of pervasive positive epistasis among alleles that confer resistance to antibiotics. Competition assays reveal that the cost of resistance to a given antibiotic is dependent on the presence of resistance alleles for other antibiotics. Surprisingly we find that a significant fraction of resistant mutations can be beneficial in certain resistant genetic backgrounds, that some double resistances entail no measurable cost, and that some allelic combinations are hotspots for rapid compensation. These results provide additional insight as to why multi-resistant bacteria are so prevalent and reveal an extra layer of complexity on epistatic patterns previously unrecognized, since it is hidden in genome-wide studies of genetic interactions using gene knockouts.

Impact of chromossomal structure on the evolution of Schizosaccharomyces pombe undergoing Mutation Accumulation

Large chromosomal rearrangements are common in natural populations and thought to be involved in speciation events. In this project, we used experimental evolution to determine how the speed of evolution and the type of accumulated mutations depend on the ancestral chromosomal structure and genotype. We utilized two Wild Type strains and a set of genetically engineered Schizosaccharomyces pombe strains, different solely in the presence of a certain type of chromosomal variant (inversions or translocations), along with respective controls. Previous research has shown that these chromosomal variants have different fitness levels in several environments, probably due to changes in the gene expression along the genome. These strains were propagated in the laboratory at very low population sizes, in which we expect natural selection to be less efficient at purging deleterious mutations. We then measured these strains’ changes in fitness throughout this accumulation of deleterious mutations, comparing the evolutionary trajectories in the different rearrangements to understand if the chromosomal structure affected the speed of evolution. We also tested these mutations for possible epistatic effects and estimated their parameters: the number of arising deleterious mutations per generation (Ud) and each one’s mean effect (sd).

Landscape genetics of the Alpine newt (Mesotriton alpestris) inferred from a strip-based approach

Habitat destruction and fragmentation are known to strongly affect dispersal by altering the quality of the environment between populations. As a consequence, lower landscape connectivity is expected to enhance extinction risks through a decrease in gene flow and the resulting negative effects of genetic drift, accumulation of deleterious mutations and inbreeding depression. Such phenomena are particularly harmful for amphibian species, characterized by disjunct breeding habitats. The dispersal behaviour of amphibians being poorly understood, it is crucial to develop new tools, allowing us to determine the influence of landscape connectivity on the persistence of populations. In this study, we developed a new landscape genetics approach that aims at identifying land-uses affecting genetic differentiation, without a priori assumptions about associated ecological costs. We surveyed genetic variation at seven microsatellite loci for 19 Alpine newt (Mesotriton alpestris) populations in western Switzerland. Using strips of varying widths that define a dispersal corridor between pairs of populations, we were able to identify land-uses that act as dispersal barriers (i.e. urban areas) and corridors (i.e. forests). Our results suggest that habitat destruction and landscape fragmentation might in the near future affect common species such as M. alpestris. In addition, by identifying relevant landscape variables influencing population structure without unrealistic assumptions about dispersal, our method offers a simple and flexible tool of investigation as an alternative to least-cost models and other approaches.

Multiple causes of the evolution of dispersal in metapopulations

Abstract This work investigates the outcome of the interaction of the multiple causes of selection acting on dispersal in metapopulations. Dispersal, defined here as the ability of individuals to move out of their natal population to reproduce in an other one, has three main causes. First, population variability, as caused by random population extinctions, induces high incentives to disperse through the probability to recolonize an empty population and thus to escape competition for space. This adds to the second cause, kin competition avoidance where individuals in a crowded patch will benefit from the release of competition with relatives caused by dispersal. Dispersal may thus be viewed as an altruistic act. Third, dispersal might evolve as a strategy of avoiding inbred matings which are expected to bear fitness costs due to the presence of a mutation load. The interaction of inbreeding avoidance and kin competition is explored in chapter 2. Conditions conducive to the establishment of a high relatedness within population are expected to induce high dispersal through both kin competition avoidance and inbreeding avoidance. However, the dynamics of inbreeding depression is bound to depend on the level of gene flow as well as on the deleterious mutation parameters. Mutations more prone to settle a high level of inbreeding depression will select for increased dispersal. Chapter 3 investigates the effect of the mating system on the joint dynamics of dispersal and inbreeding depression. Higher inbreeding rates as those found in various mating systems lead to a more efficient purge of the deleterious mutations. However, this decrease in the costs of inbreeding are usually accompanied by a higher within deme relatedness which balances the decreased effect of inbreeding avoidance on the evolution of dispersal. Finally, population turnover, as found in most natural populations has a dual effect on dispersal. Indeed, it increases dispersal by the increased probability of winning a breeding slot in extinct demes it creates but, on the other hand, it counter-selects for dispersal through the slow establishment of unsaturated demic conditions which contribute to lower the local competition for space. Résumé Ce travail se propose d'étudier les effets conjoints des multiples causes de l'évolution de la dispersion en métapopulation. La dispersion, définie ici comme étant la capacité de quitter sa population d'origine pour se reproduire dans une antre population, possède trois principales causes. Premièrement, l'extinction aléatoire de populations sélectionne pour plus de dispersion car elle augmente la Probabilité de recoloniser un patch éteint et donc d'échapper à la compétition locale. La seconde cause, l'évitement de la compétition de parentèle, sélectionne pour plus de dispersion par les bénéfices qu'elle apporte par diminution de la compétition entre individus apparentés. Troisièmement, la dispersion évolue "comme stratégie d'évitement de la dépression de consanguinité présente dans des petites populations isolées. L'interaction entre l'évitement de la consanguinité et de la compétition de parentèle est étudiée dans le chapitre 2. Les conditions conduisant à l'établissement d'un fort apparentement à l'intérieur des populations sont celles qui génèrent le plus de sélection pour la dispersion. Cependant, la dynamique de la dépression de consanguinité est dépendante de la dispersion entre populations ainsi que des paramètres des mutations délétères. Les mutations créant le plus de dépression de consanguinité sont celles qui sélectionneront le plus pour de la dispersion. Le chapitre 3 s'intéresse aux effets du système de reproduction sur la dynamique conjointe du fardeau de mutation et de la dispersion. La purge des mutations délétère étant plus sévère dans des conditions de forte consanguinité, elle diminue les coûts de la consanguinité mais est habituellement accompagné par une augmentation de l'apparentement et donc l'effet peut être neutre sur la dispersion. Finalement, le turnover de populations a un effet dual sur la dispersion. La dispersion est sélectionnée par l'augmentation de la probabilité de gagner une place de reproduction dans des patchs éteints mais elle est également contre sélectionnée par la désaturation des patchs causée par l'extinction et la diminution de la compétition pour l'espace qui intervient dans ce cas.

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.