Higher mutation rate helps to rescue genes from the elimination by selection

Directional mutation pressure associated with replication processes is the main cause of the asymmetry between the leading and lagging DNA strands in bacterial genomes. On the other hand, the asymmetry between sense and antisense strands of protein coding sequences is a result of both mutation and selection pressures. Thus, there are two different ways of superposition of the sense strand, on the leading or lagging strand. Besides many other implications of these two possible situations, one seems to be very important - because of the asymmetric replication-associated mutation pressure, the mutation rate of genes depends on their location. Using Monte Carlo methods, we have simulated, under experimentally determined directional mutation pressure, the divergence rate and the elimination rate of genes depending on their location in respect to the leading/lagging DNA strands in the asymmetric prokaryotic genome. We have found that the best survival strategy for the majority of genes is to sometimes switch between DNA strands. Paradoxically, this strategy results in higher substitution rates but remains in agreement with observations in bacterial genomes that such inversions are very frequent and divergence rate between homologs lying on different DNA strands is very high.     

Proteome composition and codon usage in spirochaetes: species-specific and DNA strand-specific mutational biases

The genomes of the spirochaetes Borrelia burgdorferi and Treponema pallidum show strong strand-specific skews in nucleotide composition, with the leading strand in replication being richer in G and T than the lagging strand in both species. This mutation bias results in codon usage and amino acid composition patterns that are significantly different between genes encoded on the two strands, in both species. There are also substantial differences between the species, with T.pallidum having a much higher G+C content than B. burgdorferi. These changes in amino acid and codon compositions represent neutral sequence change that has been caused by strong strand- and species-specific mutation pressures. Genes that have been relocated between the leading and lagging strands since B. burgdorferi and T.pallidum diverged from a common ancestor now show codon and amino acid compositions typical of their current locations. There is no evidence that translational selection operates on codon usage in highly expressed genes in these species, and the primary influence on codon usage is whether a gene is transcribed in the same direction as replication, or opposite to it. The dnaA gene in both species has codon usage patterns distinctive of a lagging strand gene, indicating that the origin of replication lies downstream of this gene, possibly within dnaN. Our findings strongly suggest that gene-finding algorithms that ignore variability within the genome may be flawed.     

Differential DNA secondary structure-mediated deletion mutation in the leading and lagging strands.

The frequencies of deletion of short sequences (mutation inserts) inserted into the chloramphenicol acetyl-transferase (CAT) gene were measured for pBR325 and pBR523, in which the orientation of the CAT gene was reversed, in Escherichia coli. Reversal of the CAT gene changes the relationship between the transcribed strand and the leading and lagging strands of the DNA replication fork in pBR325-based plasmids. Deletion of these mutation inserts may be mediated by slipped misalignment during DNA replication. Symmetrical sequences, in which the same potential DNA structural misalignment can form in both the leading and lagging strands, exhibited an approximately twofold difference in the deletion frequencies upon reversal of the CAT gene. Sequences that contained an inverted repeat that was asymmetric with respect to flanking direct repeats were designed. With asymmetric mutation inserts, different misaligned structural intermediates could form in the leading and lagging strands, depending on the orientation of the insert and/or of the CAT gene. When slippage could be stabilized by a hairpin in the lagging strand, thereby forming a three-way junction, deletion occurred by up to 50-fold more frequently than when this structure formed in the leading strand. These results support the model that slipped misalignment involving DNA secondary structure occurs preferentially in the lagging strand during DNA replication.     

Balancing eukaryotic replication asymmetry with replication fidelity

Coordinated replication of eukaryotic nuclear genomes is asymmetric, with copying of a leading strand template preceding discontinuous copying of the lagging strand template. Replication is catalyzed by DNA polymerases α, δ and ?, enzymes that are related yet differ in physical and biochemical properties, including fidelity. Recent studies suggest that Pol ? is normally the primary leading strand replicase, whereas most synthesis by Pol δ occurs during lagging strand replication. New studies show that replication asymmetry can generate strand-specific genome instability resulting from biased deoxynucleotide pools and unrepaired ribonucleotides incorporated into DNA during replication, and that the eukaryotic replication machinery has evolved to most efficiently correct those replication errors that are made at the highest rates.     

Rearrangements between differently replicating DNA strands in asymmetric bacterial genomes

