Substitution Hotspots Research Articles

Structural topological analysis of spike proteins of SARS-CoV-2 variants of concern highlight distinctive amino acid substitution patterns

Since the onset of pandemic in 2019, SARS-CoV-2 has diverged into numerous variants driven by antigenic and infectivity-oriented selection. Some variants have accumulated fitness-enhancing mutations, evaded immunity and spread despite global vaccination campaigns. The spike (S) glycoprotein of SARS-CoV-2 demonstrated the greatest immunogenicity and amino acid substitution diversity owing to its importance in the interaction with human angiotensin receptor 2 (hACE2). The S protein consistently emerges as an amino acid substitution (AAS) hotspot in all six lineages, however, in Omicron this enrichment is significantly higher. This study attempts to design and validate a method of mapping S-protein substitution profile across variants to identify the conserved and AAS regions. A substitution matrix was created based on publicly available databases, and the substitution localization was illustrated on a cryo-electron microscopy generated S-protein model. Our analyses indicated that the diversity of N-terminal (NTD) and receptor-binding (RBD) domains exceeded that of any other regions but still contained extended low substitution density regions particularly considering significantly broader substitution profiles of Omicron BA.2 and BA.4/5. Finally, the substitution matrix was compared to a random sample alignment of variant sequences, revealing discrepancies. Therefore, it was suggested to improve matrix accuracy by processing a large number of S-protein sequences using an automated algorithm. Several critical immunogenic and receptor-interacting residues were identified in the conserved regions within NTD and RBD. In conclusion, the structural and topological analysis of S proteins of SARS-CoV-2 variants highlight distinctive amino acid substitution patterns which may be foundational in predicting future variants.

Risk Assessment of Hepatitis C Virus Resistance to Ns3/4a Protease Inhibitors (Pis)

Hepatitis C virus (HCV) persistently infects more than 175-185 million people globally. Drug resistance to anti-HCV DAAs therapeutics particularly NS3/4A has been a major concern for efficacy of HCV treatment due to generation of HCV quasispecies. Risk determination of the resistance of HCV to anti-NS3/4A inhibitors was accomplished using ViPR algorithm that computes and analyzes resistance-associated amino acid substitutions in HCV NS3 protein target site and reveal altered response to NS3/4A PIs. Spot comparison and data visualization was carried out using ncbi.nlm.nih.gov protein sequence graphics and R-package for Ubuntu version 12.04. Resistance to multiple NS3/4A PIs (telaprevir, simeprevir, faldaprevir and asunaprevir) was exhibited by HCV genotype 1 variant D1194A/E/T/G-NS3 (a substitution hot spot) while the risk of resistance to telaprevir (V1062 L NS3 variant) and simeprevir (S1148R variant) was demonstrated by HCV genotype 2. Genotype 3 & 4 revealed resistance against only telaprevir. Another almost pan-genotype (genotype 2,3,4,5 and unclassified) substitution hot spot was identified as anti-telaprevir variant V1062 L. Single substitution conferring multiple anti-NS3/4A PIs resistance was observed in D1194A NS3 variant (anti-simeprevir, anti-faldaprevir, anti-asunaprevir in genotype 1), Q1106K NS3 variant (anti-simeprevir, anti-faldaprevir in genotype 5 and 6) and T1080S NS3 variant (anti-faldaprevir, anti-telaprevir in genotype 5 and unclassified NS3 sequence). In addition to that, variant S1148G exhibited increased sensitivity to simeprevir in case of genotype 1. Because of a higher degree of HCV genomic variability, NS3 variants with decreased susceptibility to NS3/4A PIs exists, and therefore the resistance profile along with potency, adverse effect of PIs is noteworthy and should be well-considered while developing different protease inhibitors for effective treatment of HCV patients.

BioAider: An efficient tool for viral genome analysis and its application in tracing SARS-CoV-2 transmission

