Estimation of DNA Sequence Context-dependent Mutation Rates Using Primate Genomic Sequences

It is understood that DNA and amino acid substitution rates are highly sequence context-dependent, e.g., C --> T substitutions in vertebrates may occur much more frequently at CpG sites and that cysteine substitution rates may depend on support of the context for participation in a disulfide bond. Furthermore, many applications rely on quantitative models of nucleotide or amino acid substitution, including phylogenetic inference and identification of amino acid sequence positions involved in functional specificity. We describe quantification of the context dependence of nucleotide substitution rates using baboon, chimpanzee, and human genomic sequence data generated by the NISC Comparative Sequencing Program. Relative mutation rates are reported for the 96 classes of mutations of the form 5' alphabetagamma 3' --> 5' alphadeltagamma 3', where alpha, beta, gamma, and delta are nucleotides and beta not equal delta, based on maximum likelihood calculations. Our results confirm that C --> T substitutions are enhanced at CpG sites compared with other transitions, relatively independent of the identity of the preceding nucleotide. While, as expected, transitions generally occur more frequently than transversions, we find that the most frequent transversions involve the C at CpG sites (CpG transversions) and that their rate is comparable to the rate of transitions at non-CpG sites. A four-class model of the rates of context-dependent evolution of primate DNA sequences, CpG transitions > non-CpG transitions approximately CpG transversions > non-CpG transversions, captures qualitative features of the mutation spectrum. We find that despite qualitative similarity of mutation rates among different genomic regions, there are statistically significant differences.     

Correlation of the Helicobacter pylori adherence factor BabA with duodenal ulcer disease in four European countries

Helicobacter pylori strains harboring the vacAs1, cagA and babA2 have been associated with ulcer disease (UD). We compared the prevalence of these different genotypes and adhesive properties in H. pylori infected patients with UD in four European countries. Genomic DNA was isolated from 314 H. pylori strains: Germany (GER; n=92), Sweden (SWE, n=74), Portugal (POR, n=91) and Finland (FIN, n=57). The frequencies of babA2 genotype varied from 35% to 60%. Triple-positive strains (vacAs1+, cagA+ and babA2+) were significantly associated with UD in GER and POR and were closely correlated with UD in FIN, but not in SWE. Classification as triple-positive strains had a higher specificity for detection of UD in GER, POR and FIN than type1 or cagA+ strains. In vitro adhesion assays revealed that Swedish strains showed high adhesion properties and were thus correlated with the diagnosis of UD, although PCR detected the babA2 gene at lower frequencies and failed to show a correlation with UD. This finding appears to reflect allelic variations of the babA2 gene in SWE, although adhesive properties of the strains are retained.     

Identification of a new segment involved in cagA 3' region variation of Helicobacter pylori

The cagA 3' region shows marked variation among Helicobacter pylori strains. Two segments of 102 bp and 57 bp are reportedly responsible for this variation. We analysed the cagA 3' region in 70 H. pylori strains using polymerase chain reaction and sequencing. We found that another segment, namely beta segment, was also involved in the variation of this region. The beta segment was 105 bp long and located between the aforementioned two segments. Six genotypes were identified based on the structure of the cagA 3' region. No relationship was found between these genotypes and the clinical outcomes or vacA genotypes. The numbers of tyrosine phosphorylation sites within the cagA 3' region varied among strains, but this was not related to the cagA genotypes. Our data suggest that the cagA 3' region is significantly variable. It appears that the variation of the cagA 3' region might contribute to the modification of virulence.     

Similar relative mutation rates in the three genetic compartments of Mesostigma and Chlamydomonas

