Sequence variation in the bovine and ovine PRNP genes

A resequencing approach was adopted to identify sequence variants in the PRNP gene that may affect susceptibility or resistance to bovine spongiform encephalopathy. The entire PRNP gene (>21 kb) was sequenced from 26 chromosomes from a group of Holstein-Friesian cows, as well as exon 3 of PRNP (>4 kb) from a further 24 chromosomes from six diverse breeds. We identified 51 variant sequences of which 42 were single nucleotide polymorphisms and nine were insertion/deletion (indel) events. The study was extended to exon 3 of the sheep PRNP gene where 23 sequence variants were observed, four of which were indels. The level of nucleotide diversity in the complete bovine PRNP gene was pi = 0.00079, which is similar to that found at the bovine T-cell receptor alpha delta joining region (pi = 0.00077), but somewhat less than that observed for the bovine leptin (pi = 0.00265). Sequence variation within exon 3 of PRNP in both cattle (pi = 0.00102) and sheep (pi = 0.00171) was greater than that for the complete PRNP gene, with sheep showing greater sequence variation in exon 3 than cattle. The level of sequence variation reported here is greater than previously thought for the bovine PRNP gene in cattle. This study highlights the contribution that recombination plays in increasing allelic diversity in this species.     

Detection of Polymorphism and Sequence Characterization of Toll-Like Receptor 7 Gene of Indian Goat Revealing Close Relationship Between Ruminant Species

In this study, approximately 3.4 kb nucleotide sequence of caprine TLR7 (Toll-like receptor 7) gene was generated from twelve different Indian goat breeds belonging to different geographical regions. Goat TLR7 gene ORF (Open Reading Frame) was found to be 3141 nucleotides long coding for 1046 amino acids similar to sheep. The sequence analysis at nucleotide level revealed goat TLR7 having 99.5% homology with sheep, followed by other livestock species. Simple Modular Architecture Research Tool (SMART) was used for the structural analysis of goat TLR7 that showed the presence of 22 leucine rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domains. TIR domain, when compared, was found to be similar in ruminant species, goat, sheep, cattle, and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to sheep followed by cattle and buffalo. A total of 22 polymorphic sites were observed in TLR7 gene of 24 goats representing 12 different breeds, out of which 19 were present within the coding region and three in 3'UTR. Out of the seven nonsynonymous SNPs, two were in ectodomains and one in TIR domain. Overall our results indicate substantial variation within goat TLR7 gene, which could be exploited for association with disease susceptibility.     

Coalescent processes and relaxation of selective constraints leading to contrasting genetic diversity at paralogs AtHVA22d and AtHVA22e in Arabidopsis thaliana

Duplicate loci offer a very powerful system for understanding the complicated genome structure and adaptive evolution of a gene family. In this study, the genetic variation at paralogs AtHVA22d and AtHVA22e, members of an ABA- and stress-inducible gene family, is examined in the selfing Arabidopsis thaliana. Population genetic analysis indicates contrasting levels of nucleotide diversity at overall exon sequence and nonsynonymous sites between AtHVA22d (pi = 0.00337, pi(rep) = 0.00158) and AtHVA22e (pi = 0.00054, pi(rep) = 0.00023). The fact of Ka/Ks ratios significantly less than 1 in all sequences indicates that both genes are functional and subjected to purifying selection. In addition, rooted at barley HVA22, accelerated evolution is detected at replacement changes in the AtHVA22d locus, indicating relaxation of purifying selection after gene duplication. However, relative rate tests reveal no deviation from the neutrality at synonymous sites between the two paralogs. Based on clock-like evolution, the rate of synonymous substitution is estimated at 1.83 x 10(-9) substitutions per site per year; and the divergence of the two paralogs is traced to 90 MYA, coinciding with a period of the diversification of angiosperms. Given no codon usage bias in both genes, natural selection alone cannot account for the 6.4-fold differences in the nucleotide variation at synonymous sites between the two paralogs. Random processes resulting in different coalescence times, 3.65 MYA at AtHVA22d vs. 1.20 MYA at AtHVA22e, may have predominantly contributed to the evident differences of the genetic diversity. Partially nonoverlapping modes of expression between the two functional paralogs suggest a subfunctionalization hypothesis for explaining the fates of duplicate loci.     

Detection of polymorphism and sequence characterization of toll-like receptor 7 gene of Indian goat revealing close relationship between ruminant species

