Nucleotide Sequence Context Research Articles

Tumor mutation burden in colorectal cancers with POLE exonuclease and non-exonuclease domain variants.

224 Background: Driver mutations in the exonuclease domain (ExoD) of DNA polymerase epsilon ( POLE) have been associated with high tumor mutation burden (TMB). High TMB is an important marker for immunotherapy, and thus understanding the mechanisms that lead to high TMB is an area of active investigation. Many tumors with POLE ExoD driver mutations have additional POLE variants of uncertain significance (VUS). Some of these VUS lie within the ExoD, while most lie outside the ExoD. This study investigates the role of such POLE VUS on TMB and POLE functionality of ExoD-mutated tumors. Methods: Retrospective colorectal cancer (CRC) genomic profiles from Caris Life Sciences (n=1,870) were separated into groups based on the presence of POLE ExoD mutations and the TMB status. Groups were defined as: ‘TMB-high (TMB-H) POLE ExoD driver without additional POLE VUS’, ‘TMB-H POLE ExoD driver plus POLE VUS’, and ‘TMB-low (TMB-L) POLE variant(s) without ExoD driver’. An additional group with ‘TMB-H POLE variant(s) without ExoD driver’ was identified as potentially novel drivers for future study. The microsatellite instability status (MSI or MSS) was also considered. The level of TMB, the presence of polymerase proofreading-associated mutational signatures and POLE stability/functionality analyzed using AlphaFold2 models and Rosetta Molecular Modeling Suite, were evaluated. Results: POLE variants were identified in 4.9% of the CRCs studied (92/1870). In all, 52.2% (48/92) of those tumors were TMB-H, and 38.0% (35/92) of them contained POLE ExoD drivers. The highest median TMB (mTMB) was found in the ‘TMB-H POLE ExoD driver plus POLE VUS’ group, and this was significantly different from the ‘TMB-H POLE ExoD driver without additional POLE VUS’ and the ‘TMB-L POLE variant without ExoD driver’ group, even when MSI tumors were excluded (p<0.001, Mann-Whitney). These results were validated in TCGA CRC tumors and in other tumor types. Within the ‘TMB-H POLE ExoD driver plus POLE VUS’ group, the mTMB increased with the number of accumulated POLE VUS (p<0.001, Mann Whitney). The nucleotide sequence context of 77.8% of the POLE VUS matched mutational spectra associated with POLE ExoD defects, suggesting that most of the POLE VUS are secondary to the ExoD driver mutation. Finally, the majority of POLE VUS within ‘TMB-H POLE ExoD driver plus POLE VUS’ group have a predicted effect on POLE stability and/or function. Conclusions: Our results suggest that the presence of secondary POLE non-driver variants increase TMB in POLE ExoD-mutated tumors. Although still preliminary, these data could impact prognosis and response to immune checkpoint inhibitors in patients. [Table: see text]

Ultra-precise single-molecule nanopore tool to watch enzymes at work.

We developed a nanopore tool to measure how enzymes move along nucleic acid chains. In this technique, which evolved from our effort in developing nanopore sequencing, the ion current through an engineered protein nanopore is used to measure the motion of the DNA or RNA. With its unprecedented spatiotemporal resolution this new single-molecule tool resolves sub-millisecond and sub-nucleotide-long enzyme steps while providing the exact nucleotide sequence context within the enzyme. In addition, the tool has high throughput and is relatively simple to set up. We find surprisingly strong sequence-dependence in the kinetics of many types of helicases such as in Hel308, PcrA, RecQ, or the SARS-CoV-2 helicase nsp13. Furthermore, we have recorded data with DNA- and RNA-polymerases and reverse transcriptases. With each enzyme we discover behavior that had not been observable with any other single-molecule technique.

Identification of 2,6-Diaminopurine As a Candidate for Hemophilia a Ribosomal Readthrough Therapy