The novel human coronavirus (SARS-CoV-2) causes the coronavirus disease 2019 (COVID-19) pandemic worldwide. Control of COVID-19 pandemic is vital for public health and is the prerequisite to maintain social stability. However, the origin and transmission route of SARS-CoV-2 is unclear, bringing huge difficult to virus control. Monitoring viral variation and screening functional mutation sites are crucial for prevention and control of infectious diseases. In this study, we developed a user-friendly software, named BioAider, for quick sequence annotation and mutation analysis on large-scale genome-sequencing data. Herein, we detected 14 substitution hotspots within 3,240 SARS-CoV-2 genome sequences, including 3 groups of potentially linked substitution. NSP13-Y541C was crucial substitution which might affect the unwinding activity of the viral helicase. In particular, we discovered a SR-rich region of SARS-CoV-2 distinct from SARS-CoV, indicating more complex replication mechanism and unique N-M interaction of SARS-CoV-2. Interestingly, the quantity of SSRX repeat fragments in SARS-CoV-2 provided further evidence of its animal origin. Overall, we developed an efficient tool for rapid identification of viral genome mutations which could facilitate viral genomic studies. Using this tool, we have found critical clues for the transmission route of SARS-CoV-2 which would provide theoretical support for the epidemic control of pathogenic coronaviruses.

Identifying recurrent mutations in cancer reveals widespread lineage diversity and mutational specificity.

Mutational hotspots indicate selective pressure across a population of tumor samples, but their prevalence within and across cancer types is incompletely characterized. An approach to detect significantly mutated residues, rather than methods that identify recurrently mutated genes, may uncover new biologically and therapeutically relevant driver mutations. Here we developed a statistical algorithm to identify recurrently mutated residues in tumour samples. We applied the algorithm to 11,119 human tumors, spanning 41 cancer types, and identified 470 hotspot somatic substitutions in 275 genes. We find that half of all human tumors possess one or more mutational hotspots with widespread lineage-, position-, and mutant allele-specific differences, many of which are likely functional. In total, 243 hotspots were novel and appeared to affect a broad spectrum of molecular function, including hotspots at paralogous residues of Ras-related small GTPases RAC1 and RRAS2. Redefining hotspots at mutant amino acid resolution will help elucidate the allele-specific differences in their function and could have important therapeutic implications.

Dynamic evolutionary pattern of α 2-macroglobulin in a model organism, the zebrafish ( Danio rerio)

Studies of alpha-2-macroglobulin (α 2M), a universal protease inhibitor, have indicated that it plays a unique and critical role in the innate immune system of vertebrate and invertebrate animals. The distinctive mechanism of pathogen inhibition – through physical entrapment of the pathogen-derived protease – makes α 2M an ideal candidate for molecular evolutionary analysis. Furthermore, recent studies revealed that the Osteichthyes are characterized by levels of α 2M diversification that exceed those recorded in other animal groups. Our study of Danio rerio (zebrafish) indicated that (1) two distinct lineages of α 2M and α 2M-like isoforms exist and (2) at least some codons in the functional domains of α 2M have been subjected to positive Darwinian selection. The findings of several hot-spots for nonsynonymous substitutions in the two functional domains such as bait region and receptor binding domain, suggest that host-immune selection have played a dominant role in these two genomic regions of α 2M. The presence of two, non-monophyletic α 2M lineages in zebrafish provides compelling evidence of an ancient gene duplication event. The accelerated rate of nucleotide substitution in the functional domains of α 2M is consistent with similar observations of other immune system components.

The Impact of Recombination on Nucleotide Substitutions in the Human Genome