In this study, approximately 3.4?kb nucleotide sequence of caprine TLR7 (Toll-like receptor 7) gene was generated from twelve different Indian goat breeds belonging to different geographical regions. Goat TLR7 gene ORF (Open Reading Frame) was found to be 3141 nucleotides long coding for 1046?amino acids similar to sheep. The sequence analysis at nucleotide level revealed goat TLR7 having 99.5% homology with sheep, followed by other livestock species. Simple Modular Architecture Research Tool (SMART) was used for the structural analysis of goat TLR7 that showed the presence of 22 leucine rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domains. TIR domain, when compared, was found to be similar in ruminant species, goat, sheep, cattle, and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to sheep followed by cattle and buffalo. A total of 22 polymorphic sites were observed in TLR7 gene of 24 goats representing 12 different breeds, out of which 19 were present within the coding region and three in 3'UTR. Out of the seven nonsynonymous SNPs, two were in ectodomains and one in TIR domain. Overall our results indicate substantial variation within goat TLR7 gene, which could be exploited for association with disease susceptibility.     

Selection in favor of nucleotides G and C diversifies evolution rates and levels of polymorphism at mammalian synonymous sites

The impact of synonymous nucleotide substitutions on fitness in mammals remains controversial. Despite some indications of selective constraint, synonymous sites are often assumed to be neutral, and the rate of their evolution is used as a proxy for mutation rate. We subdivide all sites into four classes in terms of the mutable CpG context, nonCpG, postC, preG, and postCpreG, and compare four-fold synonymous sites and intron sites residing outside transposable elements. The distribution of the rate of evolution across all synonymous sites is trimodal. Rate of evolution at nonCpG synonymous sites, not preceded by C and not followed by G, is approximately 10% below that at such intron sites. In contrast, rate of evolution at postCpreG synonymous sites is approximately 30% above that at such intron sites. Finally, synonymous and intron postC and preG sites evolve at similar rates. The relationship between the levels of polymorphism at the corresponding synonymous and intron sites is very similar to that between their rates of evolution. Within every class, synonymous sites are occupied by G or C much more often than intron sites, whose nucleotide composition is consistent with neutral mutation-drift equilibrium. These patterns suggest that synonymous sites are under weak selection in favor of G and C, with the average coefficient s approximately 0.25/Ne approximately 10(-5), where Ne is the effective population size. Such selection decelerates evolution and reduces variability at sites with symmetric mutation, but has the opposite effects at sites where the favored nucleotides are more mutable. The amino-acid composition of proteins dictates that many synonymous sites are CpGprone, which causes them, on average, to evolve faster and to be more polymorphic than intron sites. An average genotype carries approximately 10(7) suboptimal nucleotides at synonymous sites, implying synergistic epistasis in selection against them.     

Gene conversion as a source of nucleotide diversity in Plasmodium falciparum

Examination of polymorphisms in the Plasmodium falciparum gene for falcipain 2 revealed that this gene is one of two paralogs separated by 10.8 kb in chromosome 11. We designate the annotated gene denoted chr11.gen_424 as encoding falcipain 2A and the annotated gene denoted chr11.gen_427 as encoding falcipain 2B. The paralogs are 96% identical at the nucleotide level and 93% identical at the amino acid level. The consensus sequences differ in 31/309 synonymous sites and 45/1140 nonsynonymous sites, including three amino acid replacements (V393I, A400P, and Q414E) that are near the catalytic site and that may affect substrate affinity or specificity. In six reference isolates, among 36 synonymous sites and 46 nonsynonymous sites that are polymorphic in the gene for falcipain 2A, falcipain 2B, or both, significant spatial clustering is observed. All but one of the polymorphisms appear to result from gene conversion between the paralogs. The estimated rate of gene conversion between the paralogs may be as many as 1,400 to 1,700 times greater than the rate of mutation. Owing to gene conversion, one of the falcipain 2A alleles is more similar to the falcipain 2B alleles than it is to other falcipain 2A alleles. Divergence among the synonymous sites suggests that the paralogous genes last shared a common ancestor 15.2 MYA, with a range of 8.8 to 20.6 MYA. During this period, the paralogs have acquired 0.10 synonymous substitutions per synonymous site in the coding region. The 5' and 3' flanking regions differ in 47.7% and 39.8% of the nucleotide sites, respectively. Hence synonymous sites and flanking regions are not conserved in sequence in spite of their high AT content and T skew.     