Nonsense mutations introduce premature termination codons (PTCs) and are responsible for 11% hereditary genetic disorders including hemophilias. Drug-induced ribosomal readthrough over PTCs is an emerging strategy to treat the genetic disorders caused by nonsense mutations. Ribosomal readthrough therapy aims to allow translation beyond the PTC, creating a full-length protein and possible phenotypic rescue. Ataluren (PTC124), which induces readthrough by binding to the ribosomal 60S subunit, was recently granted conditional approval in the European Union for the treatment of Duchenne muscular dystrophy caused by nonsense mutations. Most nonsense mutations in FVIII and FIX lead to severe diseases. In theory, hemophilias are ideally suited for drug-induced readthrough therapy, because even a small increase in FVIII or FIX activity to over 1% could significantly improve bleeding symptoms in patients with nonsense mutations. Thus, ribosomal readthrough therapy has the potential to provide a small-molecule, personalized treatments based on patient's mutation profile. We previously investigated different readthrough drugs on different FVIII nonsense mutations and found that geneticin, but not gentamicin and Ataluren, is an effective readthrough agent for a number of such mutations. Recently, 2,6-diaminopurine (DAP), a natural compound found in extracts of the edible mushroom Lepista inversa, has shown high-efficiency correction of UGA nonsense mutations and low toxicity. DAP is thought to act as an inhibitor of a tRNA-specific 2'-O-methyltransferase and promotes the interaction of the near-cognate tRNATrp with the UGA stop codon. We introduced 22 UGA mutations, including the most common recurring mutations from hemophilia A (HA) patients into the B domain deleted FVIII-Gassia luciferase fusion (BDD-Gluc) protein. Cells stably expressing these mutant fusion proteins were treated with DAP and geneticin. DAP induced higher secreted luciferase activity, FVIII antigen and cofactor activity levels in cell culture media than geneticin treatments. A majority of DAP-treated mutants produced FVIII co-factor activities and antigens that are over 2% of wild-type level, with several mutants reaching 10-15% levels. Using an immunoprecipitation-immunoblotting assay, we detected full-length BDD-Gluc fusion proteins in culture media and cell lysates of both DAP- and geneticin-treated cells expressing nonsense mutants. Two major factors affect the outcome of readthrough therapies are the readthrough efficiency of drugs and the nature of replacement of amino acid (AA) residues. Studies had shown that a cytosine at +4 position is a favorable nucleotide sequence context. Readthrough efficiency of F8 mRNA is influenced, but not determined by nucleotide sequence contexts of the PTC. It was reported that only tryptophan was replaced with DAP-induced readthrough while three potential AAs can be replaced with geneticin treatment. We systematically substituted PTCs with potential replacement AAs and found that activities of readthrough products are determined by the compatibility of replacement AA residues with FVIII secretion and activity. Our results suggest that DAP, with its high readthrough efficiency and low toxicity, is a promising drug candidate for HA patients carrying susceptible UGA mutations. This cell-based study identifies candidate patient mutations and provides a rationale for a future clinical study of DAP or one of its derivatives as a therapeutic for these patients.

The effect of the nucleotides immediately upstream of the AUG start codon on the efficiency of translation initiation in sperm cells.

It is widely known that an optimal nucleotide sequence context immediately upstream of the AUG start codon greatly improves the efficiency of translation initiation of mRNA in mammalian and plant somatic cells, which in turn increases protein levels. However, it is still unclear whether a similar regulatory mechanism is also present in highly differentiated cells. Here, we surveyed this issue in Arabidopsis thaliana sperm cells and found that the sequence context-mediated regulation of translation initiation in sperm cells is generally similar to that in somatic cells. A simple motif of four adenine nucleotides at positions -1 to -4 greatly improved the efficiency of translation initiation, and when the motif was present there, translation was even initiated at some non-AUG codons in sperm cells. However, unlike that in mammalian cells, a mainly effective nucleotide site to regulate the efficiency of translation initiation was not present at positions -1 to -4 in sperm cells. Meanwhile, different from somatic cells, sperm cells did not use eukaryotic translation initiation factor 1 to regulate the efficiency in a poor context consisting of the lowest frequency nucleotides. All these results contribute to our understanding of the cytoplasmic event of translation initiation in highly differentiated sperm cells.

DNA hypomethylation is the plausible driver of heat stress adaptation in Linum usitatissimum.

