Nonsynonymous Nucleotide Diversity Research Articles

More effective purifying selection on RNA viruses than in DNA viruses

Analysis of the pattern of nucleotide diversity in 222 independent viral sequence datasets showed the prevalence of purifying selection. In spite of the higher mutation rate of RNA viruses, our analyses revealed stronger evidence of the action of purifying selection in RNA viruses than in DNA viruses. The ratio of nonsynonymous to synonymous nucleotide diversity was significantly lower in RNA viruses than in DNA viruses, indicating that nonsynonymous mutations have been removed at a greater rate (relative to the mutation rate) in the former than in the latter. Moreover, statistics that measure the occurrence of rare polymorphisms revealed significantly a greater excess of rare nonsynonymous polymorphisms in RNA viruses than in DNA viruses but no difference with respect to synonymous polymorphisms. Since rare nonsynonymous polymorphisms are likely to be undergoing the effects of purifying selection acting to eliminate them, this result implies a stronger signature of ongoing purifying selection in RNA viruses than in DNA viruses. Across datasets from both DNA viruses and RNA viruses, we found a negatively allometric relationship between nonsynonymous and synonymous nucleotide diversity; in other words, nonsynonymous nucleotide diversity increased with synonymous nucleotide diversity at a less than linear rate. These findings are most easily explained by the occurrence of slightly deleterious mutations. The fact that the negative allometry was more pronounced in RNA viruses than in DNA viruses provided additional evidence that purifying selection is more effective in the former than in the latter.

Micro-scale signature of purifying selection in Marburg virus genomes

In the seven protein-coding genes in the Marburg virus (MARV) genome, the synonymous nucleotide diversity substantially exceeded the nonsynonymous nucleotide diversity, indicating strong purifying selection. Likewise, there was evidence of purifying selection on 5′UTR and 3′UTR, where nucleotide diversity ( π) was significantly less than π S in the coding regions. Nonsynonymous polymorphic sites showed significantly reduced mean gene diversity in comparison to other polymorphic sites, indicating that purifying selection at certain slightly deleterious nonsynonymous polymorphisms is ongoing. Moreover, nonsynonymous polymorphic sites showed significantly reduced gene diversity in comparison to adjacent synonymous sites, even though the vast majority of such adjacent synonymous sites were in the same codon or an adjacent codon. Thus purifying selection, in conjunction with recombination and/or backward mutation, can act to break up linkage relationships at a micro-scale in the MARV genome. The ability of purifying selection to break up linkage between synonymous and nonsynonymous polymorphisms on such a fine scale has not been reported in any other genome.

Evolution and Differentiation of MSHR Gene in Different Species

Coat color offers some prospects for evolutionary studies due to its large amount of presumably adaptive coat color variation and conserved genetic mechanisms of generating different coat colors in different species. Melanocyte-stimulating hormone receptor (MSHR) gene is responsible for intraspecific and interspecific color variation in mammals and birds. A total number of 206 MSHR gene sequences belonging to 84 species, 58 genera, and 20 families were analyzed to investigate its evolution and differentiation in different species. Most of the species have 954 bp and stop codon TGA. Species in Callithrix and Callimico have a stop codon mutation from TGA to TGG and elongate 81 bp with TAG as stop codon. Species in Phasianidae, Fringillidae, and Lemuridae also use TAG as stop codon. The Sus scrofa had an insertion of AACCAGACC encoding Asn-Gln-Thr from 85 to 93 bp. In Bovidae, a brown strain of cow with 966 bp due to the 12-bp duplication of GGCATTGCCCGG from 670 to 681 bp encoding for Gly-Ile-Ala-Arg was found. Teiidae has the smallest number of total mutations (6), silent mutations (3), nonsynonymous mutations (3), average number of nucleotide differences (1.519), synonymous nucleotide diversity (pi(s) = 0.0030), and nonsynonymous nucleotide diversity (pi(a) = 0.0029), and Hominidae, Lemuridae, Canidae, and Teiidae have higher ratio of pi(a)/pi(s) (0.537-0.973). The reconstructed phylogenetic tree of MSHR gene of families is basically consistent with the taxonomy of National Center for Biotechnology Information.