Unraveling the evolutionary forces responsible for variations of neutral substitution patterns among taxa or along genomes is a major issue for detecting selection within sequences. Mammalian genomes show large-scale regional variations of GC-content (the isochores), but the substitution processes at the origin of this structure are poorly understood. We analyzed the pattern of neutral substitutions in 1 Gb of primate non-coding regions. We show that the GC-content toward which sequences are evolving is strongly negatively correlated to the distance to telomeres and positively correlated to the rate of crossovers (R2 = 47%). This demonstrates that recombination has a major impact on substitution patterns in human, driving the evolution of GC-content. The evolution of GC-content correlates much more strongly with male than with female crossover rate, which rules out selectionist models for the evolution of isochores. This effect of recombination is most probably a consequence of the neutral process of biased gene conversion (BGC) occurring within recombination hotspots. We show that the predictions of this model fit very well with the observed substitution patterns in the human genome. This model notably explains the positive correlation between substitution rate and recombination rate. Theoretical calculations indicate that variations in population size or density in recombination hotspots can have a very strong impact on the evolution of base composition. Furthermore, recombination hotspots can create strong substitution hotspots. This molecular drive affects both coding and non-coding regions. We therefore conclude that along with mutation, selection and drift, BGC is one of the major factors driving genome evolution. Our results also shed light on variations in the rate of crossover relative to non-crossover events, along chromosomes and according to sex, and also on the conservation of hotspot density between human and chimp.

Detecting molecular adaptation at individual codons in the pattern recognition protein, lipopolysaccharide- and β-1,3-glucan-binding protein of decapods

Pattern recognition proteins play an important role in the innate immune response of invertebrates. Herein we report the evolutionary relationships among Gram-negative bacteria binding proteins (GNBPs) that were previously identified and characterized from a wide array of invertebrates. Our results, together with those obtained in previous studies, indicate that decapod lipopolysaccharide- and β-1,3-glucan binding protein (LGBP/BGBP) has retained the crucial components for glucanase activity, and shares a common ancestor with GNBPs, as well as with the glucanase proteins of a wide range of invertebrates, rather than with GNBPs of some arthropods. However, experimental evidence of earlier studies suggested a lack of glucanase activity by these proteins, thus implying that during evolutionary time these proteins might have lost their glucan binding protein, but retained their glucan binding activity. The present results have also revealed that although a vast majority of the decapod LGBP/BGBP codons are constrained to purifying selection, certain codons are shown to have a higher rate of nonsynonymous substitutions per nonsynonymous site ( d N) than synonymous substitutions per synonymous site ( d S), indicating these codons have evolved adaptively ( d N/ d S > 1). Although purifying selection ( d N/ d S < 1) appears to be the major driving force in the evolution of a vast majority of LGBP/BGBP codons in decapods, the findings of several hotspots for nonsynonymous substitutions in this protein indicate host immune selection might play an important role in maintaining diversity among these ecologically diversified decapod species.

Metagenomic DNA fragments that affect Escherichia coli mutational pathways

A multicopy cloning approach was used to search for metagenomic DNA fragments that affect Escherichia coli mutational pathways. Soil metagenomic expression libraries were constructed with DNA samples prepared directly from soil samples collected from the UCLA Botanical Garden. Using frameshift mutator screening, we obtained a total of 26 unique metagenomic fragments that stimulate frameshift rates in an E. coli wild-type host. Mutational enhancer strains such as an ndk-deficient strain and a temperature sensitive mutS strain (mutS60) were used to further verify the mutator phenotype. We found that the presence of multiple copies of certain types of metagenomic DNA sequence repeats cause general genome instability in the wild-type E. coli host and the effect can be suppressed by overproducing a DNA mismatch component MutL. In addition, we identified nine metagenomic mutator genes (designated as smu genes) that encode proteins that have not been linked to mutator phenotypes prior to this study including a putative RNA methyltransferase Smu10A. The strain overproducing Smu10A displays one prominent base substitution hotspot in the rpoB gene, which coincides with the base substitution hotspot we have observed in cells that are partially deficient in the proofreading function carried out by the DNA polymerase III epsilon subunit. Based on the structural conservation of DNA replication/recombination/repair machineries among microorganisms, this approach would allow us to both identify new mutational pathways in E. coli and to find genes involved in DNA replication, recombination or DNA repair from vast unculturable microbes.

Developing a genetic system in Deinococcus radiodurans for analyzing mutations.

