Fixation Of Mutations Research Articles

The SARS-CoV nsp12 Polymerase Active Site Is Tuned for Large-Genome Replication.

ABSTRACTPositive-strand RNA viruses replicate their genomes using virally encoded RNA-dependent RNA polymerases (RdRP) with a common active-site structure and closure mechanism upon which replication speed and fidelity can evolve to optimize virus fitness. Coronaviruses (CoV) form large multicomponent RNA replication-transcription complexes containing a core RNA synthesis machine made of the nsp12 RdRP protein with one nsp7 and two nsp8 proteins as essential subunits required for activity. We show that assembly of this complex can be accelerated 5-fold by preincubation of nsp12 with nsp8 and further optimized with the use of a novel nsp8L7 heterodimer fusion protein construct. Using rapid kinetics methods, we measure elongation rates of up to 260 nucleotides (nt)/s for the core replicase, a rate that is unusually fast for a viral polymerase. To address the origin of this fast rate, we examined the roles of two CoV-specific residues in the RdRP active site: Ala547, which replaces a conserved glutamate above the bound NTP, and Ser759, which mutates the palm domain GDD sequence to SDD. Our data show that Ala547 allows for a doubling of replication rate, but this comes at a fidelity cost that is mitigated by using a SDD sequence in the palm domain. Our biochemical data suggest that fixation of mutations in polymerase motifs F and C played a key role in nidovirus evolution by tuning replication rate and fidelity to accommodate their large genomes.IMPORTANCE Replicating large genomes represents a challenge for RNA viruses because fast RNA synthesis is needed to escape innate immunity defenses, but faster polymerases are inherently low-fidelity enzymes. Nonetheless, the coronaviruses replicate their ≈30-kb genomes using the core polymerase structure and mechanism common to all positive-strand RNA viruses. The classic explanation for their success is that the large-genome nidoviruses have acquired an exonuclease-based repair system that compensates for the high polymerase mutation rate. In this work, we establish that the nidoviral polymerases themselves also play a key role in maintaining genome integrity via mutations at two key active-site residues that enable very fast replication rates while maintaining typical mutation rates. Our findings further demonstrate the evolutionary plasticity of the core polymerase platform by showing how it has adapted during the expansion from short-genome picornaviruses to long-genome nidoviruses.

Fitness cost of target-site and metabolic resistance to pyrethroids drives restoration of susceptibility in a highly resistant Anopheles gambiae population from Uganda.

BackgroundInsecticide resistance threatens the effectiveness of malaria vector control, calling for an urgent need to design suitable resistance management strategies. Here, we established the resistance profiling of an Ugandan Anopheles gambiae population to insecticides using WHO procedures and assessed the potential restoration of susceptibility in the hybrid line Mayuge/KISUMU in an insecticide-free environment for eighteen (18) generations.ResultsThis An gambiae population exhibited a very high intensity of resistance to permethrin, deltamethrin, and alphacypermethrin with a consistent loss of efficacy of all long-lasting insecticidal nets (LLINs) tested including PBO-based and new generation nets Interceptor G2 (IG2) and Royal guard. Molecular analysis revealed a fixation of the L1014S-kdr mutation together with the overexpression of some P450 metabolic genes (CYP6Z1, CYP9K1, CYP6P1, 3 & 4) besides the cuticular resistance-related genes (CYP4G16) and sensorial appendage proteins (SAP1, SAP2, and SAP3) but no GSTe2 overexpression. In the absence of selection pressure, the mortality rate after exposure to insecticides increased significantly over generations, and restoration of susceptibility was observed for most of the insecticides in less than 10 generations. Accordingly, a significant reduction in the frequency of KdrE was observed after 13 generations coupled with reduced expression of most metabolic resistance genes.ConclusionsThe results of this study show that the high intensity of pyrethroid resistance observed in An gambiae from Uganda associated with the loss of efficacy of LLINs could compromise vector control efforts. The study also highlights that an early rotation of insecticides could help manage resistance to insecticides by restoring the susceptibility. However, the persistence of Kdr mutation together with overexpression of some metabolic genes after many generations in the absence of selection pressure indicates the potential implication of modifiers alleviating the cost of resistance which needs to be further investigated.

Grey wolf genomic history reveals a dual ancestry of dogs

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived1–8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000–30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