Synonymous Nucleotide Divergence and Saturation: Effects of Site-Specific Variations in Codon Bias and Mutation Rates

The synonymous divergence between Escherichia coli and Salmonella typhimurium is explained in a model where there is a large variation between mutation rates at different nucleotide sites in the genome. The model is based on the experimental observation that spontaneous mutation rates can vary over several orders of magnitude at different sites in a gene. Such site-specific variation must be taken into account when studying synonymous divergence and will result in an apparent saturation below the level expected from an assumption of uniform rates. Recently, it has been suggested that codon preference in enterobacteria has a very large site-specific variation and that the synonymous divergence between different species, e.g., E. coli and Salmonella, is saturated. In the present communication it is shown that when site-specific variation in mutation rates is introduced, there is no need to invoke assumptions of saturation and a large variability in codon preference. The same rate variation will also bring average mutation rates as estimated from synonymous sequence divergence into numerical agreement with experimental values. Received: 10 July 1998 / Accepted: 20 August 1998     

Evolution of the autosomal chorion cluster inDrosophila. IV. The HawaiianDrosophila: Rapid protein evolution and constancy in the rate of DNA divergence

Summary Autosomal chorion geness18, s15, ands19 are shown to diverge at extremely rapid rates in closely related taxa of HawaiianDrosophila. Their nucleotide divergence rates are at least as fast as those of intergenic regions that are known to evolve more extensively between distantly related species. Their amino acid divergence rates are the fastest known to date. There are two nucleotide replacement substitutions for every synonymous one. The molecular basis for observed length and substitution mutations is analyzed. Length mutations are strongly associated with direct repeats in general, and with tandem repeats in particular, whereas the rate for an average transition is twice that for an average transversion.The DNA sequence of the cluster was used to construct a phylogenetic tree for five taxa of the Hawaiian picture-winged species group ofDrosophila. Assignment of observed base substitutions occurring in various branches of the tree reveals an excess of would-be homoplasies in a centrally localized 1.8-kb segment containing thes15 gene. This observation may be a reflection of ancestral excess polymorphisms in the segment. The chorion cluster appears to evolve at a constant rate regardless of whether the central 1.8-kb segment is included or not in the analysis. Assuming that the time of divergence ofDrosophila grimshawi and theplanitibia subgroup coincides with the emergence of the island of Kauai, the overall rate of base substitution in the cluster is estimated to be 0.8% million years, whereas synonymous sites are substituted at a rate of 1.2%/million years.     

Sequence of specific mitochondrial 16S rRNA gene fragment from Egyptian buffalo is used as a pattern for discrimination between river buffaloes,cattle, sheep and goats

Characterization of molecular markers and the development of better assays for precise and rapid detection of domestic species are always in demand. This is particularly due to recent food scares and the crisis of biodiversity resulting from the huge ongoing illegal traffic of endangered species. The aim of this study was to develop a new and easy method for domestic species identification (river buffalo, cattle, sheep and goat) based on the analysis of a specific mitochondrial nucleotide sequence. For this reason, a specific fragment of Egyptian buffalo mitochondrial 16S rRNA gene (422 bp) was amplified by PCR using two universal primers. The sequence of this specific fragment is completely conserved between all tested Egyptian buffaloes and other river buffaloes in different places in the world. Also, the lengths of the homologous fragments were less by one nucleotide (421 bp) in case of goats and two nucleotides (420 bp) in case of both cattle and sheep. The detection of specific variable sites between investigated species within this fragment was sufficient to identify the biological origin of the samples. This was achieved by alignment between the unknown homologous sequence and the reference sequences deposited in GenBank database (accession numbers, FJ748599–FJ748607). Considering multiple alignment results between 16S rRNA homologous sequences obtained from GenBank database with the reference sequence, it was shown that definite nucleotides are specific for each of the four studied species of the family Bovidae. In addition, other nucleotides are detected which can allow discrimination between two groups of animals belonging to two subfamilies of family Bovidae, Group one (closely related species like cattle and buffalo, Subfamily Bovinae) and Group two (closely related species like sheep and goat, Subfamily Caprinae). This 16S DNA barcode character-based approach could be used to complement cytochrome c oxidase I (COI) in DNA barcoding. Also, it is a good tool for identification of unknown sample belonging to one of the four domestic animal species of family Bovidae quickly and easily.     

