The mutational profile of the yeast genome is shaped by replication

Despite the scrutiny that has been directed for years at the yeast genome, relatively little is known about the impact of replication on the substitution dynamics in Saccharomyces cerevisiae. Here, we show that the mutation rate increases with the replication timing by more than 30% between the earliest and the latest replicating regions. In addition, we found a mutational asymmetry associated with the polarity of replication resulting in higher rates of substitutions toward C and A than toward G and T in leading strands (reciprocally more substitutions toward G and T in lagging strands). Such mutational asymmetries applied over long evolutionary periods should generate compositional skews between the two DNA strands. Thus, we show that the leading replicating strands present an excess of C over G and of A over T in the genome of S. cerevisiae (reciprocally an excess of G + T over C + A in lagging strands). We also show that the nucleotide frequencies at mutational equilibrium predict a compositional skew at equilibrium very close to the observed skew between leading and lagging strands, suggesting that compositional equilibrium has been nearly attained in the present day genome of S. cerevisiae. Surprisingly, the direction of this skew is inverted compared with the one in the human genome.     

DNA asymmetric strand bias affects the amino acid composition of mitochondrial proteins.

Variations in GC content between genomes have been extensively documented. Genomes with comparable GC contents can, however, still differ in the apportionment of the G and C nucleotides between the two DNA strands. This asymmetric strand bias is known as GC skew. Here, we have investigated the impact of differences in nucleotide skew on the amino acid composition of the encoded proteins. We compared orthologous genes between animal mitochondrial genomes that show large differences in GC and AT skews. Specifically, we compared the mitochondrial genomes of mammals, which are characterized by a negative GC skew and a positive AT skew, to those of flatworms, which show the opposite skews for both GC and AT base pairs. We found that the mammalian proteins are highly enriched in amino acids encoded by CA-rich codons (as predicted by their negative GC and positive AT skews), whereas their flatworm orthologs were enriched in amino acids encoded by GT-rich codons (also as predicted from their skews). We found that these differences in mitochondrial strand asymmetry (measured as GC and AT skews) can have very large, predictable effects on the composition of the encoded proteins.     

Evidence from simulation studies for selective constraints on the codon usage of the Angiosperm psbA gene

The codon usage of the Angiosperm psbA gene is atypical for flowering plant chloroplast genes but similar to the codon usage observed in highly expressed plastid genes from some other Plantae, particularly Chlorobionta, lineages. The pattern of codon bias in these genes is suggestive of selection for a set of translationally optimal codons but the degree of bias towards these optimal codons is much weaker in the flowering plant psbA gene than in high expression plastid genes from lineages such as certain green algal groups. Two scenarios have been proposed to explain these observations. One is that the flowering plant psbA gene is currently under weak selective constraints for translation efficiency, the other is that there are no current selective constraints and we are observing the remnants of an ancestral codon adaptation that is decaying under mutational pressure. We test these two models using simulations studies that incorporate the context-dependent mutational properties of plant chloroplast DNA. We first reconstruct ancestral sequences and then simulate their evolution in the absence of selection on codon usage by using mutation dynamics estimated from intergenic regions. The results show that psbA has a significantly higher level of codon adaptation than expected while other chloroplast genes are within the range predicted by the simulations. These results suggest that there have been selective constraints on the codon usage of the flowering plant psbA gene during Angiosperm evolution.     

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.     

Asymmetric directional mutation pressures in bacteria

Background

When there are no strand-specific biases in mutation and selection rates (that is, in the substitution rates) between the two strands of DNA, the average nucleotide composition is theoretically expected to be A = T and G = C within each strand. Deviations from these equalities are therefore evidence for an asymmetry in selection and/or mutation between the two strands. By focusing on weakly selected regions that could be oriented with respect to replication in 43 out of 51 completely sequenced bacterial chromosomes, we have been able to detect asymmetric directional mutation pressures.

Results

Most of the 43 chromosomes were found to be relatively enriched in G over C and T over A, and slightly depleted in G+C, in their weakly selected positions (intergenic regions and third codon positions) in the leading strand compared with the lagging strand. Deviations from A = T and G = C were highly correlated between third codon positions and intergenic regions, with a lower degree of deviation in intergenic regions, and were not correlated with overall genomic G+C content.

Conclusions

During the course of bacterial chromosome evolution, the effects of asymmetric directional mutation pressures are commonly observed in weakly selected positions. The degree of deviation from equality is highly variable among species, and within species is higher in third codon positions than in intergenic regions. The orientation of these effects is almost universal and is compatible in most cases with the hypothesis of an excess of cytosine deamination in the single-stranded state during DNA replication. However, the variation in G+C content between species is influenced by factors other than asymmetric mutation pressure.

     

The mitochondrial genome of red-necked phalarope <Emphasis Type="Italic">Phalaropus lobatus</Emphasis> (Charadriiformes: Scolopacidae) and phylogeny analysis among Scolopacidae

