The influence of deleterious mutations on adaptation in asexual populations

We study the dynamics of adaptation in asexual populations that undergo both beneficial and deleterious mutations. In particular, how the deleterious mutations affect the fixation of beneficial mutations was investigated. Using extensive Monte Carlo simulations, we find that in the "strong-selection weak mutation (SSWM)" regime or in the "clonal interference (CI)" regime, deleterious mutations rarely influence the distribution of "selection coefficients of the fixed mutations (SCFM)"; while in the "multiple mutations" regime, the accumulation of deleterious mutations would lead to a decrease in fitness significantly. We conclude that the effects of deleterious mutations on adaptation depend largely on the supply of beneficial mutations. And interestingly, the lowest adaptation rate occurs for a moderate value of selection coefficient of deleterious mutations.     

扩增子测序法筛选结核分枝杆菌的耐药突变

结核病是严重的公共健康问题之一,而耐药结核病的增加是控制结核病流行的难点之一。快速、准确的诊断是提高结核患者治愈率和降低死亡率的关键因素。本研究建立了基于二代测序技术的扩增子测序方法,对5种一线抗结核药物的17个耐药基因进行检测。在26个临床耐药结核菌株中共鉴定出65个突变,包括33个热点突变,9个稀有突变和23个新突变。对18个新发现的错义突变进行了蛋白质序列保守性和蛋白质局部结构的分析。结果表明,14个新的错义突变在9种分枝杆菌中显示出高度保守性,并且导致了该蛋白质局部结构的改变。根据本研究检测和分析结果,推测这些新发现的突变可能是潜在的耐药突变。在本研究中,构建了扩增子测序的检测方法,可同时检测10株临床结核菌株的17个耐药基因,是一种快速、准确并且全面的检测耐药结核分枝杆菌一线治疗药物耐药突变的方法,该方法不仅能检测热点突变和稀有突变,还能发现一些未报道过的新突变。该检测方法或可用于临床诊断和基础研究。     

Genetic Constraints on Protein Evolution

ABSTRACT

Evolution requires the generation and optimization of new traits (“adaptation”) and involves the selection of mutations that improve cellular function. These mutations were assumed to arise by selection of neutral mutations present at all times in the population. Here we review recent evidence that indicates that deleterious mutations are more frequent in the population than previously recognized and that these mutations play a significant role in protein evolution through continuous positive selection. Positively selected mutations include adaptive mutations, i.e. mutations that directly affect enzymatic function, and compensatory mutations, which suppress the pleiotropic effects of adaptive mutations. Compensatory mutations are by far the most frequent of the two and would allow potentially adaptive but deleterious mutations to persist long enough in the population to be positively selected during episodes of adaptation. Compensatory mutations are, by definition, context-dependent and thus constrain the paths available for evolution. This provides a mechanistic basis for the examples of highly constrained evolutionary landscapes and parallel evolution reported in natural and experimental populations. The present review article describes these recent advances in the field of protein evolution and discusses their implications for understanding the genetic basis of disease and for protein engineering in vitro.     

Genetic constraints on protein evolution

Evolution requires the generation and optimization of new traits ("adaptation") and involves the selection of mutations that improve cellular function. These mutations were assumed to arise by selection of neutral mutations present at all times in the population. Here we review recent evidence that indicates that deleterious mutations are more frequent in the population than previously recognized and that these mutations play a significant role in protein evolution through continuous positive selection. Positively selected mutations include adaptive mutations, i.e. mutations that directly affect enzymatic function, and compensatory mutations, which suppress the pleiotropic effects of adaptive mutations. Compensatory mutations are by far the most frequent of the two and would allow potentially adaptive but deleterious mutations to persist long enough in the population to be positively selected during episodes of adaptation. Compensatory mutations are, by definition, context-dependent and thus constrain the paths available for evolution. This provides a mechanistic basis for the examples of highly constrained evolutionary landscapes and parallel evolution reported in natural and experimental populations. The present review article describes these recent advances in the field of protein evolution and discusses their implications for understanding the genetic basis of disease and for protein engineering in vitro.     

Mutations that eliminate the requirement for the vertex protein in bacteriophage T4 capsid assembly.

The capsid of bacteriophage T4 is composed of two essential structural proteins, gp23, the major constituent of the capsid, and gp24, a less prevalent protein that is located in the pentameric vertices of the capsid. gp24 is required both to stabilize the capsid and to allow it to be further matured. This requirement can be eliminated by bypass-24 (byp24) mutations within g23. We have isolated, cloned and sequenced several new byp24 mutations. These mutations are cold-sensitive in the absence of gp24, and are located in regions of g23 not known to contain any other mutations affecting capsid assembly. The cold-sensitivity of the byp24 mutations can be reduced by further mutations within g23 (trb mutations). Cloning and sequencing of these trb mutations has revealed that they lie in regions of g23 that contain clusters of mutations that cause the production of high levels of petite and giant phage (ptg mutations). Despite the proximity of the trb mutations to the ptg mutations, none of the ptg mutations has a Trb phenotype. The mutation ptE920g, which is also located near one of the ptg clusters, and which produces only petite and wild-type phage, has been shown to confer a Trb but not a Byp24 phenotype. The relevance of these observations to our understanding of capsid assembly is discussed.     