Narrow transmission bottlenecks and limited within-host viral diversity during a SARS-CoV-2 outbreak on a fishing boat.

The long-term evolution of viruses is ultimately due to viral mutants that arise within infected individuals and transmit to other individuals. Here, we use deep sequencing to investigate the transmission of viral genetic variation among individuals during a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak that infected the vast majority of crew members on a fishing boat. We deep-sequenced nasal swabs to characterize the within-host viral population of infected crew members, using experimental duplicates and strict computational filters to ensure accurate variant calling. We find that within-host viral diversity is low in infected crew members. The mutations that did fix in some crew members during the outbreak are not observed at detectable frequencies in any of the sampled crew members in which they are not fixed, suggesting that viral evolution involves occasional fixation of low-frequency mutations during transmission rather than persistent maintenance of within-host viral diversity. Overall, our results show that strong transmission bottlenecks dominate viral evolution even during a superspreading event with a very high attack rate.

Abstract 718: Significantly strong protective efficacy of silibinin against basal cell carcinoma growth and progression

Abstract Basal cell carcinoma (BCC) is the most common skin malignancy accounting for ∼80% of non-melanoma skin cancers (NMSCs). It arises in basal cells of skin that line the deepest layer of the epidermis. The most important risk factor for BCC is solar ultraviolet B (UVB) radiation that results in DNA damage (CPDs formation) in skin epidermis, which, if not repaired, leads to fixation of mutations and initiation of skin carcinogenesis. One of the key molecular features of BCC is sustained activation of Hedgehog signaling (Hh) pathway through inactivating mutations in tumor suppressor gene Patched (Ptch) or activating mutations in Smoothened (Smo). Therefore, extensive efforts have been made to target activated Hh pathway for the treatment of BCC, though with toxic side effects. With regards to chemoprevention of BCC, the agents which target the events associated with UVB-induced DNA damage repair, together with targeting promotion/progression stages, might provide more effective broad-spectrum opportunities to intervene at the earliest. Utilizing the well-established patched (Ptch)+/- mouse model of UVB radiation-induced BCC formation, our approach in this study was to target BCC development and its prevention by employing silibinin (a natural flavonolignan from milk thistle seeds). Previously, we have extensively reported the efficacy of silibinin against UVB-induced photodamage and photo-carcinogenesis, and also recently reported silibinin efficacy against BCC growth in cell culture (including anti-BCC drug-resistant cell lines) and mouse BCC allograft tumors. In this study, both male and female (Ptch)+/- mice were irradiated with 240 mJ/cm2 UVB dose, 3 times per week (M, W, F) for 26 and 46 weeks with or without silibinin. Silibinin (9 mg/200μl of acetone) was applied topically, 30 min post and pre-UVB exposure. Our study shows that chronic UVB exposure induced BCCs in Ptch+/- mice. Treatment with silibinin post and pre-UVB exposure for 26 weeks decreased BCC lesion numbers (39-65%), and cross-sectional area (45-72%) p<0.001, respectively compared to UVB alone. Furthermore, continuous UVB exposure up to 46 weeks increased the BCC lesion number by ∼6 folds and cross-sectional area by ∼3.4 folds (p<0.001) respectively. Notably, silibinin (irrespective of post and pre-UVB treatment) significantly halted the progression of BCC (81-94%, p<0.001), even upon prolonged UVB exposure. Histological analysis showed increased dysplasia, fibro-sarcoma, and squamous cell carcinoma upon UVB exposure, which significantly decreased upon silibinin treatment. In addition, mechanistic studies revealed that silibinin significantly decreased basal cell proliferation (Ki-67) and the expression of cytokeratins (14 and 15), and Hh signaling mediators Smo and Gli1 in the BCC lesions. Together, our findings demonstrate that silibinin has the potential to prevent the growth and progression of UVB-induced BCC. Citation Format: Sandeep Paudel, Komal Raina, Vasundhara R. Tiku, Akhilendra Maurya, David J. Orlicky, Zhiying You, Cindy M. Rigby, Gagan Deep, Rama Kant, Bupinder Raina, Chapla Agarwal, Rajesh Agarwal. Significantly strong protective efficacy of silibinin against basal cell carcinoma growth and progression [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 718.

A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication.

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, “marking” one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.

Optimal pathways control fixation of multiple mutations during cancer initiation

Cancer starts after initially healthy tissue cells accumulate several specific mutations or other genetic alterations. The dynamics of tumor formation is a very complex phenomenon due to multiple involved biochemical and biophysical processes. It leads to a very large number of possible pathways on the road to final fixation of all mutations that marks the beginning of the cancer, complicating the understanding of microscopic mechanisms of tumor formation. We present a new theoretical framework of analyzing the cancer initiation dynamics by exploring the properties of effective free-energy landscape of the process. It is argued that although there are many possible pathways for the fixation of all mutations in the system, there are only a few dominating pathways on the road to tumor formation. The theoretical approach is explicitly tested in the system with only two mutations using analytical calculations and Monte Carlo computer simulations. Excellent agreement with theoretical predictions is found for a large range of parameters, supporting our hypothesis and allowing us to better understand the mechanisms of cancer initiation. Our theoretical approach clarifies some important aspects of microscopic processes that lead to tumor formation.

Superinfection exclusion: A viral strategy with short-term benefits and long-term drawbacks.

Viral superinfection occurs when multiple viral particles subsequently infect the same host. In nature, several viral species are found to have evolved diverse mechanisms to prevent superinfection (superinfection exclusion) but how this strategic choice impacts the fate of mutations in the viral population remains unclear. Using stochastic simulations, we find that genetic drift is suppressed when superinfection occurs, thus facilitating the fixation of beneficial mutations and the removal of deleterious ones. Interestingly, we also find that the competitive (dis)advantage associated with variations in life history parameters is not necessarily captured by the viral growth rate for either infection strategy. Putting these together, we then show that a mutant with superinfection exclusion will easily overtake a superinfecting population even if the latter has a much higher growth rate. Our findings suggest that while superinfection exclusion can negatively impact the long-term adaptation of a viral population, in the short-term it is ultimately a winning strategy.

Trends of Transmitted and Acquired Drug Resistance in Europe From 1981 to 2019: A Comparison Between the Populations of Late Presenters and Non-late Presenters.

BackgroundThe increased use of antiretroviral therapy (ART) has decreased mortality and morbidity of HIV-1 infected people but increasing levels of HIV drug resistance threatens the success of ART regimens. Conversely, late presentation can impact treatment outcomes, health costs, and potential transmission of HIV.ObjectiveTo describe the patterns of transmitted drug resistance (TDR) and acquired drug resistance (ADR) in HIV-1 infected patients followed in Europe, to compare its patterns in late presenters (LP) vs non-late presenters (NLP), and to analyze the most prevalent drug resistance mutations among HIV-1 subtypes.MethodsOur study included clinical, socio-demographic, and genotypic information from 26,973 HIV-1 infected patients from the EuResist Integrated Database (EIDB) between 1981 and 2019.ResultsAmong the 26,973 HIV-1 infected patients in the analysis, 11,581 (42.9%) were ART-naïve patients and 15,392 (57.1%) were ART-experienced. The median age was 37 (IQR: 27.0–45.0) years old and 72.6% were males. The main transmission route was through heterosexual contact (34.9%) and 81.7% of patients originated from Western Europe. 71.9% of patients were infected by subtype B and 54.8% of patients were classified as LP. The overall prevalence of TDR was 12.8% and presented an overall decreasing trend (p for trend < 0.001), the ADR prevalence was 68.5% also with a decreasing trend (p for trend < 0.001). For LP and NLP, the TDR prevalence was 12.3 and 12.6%, respectively, while for ADR, 69.9 and 68.2%, respectively. The most prevalent TDR drug resistance mutations, in both LP and NLP, were K103N/S, T215rev, T215FY, M184I/V, M41I/L, M46I/L, and L90M.ConclusionOur study showed that the overall TDR (12.8%) and ADR (68.5%) presented decreasing trends during the study time period. For LP, the overall TDR was slightly lower than for NLP (12.3 vs 12.6%, respectively); while this pattern was opposite for ADR (LP slightly higher than NLP). We suggest that these differences, in the case of TDR, can be related to the dynamics of fixation of drug resistance mutations; and in the case of ADR with the more frequent therapeutic failure in LPs.

Selection on the gametophyte: Modeling alternation of generations in plants.

PremiseThe degree of gametophyte dependence on the sporophyte life stage is a major feature that differentiates the life cycles of land plants, yet the evolutionary consequences of this difference remain poorly understood. Most evolutionary models assume organisms are either haploid or diploid for their entire lifespan, which is not appropriate for simulating plant life cycles. Here, we introduce shadie (Simulating Haploid–Diploid Evolution), a new, simple Python program for implementing simulations with biphasic life cycles and analyzing their results, using SLiM 3 as a simulation back end.MethodsWe implemented evolutionary simulations under three realistic plant life cycle models supported in shadie, using either standardized or biologically realistic parameter settings to test how variation in plant life cycles and sexual systems affects patterns of genome diversity.ResultsThe dynamics of single beneficial mutation fixation did not vary dramatically between different models, but the patterns of spatial variation did differ, demonstrating that different life histories and model parameters affect both genetic diversity and linkage disequilibrium. The rate of linkage disequilibrium decay away from selected sites varied depending on model parameters such as cloning and selfing rates, through their impact on effective population sizes.DiscussionEvolutionary simulations are an exciting, underutilized approach in evolutionary research and education. shadie can aid plant researchers in developing null hypotheses, examining theory, and designing empirical studies, in order to investigate the role of the gametophyte life stage, and the effects of variation in plant life cycles, on plant genome evolution.

The effects of weak selection on neutral diversity at linked sites.

The effects of selection on variability at linked sites have an important influence on levels and patterns of within-population variation across the genome. Most theoretical models of these effects have assumed that selection is sufficiently strong that allele frequency changes at the loci concerned are largely deterministic. These models have led to the conclusion that directional selection for selectively favorable mutations, or against recurrent deleterious mutations, reduces nucleotide site diversity at linked neutral sites. Recent work has shown, however, that fixations of weakly selected mutations, accompanied by significant stochastic changes in allele frequencies, can sometimes cause higher diversity at linked sites when compared with the effects of fixations of neutral mutations. This study extends this work by deriving approximate expressions for the mean conditional times to fixation and loss of mutations subject to selection, and analyzing the conditions under which selection increases rather than reduces these times. Simulations are used to examine the relations between diversity at a neutral site and the fixation and loss times of mutations at a linked site that is subject to selection. It is shown that the long-term level of neutral diversity can be increased over the purely neutral value by recurrent fixations and losses of linked, weakly selected dominant or partially dominant favorable mutations, or linked recessive or partially recessive deleterious mutations. The results are used to examine the conditions under which associative overdominance, as opposed to background selection, is likely to operate.

The role of mitonuclear incompatibilities in allopatric speciation.

Aerobic metabolism in eukaryotic cells requires extensive interactions between products of the nuclear and mitochondrial genomes. Rapid evolution of the mitochondrial genome, including fixation of both adaptive and deleterious mutations, creates intrinsic selection pressures favoring nuclear gene mutations that maintain mitochondrial function. As this process occurs independently in allopatry, the resulting divergence between conspecific populations can subsequently be manifest in mitonuclear incompatibilities in inter-population hybrids. Such incompatibilities, mitonuclear versions of Bateson-Dobzhansky-Muller incompatibilities that form the standard model for allopatric speciation, can potentially restrict gene flow between populations, ultimately resulting in varying degrees of reproductive isolation. The potential role of mitonuclear incompatibilities in speciation is further enhanced where mtDNA substitution rates are elevated compared to the nuclear genome and where population structure maintains allopatry for adequate time to evolve multiple mitonuclear incompatibilities. However, the fact that mitochondrial introgression occurs across species boundaries has raised questions regarding the efficacy of mitonuclear incompatibilities in reducing gene flow. Several scenarios now appear to satisfactorily explain this phenomenon, including cases where differences in mtDNA genetic load may drive introgression or where co-introgression of coadapted nuclear genes may support the function of introgressed mtDNA. Although asymmetries in reproductive isolation between taxa are consistent with mitonuclear incompatibilities, interactions between autosomes and sex chromosomes yield similar predictions that are difficult to disentangle. With regard to establishing reproductive isolation while in allopatry, existing studies clearly suggest that mitonuclear incompatibilities can contribute to the evolution of barriers to gene flow. However, there is to date relatively little definitive evidence supporting a primary role for mitonuclear incompatibilities in the speciation process.

Experimental Evolution Expands the Breadth of Adaptation to an Environmental Gradient Correlated With Genome Reduction.

Whether and how adaptive evolution adjusts the breadth of adaptation in coordination with the genome are essential issues for connecting evolution with ecology. To address these questions, experimental evolution in five Escherichia coli strains carrying either the wild-type genome or a reduced genome was performed in a defined minimal medium (C0). The ancestral and evolved populations were subsequently subjected to fitness and chemical niche analyses across an environmental gradient with 29 combinations of eight chemical components of the minimal medium. The results showed that adaptation was achieved not only specific to the evolutionary condition (C0), but also generally, to the environmental gradient; that is, the breadth of adaptation to the eight chemical niches was expanded. The magnitudes of the adaptive improvement and the breadth increase were both correlated with genome reduction and were highly significant in two out of eight niches (i.e., glucose and sulfate). The direct adaptation-induced correlated adaptation to the environmental gradient was determined by only a few genome mutations. An additive increase in fitness associated with the stepwise fixation of mutations was consistently observed in the reduced genomes. In summary, this preliminary survey demonstrated that evolution finely tuned the breadth of adaptation correlated with genome reduction.

Mutation accumulation in inbreeding populations under evolution of the selfing rate.

It is theoretically established that self-fertilization can facilitate mutation accumulation, thus increasing extinction risk. However, in previous studies, selfing rates are often set as fixed parameters, but in natural systems, evolution of selfing rates and deleterious mutations may mutually affect each other. I carried out simulations to investigate the dynamics of selfing rates and mutation accumulation, by allowing deleterious mutations to coevolve with alleles that modify the selfing rate (selfing modifiers). I found that selfing rates will often fluctuate over time, due to successive invasion of alleles that increase selfing and outcrossing. Since mutation fixation is mainly caused by Muller's ratchet, its rate is sensitive to the change of the selfing rate mutations will accumulate in a punctuated pattern. The dynamics are influenced by several factors, such as recombination and the selfing rate effects of selfing modifier loci. Also, such temporal variation produces variation of selfing rates and mutation accumulation rates between multiple conspecific populations, which can increase the average fitness across populations. As factors, such as the genomic mutation rate of deleterious mutations, can simultaneously influence the selfing rate and mutation fixation, effects of these factors on mutation accumulation rates can be complicated and non-monotonic.

Slow Recovery from Inbreeding Depression Generated by the Complex Genetic Architecture of Segregating Deleterious Mutations.

The deleterious effects of inbreeding have been of extreme importance to evolutionary biology, but it has been difficult to characterize the complex interactions between genetic constraints and selection that lead to fitness loss and recovery after inbreeding. Haploid organisms and selfing organisms like the nematode Caenorhabditis elegans are capable of rapid recovery from the fixation of novel deleterious mutation; however, the potential for recovery and genomic consequences of inbreeding in diploid, outcrossing organisms are not well understood. We sought to answer two questions: 1) Can a diploid, outcrossing population recover from inbreeding via standing genetic variation and new mutation? and 2) How does allelic diversity change during recovery? We inbred C. remanei, an outcrossing relative of C. elegans, through brother-sister mating for 30 generations followed by recovery at large population size. Inbreeding reduced fitness but, surprisingly, recovery from inbreeding at large populations sizes generated only very moderate fitness recovery after 300 generations. We found that 65% of ancestral single nucleotide polymorphisms (SNPs) were fixed in the inbred population, far fewer than the theoretical expectation of ∼99%. Under recovery, 36 SNPs across 30 genes involved in alimentary, muscular, nervous, and reproductive systems changed reproducibly across replicates, indicating that strong selection for fitness recovery does exist. Our results indicate that recovery from inbreeding depression via standing genetic variation and mutation is likely to be constrained by the large number of segregating deleterious variants present in natural populations, limiting the capacity for recovery of small populations.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing.

Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnPGGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