Levels of nucleotide substitution at silent sites in organelle versus nuclear DNAs have been used to estimate relative mutation rates among these compartments and explain lineage-specific features of genome evolution. Synonymous substitution divergence values in animals suggest that the rate of mutation in the mitochondrial DNA is 10-50 times higher than that of the nuclear DNA, whereas overall data for most seed plants support relative mutation rates in mitochondrial, plastid, and nuclear DNAs of 1:3:10. Little is known about relative mutation rates in green algae, as substitution rate data is limited to only the mitochondrial and nuclear genomes of the chlorophyte Chlamydomonas. Here, we measure silent-site substitution rates in the plastid DNA of Chlamydomonas and the three genetic compartments of the streptophyte green alga Mesostigma. In contrast to the situation in animals and land plants, our results support similar relative mutation rates among the three genetic compartments of both Chlamydomonas and Mesostigma. These data are discussed in relation to published intra-species genetic diversity data for the three genetic compartments of Chlamydomonas and are ultimately used to address contemporary hypotheses on the organelle genome evolution. To guide future work, we describe evolutionary divergence data of all publically available Mesostigma viride strains and identify, for the first time, three distinct lineages of Mesostigma.     

High frequency of spontaneous Minute mutations detected in the F1 progeny of interstrain matings between a recombination-defective mei-9L1 and the y w m f/sc8(y+) Y BS; dp strain of Drosophila melanogaster

The yield of spontaneous Minute mutations was recorded in the F1 progeny of interstrain (reciprocal) and intrastrain matings between a recombination- and excision repair-defective mei-9L1 (mei-9) strain and the y w m f/sc8(y+) Y BS; dp (ywmf-2) strain of Drosophila melanogaster. As a comparison, interstrain matings between a postreplication repair-defective st mus(3)302D1 (mus(3)) strain and the ywmf-2 strain were also studied for Minute mutations. The results show that: (1) a strikingly high frequency of Minute mutations is observed in the progeny of mei-9 female X ywmf-2 male crosses, but not in that of ywmf-2 females X mei-9 males; (2) no such difference exists in the progeny of intrastrain matings; and (3) there exists no marked inequality of Minute frequencies in the progeny of reciprocal crosses of mus(3) and ywmf-2 strains.     

Whole-Genome Sequencing of Theileria parva Strains Provides Insight into Parasite Migration and Diversification in the African Continent

The disease caused by the apicomplexan protozoan parasite Theileria parva, known as East Coast fever or Corridor disease, is one of the most serious cattle diseases in Eastern, Central, and Southern Africa. We performed whole-genome sequencing of nine T. parva strains, including one of the vaccine strains (Kiambu 5), field isolates from Zambia, Uganda, Tanzania, or Rwanda, and two buffalo-derived strains. Comparison with the reference Muguga genome sequence revealed 34 814–121 545 single nucleotide polymorphisms (SNPs) that were more abundant in buffalo-derived strains. High-resolution phylogenetic trees were constructed with selected informative SNPs that allowed the investigation of possible complex recombination events among ancestors of the extant strains. We further analysed the dN/dS ratio (non-synonymous substitutions per non-synonymous site divided by synonymous substitutions per synonymous site) for 4011 coding genes to estimate potential selective pressure. Genes under possible positive selection were identified that may, in turn, assist in the identification of immunogenic proteins or vaccine candidates. This study elucidated the phylogeny of T. parva strains based on genome-wide SNPs analysis with prediction of possible past recombination events, providing insight into the migration, diversification, and evolution of this parasite species in the African continent.     

The Evolutionary Dynamics of Bluetongue Virus

Bluetongue virus (BTV) is a midge-borne member of the genus Orbivirus that causes an eponymous debilitating livestock disease of great agricultural impact and which has expanded into Europe in recent decades. Reassortment among the ten segments comprising the double-stranded (ds) RNA genome of BTV has played an important role in generating the epidemic strains of this virus in Europe. In this study, we investigated the dynamics of BTV genome segment evolution utilizing time-structured data sets of complete sequences from four segments, totalling 290 sequences largely sampled from ruminant hosts. Our analysis revealed that BTV genome segments generally evolve under strong purifying selection and at substitution rates that are generally lower (mean rates of ~0.5–7 × 10⁻⁴ nucleotide substitutions per site, per year) than vector-borne positive-sense viruses with single-strand (ss) RNA genomes. These also represent the most robust estimates of the nucleotide substitution rate in a dsRNA virus generated to date. Additionally, we determined that patterns of geographic structure and times to most recent common ancestor differ substantially between each segment, including a relatively recent origin for the diversity of segment 10 within the past millennium. Together, these findings demonstrate the effect of reassortment to decouple the evolutionary dynamics of BTV genome segments.     