Uncorrected nucleotide bias in mtDNA can mimic the effects of positive Darwinian selection

The relative rates of nucleotide substitution at synonymous and nonsynonymous sites within protein-coding regions have been widely used to infer the action of natural selection from comparative sequence data. It is known, however, that mutational and repair biases can affect rates of evolution at both synonymous and nonsynonymous sites. More importantly, it is also known that synonymous sites are particularly prone to the effects of nucleotide bias. This means that nucleotide biases may affect the calculated ratio of substitution rates at synonymous and nonsynonymous sites. Using a large data set of animal mitochondrial sequences, we demonstrate that this is, in fact, the case. Highly biased nucleotide sequences are characterized by significantly elevated dN/dS ratios, but only when the nucleotide frequencies are not taken into account. When the analysis is repeated taking the nucleotide frequencies at each codon position into account, such elevated ratios disappear. These results suggest that the recently reported differences in dN/dS ratios between vertebrate and invertebrate mitochondrial sequences could be explained by variations in mitochondrial nucleotide frequencies rather than the effects of positive Darwinian selection.     

Nucleotide sequence of the Xdh region in Drosophila pseudoobscura and an analysis of the evolution of synonymous codons

The nucleotide sequence of the Xdh region of Drosophila pseudoobscura is presented. The Xdh gene structure and organization are compared with the homologous region in D. melanogaster. This locus is shown to have similar organization in the two species, although an additional intron and three insertion/deletion events are described for the D. pseudoobscura coding region. The encoded proteins are predicted to have very similar charges and hydrophobic/hydrophilic domains even though 11% of the amino acids are different. A gene 5' to Xdh, putative l(3)s12, is suggested from sequence similarity between the species. Synonymous differences at the Xdh locus between the two species are analyzed using a new method described in the preceding paper by Lewontin. This analysis shows that synonymous positions within the Xdh locus are evolving at very different rates, being dependent on level of codon redundancy. A comparison of synonymous divergence between D. melanogaster and D. pseudoobscura in five additional genes reveals variation in the level of synonymous substitution.      

Testing for Selection on Synonymous Sites in Plant Mitochondrial DNA: The Role of Codon Bias and RNA Editing

Since plant mitochondrial genomes exhibit some of the slowest known synonymous substitution rates, it is generally believed that they experience exceptionally low mutation rates. However, the use of synonymous substitution rates to infer mutation rates depends on the implicit assumption that synonymous sites are evolving neutrally (or nearly so). To assess the validity of this assumption in plant mitochondrial genomes, we examined coding sequence for footprints of selection acting at synonymous sites. We found that synonymous sites exhibit an AT rich and pyrimidine skewed nucleotide composition compared to both non-synonymous sites and non-coding regions. We also found some evidence for selection associated with both biased codon usage and conservation of regulatory sequences involved in mRNA processing, although some of these findings are subject to alternative non-adaptive interpretations. Regardless, the inferred strength of selection appears too weak to account for the variation in substitution rates between the mitochondrial genomes of plants and other multicellular eukaryotes. Therefore, these results are consistent with the interpretation that plant mitochondrial genomes experience a substantially lower mutation rate rather than increased functional constraints acting on synonymous sites. Nevertheless, there are important nucleotide composition patterns (particularly the differences between synonymous sites and non-coding DNA) that remain largely unexplained.     

Prion gene sequence variation within diverse groups of U.S. sheep,beef cattle,and deer

Prions are proteins that play a central role in transmissible spongiform encephalopathies in a variety of mammals. Among the most notable prion disorders in ungulates are scrapie in sheep, bovine spongiform encephalopathy in cattle, and chronic wasting disease in deer. Single nucleotide polymorphisms in the sheep prion gene (PRNP) have been correlated with susceptibility to natural scrapie in some populations. Similar correlations have not been reported in cattle or deer; however, characterization of PRNP nucleotide diversity in those species is incomplete. This report describes nucleotide sequence variation and frequency estimates for the PRNP locus within diverse groups of U.S. sheep, U.S. beef cattle, and free-ranging deer (Odocoileus virginianus and O. hemionus from Wyoming). DNA segments corresponding to the complete prion coding sequence and a 596-bp portion of the PRNP promoter region were amplified and sequenced from DNA panels with 90 sheep, 96 cattle, and 94 deer. Each panel was designed to contain the most diverse germplasm available from their respective populations to facilitate polymorphism detection. Sequence comparisons identified a total of 86 polymorphisms. Previously unreported polymorphisms were identified in sheep (9), cattle (13), and deer (32). The number of individuals sampled within each population was sufficient to detect more than 95% of all alleles present at a frequency greater than 0.02. The estimation of PRNP allele and genotype frequencies within these diverse groups of sheep, cattle, and deer provides a framework for designing accurate genotype assays for use in genetic epidemiology, allele management, and disease control.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Sequence analysis of UTR and coding region of kappa-casein gene of Indian riverine buffalo (Bubalus bubalis).