Linked selection shapes the landscape of genomic variation in three oak species.

Natural selection shapes genome-wide patterns of diversity within species and divergence between species. However, quantifying the efficacy of selection and elucidating the relative importance of different types of selection in shaping genomic variation remain challenging. We sequenced whole genomes of 101 individuals of three closely related oak species to track the divergence history, and to dissect the impacts of selective sweeps and background selection on patterns of genomic variation. We estimated that the three species diverged around the late Neogene and experienced a bottleneck during the Pleistocene. We detected genomic regions with elevated relative differentiation ('FST -islands'). Population genetic inferences from the site frequency spectrum and ancestral recombination graph indicated that FST -islands were formed by selective sweeps. We also found extensive positive selection; the fixation of adaptive mutations and reduction neutral diversity around substitutions generated a signature of selective sweeps. Prevalent negative selection and background selection have reduced genetic diversity in both genic and intergenic regions, and contributed substantially to the baseline variation in genetic diversity. Our results demonstrate the importance of linked selection in shaping genomic variation, and illustrate how the extent and strength of different selection models vary across the genome.

Hidden role of mutations in the evolutionary process.

Mutations not only alter allele frequencies in a genetic pool but may also determine the fate of an evolutionary process. Here we study which allele fixes in a one-step, one-way model including the wild type and two adaptive mutations. We study the effect of the four basic evolutionary mechanisms-genetic drift, natural selection, mutation, and gene flow-on mutant fixation and its kinetics. Determining which allele is more likely to fix is not simply a question of comparing fitnesses and mutation rates. For instance, if the allele of interest is less fit than the other, then not only must it have a greater mutation rate, but also its mutation rate must exceed a specific threshold for it to prevail. We find exact expressions for such conditions. Our conclusions are based on the mathematical description of two extreme but important regimes, as well as on simulations.