Heat stress has a significant impact on the climatic adaptation of flax, a cool-season economic crop. Genome-wide DNA methylation patterns are crucial for understanding how flax cultivars respond to heat adversities. It is worth noting that the DNA methylome in flax has yet to be investigated at the nucleotide level. Although heat stress above 40°C caused oxidative damage in flax leaves, 5-azacytidine, a hypomethylating agent, reduced this effect by 15%-24%. Differences in the expression of the LuMET1 (DNA methyltransferase) gene suggested that DNA methylation/demethylation may play a major role in the flax heat stress response. Thus, whole-genome bisulfite sequencing-derived DNA methylation profiles in flax, with or without heat stress and 5-azaC, were developed and analyzed here. In response to heat stress, a high percentage of significant differentially methylated regions (DMRs), particularly hypomethylated DMRs, were identified in the CHH nucleotide sequence context (H=A/T/C). Some of these DMRs overlapped with transposable element insertions. The majority of DMRs were discovered in intergenic regions, but several DMR loci were also found near genes relevant to heat stress response and epigenetic processes. These DMRs, in particular, are linked to CpG islands, implying a possible role in promoter methylation and gene silencing. The DMRs discovered in this study are crucial for understanding and identifying the key players in heat stress response in flax, which will help in developing climate-smart flax varieties.

Ribose-Map: A bioinformatics toolkit for ribonucleotide sequencing experiments

<h2>Abstract</h2> Current ribonucleotide mapping software is customized for data generated using only one type of ribonucleotide sequencing technique. To standardize the analysis of ribonucleotide sequencing experiments, we developed the bioinformatics toolkit, Ribose-Map. Through a series of analytical modules, Ribose-Map transforms raw sequencing data into summary datasets and publication-ready visualizations of results, allowing researchers to locate the genomic coordinates of ribonucleotides to single-nucleotide resolution, study the nucleotide sequence context of ribonucleotides, explore their genome-wide distribution, and identify hotspots of ribonucleotide incorporation. Ribose-Map does not require background knowledge of ribonucleotide sequencing analysis and assumes only basic command-line skills. Ribose-Map is available at https://github.com/agombolay/ribose-map.

Abstract 244: A computation method for noise reduction based on ultra-deep targeted sequencing data