We have applied a genetic system for analyzing mutations in Escherichia coli to Deinococcus radiodurans, an extremeophile with an astonishingly high resistance to UV- and ionizing-radiation-induced mutagenesis. Taking advantage of the conservation of the beta-subunit of RNA polymerase among most prokaryotes, we derived again in D. radiodurans the rpoB/Rif(r) system that we developed in E. coli to monitor base substitutions, defining 33 base change substitutions at 22 different base pairs. We sequenced >250 mutations leading to Rif(r) in D. radiodurans derived spontaneously in wild-type and uvrD (mismatch-repair-deficient) backgrounds and after treatment with N-methyl-N'-nitro-N-nitrosoguanidine (NTG) and 5-azacytidine (5AZ). The specificities of NTG and 5AZ in D. radiodurans are the same as those found for E. coli and other organisms. There are prominent base substitution hotspots in rpoB in both D. radiodurans and E. coli. In several cases these are at different points in each organism, even though the DNA sequences surrounding the hotspots and their corresponding sites are very similar in both D. radiodurans and E. coli. In one case the hotspots occur at the same site in both organisms.

Developing a Genetic System in Deinococcus radiodurans for Analyzing Mutations

Abstract We have applied a genetic system for analyzing mutations in Escherichia coli to Deinococcus radiodurans, an extremeophile with an astonishingly high resistance to UV- and ionizing-radiation-induced mutagenesis. Taking advantage of the conservation of the β-subunit of RNA polymerase among most prokaryotes, we derived again in D. radiodurans the rpoB/Rif r system that we developed in E. coli to monitor base substitutions, defining 33 base change substitutions at 22 different base pairs. We sequenced &amp;gt;250 mutations leading to Rif r in D. radiodurans derived spontaneously in wild-type and uvrD (mismatch-repair-deficient) backgrounds and after treatment with N-methyl-N′-nitro-N-nitrosoguanidine (NTG) and 5-azacytidine (5AZ). The specificities of NTG and 5AZ in D. radiodurans are the same as those found for E. coli and other organisms. There are prominent base substitution hotspots in rpoB in both D. radiodurans and E. coli. In several cases these are at different points in each organism, even though the DNA sequences surrounding the hotspots and their corresponding sites are very similar in both D. radiodurans and E. coli. In one case the hotspots occur at the same site in both organisms.

Differential Influence of Nucleoside Analog-resistance Mutations K65R and L74V on the Overall Mutation Rate and Error Specificity of Human Immunodeficiency Virus Type 1 Reverse Transcriptase

Differential Influence of Nucleoside Analog-resistance Mutations K65R and L74V on the Overall Mutation Rate and Error Specificity of Human Immunodeficiency Virus Type 1 Reverse Transcriptase

Evaluation of the database on mutant frequencies and DNA sequence alterations of vermilion mutations induced in germ cells of Drosophila shows the importance of a neutral mutation detection system

The vermilion gene in Drosophila has extensively been used for the molecular analysis of mutations induced by chemicals in germ cells in vivo. The gene is located on the X-chromosome and is a useful target for the study of mutagenesis since all types of mutations are generated. We have critically evaluated this system with respect to sensitivity for mutation induction and selectivity for different types of mutations, using a database of more than 600 vermilion mutants induced in postmeiotic male germ cells by 18 mutagens. From most of these mutants the mutation has been analysed. These data showed 336 base substitutions, 96 intra-locus DNA rearrangements and 78 multi-locus deletions (MLD). Mutants containing a MLD were either heterozygous sterile or homozygous and hemizygous lethal. The distribution of both basepair (bp) changes and intra-locus rearrangements over the coding region of the vermilion gene was uniform with no preferences concerning 5′ or 3′ regions, certain exons, splice sites, specific amino acid changes or nonsense mutations. Possible hotspots for base substitutions seem to be related to the type of DNA damage rather than to the vermilion system. Gene mutations other than bp changes were examined on sequence characteristics flanking the deletion breakpoints. Induction frequencies of vermilion mosaic mutants were, in general, higher than those of vermilion complete mutants, suggesting that persistent lesions are the main contributors to the molecular spectra. Comparison of induction frequencies of vermilion mutants and sex-linked recessive lethal (SLRL) mutants for the 18 mutagens showed that the sensitivity of the vermilion gene against a mutagenic insult is representative for genes located on the X-chromosome. The effect of nucleotide excision repair (NER) on the formation of SLRL mutants correlated with an increase of transversions in the vermilion spectra under NER deficient conditions. Furthermore, the clastogenic potency of the mutagens, i.e., the efficiency to induce chromosomal-losses vs. SLRL forward mutations, shows a positive correlation with the percentage of DNA deletions in the molecular spectra of vermilion mutants.