Many bacterial genomes are under asymmetric mutational pressure which introduces compositional asymmetry into DNA molecule resulting in many biases in coding structure of chromosomes. One of the processes affected by the asymmetry is translocation changing the position of the coding sequence on chromosome in respect to the orientation on the leading and lagging DNA strand. When analysing sets of paralogs in 50 genomes, we found that the number of observed genes which switched their positions on DNA strand is lowest for genomes with the highest DNA asymmetry. However, the number of orthologs which changed DNA strand increases with the phylogenetic distance between the compared genomes. Nevertheless, there is a fraction of coding sequences that stay on the leading strand in all analysed genomes, whereas there are no sequences that stay always on the lagging strand. Since sequences diverge very fast after switching the DNA strand, this bias in mobility of sequences is responsible, in part, for higher divergence rates among some of coding sequences located on the lagging DNA strand.     

强烈的复制链组成偏差——某些专性寄生(共生)菌特有的基因组学特征

DNA复制是一个不对称的过程。不对称的一个体现是复制链分为前导链和滞后链,前者连续复制而后者的复制却不连续。这种不对称最终导致两条链上核酸组成的不对称。链特异的核酸组成偏差最先发现于棘皮类动物和脊椎动物的线粒体DNA上。随着全基因组的大量测序,越来越多的细菌被发现具有类似的链特异的核酸组成偏差,甚至很多真核生物的染色体基因组也被发现具有类似偏差。在某些细菌中,链特异的组成偏差强烈到足以使前导链和滞后链的基因间具有分离的密码子使用。至今,共有11种细菌被发现具有和复制相关的分离的密码子使用。这11种细菌无一例外都属于专性寄生(共生)菌。对于链组成偏差产生的内在机制以及特定细菌具有强烈链组成偏差的成因,目前学术界尚没有统一的理论解释。文章对这一遗传学及基因组学的重要问题进行了综述和展望。     

基于复制-转录碰撞的突变和进化意义

复制和转录机器会同时使用相同的DNA区域作为模板,因此复制和转录不可避免地以头对头或追尾方式相互碰撞。头对头碰撞和追尾碰撞均会导致复制机器停留,从而造成DNA损伤和基因组不稳定。就基因组完整性而言,头对头碰撞比追尾碰撞的后果更严重。本文回顾总结了复制-转录冲突的解决机制和进化影响。相对于前导链,滞后链上非同义（氨基酸改变）突变的发生率更高,并且滞后链上基因的高频诱变取决于转录本和基因大小,因此,较快的适应性突变发生在滞后链上。头对头基因的高度转录增加了复制过程中响应压力的突变率。无论是头对头还是追尾模式,复制-转录冲突都可能是适应性进化的驱动力。     

GC skew in protein-coding genes between the leading and lagging strands in bacterial genomes: new substitution models incorporating strand bias

The DNA strands in most prokaryotic genomes experience strand-biased spontaneous mutation, especially C→T mutations produced by deamination that occur preferentially in the leading strand. This has often been invoked to account for the asymmetry in nucleotide composition, typically measured by GC skew, between the leading and the lagging strand. Casting such strand asymmetry in the framework of a nucleotide substitution model is important for understanding genomic evolution and phylogenetic reconstruction. We present a substitution model showing that the increased C→T mutation will lead to positive GC skew in one strand but negative GC skew in the other, with greater C→T mutation pressure associated with greater differences in GC skew between the leading and the lagging strand. However, the model based on mutation bias alone does not predict any positive correlation in GC skew between the leading and lagging strands. We computed GC skew for coding sequences collinear with the leading and lagging strands across 339 prokaryotic genomes and found a strong and positive correlation in GC skew between the two strands. We show that the observed positive correlation can be satisfactorily explained by an improved substitution model with one additional parameter incorporating a general trend of C avoidance.     

Duplications between direct repeats stabilized by DNA secondary structure occur preferentially in the leading strand during DNA replication