Abstract Detection of circulating tumor DNA is a promising approach for the diagnosis and monitor of solid tumors. However, sensitive detection of ctDNA remains challenging when minor allele frequency (MAF) is ultra-low. Numerous noise reduction data algorithms have been developed to reduce next-generation sequencing (NGS) errors, but variant callers that can take advantage of ultra-deep sequencing data are lacking. Here we develop a bin-based genotyping (BinGo) method for canceling sequencing noise using ultra-deep (&amp;gt; 100,000 X depth) targeted sequencing data. A mutation reference sample with 44 mutations was diluted to different MAFs at 1% and 0.5%. One blood sample for proficiency test was included. We constructed targeted sequencing libraries by anchored multiplex enrichment for a panel of targets that cover the 44 mutations, and sequenced to a very high depth (Mean X ranged: 71K to 350K).After raw NGS reads were mapped by BWA MEM, BinGo tags unique molecular index (UMI) to each alignment (BAM) to preserve alignment signal of all bases in all raw reads. BinGo then bins all reads with the same UMI into one bin, generating millions of bins. Millions of parallel Samtools mpileup was enabled by the GNU parallel tool. Typical mpileup analyses use the human genome (requires ~4GB RAM) as a reference, making millions of such analyses essentially impossible in a regular computing server. To overcome this hurdle, BinGo uses tailor-made reference sequences comprising the targets only (&amp;lt; 100KB), drastically reduces computation memory consumption and enables millions of mpileup analyses. BinGo recompiles Samtools to increase its default max depth from 8,000 to 8,000,000 to accommodate the ultra-deep data. Variants derived from the same bin are effectively phased, thereby enabling calling of di-nuclei mutations such as BRAF V600K (GT&amp;gt;TT) and V600R (GT&amp;gt;AG). This binning also allows for sensitive calling of insertions and deletions even down to a single read. To profile noise patterns, BinGo uses Samtools tview to acquire every base in every alignment position of the target sequences in all raw reads. Finally, BinGo uses multiple statistical testing to call variants, including calling nucleotide changes over the background noise model that is defined by nucleotide change and sequence context; filtering out strand imbalance likely due to DNA template damage, and a low number of supporting ligation starting sites. With 10 ng DNA as input and at an expected MAF of 0.5%, BinGo demonstrates a 100% (44/44) analytical sensitivity and 100% specificity. When using 1 ng DNA as input and at an expected MAF of 0.5%, only one mutation was missed, leading to a sensitivity of 97.73% (43/44) and 100% specificity. Furthermore, BinGo detected in the blood DNA sample a 15-bp EGFR exon 19 deletion at 0.05%.Bingo integrates UMI-bin derived variants and raw sequencing background noise modeling, enabling analysis of ultra-deep targeted NGS data for rare variant detection. Citation Format: Weiwei Bian, Firaol Kebede, Zongli Zheng. A computation method for noise reduction based on ultra-deep targeted sequencing data [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 244.

Mapping ribonucleotides embedded in genomic DNA to single-nucleotide resolution using Ribose-Map.

The most common nonstandard nucleotides found in genomic DNA are ribonucleotides. Although ribonucleotides are frequently incorporated into DNA by replicative DNA polymerases, very little is known about the distribution and signatures of ribonucleotides incorporated into DNA. Recent advances in high-throughput ribonucleotide sequencing can capture the exact locations of ribonucleotides in genomic DNA. Ribose-Map is a user-friendly, standardized bioinformatics toolkit for the comprehensive analysis of ribonucleotide sequencing experiments. It allows researchers to map the locations of ribonucleotides in DNA to single-nucleotide resolution and identify biological signatures of ribonucleotide incorporation. In addition, it can be applied to data generated using any currently available high-throughput ribonucleotide sequencing technique, thus standardizing the analysis of ribonucleotide sequencing experiments and allowing direct comparisons of results. This protocol describes in detail how to use Ribose-Map to analyze ribonucleotide sequencing data, including preparing the reads for analysis, locating the genomic coordinates of ribonucleotides, exploring the genome-wide distribution of ribonucleotides, determining the nucleotide sequence context of ribonucleotides and identifying hotspots of ribonucleotide incorporation. Ribose-Map does not require background knowledge of ribonucleotide sequencing analysis and assumes only basic command-line skills. The protocol requires less than 3 h of computing time for most datasets and ~30 min of hands-on time. Ribose-Map is available at https://github.com/agombolay/ribose-map .

Translesion Synthesis Past 5-Formylcytosine-Mediated DNA-Peptide Cross-Links by hPolη Is Dependent on the Local DNA Sequence.

DNA-protein cross-links (DPCs) are unusually bulky DNA lesions that form when cellular proteins become trapped on DNA following exposure to ultraviolet light, free radicals, aldehydes, and transition metals. DPCs can also form endogenously when naturally occurring epigenetic marks [5-formyl cytosine (5fC)] in DNA react with lysine and arginine residues of histones to form Schiff base conjugates. Our previous studies revealed that DPCs inhibit DNA replication and transcription but can undergo proteolytic cleavage to produce smaller DNA-peptide conjugates. We have shown that 5fC-conjugated DNA-peptide cross-links (DpCs) placed within the CXA sequence (X = DpC) can be bypassed by human translesion synthesis (TLS) polymerases η and κ in an error-prone manner. However, the local nucleotide sequence context can have a strong effect on replication bypass of bulky lesions by influencing the geometry of the ternary complex among the DNA template, polymerase, and the incoming dNTP. In this work, we investigated polymerase bypass of 5fC-DNA-11-mer peptide cross-links placed in seven different sequence contexts (CXC, CXG, CXT, CXA, AXA, GXA, and TXA) in the presence of human TLS polymerase η. Primer extension products were analyzed by gel electrophoresis, and steady-state kinetics of the misincorporation of dAMP opposite the DpC lesion in different base sequence contexts was investigated. Our results revealed a strong impact of nearest neighbor base identity on polymerase η activity in the absence and presence of a DpC lesion. Molecular dynamics simulations were used to structurally explain the experimental findings. Our results suggest a possible role of local DNA sequence in promoting TLS-related mutational hot spots in the presence and absence of DpC lesions.

Analysis of the epitranscriptome with ion-pairing reagent free oligonucleotide mass spectrometry.

RNA modifications gain growing attention as a new frontier in the life sciences but with the rise of RNA vaccines also in biomedical drug development. Impeccable characterization of RNA modifications within their sequence context remains an analytical challenge. Oligonucleotide mass spectrometry (ON-MS), an approach similar to bottom-up proteome analysis, is capable of defining a short 5-15 nucleotide sequence context of an RNA modification while delivering information on the chemical character of the modified nucleotide. Commonly, ON-MS requires the use of ion pairing reagents for ON separation which is not compatible with most other MS-based applications and only few laboratories run dedicated MS instruments for the task. Here, we present an ON-MS technique which is independent of ion pairing reagents and can be used on any available mass spectrometer without risking its sensitivity for other analytes. In this chapter, we describe the experiments necessary for ON-MS method development, ON-MS application to native and synthetic RNAs and finally a guideline for data analysis.

Ribose-Map: a bioinformatics toolkit to map ribonucleotides embedded in genomic DNA.

Recent advances in high-throughput sequencing techniques have made it possible to tag ribonucleoside monophosphates (rNMPs) embedded in genomic DNA for sequencing. rNMP sequencing experiments generate large, complex datasets that require efficient, scalable software that can accurately map embedded rNMPs independently of the particular sequencing technique used. Current computational pipelines designed to map rNMPs embedded in genomic DNA are customized for data generated using only one type of rNMP sequencing technique. To standardize the processing and analysis of rNMP sequencing experiments, we developed Ribose-Map. Through a series of analytical modules, Ribose-Map transforms raw sequencing data into summary datasets and publication-ready visualizations of results, allowing biologists to identify sites of embedded rNMPs, study the nucleotide sequence context of these rNMPs and explore their genome-wide distribution. By accommodating data from any of the available rNMP sequencing techniques, Ribose-Map can increase the reproducibility of rNMP sequencing experiments and enable a head-to-head comparison of these experiments.

Genome-wide chemical mutagenesis screens allow unbiased saturation of the cancer genome and identification of drug resistance mutations

Drug resistance is an almost inevitable consequence of cancer therapy and ultimately proves fatal for the majority of patients. In many cases, this is the consequence of specific gene mutations that have the potential to be targeted to resensitize the tumor. The ability to uniformly saturate the genome with point mutations without chromosome or nucleotide sequence context bias would open the door to identify all putative drug resistance mutations in cancer models. Here, we describe such a method for elucidating drug resistance mechanisms using genome-wide chemical mutagenesis allied to next-generation sequencing. We show that chemically mutagenizing the genome of cancer cells dramatically increases the number of drug-resistant clones and allows the detection of both known and novel drug resistance mutations. We used an efficient computational process that allows for the rapid identification of involved pathways and druggable targets. Such a priori knowledge would greatly empower serial monitoring strategies for drug resistance in the clinic as well as the development of trials for drug-resistant patients.

Guanine holes are prominent targets for mutation in cancer and inherited disease.

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G•C bp in the context of all 64 5′-NGNN-3′ motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

Viral Genome Segmentation Can Result from a Trade-Off between Genetic Content and Particle Stability

The evolutionary benefit of viral genome segmentation is a classical, yet unsolved question in evolutionary biology and RNA genetics. Theoretical studies anticipated that replication of shorter RNA segments could provide a replicative advantage over standard size genomes. However, this question has remained elusive to experimentalists because of the lack of a proper viral model system. Here we present a study with a stable segmented bipartite RNA virus and its ancestor non-segmented counterpart, in an identical genomic nucleotide sequence context. Results of RNA replication, protein expression, competition experiments, and inactivation of infectious particles point to a non-replicative trait, the particle stability, as the main driver of fitness gain of segmented genomes. Accordingly, measurements of the volume occupation of the genome inside viral capsids indicate that packaging shorter genomes involves a relaxation of the packaging density that is energetically favourable. The empirical observations are used to design a computational model that predicts the existence of a critical multiplicity of infection for domination of segmented over standard types. Our experiments suggest that viral segmented genomes may have arisen as a molecular solution for the trade-off between genome length and particle stability. Genome segmentation allows maximizing the genetic content without the detrimental effect in stability derived from incresing genome length.

Suppression of Nonsense Mutations in Rett Syndrome by Aminoglycoside Antibiotics

Rett Syndrome (RTT) is caused in more than 60% of cases by nonsense mutations in the MECP2 gene. So far, no curative therapy for RTT has become available. In other genetic disorders, it has been shown that aminoglycosides can cause a read-through of nonsense mutations with an efficiency of up to 20%. The aim of this study was to evaluate if this therapeutic concept is applicable to RTT. HeLa cells were transfected with eukaryotic expression vectors carrying mutant alleles of frequently occurring MECP2 nonsense mutations that were N-terminally fused to a FLAG tag. Transfected cells were incubated 24 h in the presence of gentamicin. The expression of full-length protein was analyzed by Western blotting and immunofluorescent cell staining. In the presence of gentamicin a read-through varying between 10 and 21.8% was found, depending on the nucleotide sequence context of the nonsense mutations. The full-length protein was located correctly in the nucleus. We have shown that aminoglycoside-mediated read-through of nonsense mutations in the MECP2 gene can be achieved in vitro with efficiency comparable with that seen in other disorders.

The intraleader AUG nucleotide sequence context is important for equine arteritis virus replication.

The 5(-terminal leader sequence of the equine arteritis virus (EAV) genome contains an open reading frame (ORF) with an AUG codon in a suboptimal context for initiation of protein synthesis. To investigate the significance of this intraleader ORF (ILO), an expression plasmid was generated carrying a DNA copy of the subgenomic mRNA7 behind a T7 promoter. Capped RNA transcribed from this construct was shown to direct, in an in vitro translation system, the synthesis of leader peptide as well as N protein. Site-directed mutations aimed to either optimize or weaken the sequence context of the ILO start codon affected leader peptide synthesis as predicted; no peptide was detected when the initiation codon was incapacitated. Translation of the downstream N gene was inversely affected by leader peptide production, consistent with a ribosomal leaky scanning mechanism. To investigate the role of the leader peptide in the EAV replication life cycle we generated, using an infectious EAV cDNA clone, two mutant viruses in one of which the ILO start codon was in an optimal Kozak context for translation initiation while in the other the codon was again incapacitated. Surprizingly, both mutant viruses were equally viable and exhibited similar phenotypes in BHK-21 cells. However, their replication kinetics and viral yields were reduced relative to that of the wild-type parental virus, as were their plaque sizes. Importantly, the mutations introduced into the viruses appeared to be rapidly and precisely repaired upon passaging. Already after one viral passage a significant fraction of the viruses had regained the wild-type sequence as well as its phenotype. The results demonstrate that EAV replication is not dependent on the synthesis of the intraleader peptide. Rather, the leader peptide does not seem to have any function in the EAV life cycle. As we discuss, the available data indicate that the ILO 5( nucleotide sequence per se, not its functioning in translation initiation, is of critical importance for EAV replication.

Cumulative Author Index for MUTREV 2003

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The hotspots might also reflect structural and functional features of the respective DNA sequences. When mutations in a gene are identified using a particular experimental system, resulting hotspots could reflect the properties of the gene product and the mutant selection scheme. Analysis of the nucleotide sequence context of hotspots can provide information on the molecular mechanisms of mutagenesis. However, the determinants of mutation frequency and specificity are complex, and there are many analytical methods for their study. Here we review computational approaches for analyzing mutation spectra (distribution of mutations along the target genes) that include many mutable (detectable) positions. The following methods are reviewed: derivation of a consensus sequence, application of regression approaches to correlate nucleotide sequence features with mutation frequency, mutation hotspot prediction, analysis of oligonucleotide composition of regions containing mutations, pairwise comparison of mutation spectra, analysis of multiple spectra, and analysis of “context-free” characteristics. The advantages and pitfalls of these methods are discussed and illustrated by examples from the literature. The most reliable analyses were obtained when several methods were combined and information from theoretical analysis and experimental observations was considered simultaneously. Simple, robust approaches should be used with small samples of mutations, whereas combinations of simple and complex approaches may be required for large samples. We discuss several well-documented studies where analysis of mutation spectra has substantially contributed to the current understanding of molecular mechanisms of mutagenesis. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and DNA repair, replication, and modification enzymes, and the analysis of hotspot context provides evidence of such interactions.

Synchronous Loss of Quasispecies Memory in Parallel Viral Lineages: A Deterministic Feature of Viral Quasispecies

Viral quasispecies are endowed with a memory of their past evolutionary history in the form of minority genomes of their mutant spectra. To determine the fate of memory genomes in evolving viral quasispecies, we have measured memory levels of antigenic variant of foot-and-mouth disease virus (FMDV) RED, which includes an Arg-Glu-Asp (RED) at a surface antigenic loop of the viral capsid. The RED reverted to the standard Arg-Gly-Asp (RGD), and the RED remained as memory in the evolving quasispecies. In four parallel evolutionary lineages, memory reduction followed a strikingly similar pattern, and at passage 60 memory levels were indistinguishable from those of control populations (devoid of memory). Nucleotide sequence analyses indicated that memory loss occurred synchronously despite its ultimate molecular basis being the stochastic occurrence of mutations in the evolving quasispecies. These results on the kinetics of memory levels have unveiled a deterministic feature of viral quasispecies. Molecular mechanisms that may underlie synchronous memory loss are the averaging of noise signals derived from mutational input, and constraints to genome diversification imposed by a nucleotide sequence context in the viral genome. Possible implications of the behaviour of complex, adaptive viral systems as experimental models to address primary mechanisms of neurological memory are discussed.

Concerted changes in the nucleotide sequences of the intragenic promoter regions of eukaryotic genes for tRNAs of all specificities.

RNA polymerase III promoter is located within the coding region in all eukaryotic tRNA genes, whereas in prokaryotic tRNA genes, the promoter is located upstream of the transcription initiation site. We analyzed the nucleotide sequence context of the A and B boxes of RNA polymerase III promoters from different unicellular eukaryotes, plants, and animals and the homologous sequences in the tRNA genes of prokaryotic species (Archaea, Eubacteria). The long and short sequence variants of the A box are nonrandomly distributed across different types of the eukaryotic tRNA genes. In contrast, the sequences from the prokaryotic counterparts of the long and short variants are randomly scattered among prokaryotic tRNA genes. The sequence diversity of the 3' portion of the A box in the eukaryotic tRNA genes was sharply reduced compared to that of the homologous segments of the prokaryotic tRNA genes. Analysis of the frequencies of oligonucleotide variants from different parts of the A box suggested that the prominent changes in the structure of the A box occurred at about the same time in all tRNA genes during the prokaryote eukaryote evolutionary transition.

Theoretical analysis of mutation hotspots and their DNA sequence context specificity.

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The hotspots might also reflect structural and functional features of the respective DNA sequences. When mutations in a gene are identified using a particular experimental system, resulting hotspots could reflect the properties of the gene product and the mutant selection scheme. Analysis of the nucleotide sequence context of hotspots can provide information on the molecular mechanisms of mutagenesis. However, the determinants of mutation frequency and specificity are complex, and there are many analytical methods for their study. Here we review computational approaches for analyzing mutation spectra (distribution of mutations along the target genes) that include many mutable (detectable) positions. The following methods are reviewed: derivation of a consensus sequence, application of regression approaches to correlate nucleotide sequence features with mutation frequency, mutation hotspot prediction, analysis of oligonucleotide composition of regions containing mutations, pairwise comparison of mutation spectra, analysis of multiple spectra, and analysis of "context-free" characteristics. The advantages and pitfalls of these methods are discussed and illustrated by examples from the literature. The most reliable analyses were obtained when several methods were combined and information from theoretical analysis and experimental observations was considered simultaneously. Simple, robust approaches should be used with small samples of mutations, whereas combinations of simple and complex approaches may be required for large samples. We discuss several well-documented studies where analysis of mutation spectra has substantially contributed to the current understanding of molecular mechanisms of mutagenesis. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and DNA repair, replication, and modification enzymes, and the analysis of hotspot context provides evidence of such interactions.