Molecular Evolution of the Duplicated Amy Locus in the Drosophila Melanogaster Species Subgroup: Concerted Evolution Only in the Coding Region and an Excess of Nonsynonymous Substitutions in Speciation

From the analysis of restriction maps of the Amy region in eight sibling species belonging to the Drosophila melanogaster species subgroup, we herein show that the patterns of duplication of the Amy gene are almost the same in all species. This indicates that duplication occurred before speciation within this species subgroup. From the nucleotide sequence data, we show a strong within-species similarity between the duplicated loci in the Amy coding region. This is in contrast to a strong similarity in the 5' and 3' flanking regions within each locus (proximal or distal) throughout the species subgroup. This means that concerted evolution occurred only in the Amy coding region and that differentiated evolution between the duplication occurred in the flanking regions. Moreover, when comparing the species, we also found a significant excess of nonsynonymous substitutions. In particular, all the fixed substitutions specific to D. erecta were found to be nonsynonymous. We thus conclude that adaptive protein evolution occurred in the lineage of D. erecta that is a ``specialist' species for host plants and probably also occurs in the process of speciation in general.     

Molecular evolution of the histone 3 multigene family in the Drosophila melanogaster species subgroup

Molecular evolution of the histone multigene family was studied by cloning and sequencing regions of the histone 3 gene in the Drosophila melanogaster species subgroup. Analysis of the nucleotide substitution pattern showed that in the coding region synonymous changes occurred more frequently to A or T in contrast to the GC-rich base composition, while in the 3' region the nucleotide substitutions were most likely in equilibrium. These results suggested that the base composition at the third codon position of the H3 gene, i.e., codon usage, has been changing to A or T in the Drosophila melanogaster species subgroup.     

Ectopic Gene Conversions in Four <Emphasis Type="BoldItalic">Escherichia coli</Emphasis> Genomes: Increased Recombination in Pathogenic Strains

We characterized the ectopic gene conversions in the genomes of the K-12 MG1655, O157:H7 Sakai, O157:H7 EDL933, and CFT073 strains of E coli. Compared to the three pathogenic strains, the K-12 strain has a much smaller number of gene families, its gene families contain fewer genes, and gene conversions are less frequent. Whereas the three pathogenic strains have gene conversions covering hundreds of nucleotides when their flanking regions have as little as 50% similarity, flanking region similarity of at least 94% on both sides of the converted region is required to observe conversions of more than 87 nucleotides in the K-12 strain. Recombination is therefore more frequent and requires less sequence similarity in the three pathogenic strains than in K-12. This higher recombination level might be due to mutations in some of their mismatch-repair genes. In contrast with the gene conversions present in the yeast genome, the gene conversions found in the E. coli genomes do not occur more frequently between duplicated genes that are close to one another than between duplicated genes that are far apart and are randomly distributed along the length of the genes. In E. coli, gene conversions are not more frequent near the origin of replication. However, they do occur more frequently near the terminus of replication of the Sakai genome, where multigene family members are more abundant. This suggests that, in E. coli, gene conversions occur randomly between genes located in different chromosomal locations or located on different copies of the multiple chromosomes found in E. coli cells.     

Genetic battle between Helicobacter pylori and humans. The mechanism underlying homologous recombination in bacteria,which can infect human cells

Helicobacter pylori is a gram-negative pathogenic bacterium that colonises the human stomach. The chronic infection it causes results in peptic ulcers and gastric cancers. H. pylori can easily establish a chronic infection even if the immune system attacks this pathogen with oxidative stress agents and immunoglobulins. This is attributed to bacterial defence mechanisms against these stresses. As a defence mechanism against oxidative stresses, in bacterial genomes, homologous recombination can act as a repair pathway of DNA's double-strand breaks (DSBs). Moreover, homologous recombination is also involved in the antigenic variation in H. pylori. Gene conversion alters genomic structures of babA and babB (encoding outer membrane proteins), resulting in escape from immunoglobulin attacks. Thus, homologous recombination in bacteria plays an important role in the maintenance of a chronic infection. In addition, H. pylori infection causes DSBs in human cells. Homologous recombination is also involved in the repair of DSBs in human cells. In this review, we describe the roles of homologous recombination with an emphasis on the maintenance of a chronic infection.     