To ascertain a leading or lagging strand preference for duplication mutations, several short DNA sequences, i.e. mutation inserts, were designed that should demonstrate an asymmetric propensity for duplication mutations in the two complementary DNA strands during replication. The design of the mutation insert involved a 7-bp quasi inverted repeat that forms a remarkably stable hairpin in one DNA strand, but not the other. The inverted repeat is asymmetrically placed between flanking direct repeats. This sequence was cloned into a modified chloramphenicol acetyltransferase (CAT) gene containing a −1 frameshift mutation. Duplication of the mutation insert restores the reading frame of the CAT gene resulting in a chloramphenicol resistant phenotype. The mutation insert showed greater than a 200-fold preference for duplication mutations during leading strand, compared with lagging strand, replication. This result suggests that misalignment stabilized by DNA secondary structure, leading to duplication between direct repeats, occurred preferentially during leading strand synthesis.     

Fidelity of replication of the leading and the lagging DNA strands opposite N-methyl-N-nitrosourea-induced DNA damage in human cells.

Semi-conservative replication of double-stranded DNA in eukaryotic cells is an asymmetric process involving leading and lagging strand synthesis and different DNA polymerases. We report a study to analyze the effect of these asymmetries when the replication machinery encounters alkylation-induced DNA adducts. The model system is an EBV-derived shuttle vector which replicates in synchrony with the host human cells and carries as marker gene the bacterial gpt gene. A preferential distribution of N-methyl-N-nitrosourea (MNU)-induced mutations in the non transcribed DNA strand of the shuttle vector pF1-EBV was previously reported. The hypermutated strand was the leading strand. To test whether the different fidelity of DNA polymerases synthesizing the leading and the lagging strands might contribute to MNU-induced mutation distribution the mutagenesis study was repeated on the shuttle vector pTF-EBV which contains the gpt gene in the inverted orientation. We show that the base substitution error rates on an alkylated substrate are similar for the replication of the leading and lagging strands. Moreover, we present evidence that the fidelity of replication opposite O6-methylguanine adducts of both the leading and lagging strands is not affected by the 3' flanking base. The preferential targeting of mutations after replication of alkylated DNA is mainly driven by the base at the 5' side of the G residues.     

76种细菌DNA双链碱基使用频率的比较及其意义

应用生物信息学方法,对已完成测序的76种细菌基因组进行比较,分析细菌基因组中编码区及密码子上碱基使用频率情况,结果显示：1.先导链与滞后链上在编码区的碱基使用频率无明显差异且显著正相关;2.先地链与滞后链在第一,第二,第三密码子碱基使用频率基本一致且显著正相关,结果表明,选择压力及自然突变对DNA双链总体碱基分布的影响相等。     

Evidence for preferential mismatch repair of lagging strand DNA replication errors in yeast

Duplex DNA is replicated in the 5'-3' direction by coordinated copying of leading and lagging strand templates with somewhat different proteins and mechanics, providing the potential for differences in the fidelity of replication of the two strands. We previously showed that in Saccharomyces cerevisiae, active replication origins establish a strand bias in the rate of base substitutions resulting from replication of unrepaired 8-oxo-guanine (GO) in DNA. Lower mutagenesis was associated with replicating lagging strand templates. Here, we test the hypothesis that this bias is due to more efficient repair of lagging stand mismatches by measuring mutation rates in ogg1 strains with a reporter allele in two orientations at loci on opposite sides of a replication origin on chromosome III. We compare a MMR-proficient strain to strains deleted for the MMR genes MSH2, MSH6, MLH1, or EXOI. Loss of MMR reduces the strand bias by preferentially increasing mutagenesis for lagging strand replication. We conclude that GO-A mismatches generated during lagging strand replication are more efficiently repaired. This is consistent with the hypothesis that 5' ends of Okazaki fragments and PCNA, present at high density during lagging strand replication, are used as strand discrimination signals for mismatch repair in vivo.     

Strong asymmetric mutation bias in endosymbiont genomes coincide with loss of genes for replication restart pathways

A large majority of bacterial genomes show strand asymmetry, such that G and T preferentially accumulate on the leading strand. The mechanisms are unknown, but cytosine deaminations are thought to play an important role. Here, we have examined DNA strand asymmetry in three strains of the aphid endosymbiont Buchnera aphidicola. These are phylogenetically related, have similar genomic GC contents, and conserved gene order structures, yet B. aphidicola (Bp) shows a fourfold higher replication-induced strand bias than B. aphidicola (Sg) and (Ap). We rule out an increase in the overall substitution frequency as the major cause of the stronger strand bias in B. aphidicola (Bp). Instead, the results suggest that the higher GC skew in this species is caused by a different spectrum of mutations, including a relatively higher frequency of C to T mutations on the leading strand and/or of G to A mutations on the lagging strand. A comparative analysis of 20 gamma-proteobacterial genomes revealed that endosymbiont genomes lacking recA and other genes involved in replication restart processes, such as priA, which codes for primosomal helicase PriA, displayed the strongest strand bias. We hypothesize that cytosine deaminations accumulate during single-strand exposure at arrested replication forks and that inefficient restart mechanisms may lead to high DNA strand asymmetry in bacterial genomes.     