Accessory mutations balance the marginal stability of the HIV-1 protease in drug resistance

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations.     

Role in tumorigenesis of silent mutations in the TP53 gene

Over 10,000 mutations in the TP53 suppressor gene have been recorded in the International Agency for Research on Cancer (IARC) tumor data base. About 4% of these mutations are silent. It is a question whether these mutations play a role in tumor development. In order to approach this question, we asked whether the reported silent mutations are randomly distributed throughout the TP53 gene. The p53 data base was searched exon by exon. From the frequency of codons with no silent mutations, the average number of silent mutations per codon for each exon was calculated using the Poisson distribution. The results indicate the distribution to be non-random. About one-third of all silent mutations occur in "hot-spots" and after subtraction of these hot-spots, the remaining silent mutations are randomly distributed. In addition, the percentage of silent mutations among the total in the silent mutation hot-spots is close to that expected for random mutation. We conclude that most of the silent mutations recorded in tumors play no role in tumor development and that the percentage of silent mutation is an indication of the amount of random mutation during tumorigenesis. Silent mutations occur to a significantly different extent in different tumor types. Tumors of the esophagus and colon have a low frequency of silent mutations, tumors of the prostate have a high frequency.     

Analysis of ENU-induced mutations at the Adh locus in Drosophila melanogaster

N-Ethyl-N-nitrosourea (ENU) was used to induce mutations in the Drosophila melanogaster, alcohol dehydrogenase (Adh) gene. Flies were treated with ENU and mated to homozygous intragenic Adh null mutants; Adh null mutations were selected by exposure of the F1 generation to 1-penten-3-ol. Fourteen Adh null mutations were recovered which included 11 from spermatozoa, 2 from oocytes and 1 from a premeiotic spermatocyte. 2 mutations from spermatozoa and 1 of the mutations from oocytes were multilocus deficiencies which included the Adh locus as determined by complementation tests. The remaining 11 intragenic Adh null mutations were sequenced using the Sanger dideoxy method. One Adh null mutation induced in an oocyte was an AT to TA transversion and the mutation induced in a premeiotic spermatocyte was a GC to AT transition, both of which resulted in a single amino acid substitution. The 11 null mutations induced in spermatozoa were a data set in which both the dose of ENU and the treated germ-cell stage were held constant; therefore, only these 11 mutations were used to calculate the mutation frequency and compare the mutations at the Adh locus with those recovered in other studies. The dose of ENU induced a sex-linked recessive lethal frequency approximately 300 times that of the spontaneous frequency; therefore, these mutations were assumed to have been induced by ENU. 2 of the 11 mutations induced in spermatozoa were multilocus deficiencies and 9 were intragenic mutations. 7 of the 9 intragenic mutations were GC to AT transitions which resulted in 5 single amino acid substitutions, 1 premature translation termination codon, and 1 splice site mutation.(ABSTRACT TRUNCATED AT 400 WORDS)     

Adaptive evolution of asexual populations under Muller's ratchet

We study the population genetics of adaptation in nonequilibrium haploid asexual populations. We find that the accumulation of deleterious mutations, due to the operation of Muller's ratchet, can considerably reduce the rate of fixation of advantageous alleles. Such reduction can be approximated reasonably well by a reduction in the effective population size. In the absence of Muller's ratchet, a beneficial mutation can only become fixed if it creates the best possible genotype; if Muller's ratchet operates, however, mutations initially arising in a nonoptimal genotype can also become fixed in the population, since the loss of the least-loaded class implies that an initially nonoptimal background can become optimal. We show that, while the rate at which adaptive mutations become fixed is reduced, the rate of fixation of deleterious mutations due to the ratchet is not changed by the presence of beneficial mutations as long as the rate of their occurrence is low and the deleterious effects of mutations (s(d)) are higher than the beneficial effects (s(a)). When s(a) > s(d), the advantage of a beneficial mutation can outweigh the deleterious effects of associated mutations. Under these conditions, a beneficial allele can drag to fixation deleterious mutations initially associated with it at a higher rate than in the absence of advantageous alleles. We propose analytical approximations for the rates of accumulation of deleterious and beneficial mutations. Furthermore, when allowing for the possible occurrence of interference between beneficial alleles, we find that the presence of deleterious mutations of either very weak or very strong effect can marginally increase the rate of accumulation of beneficial mutations over that observed in the absence of such deleterious mutations.     