In this study, complete nucleotide as well as derived amino acid sequence characterization of water buffalo (Bubalus bubalis) kappa-casein gene has been presented. Kappa-casein cDNA clones were identified and isolated from a buffalo lactating mammary gland cDNA library. Sequence analysis of kappa-casein cDNA revealed 850 nucleotides with an open reading frame (ORF) of 573 nucleotides, encoding mature peptide of 169 amino acids. The 5' untranslated region (UTR) comprised 71 nucleotides, while 3' UTR was of 206 nucleotides. A total of 11 nucleotide and seven amino acid changes were observed in, buffalo (Bubalus bubalis) as compared to cattle (Bos taurus), sheep (Ovis aries) and goat (Capra hircus). Among these nucleotide changes, eight were unique in buffalo as they were fully conserved in cattle, sheep and goat. Majority of the nucleotide changes and all the amino acid changes; 14 (Asp-Glu), 19(Asp/Ser-Asn), 96(Ala-Thr), 126(Ala-Val), 128(Ala/Gly-Val), 156(Ala/Pro-Val) and 168(Ala/Glu-Val) were limited to exon IV. Three glycosylation sites, Thr 131, Thr 133 and Thr 142 reported in cattle and goat kappa-casein gene were also conserved in buffalo, however, in sheep Thr 142 was replaced by Ala. Chymosin hydrolysis site, between amino acids Phe 105 and Met 106, important for rennet coagulation process, were found to be conserved across four bovid species. Buffalo kappa-casein with the presence of amino acids Thr 136 and Ala 148 seems to be an intermediate of "A" and "B" variants of cattle. Comparison with other livestock species revealed buffalo kappa-casein sharing maximum nucleotide (95.5%) and amino acid (92.6%) similarity with cattle, whereas with pig it showed least sequence similarity of 76.0% and 53.2%, respectively. Phylogenetic analysis based on both nucleotide and amino acid sequence indicated buffalo kappa-casein grouping with cattle, while sheep and goat forming a separate cluster close to them. The non-ruminant species viz. camel, horse and pig were distantly placed, in separate lineages.     

Chromosomal location effects on gene sequence evolution in mammals.

BACKGROUND: Nucleotide substitution rates and G + C content vary considerably among mammalian genes. It has been proposed that the mammalian genome comprises a mosaic of regions - termed isochores - with differing G + C content. The regional variation in gene G + C content might therefore be a reflection of the isochore structure of chromosomes, but the factors influencing the variation of nucleotide substitution rate are still open to question. RESULTS: To examine whether nucleotide substitution rates and gene G + C content are influenced by the chromosomal location of genes, we compared human and murid (mouse or rat) orthologues known to belong to one of the chromosomal (autosomal) segments conserved between these species. Multiple members of gene families were excluded from the dataset. Sets of neighbouring genes were defined as those lying within 1 centiMorgan (cM) of each other on the mouse genetic map. For both synonymous substitution rates and G + C content at silent sites, neighbouring genes were found to be significantly more similar to each other than sets of genes randomly drawn from the dataset. Moreover, we demonstrated that the regional similarities in G + C content (isochores) and synonymous substitution rate were independent of each other. CONCLUSIONS: Our results provide the first substantial statistical evidence for the existence of a regional variation in the synonymous substitution rate within the mammalian genome, indicating that different chromosomal regions evolve at different rates. This regional phenomenon which shapes gene evolution could reflect the existence of 'evolutionary rate units' along the chromosome.     

Low sequence variation in the gene encoding the human beta-myosin heavy chain