The major roles of DNA polymerases epsilon and delta at the eukaryotic replication fork are evolutionarily conserved

Coordinated replication of eukaryotic genomes is intrinsically asymmetric, with continuous leading strand synthesis preceding discontinuous lagging strand synthesis. Here we provide two types of evidence indicating that, in fission yeast, these two biosynthetic tasks are performed by two different replicases. First, in Schizosaccharomyces pombe strains encoding a polδ-L591M mutator allele, base substitutions in reporter genes placed in opposite orientations relative to a well-characterized replication origin are strand-specific and distributed in patterns implying that Polδ is primarily involved in lagging strand replication. Second, in strains encoding a polε-M630F allele and lacking the ability to repair rNMPs in DNA due to a defect in RNase H2, rNMPs are selectively observed in nascent leading strand DNA. The latter observation demonstrates that abundant rNMP incorporation during replication can be tolerated and that they are normally removed in an RNase H2-dependent manner. This provides strong physical evidence that Polε is the primary leading strand replicase. Collectively, these data and earlier results in budding yeast indicate that the major roles of Polδ and Polε at the eukaryotic replication fork are evolutionarily conserved.     

Strand bias influences the mechanism of gene editing directed by single-stranded DNA oligonucleotides

Gene editing directed by modified single-stranded DNA oligonucleotides has been used to alter a single base pair in a variety of biological systems. It is likely that gene editing is facilitated by the direct incorporation of the oligonucleotides via replication and/or by direct conversion, most likely through the DNA mismatch repair pathway. The phenomenon of strand bias, however, as well as its importance to the gene editing reaction itself, has yet to be elucidated in terms of mechanism. We have taken a reductionist approach by using a genetic readout in Eschericha coli and a plasmid-based selectable system to evaluate the influence of strand bias on the mechanism of gene editing. We show that oligonucleotides (ODNs) designed to anneal to the lagging strand generate 100-fold greater 'editing' efficiency than 'those that anneal to' the leading strand. The majority of editing events (～70%) occur by the incorporation of the ODN during replication within the lagging strand. Conversely, ODNs that anneal to the leading strand generate fewer editing events although this event may follow either the incorporation or direct conversion pathway. In general, the influence of DNA replication is independent of which ODN is used suggesting that the importance of strand bias is a reflection of the underlying mechanism used to carry out gene editing.     

The Differential Killing of Genes by Inversions in Prokaryotic Genomes

We have elaborated a method which has allowed us to estimate the direction of translocation of orthologs which have changed, during the phylogeny, their positions on chromosome in respect to the leading or lagging role of DNA strands. We have shown that the relative number of translocations which have switched positions of genes from the leading to the lagging DNA strand is lower than the number of translocations which have transferred genes from the lagging strand to the leading strand of prokaryotic genomes. This paradox could be explained by assuming that the stronger mutation pressure and selection after inversion preferentially eliminate genes transferred from the leading to the lagging DNA strand. Received: 12 December 2000 / Accepted: 20 April 2001     

Division of labor at the eukaryotic replication fork

DNA polymerase delta (Pol delta) and DNA polymerase epsilon (Pol epsilon) are both required for efficient replication of the nuclear genome, yet the division of labor between these enzymes has remained unclear for many years. Here we investigate the contribution of Pol delta to replication of the leading and lagging strand templates in Saccharomyces cerevisiae using a mutant Pol delta allele (pol3-L612M) whose error rate is higher for one mismatch (e.g., T x dGTP) than for its complement (A x dCTP). We find that strand-specific mutation rates strongly depend on the orientation of a reporter gene relative to an adjacent replication origin, in a manner implying that >90% of Pol delta replication is performed using the lagging strand template. When combined with recent evidence implicating Pol epsilon in leading strand replication, these data support a model of the replication fork wherein the leading and lagging strand templates are primarily copied by Pol epsilon and Pol delta, respectively.