Excess of nonsynonymous polymorphism at Acph-1 in different gene arrangements of Drosophila subobscura.

Nucleotide variation in the Acph-1 gene region was analyzed in a natural population of Drosophila subobscura from Bizerte (Tunisia). The lines studied differed in their gene arrangement for segment I of the O chromosome: 21 lines were O3+4+8, 21 were O3+4+23, and 3 were O3+4. According to chromosomal phylogenies, O3+4 is a central arrangement from which O3+4+8 and O3+4+23 originated. Strong genetic differentiation at Acph-1 was detected among the different arrangements, which is reflected in strong linkage disequilibrium between the variants at informative polymorphic sites and the type of arrangement. Estimates of silent nucleotide diversity are slightly lower within O3+4+23 (pisilent = 0.0166) than within O3+4+8 (pisilent = 0.0228) or O3+4 (pisilent = 0.0234). In contrast, nonsynonymous nucleotide diversity estimates (around 0.1%) are similar in the three arrangements. Most nonsynonymous rare variants are singletons, which results in highly significant Tajima's neutrality tests within the young O3+4+8 and O3+4+23 arrangements. This test is not significant for nonsynonymous mutations within a large Spanish O3+4 sample. In addition, a significant and marginally significant excess of nonsynonymous polymorphism was detected by the McDonald and Kreitman test within O3+4+23 and O3+4+8, respectively. This excess results in a rather high neutrality index (NI = 5.25) when both arrangements are jointly analyzed, in contrast to its value within the old O3+4 arrangement (NI = 1.74). The pattern of variation at Acph-1 within the young arrangements is unusual for nuclear genes and has the same characteristics previously detected in most genes of the mitochondrial genome. Assuming that most nonsynonymous mutations at Acph-1 are under weak negative selection, a smaller effective size of the young arrangements relative to O3+4 might explain the observed results. The relatively low frequency of O3+4+8 and O3+4+23 in the distribution area of D. subobscura, the more recent origin of these arrangements relative to O3+4 and the suppression of recombination in heterokaryotypes might contribute to the relatively small effective size of the young arrangements. Therefore, present results indicate that the differences in effective size and recombination caused by chromosomal arrangements are modulating nonsynonymous variation at Acph-1.

DNA sequence variation in BpMADS2 gene in two populations of Betula pendula

The PISTILLATA (PI) homologue, BpMADS2, was isolated from silver birch (Betula pendula Roth) and used to study nucleotide polymorphism. Two regions (together about 2450 bp) comprising mainly untranslated sequences were sequenced from 10 individuals from each of two populations in Finland. The nucleotide polymorphism was low in the BpMADS2 locus, especially in the coding region. The synonymous site overall nucleotide diversity (pis) was 0.0043 and the nonsynonymous nucleotide diversity (pia) was only 0.000052. For the whole region, the pi values for the two populations were 0.0039 and 0.0045, and for the coding regions, the pi values were only 0 and 0.00066 (for the corresponding coding regions of Arabidopsis thaliana PI world-wide pi was 0.0021). Estimates of pi or theta did not differ significantly between the two populations, and the two populations were not diverged from each other. Two classes of BpMADS2 alleles were present in both populations, suggesting that this gene exhibits allelic dimorphism. In addition to the nucleotide site variation, two microsatellites were also associated within the haplotypes. This allelic dimorphism might be the result of postglacial re-colonization partly from northwestern, partly from southeastern/eastern refugia. The sequence comparison detected five recombination events in the regions studied. The large number of microsatellites in all of the three introns studied suggests that BpMADS2 is a hotspot for microsatellite formation.