Heterogeneous Nature and Distribution of Interruptions in Dinucleotides May Indicate the Existence of Biased Substitutions Underlying Microsatellite Evolution

Some aspects of microsatellite evolution, such as the role of base substitutions, are far from being fully understood. To examine the significance of base substitutions underlying the evolution of microsatellites we explored the nature and the distribution of interruptions in dinucleotide repeats from the human genome. The frequencies that we inferred in the repetitive sequences were statistically different from the frequencies observed in other noncoding sequences. Additionally, we detected that the interruptions tended to be towards the ends of the microsatellites and 5'-3' asymmetry. In all the estimates nucleotides forming the same repetitive motif seem to be affected by different base substitution rates in AC and AG. This tendency itself could generate patterning and similarity in flanking sequences and reconcile these phenomena with the high mutation rate found in flanking sequences without invoking convergent evolution. Nevertheless, our data suggest that there is a regional bias in the substitution pattern of microsatellites. The accumulation of random substitutions alone cannot explain the heterogeneity and the asymmetry of interruptions found in this study or the relative frequency of different compound microsatellites in the human genome. Therefore, we cannot rule out the possibility of a mutational bias leading to convergent or parallel evolution in flanking sequences.     

Gene conversions in genes encoding outer-membrane proteins in H. pylori and C. pneumoniae

Helicobacter pylori and Chlamydia pneumoniae are both pathogenic to humans. Their genomes have recently been completed, allowing detailed study of their evolution and organization. Here we describe an evolutionary analysis of the H. pylori and C. pneumoniae genes that encode their outer-membrane proteins. By comparing complete genome sequences of two H. pylori strains and two C. pneumoniae strains, we identify multiple independent conversions among these genes. Such recombination events might provide a selective advantage for these bacterial pathogens.     

Zea ribosomal repeat evolution and substitution patterns

Zea and Tripsacum nuclear ribosomal internal transcribed spacer (ITS) sequences were used to evaluate patterns of concerted evolution, rates of substitutions, patterns of methylation-induced deamination, and structural constraints of the ITS. ITS pseudogenes were identified by their phylogenetic position, differences in nucleotide composition, extensive deamination at ancestral methylation sites, and substitutions resulting in low-stability secondary RNA structures. Selection was important in shaping the kinds of polymorphisms and substitutions observed in the ITS. ITS substitution rates were significantly different among the Zea taxa. Deamination of cytosines at methylation sites was a potent mutation source, but selection appeared to maintain high methylation site density throughout the ribosomal repeat except for the gene promoter. Nucleotide divergence statistics identified selectively constrained regions at the 5' ends of the ITS1 and ITS2.      

PCR-based molecular characterization and insilico analysis of food-borne trematode parasites Paragonimus westermani, Fasciolopsis buski and Fasciola gigantica from Northeast India using ITS2 rDNA

In early 2009, new swine-origin influenza A (H1N1) virus emerged in Mexico and the United States. The emerging influenza virus had made global influenza pandemic for nearly one year. To every emerging pathogen, exploring the origin sources is vital for viral control and clearance. Influenza virus is different from other virus in that it has 8 segments, making the segment reassortment a main drive in virus evolution. In exploring reassortment evolution origins of a newly emerging influenza virus, integrated comparing of the origin sources of all the segments is necessary. If some segments have high homologous with one parental strain, lower homologous with another parental strain, while other segments are reverse, can we proposed that this emerging influenza virus may re-assort from the two parental strains. Here we try to explore the multilevel reassortment evolution origins of 2009 H1N1 influenza virus using this method. By further validating the fidelity of this strategy, this method might be useful in judging the reassortment origins of newly emerging influenza virus.     

The Determinants of Stability and Folding in Evolutionarily Diverged Cytochromes c