Lessons learned from xeroderma pigmentosum.

Lessons learned from xeroderma pigmentosum.

The Influence of Single Base Changes on UV Mutational Activity at Two Translocated Hotspots

Mutation hotspots have been a staple of mutation spectra since the introduction of fine structure mutation mapping almost 40 years ago. It has been well established that sequence context is an important determinant of mutational activity at mutagen induced hotspots and coldspots. However, our understanding of the sequence effectors of base substitution hotspots is quite limited. This is because manipulation of the sequence about a hotspot site in a marker gene is restricted by the need to maintain a functional marker. In this work, we describe a generalizable system for studying sequence context effects on mutagenesis. We have prepared a variant of the supF tRNA gene (a marker used by us in previous studies) in which an eight-base palindrome, the site of two UV hotspots in the interior of the gene, was copied into the acceptor stem and pre-tRNA region. The variant tRNA was active. The UV mutation spectrum of this variant showed that the new copy of the palindrome generated two hotspots which were as intense as the original sites in the interior of the gene. Variant genes were constructed with all possible bases at the first position in the palindrome in the pre-tRNA sequence, which does not affect tRNA function. The mutation analysis showed that activity at one of the hotspots could be reduced or enhanced by the changes, while activity at the other site was not significantly affected. The base changes did not influence the frequency of cyclobutane dimer or (6-4) photoproduct formation at the two hotspot sites. Thus, the changes in mutational activity were due to the influence of sequence context on the efficiency of mutation formation at the sites of UV lesions.

The major, N2-Gua adduct of the (+)-anti-benzo[a]pyrene diol epoxide can be unstable in double-stranded DNA.

The mechanisms of mutagenesis by the (+)-anti diol epoxide of benzo[a]pyrene [(+)-anti-B[a]PDE] was investigated in supF of the Escherichia coli plasmid pUB3 [Rodriguez & Loechler (1993) Biochemistry 32, 1759]. pUB3 was reacted with (+)-anti-B[a]PDE, then either (1) transformed immediately into E. coli or (2) heated at 80 degrees C for 10 min prior to transformation--the latter to probe mechanism. Qualitatively, heating did not have a statistically significant effect on the mutagenic pattern, except at the major base substitution hot spot, G115, in supF; principally, G115-->T mutations were obtained prior to heating, while after heating, G115-->A and G115-->C mutations became more prevalent. Quantitatively, heating caused an approximately 2-fold decrease in mutation frequency. Heating released a small fraction of adducts (approximately 5%), and the chemistry and implications of this reaction are investigated herein. Although the major adduct of (+)-anti-B[a]PDE (formed at N2-Gua) is generally regarded to be heat stable, it can be quite unstable in double-stranded (but not single-stranded) DNA at low [Mg2+]. Heating releases the corresponding tetraols from (+)-anti-B[a]P-N2-Gua in approximately the same ratio as for simple hydrolysis of (+)-anti-B[a]PDE itself. This and other results suggest that guanine remains in DNA when (+)-anti-B[a]P-N2-Gua adducts are hydrolyzed. [No evidence for any other chemical change in (+)-anti-B[a]PDE adducts was found.] While no general acid/base or nucleophilic catalysis was observed, adduct hydrolysis was specific acid catalyzed down to pH approximately 5.6, where the pH-rate profile showed a break to a slope of approximately 0.0. This break probably indicates the pKa of (+)-anti-B[a]P-N2-Gua protonated at the N2-position, which is higher than expected.(ABSTRACT TRUNCATED AT 250 WORDS)

DNA sequence analysis of spontaneous hprt mutations arising in vivo in cynomolgus monkey T-lymphocytes.