Over 40 different mutations in the cardiac myosin heavy chain gene (MYH7) have been associated with familial hypertrophic cardiomyopathy (FHC), but no study has analyzed variation at this locus within the normal human population. Here we determine the extent and distribution of nucleotide variation in the 5808-bp MYH7 coding sequence in 25 normal individuals without FHC. We identified six single-nucleotide polymorphisms, none of which changes the encoded amino acid. At one of these sites, the frequencies of both alleles are equal; at the other five sites, the frequency of the rarer allele varies from 0.02 to 0.08. The nucleotide diversity (pi) calculated from these data is 1.73x10(-4)+/-0.49x10(-4), which is lower than the nucleotide diversity found in most other human autosomal genes. Substitution analysis of homologous genes between human and rodent also indicates that the MYH7 sequence has evolved at a very slow rate. The rate of both synonymous and nonsynonymous substitutions, especially in the portion of the sequence that encodes the alpha-helical myosin rod, is extremely low. The low level of even silent sequence variation in MYH7 in comparisons between human sequences and between human and rodent sequences may be a consequence of strong selective pressure against mutations that cause cardiomyopathy.     

A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome

A model of DNA sequence evolution applicable to coding regions is presented. This represents the first evolutionary model that accounts for dependencies among nucleotides within a codon. The model uses the codon, as opposed to the nucleotide, as the unit of evolution, and is parameterized in terms of synonymous and nonsynonymous nucleotide substitution rates. One of the model's advantages over those used in methods for estimating synonymous and nonsynonymous substitution rates is that it completely corrects for multiple hits at a codon, rather than taking a parsimony approach and considering only pathways of minimum change between homologous codons. Likelihood-ratio versions of the relative-rate test are constructed and applied to data from the complete chloroplast DNA sequences of Oryza sativa, Nicotiana tabacum, and Marchantia polymorpha. Results of these tests confirm previous findings that substitution rates in the chloroplast genome are subject to both lineage-specific and locus-specific effects. Additionally, the new tests suggest tha the rate heterogeneity is due primarily to differences in nonsynonymous substitution rates. Simulations help confirm previous suggestions that silent sites are saturated, leaving no evidence of heterogeneity in synonymous substitution rates.      

The cDNA and protein sequences of mouse lactate dehydrogenase B. Molecular evolution of vertebrate lactate dehydrogenase genes A (muscle), B (heart) and C (testis)

Mouse lactate dehydrogenase-B cDNAs were isolated from cDNA libraries of macrophage (ICR strain) and thymus (F1 hybrid of C57BL/6 and CBA strains), and their nucleotide sequences determined. The lactate dehydrogenase-B cDNA insert of thymus clone mB188 consists of the protein-coding sequence (1002 nucleotides), the 5' (46 nucleotides) and 3' (190 nucleotides) non-coding regions, and poly(A) tail (19 nucleotides), while macrophage clone mB168 contains a partial lactate dehydrogenase cDNA insert from codon no. 55 to the poly(A) tail. Seven silent nucleotide substitutions at codon no. 142, 143, 186, 187, 241, 285 and 292, as well as a single nucleotide change in the 3' non-coding region, were found between these different strains of mice. The predicted sequence of 333 amino acids, excluding initiation methionine, was confirmed by sequencing and/or compositional analyses of a total of 103 (31%) amino acids from tryptic peptides of mouse lactate dehydrogenase-B protein. The nucleotide sequence of the mouse coding region for lactate dehydrogenase B shows 86% identity with that of the human isoenzyme, and only eight of the 139 nucleotide differences resulted in amino acid substitutions at residues 10, 13, 14, 17, 52, 132, 236 and 317. The rates of nucleotide substitutions at synonymous and nonsynonymous sites in the mammalian lactate dehydrogenase genes are calculated. The rates of synonymous substitutions for lactate dehydrogenase genes A (muscle) and B (heart) are considerably higher than the average rate computed from human and rodent genes. The rates of nonsynonymous substitutions for lactate dehydrogenase genes A (muscle) and B (heart), particularly the latter, are highly conservative. The rates of synonymous and nonsynonymous substitutions for the lactate dehydrogenase-C gene are about the same as the average rates for mammalian genes. A phylogenetic tree of vertebrate lactate dehydrogenase protein sequences is constructed. In agreement with the previous results, this analysis further indicates that lactate dehydrogenase-C gene branched off earlier than did lactate dehydrogenase-A and lactate dehydrogenase-B genes.