The red-necked phalarope is a wonderful species with specific morphological characters and lifestyles. Mitochondrial genomes, encoding necessary proteins involved in the system of energy metabolism, are important for the evolution and adaption of species. In this study, we determined the complete mitogenome sequence of Phalaropus lobatus (Charadriiformes: Scolopacidae). The circular genome is 16714 bp in size, containing 13 PCGs, two ribosomal RNAs and 22 tRNAs and a high AT-rich control region. The AT skew and GC skew of major strand is positive and negative respectively. Most of PCGs are biased towards A-rich except ND1. A codon usage analysis shows that 3 start codons (ATG, GTG and ATA), 4 stop codons (TAA, TAG, AGG, AGA) and two incomplete terminate codons (T–). Twenty two transfer RNAs have the typical cloverleaf structure, and a total of ten base pairs are mismatched throughout the nine tRNA genes. The phylogenetic tree based on 13 PCGs and 2 rRNA genes indicates that monophyly of the family and genus Phalaropus is close to genus Xenus plus Tringa. The analysis of selective pressure shows 13 protein-coding genes are evolving under the purifying selection and P. lobatus is different from other Scolopacidae species on the selective pressure of gene ND4. This study helps us know the inherent mechanism of mitochondrial structure and natural selection.     

Mitochondrial Gene Rearrangements and Partial Genome Duplications Detected by Multigene Asymmetric Compositional Bias Analysis

Asymmetric compositional and mutation bias between the two strands occurs in mitochondrial genomes, and an asymmetric mechanism of mtDNA replication is a potential source of this bias. Some evidence indicates that during replication the heavy strand is subject to a gradient of time spent in a single-stranded state (D _ssH) and a gradient of mutational damage. The nucleotide composition bias among genes varies with D _ssH. Consequently, partial genome duplications (PGD) will alter the skew for genes located downstream of the duplication, relatively to nascent light strand synthesis, and in the same way, gene rearrangements (GRr) will affect genes by changing their skews. We examined cases where there had been PGD or GRr and determined whether this left a trace in the form of unusual patterns of base composition. We compared the skew of genes differently located on the mtDNA genome of previously published whole mtDNA genomes from amphibians, a group that shows considerable levels of both GRr and PGD. After observing a significant correlation between AT and GC skew with D _ssH at fourfold redundant sites, we ran our analysis and detected 31.3% of the species with GRr and/or PGD. By comparing the nucleotide composition at fourfold redundant sites in normal and “abnormal” species, we found that A/C variation occurs and is associated with GRr/PGD. These results show that by analyzing the nucleotide skews of only three genes, it may be possible to predict some mitochondrial GRr and/or PGD without knowing the complete mtDNA genome sequence. [Reviewing Editor: Dr. David Pollock]     

斐豹蛱蝶线粒体基因组全序列的测定和分析(英文)

该研究对斐豹蛱蝶(Argyreus hyperbius)(鳞翅目:蛱蝶科)线粒体基因组全序列进行了测定和初步分析。结果表明:斐豹蛱蝶线粒体基因全序列全长为15156bp,包含13个蛋白质编码基因、22个tRNA和2个rRNA基因以及1个非编码的A+T富集区,基因排列顺序与其它鳞翅目种类一致;线粒体全序列核苷酸组成和密码子使用显示出明显的A+T偏好(80.8%)和轻微的AT偏移(AT skew,?0.019)。基因组中共存在11个2～52bp不等的基因间隔区,总长96bp;以及14个1～8bp不等的基因重叠区,总长34bp。除COI以CGA作为起始密码子外,13个蛋白质编码基因中的其余12个基因是以ATN作为起始密码子。除COI和COII基因是以单独的一个T为终止密码子,其余11个蛋白质编码基因都是以TAA结尾的。除了缺少DHU臂的tRNASer(AGN),其余的tRNA基因都显示典型的三叶草结构。tRNA(AGN)和ND1之间的基因间隔区包含一个ATACTAA结构域,这个结构域在鳞翅目中是保守的。A+T富集区没有较大的多拷贝重复序列,但是包含一些微小重复结构:ATAGA结构域下游的20bp poly-T结构,ATTTA结构域后的(AT)9重复,以及位于tRNAMet上游的5bp poly-A结构等。这项研究所揭示的斐豹蛱蝶的线粒体基因组特征,不仅为认识蛱蝶科的遗传多样性贡献数据,而且对于该物种的保护生物学、群体遗传学、谱系地理及演化研究等具有重要意义。     

Evidence of Selection upon Genomic GC-Content in Bacteria

The genomic GC-content of bacteria varies dramatically, from less than 20% to more than 70%. This variation is generally ascribed to differences in the pattern of mutation between bacteria. Here we test this hypothesis by examining patterns of synonymous polymorphism using datasets from 149 bacterial species. We find a large excess of synonymous GC→AT mutations over AT→GC mutations segregating in all but the most AT-rich bacteria, across a broad range of phylogenetically diverse species. We show that the excess of GC→AT mutations is inconsistent with mutation bias, since it would imply that most GC-rich bacteria are declining in GC-content; such a pattern would be unsustainable. We also show that the patterns are probably not due to translational selection or biased gene conversion, because optimal codons tend to be AT-rich, and the excess of GC→AT SNPs is observed in datasets with no evidence of recombination. We therefore conclude that there is selection to increase synonymous GC-content in many species. Since synonymous GC-content is highly correlated to genomic GC-content, we further conclude that there is selection on genomic base composition in many bacteria.     