The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels.

Selective breeding for desirable traits in strictly controlled populations has generated an extraordinary diversity in canine morphology and behaviour, but has also led to loss of genetic variation and random entrapment of disease alleles. As a consequence, specific diseases are now prevalent in certain breeds, but whether the recent breeding practice led to an overall increase in genetic load remains unclear. Here we generate whole genome sequencing (WGS) data from 20 dogs per breed from eight breeds and document a ~10% rise in the number of derived alleles per genome at evolutionarily conserved sites in the heavily bottlenecked cavalier King Charles spaniel breed (cKCs) relative to in most breeds studied here. Our finding represents the first clear indication of a relative increase in levels of deleterious genetic variation in a specific breed, arguing that recent breeding practices probably were associated with an accumulation of genetic load in dogs. We then use the WGS data to identify candidate risk alleles for the most common cause for veterinary care in cKCs-the heart disease myxomatous mitral valve disease (MMVD). We verify a potential link to MMVD for candidate variants near the heart specific NEBL gene in a dachshund population and show that two of the NEBL candidate variants have regulatory potential in heart-derived cell lines and are associated with reduced NEBL isoform nebulette expression in papillary muscle (but not in mitral valve, nor in left ventricular wall). Alleles linked to reduced nebulette expression may hence predispose cKCs and other breeds to MMVD via loss of papillary muscle integrity.