Epistasis between beneficial mutations and the phenotype-to-fitness Map for a ssDNA virus

Epistatic interactions between genes and individual mutations are major determinants of the evolutionary properties of genetic systems and have therefore been well documented, but few quantitative data exist on epistatic interactions between beneficial mutations, presumably because such mutations are so much rarer than deleterious ones. We explored epistasis for beneficial mutations by constructing genotypes with pairs of mutations that had been previously identified as beneficial to the ssDNA bacteriophage ID11 and by measuring the effects of these mutations alone and in combination. We constructed 18 of the 36 possible double mutants for the nine available beneficial mutations. We found that epistatic interactions between beneficial mutations were all antagonistic-the effects of the double mutations were less than the sums of the effects of their component single mutations. We found a number of cases of decompensatory interactions, an extreme form of antagonistic epistasis in which the second mutation is actually deleterious in the presence of the first. In the vast majority of cases, recombination uniting two beneficial mutations into the same genome would not be favored by selection, as the recombinant could not outcompete its constituent single mutations. In an attempt to understand these results, we developed a simple model in which the phenotypic effects of mutations are completely additive and epistatic interactions arise as a result of the form of the phenotype-to-fitness mapping. We found that a model with an intermediate phenotypic optimum and additive phenotypic effects provided a good explanation for our data and the observed patterns of epistatic interactions.     

X-Linked chronic granulomatous disease: mutations in the CYBB gene encoding the gp91-phox component of respiratory-burst oxidase.

Chronic granulomatous disease (CGD) is a hereditary disorder of host defense due to absent or decreased activity of phagocyte NADPH oxidase. The X-linked form of the disease derives from defects in the CYBB gene, which encodes the 91-kD glycoprotein component (termed "gp91-phox") of the oxidase. We have identified the mutations in the CYBB gene responsible for X-linked CGD in 131 consecutive independent kindreds. Screening by SSCP analysis identified mutations in 124 of the kindreds, and sequencing of all exons and intron boundary regions revealed the other seven mutations. We detected 103 different specific mutations; no single mutation appeared in more than seven independent kindreds. The types of mutations included large and small deletions (11%), frameshifts (24%), nonsense mutations (23%), missense mutations (23%), splice-region mutations (17%), and regulatory-region mutations (2%). The distribution of mutations within the CYBB gene exhibited great heterogeneity, with no apparent mutational hot spots. Evaluation of 87 available mothers revealed X-linked carrier status in all but 10. The heterogeneity of mutations and the lack of any predominant genotype indicate that the disease represents many different mutational events, without a founder effect, as is expected for a disorder with a previously lethal phenotype.     

Mitochondrial DNA somatic mutation burden and heteroplasmy are associated with chronological age,smoking, and HIV infection

The gradual accumulation of mitochondrial DNA (mtDNA) mutations is implicated in aging and may contribute to the accelerated aging phenotype seen with tobacco smoking and HIV infection. mtDNA mutations are thought to arise from oxidative damage; however, recent reports implicate polymerase γ errors during mtDNA replication. Investigations of somatic mtDNA mutations have been hampered by technical challenges in measuring low‐frequency mutations. We use primer ID‐based next‐generation sequencing to quantify both somatic and heteroplasmic blood mtDNA point mutations within the D‐loop, in 164 women and girls aged 2–72 years, of whom 35% were smokers and 56% were HIV‐positive. Somatic mutations and the occurrence of heteroplasmic mutations increased with age. While transitions are theorized to result from polymerase γ errors, transversions are believed to arise from DNA oxidative damage. In our study, both transition and transversion mutations were associated with age. However, transition somatic mutations were more prevalent than transversions, and no heteroplasmic transversions were observed. We also measured elevated somatic mutations, but not heteroplasmy, in association with high peak HIV viremia. Conversely, heteroplasmy was higher among smokers, but somatic mutations were not, suggesting that smoking promotes the expansion of preexisting mutations rather than de novo mutations. Taken together, our results are consistent with blood mtDNA mutations increasing with age, inferring a greater contribution of polymerase γ errors in mtDNA mutagenesis. We further suggest that smoking and HIV infection both contribute to the accumulation of mtDNA mutations, though in different ways.     

Role in tumorigenesis of silent mutations in the TP53 gene