To study the mechanisms of mutagenesis in vivo, we analyzed mutations at the hypoxanthine phosphoribosyl transferase (hprt) locus using cDNA from cynomolgus monkey T-lymphocytes. In the present study, the spectrum of spontaneous hprt mutations arising in vivo in wild-caught cynomolgus monkey peripheral T-lymphocytes is described. Cells were isolated from peripheral blood, and mutant clones were selected in 6-thioguanine, propagated, and stored frozen. cDNA was copied from hprt mRNA from a lysate of 7,000 to 20,000 cells. A 780-base-pairs (bp) region including the coding region was amplified by polymerase chain reaction and directly sequenced. We sequenced 40 spontaneous mutants from 11 monkeys. Of these 40 clones, 23 (57%) had base-pair substitutions, 11 (28%) had small (< 20 bp) deletions and/or insertions, and 6 (15%) had large (> 20 bp) deletions and/or insertions. Of the 23 base substitutions, 13 were transitions (11 G:C-->A:T, 1 A:T-->G:C, and 1 tandem TT-->CC) and 10 were transversions (3 G:C-->T:A, 3 G:C-->C:G, 2 A:T-->T:A, 2 A:T-->C:G). Bases 209 and 617 were apparent substitution hotspots, which have also been observed as hotspots in human hprt. In 2 clones with large insertions, the inserted bases were of intronic origin. One of these lost 272 bp from exons 2-3 and contained a 93-bp insertion from the middle of intron 3. Two clones with small deletions and 5 clones with large deletions or insertions (7/40 or 17.5%) could be splice mutants.

Spontaneous mutation in Escherichia coli containing the dnaE911 DNA polymerase antimutator allele.

We have previously isolated mutants of Escherichia coli that replicate their DNA with increased fidelity. These mutants have a mutation in the dnaE gene, encoding the alpha subunit of DNA polymerase III. They were isolated in a mismatch-repair-defective mutL background, in which mutations can be considered to represent uncorrected DNA replication errors. In the present study we analyze the effect of one of these alleles, dnaE911, on spontaneous mutagenesis in a mismatch-repair-proficient background. In this background, spontaneous mutations may be the sum of uncorrected replication errors and mutations resulting from other pathways. Hence, the effect of the dnaE allele may provide insights into the contribution of uncorrected DNA replication errors to spontaneous mutation. The data show that dnaE911 decreases the level of Rifr, lacI and galK mutations in this background by 1.5-2-fold. DNA sequencing of 748 forward mutants in the lacI gene reveals that this effect has a clear specificity. Transversions are decreased by approximately 3-fold, whereas transitions, frameshifts, deletions and duplications remain essentially unchanged. Among the transversions, A.T-->T.A are affected most strongly (approximately 6-fold). In addition to this effect on transversions within the lacI gene, one previously recognized A.T-->G.C base-pair substitution hotspot in the lac operator is also reduced (approximately 5-fold). The data are discussed in the light of the role of DNA replication errors in spontaneous mutation, as well as other possible explanations for the observed antimutator effects.

Hydrogen peroxide induces G:C to T:A and G:C to C:G transversions in the supF gene of Escherichia coli.

A vector plasmid, pZ189, carrying an Escherichia coli supF gene as a target for mutations, was treated with a combination of hydrogen peroxide and Fe3+/EDTA complex and propagated in E. coli host cells that had been induced for SOS functions by ultraviolet irradiation. The mutations frequency increased by up to 30-fold over spontaneous background levels with increasing concentrations of hydrogen peroxide. The increase in mutation frequency correlated with an increase in the formation of 8-hydroxydeoxyguanosine in the pZ189 DNA. Sequence analysis of 82 independent supF mutant plasmids revealed that 70 mutants contained base substitutions, with 63 of the 70 involving a G:C base pair, and with G:C-->C:G (28 cases) and G:C-->T:A (26 cases) transversions predominating. Investigation of the influence of the local DNA sequence on the transversions revealed that the guanine at the center of the triplet 5'-PuGA-3' was five times more likely to mutate after treatment with hydrogen peroxide than that at the center of 5'PyGN3'. G:C-->T:A transversions presumably resulted from mispairing of an altered G (probably 8-hydroxydeoxyguanosine) with deoxyadenosine. The origin of the G:C-->C:G transversions may be an as yet unidentified lesion generated by hydrogen peroxide. Mutagenic hotspots for base substitutions were found at positions 133, 160 and 168. Mutation spectra and the positions of mutagenic hotspots, when compared with a previously determined spontaneous mutagenesis spectrum, also provide information on the mechanism of spontaneous mutagenesis.