Strong within-host selection in a maternally inherited obligate symbiont: Buchnera and aphids

Numerous animal lineages have maternally inherited symbionts that are required for host reproduction and growth. Endosymbionts also pose a risk to their hosts because of the mutational decay of their genomes through genetic drift or to selfish mutations that favor symbiont fitness over host fitness. One model for heritable endosymbiosis is the association of aphids with their obligate bacterial symbiont, Buchnera We experimentally established heteroplasmic pea aphid matrilines containing pairs of closely related Buchnera haplotypes and used deep sequencing of diagnostic markers to measure haplotype frequencies in successive host generations. These frequencies were used to estimate the effective population size of Buchnera within hosts (i.e., the transmission bottleneck size) and the extent of within-host selection. The within-host effective population size was in the range of 10 to 20, indicating a strong potential for genetic drift and fixation of deleterious mutations. Remarkably, closely related haplotypes were subject to strong within-host selection, with selection coefficients as high as 0.5 per aphid generation. In one case, the direction of selection depended on the thermal environment and went in the same direction as between-host selection. In another, a new mutant haplotype had a strong within-host advantage under both environments but had no discernible effect on host-level fitness under laboratory conditions. Thus, within-host selection can be strong, resulting in a rapid fixation of mutations with little impact on host-level fitness. Together, these results show that within-host selection can drive evolution of an obligate symbiont, accelerating sequence evolution.