Over 10,000 mutations in the TP53 suppressor gene have been recorded in the International Agency for Research on Cancer (IARC) tumor data base. About 4% of these mutations are silent. It is a question whether these mutations play a role in tumor development. In order to approach this question, we asked whether the reported silent mutations are randomly distributed throughout the TP53 gene. The p53 data base was searched exon by exon. From the frequency of codons with no silent mutations, the average number of silent mutations per codon for each exon was calculated using the Poisson distribution. The results indicate the distribution to be non-random. About one-third of all silent mutations occur in “hot-spots” and after subtraction of these hot-spots, the remaining silent mutations are randomly distributed. In addition, the percentage of silent mutations among the total in the silent mutation hot-spots is close to that expected for random mutation. We conclude that most of the silent mutations recorded in tumors play no role in tumor development and that the percentage of silent mutation is an indication of the amount of random mutation during tumorigenesis. Silent mutations occur to a significantly different extent in different tumor types. Tumors of the esophagus and colon have a low frequency of silent mutations, tumors of the prostate have a high frequency.     

Mutational effects on the clonal interference phenomenon

We study the process of fixation of beneficial mutations in an asexual population by means of a theoretical model. Particularly, we wish to investigate how the supply of deleterious and beneficial mutations influences the dynamics of the adaptive process of an evolving population. It is well known that the deleterious mutations drastically affect the fate of beneficial mutations. In addition, an increasing supply of favorable mutations, to compensate the decay of the fitness due to the accumulation of deleterious mutations, produces the clonal interference phenomenon where advantageous mutations in distinct lineages compete to reach fixation. This competition imposes a limit to the speed of adaptation of the population. Intuitively, we would expect that the interplay of the two mechanisms would conspire to ensure fixation of only large-effect beneficial mutations. Our results, however, show that beneficial mutations of small effect have an increased probability of fixation when both beneficial and deleterious mutations rates are increased.     

Microsatellite mutations in litters of the Australian lizard Egernia stokesii

Seven tetranucleotide (AAAG) loci were analysed in a population (including 19 litters) of the Australian lizard Egernia stokesii. In an examination of 76 offspring we observed 13 mutations involving five loci. Two of our loci were highly mutable, with observed mutation rates of 2.7% and 4.2%, representing some of the highest mutations rates reported so far. A high proportion of mutations (46.2%) could not be assigned to changes involving only a single repeat, suggesting that mutations in at least two of the loci follow a multistep model. There was no significant bias of mutations leading to an increase or decrease in allele size; however, all multistep mutations involved a loss of repeats. These results add to increasing evidence casting doubt on microsatellite mutations being primarily single step mutations.     

Generation of recombinant human parainfluenza virus type 1 vaccine candidates by importation of temperature-sensitive and attenuating mutations from heterologous paramyxoviruses

Human parainfluenza virus type 1 (HPIV1) is a significant cause of respiratory tract disease in infants and young children for which a vaccine is needed. In the present study, we sought to attenuate HPIV1 by the importation of one or more known attenuating point mutations from heterologous paramyxoviruses into homologous sites in HPIV1. The introduced mutations were derived from three attenuated paramyxoviruses: (i) HPIV3cp45, a live-attenuated HPIV3 vaccine candidate containing multiple attenuating mutations; (ii) the respiratory syncytial virus cpts530 with an attenuating mutation in the L polymerase protein; and (iii) a murine PIV1 (MPIV1) attenuated by a mutation in the accessory C protein. Recombinant HPIV1 (rHPIV1) mutants bearing a single imported mutation in C, any of three different mutations in L, or a pair of mutations in F exhibited a 100-fold or greater reduction in replication in the upper or lower respiratory tract of hamsters. Both temperature-sensitive (ts) (mutations in the L and F proteins) and non-ts (the mutation in the C protein) attenuating mutations were identified. rHPIV1 mutants containing a combination of mutations in L were generated that were more attenuated than viruses bearing the individual mutations, showing that the systematic accretion of mutations can yield progressive increases in attenuation. Hamsters immunized with rHPIV1 mutants bearing one or two mutations developed neutralizing antibodies and were resistant to challenge with wild-type HPIV1. Thus, importation of attenuating mutations from heterologous viruses is an effective means for rapidly identifying mutations that attenuate HPIV1 and for generating live-attenuated HPIV1 vaccine candidates.     

Dominant mutations of rifampicin resistance in Escherichia coli K-12

Significant portion (up to 20%) of dominant mutations (rifd mutations) was observed among spontaneous mutations of rifampicin resistance picked up in cells of haploid Escherichia coli strain. These mutations are similar to rifd mutations obtained earlier when selecting them in rif-s/rif-s merodiploids. On the basis of analysis of nucleotide substitutions taking place in formation of spontaneous and induced mutations, it is established that rifd mutations are caused by single nucleotide substitution. The majority of rifd mutations are localized in a small region of the central part of RNA polymerase beta-subunit gene covering less than 200 base pairs. A rifd mutant has been described which occurred as a result of micro-deletion in one of the "hot" spots of the central region of beta-subunit gene.