AP sites are not significantly involved in mutagenesis by the (+)-anti diol epoxide of benzo[a]pyrene: the complexity of its mutagenic specificity is likely to arise from adduct conformational polymorphism.

In previous work, mutations induced by the (+)-anti diol epoxide of benzo[a]pyrene [(+)-anti-B[a]PDE] were scored in the supF gene of the Escherichia coli plasmid pUB3 [Rodriguez & Loechler (1993) Biochemistry 32, 1759]. pUB3 was reacted with (+)-anti-B[a]PDE and then either (1) transformed immediately into E. coli or (2) heated at 80 degrees C for 10 min prior to transformation. Heating only released a small fraction of adducts (approximately 5%) and did not significantly affect the mutagenic pattern at most sites in supF. However, at the major base substitution hotspot, G115, principally G-->T mutations (87%) were obtained prior to heating, while after heating, G-->T mutations decreased (45%) and G-->A (21%) and G-->C (33%) mutations became more prevalent. One model for this result is that prior to heating a heat-labile adduct at G115 causes one pattern of mutagenesis, but after heating the labile adduct is hydrolyzed to an apurinic site (AP site), which causes a second mutational pattern. To test this, a role for AP sites generated from labile adducts by heating at 80 degrees C for 10 min is investigated. It is shown that when plasmid pUB3 contains 22 (+)-anti-B[a]PDE adducts, 0.6% (or fewer) are converted to AP sites as determined in an assay based upon the action of an AP-endonuclease. In a separate line of investigation not involving (+)-anti-B[a]PDE adducts, mutation frequency (MF) per AP site is estimated. (In these experiments, AP sites were introduced into pUB3 by the classic procedure of heating at 70 degrees C/pH 5.0 to hydrolyze purines.(ABSTRACT TRUNCATED AT 250 WORDS)

Error-prone polymerization by HIV-1 reverse transcriptase. Contribution of template-primer misalignment, miscoding, and termination probability to mutational hot spots

We have observed previously that DNA template-directed polymerization by the type 1 human immunodeficiency virus reverse transcriptase is error-prone for single-nucleotide substitution, addition and deletion errors at homopolymeric sequences. We have also noted strong termination of processive synthesis at these positions (Bebenek, K., Abbotts, J., Roberts, J. D., Wilson, S. H., and Kunkel, T. A. (1989) J. Biol. Chem. 264, 16948-16956). Here we have tested three models to explain errors at these hot spots: template-primer misalignment for deletion errors, and dislocation and direct miscoding for substitution errors. The approach involves introducing single-nucleotide changes within or flanking the homopolymeric hot spots and examining the effects that these changes have on human immunodeficiency virus type 1 (HIV-1) reverse transcriptase error rate, error specificity, and termination probability. The results obtained suggest that single-nucleotide deletion errors in homopolymeric runs result from template-primer misalignment and that both direct miscoding and template-primer dislocation contribute to the base substitution hot spots. The data also suggest that base substitution errors at one position can be templated by the preceding nucleotide or either of the next two nucleotides. Frameshift error rates at homopolymeric sites were affected by changes in the sequences flanking the runs, including single-nucleotide differences in the single-stranded template strand and in the double-stranded primer region as many as six nucleotides distant from the hot spot. Both increases and decreases in frameshift fidelity were observed, and most of these correlated with concomitant increases or decreases in the probability that HIV-1 reverse transcriptase terminated processive synthesis within the run. These data provide further support for a relationship between the frameshift fidelity and the processivity of DNA-dependent DNA synthesis by HIV-1 reverse transcriptase.