The complete mitochondrial genome of the Chinese oak silkmoth, Antheraea pernyi (Lepidoptera: Saturniidae)

We determined the complete nucleotide sequence of the mitogenome from Chinese oak silkmoth, Antheraea pernyi (Lepidoptera: Saturniidae). The 15,566 bp circular genome contains atypical gene organization and order for lepidopteran mitogenomes. The mitogenome contains the lowest A+T content (80.16%) among the known lepidopteran mitogenome sequences. An unusual feature is the occurrence of more Ts than As, with a slightly negative AT skewness (−0.021), in the composition of the major genome strand. All protein-coding genes are initiated by ATN codons, except for cytochrome oxidase subunit I, which is proposed by the TTAG sequence as observed in other lepidopterans. All transfer RNAs (tRNAs) have a typical clover-leaf structure of mitochondrial tRNA, except for tRNA ^Ser (AGN) , the DHU arm of which could not form a stable stem-loop structure. Two aligned sequence blocks with a length of more than 50 bp and 90% of the sequence identity were identified in the A+T-rich region of the Saturniidae and Bombycoidae species.     

Origin of replication in circular prokaryotic chromosomes

To predict origins of replication in prokaryotic chromosomes, we analyse the leading and lagging strands of 200 chromosomes for differences in oligomer composition and show that these correlate strongly with taxonomic grouping, lifestyle and molecular details of the replication process. While all bacteria have a preference for Gs over Cs on the leading strand, we discover that the direction of the A/T skew is determined by the polymerase-alpha subunit that replicates the leading strand. The strength of the strand bias varies greatly between both phyla and environments and appears to correlate with growth rate. Finally we observe much greater diversity of skew among archaea than among bacteria. We have developed a program that accurately locates the origins of replication by measuring the differences between leading and lagging strand of all oligonucleotides up to 8 bp in length. The program and results for all publicly available genomes are available from http://www.cbs.dtu.dk/services/GenomeAtlas/suppl/origin.     

Assessing Substitution Variation Across Sites in Grass Chloroplast DNA

We assess the similarity of base substitution processes, described by empirically derived 4 × 4 matrices, using chi-square homogeneity tests. Such significance analyses allow us to assess variation in sequence evolution across sites and we apply them to matrices derived from noncoding sites in different contexts in grass chloroplast DNA. We show that there is statistically significant variation in rates and patterns of mutation among noncoding sites in different contexts and then demonstrate a similar and significant influence of context on substitutions at fourfold degenerate sites of coding regions from grass chloroplast DNA. These results show that context has the same general effect on substitution bias in coding and noncoding DNA: the A+T content of flanking bases is correlated with rate of substitution, transition bias, and GC → AT pressure, while the number of flanking pyrimidines on a single strand is correlated with a mutational bias, or skew, toward pyrimidines. Despite the similarity in general trends, however, when we compare coding and noncoding matrices we find that there is a statistically significant difference between them even when we control for context. Most noticeably, fourfold degenerate sites in coding sequences are undergoing substitution at a higher rate and there are also significant differences in the relationship between pyrimidines skew and the number of flanking pyrimidines. Possible reasons for the differences between coding and noncoding sites are discussed. Furthermore, our analysis illustrates a simple statistical way for comparing substitution processes across sites allowing us to better study variation in evolutionary processes across a genome. [Reviewing Editor: Dr. Martin Kreitman]     

黄脸油葫芦线粒体基因组:一种新的基因排列方式

采用长距PCR扩增及保守引物步移法测定并注释了黄脸油葫芦（Teleogryllus emma）线粒体基因组全序列。结果表明,黄脸油葫芦的线粒体基因组全长15 660 bp,A+T含量为73.1%。谷氨酸、色氨酸及天冬酰胺的转运RNA基因由N链编码,形成了直翅目中的第三种基因排列顺序,其余结构与其它螽亚目昆虫的线粒体结构一致。基因间隔序列共计73 bp,间隔长度从1—24 bp不等;有14对基因间存在共54 bp重叠,重叠碱基数在1—11 bp之间。13个蛋白质编码基因中12个基因（除COⅠ基因外）的起始密码为标准的ATN组成,COI基因的起始密码子为TTA。有10个基因在基因3＇端能找到完全的TAA或TAG终止密码子,而有三个基因（COII,ND5和ND4）终止密码子为不完整的T。除tRNASer（AGN）外,其余21个tRNA基因的二级结构均属典型的三叶草结构。黄脸油葫芦940bp的A+T富集区中存在一个被认为与复制起始有关的保守的二级结构,该结构不仅存在于直翅目昆虫中,而且也存在于双翅目、鳞翅目和膜翅目中,但是未见于昆虫纲的早期分化类群——弹尾目中。