Cytochrome c has served as a paradigm for the study of protein stability, folding, and molecular evolution, but it remains unclear how these aspects of the protein are related. For example, while the bovine and equine cytochromes c are known to have different stabilities, and possibly different folding mechanisms, it is not known how these differences arise from just three amino acid substitutions introduced during divergence. Using site-selectively incorporated carbon-deuterium bonds, we show that like the equine protein, bovine cytochrome c is induced to unfold by guanidine hydrochloride via a stepwise mechanism, but it does not populate an intermediate as is observed with the equine protein. The increased stability also results in more similar free energies of unfolding observed at different sites within the protein, giving the appearance of a more concerted mechanism. Furthermore, we show that the differences in stability and folding appear to result from a single amino acid substitution that stabilizes a helix by allowing for increased solvation of its N-terminus.     

Mosaic Structure and Molecular Evolution of the Leukotoxin Operon (lktCABD) in Mannheimia (Pasteurella) haemolytica, Mannheimia glucosida, and Pasteurella trehalosi

The mosaic structure and molecular evolution of the leukotoxin operon (lktCABD) was investigated by nucleotide sequence comparison of the lktC, lktB, and lktD genes in 23 Mannheimia (Pasteurella) haemolytica, 6 Mannheimia glucosida, and 4 Pasteurella trehalosi strains. Sequence variation in the lktA gene has been described previously (R. L. Davies et al., J. Bacteriol. 183:1394-1404, 2001). The leukotoxin operon of M. haemolytica has a complex mosaic structure and has been derived by extensive inter- and intraspecies horizontal DNA transfer and intragenic recombination events. However, the pattern of recombination varies throughout the operon and among the different evolutionary lineages of M. haemolytica. The lktA and lktB genes have the most complex mosaic structures with segments derived from up to four different sources, including M. glucosida and P. trehalosi. In contrast, the lktD gene is highly conserved in M. haemolytica. The lktC, lktA, and lktB genes of strains representing the major ovine lineages contain recombinant segments derived from bovine or bovine-like serotype A2 strains. These findings support the previous conclusion that host switching of bovine A2 strains from cattle to sheep has played a major role in the evolution of the leukotoxin operon in ovine strains of M. haemolytica. Homologous segments of donor and recipient alleles are identical, or nearly identical, indicating that the recombinational exchanges occurred relatively recent in evolutionary terms. The 5' and 3' ends of the operon are highly conserved in M. haemolytica, which suggests that multiple horizontal exchanges of the complete operon have occurred by a common mechanism such as transduction. Although the lktA and lktB genes both have complex mosaic structures and high nucleotide substitution rates, the amino acid diversity of LktB is significantly lower than that of LktA due to a higher degree of evolutionary constraint against amino acid replacement. The recombinational exchanges within the leukotoxin operon have had greatest effect on LktA and probably provide an adaptive advantage against the host antibody response by generating novel antigenic variation at surface-exposed sites.     

An RNA conformational shift in recent H5N1 influenza A viruses

Recent outbreaks of avian influenza are being caused by unusually virulent H5N1 strains. It is unknown what makes these recent H5N1 strains more aggressive than previously circulating strains. Here, we have compared more than 3000 RNA sequences of segment 8 of type A influenza viruses and found a unique single nucleotide substitution typically associated with recent H5N1 strains. By phylogenetic analysis, biochemical and biophysical experiments, we demonstrate that this substitution dramatically affects the equilibrium between a hairpin and a pseudoknot conformation near the 3' splice-site of the NS gene. This conformational shift may have consequences for splicing regulation of segment 8 mRNA. Our data suggest that besides changes at the protein level, changes in RNA secondary structure should be seriously considered when attempting to explain influenza virus evolution. Supplementary information: Supplementary data are available at Bioinformatics online.     

Mechanisms of molecular evolution

Both drift and selection are important for nucleotide substitutions in evolution. The nearly neutral theory was developed to clarify the effects of these processes. In this article, the nearly neutral theory is presented with special reference to the nature of weak selection. The mean selection coefficient is negative, and the variance is dependent on the environmental diversity. Some facts relating to the theory are reviewed. As well as nucleotide substitutions, illegitimate recombination events such as duplications, deletions and gene conversions leave indelible marks on molecular evolution. Gene duplication and conversion are sources of the evolution of new gene functions. Positive selection is necessary for the evolution of novel functions. However, many examples of current gene families suggest that both drift and selection are at work on their evolution.