The cDNA and protein sequences of human lactate dehydrogenase B.

Human lactate dehydrogenase B (LDH-B) cDNA was isolated and sequenced. The LDH-B cDNA insert consists of the protein-coding sequence (999 bp), the 5' (54 bp) and 3' (203 bp) non-coding regions, and the poly(A) tail (50 bp). The predicted sequence of 333 amino acid residues was confirmed by amino acid composition and/or sequence analyses of a total of 185 (56%) residues from tryptic peptides of human LDH-B protein. The nucleotide and amino acid sequences of the human LDH-B coding region show 68% and 75% homologies respectively with those of the human LDH-A. The peptide map and amino acid composition data have been deposited as Supplementary Publication SUP 50139 (7 pages) at the British Library Lending Division, Boston Spa, Wetherby, West Yorkshire LS23 7BQ, U.K., from whom copies are available on prepayment [see Biochem. J. (1987) 241, 5].     

Molecular cloning and nucleotide sequence of the cDNA for sperm-specific lactate dehydrogenase-C from mouse.

Mouse sperm-specific lactate dehydrogenase-C (LDH-C) cDNA was cloned and sequenced from lambda gt11 expression library. The LDH-C cDNA insert of 1236 bp consists of the protein-coding sequence (999 bp), the 5' (54 bp) and 3' (113 bp) non-coding regions, and the poly(A) tail (70 bp). The Northern blot analysis of poly(A)-containing RNAs from mouse testes and liver indicates that the LDH-C gene is expressed in testes but not in liver, and that its mRNA is approx. 1400 nucleotides in length. The nucleotide and amino acid sequences of the mouse LDH-C cDNA show 73% and 72% homologies, respectively, with those of the mouse LDH-A. The Southern blot analysis of genomic DNAs from mouse liver and human placenta indicates the presence of multiple LDH-C gene-related sequences.     

Evolution of nucleic acids coding for ribonucleases: the mRNA sequence of mouse pancreatic ribonuclease

The cDNA of mouse pancreatic mRNA has been cloned. After the library was screened with a rat ribonuclease cDNA probe, the positive clones were isolated and sequenced. There were no differences from the previously determined protein sequence. The mRNA codes for a preribonuclease of 149 amino acid residues including a signal peptide of 25 amino acids. The 3' noncoding region has a length of 260 bp, and the total mRNA length is approximately 940 bp. Comparison with the rat pancreatic ribonuclease sequence showed a high rate of nucleotide substitution. Within the coding region, nonsynonymous and synonymous substitution rates are 4.3 X 10(-9) and 15 X 10(-9) nucleotide substitutions/site/year, respectively. The latter value is one of the highest rates observed in the molecular evolution of mammalian nuclear genes. In the signal sequences the synonymous substitution rate is much lower and about the same as the nonsynonymous rate. Signal sequences of other mouse and rat proteins also exhibit little difference between synonymous and nonsynonymous rates. The sequences of rat and mouse pancreatic ribonuclease messengers were compared with those of bovine pancreatic, seminal, and brain ribonuclease. While the 3' noncoding regions of rat and mouse are very similar, as are those of the three bovine messengers, there is no significant similarity between both rodent and the three bovine messengers for the greater part of these regions. There is a duplication of approximately 50 nucleotides in the 3' noncoding region of the bovine messengers, with a region rich in A and C in between. The presence of this structural feature may be correlated with recent gene duplications that have occurred in the bovine genome.     

Molecular evolution of mammalian lactate dehydrogenase-A genes and pseudogenes: association of a mouse processed pseudogene with a B1 repetitive sequence

A mouse genomic clone containing a lactate dehydrogenase-A (LDH-A) processed pseudogene and a B1 repetitive element was isolated, and a nucleotide sequence of approximately 3 kb was determined. The pseudogene and B1 element are flanked by perfect 13-bp repeats, and the B1 sequence starts at 14 nucleotides 3' to the presumptive polyadenylation signal of the pseudogene. The nucleotide sequences of the LDH-A genes and processed pseudogenes from mouse, rat, and human were compared, and a phylogenetic tree was constructed. The rate and pattern of nucleotide substitutions in the LDH-A pseudogenes are similar to previously reported results (Li et al. 1984). The average rate of nucleotide substitutions in the LDH-A pseudogenes is 4.3 X 10(- 9)/site/year. The substitutions of C----T and G----A are most frequent, and A----G substitutions are relatively high. The rate of synonymous substitutions in the LDH-A genes is 5.3 X 10(-9), which is not significantly higher than the average rate of 4.7 X 10(-9) for 35 mammalian genes. The rate of nonsynonymous substitutions in the LDH-A genes is 0.20 X 10(-9), which is considerably lower than the average rate of 0.88 X 10(-9) for 35 mammalian genes. Thus, the mammalian LDH-A gene appears to be highly conserved in evolution.      

Human lactate dehydrogenase-B processed pseudogene: nucleotide sequence analysis and assignment to the X-chromosome

Two human genomic clones containing the lactate dehydrogenase-B processed pseudogene were isolated from two patients deficient in lactate dehydrogenase-B isozyme. The sequences of 3,287 nucleotides, including the pseudogenes and its flanking regions, from both clones were found to be identical except for three differences in the pseudogenes. The sequences of 1,286 nucleotides from these two pseudogenes exhibited 93% homology with the cDNA sequence of the lactate dehydrogenase-B functional gene, and the pseudogene contained 75/76 base substitutions, 11/12 single-base deletions, and 5 single-base insertions. This pseudogene was mapped to the x-chromosome by dot-blot analysis using a probe for the pseudogene or its 5' flanking sequence.     

The nucleotide sequence of myosin light chain (L-2A) mRNA from embryonic chicken cardiac muscle tissue.

The nucleotide sequence of a cDNA clone (pML10) for chicken cardiac myosin light chain is described. The cDNA insert contains 613 nucleotides representing the entire coding sequence, with the exception of nine NH2-terminal amino acids, and the full 3'-non-coding region of 146 nucleotides. The missing 5' terminus of the mRNA, not represented in the clone pML10, was obtained by extension of the cDNA using a 43 nucleotide long internal EcoR1 fragment as a primer. The non-coding region contains several direct and inverted repeated sequences and the polyadenylation signal sequence AATAAA. The coding portion exhibits non-random usage of synonymous codons with a strong bias for codons ending in G and C.     

Codon Substitution in Evolution and the "Saturation" of Synonymous Changes

A mathematical model for codon substitution is presented, taking into account unequal mutation rates among different nucleotides and purifying selection. This model is constructed by using a 61 X 61 transition probability matrix for the 61 nonterminating codons. Under this model, a computer simulation is conducted to study the numbers of silent (synonymous) and amino acid-altering (nonsynonymous) nucleotide substitutions when the underlying mutation rates among the four kinds of nucleotides are not equal. It is assumed that the substitution rates are constant over evolutionary time, the codon frequencies being in equilibrium, and, thus, the numbers of synonymous and nonsynonymous substitutions both increase linearly with evolutionary time. It is shown that, when the mutation rates are not equal, the estimate of synonymous substitutions obtained by F. Perler, A. Efstratiadis, P. Lomedico, W. Gilbert, R. Kolodner and J. Dodgson's "Percent Corrected Divergence" method increases nonlinearly, although the true number of synonymous substitutions increases linearly. It is, therefore, possible that the "saturation" of synonymous substitutions observed by Perler et al. is due to the inefficiency of their method to detect all synonymous substitutions.     

Nucleotide and deduced amino acid sequence of the alpha subunit of yeast pyruvate dehydrogenase

A 413-base cDNA insert encoding a portion of the alpha subunit of pyruvate dehydrogenase (E1 alpha; EC 1.2.4.1) from Saccharomyces cerevisiae was isolated from a lambda gt11 cDNA library by immunoscreening and by hybridization with an oligonucleotide probe which corresponded to the amino acid sequence around the phosphorylation site of E1 alpha. This cDNA was subcloned, sequenced and used as a probe to isolate two additional cDNA inserts which were subcloned and sequenced. These overlapping clones comprised the carboxyl-terminal part of E1 alpha. To identify the missing nucleotide sequence, the polymerase chain reaction was used to amplify yeast genomic DNA with synthetic oligonucleotide primers based on the amino-terminal sequence of E1 alpha and the 5' end of one of the cDNA clones. Three DNA fragments were isolated and sequenced. The composite nucleotide sequence has an open reading frame of 1260 nucleotides encoding a putative presequence of 33 amino acids and a mature protein of 387 amino acids (Mr = 42,703). Hybridization analysis showed that the size of the mRNA is about 1.4 kilobases.     

Isolation and characterization of a complementary DNA clone coding for the E1 beta subunit of the bovine branched-chain alpha-ketoacid dehydrogenase complex: complete amino acid sequence of the precursor protein and its proteolytic processing

A 1.7-kb cDNA clone encoding the entire precursor of the E1 beta subunit of the branched-chain alpha-ketoacid dehydrogenase (BCKDH) complex was isolated from a bovine liver cDNA library by screening with a mixture of synthetic oligonucleotide probes corresponding to the C-terminal five-residue sequence of the mature E1 beta subunit. A partial amino acid sequence was determined by Edman degradation of the intact subunit and the peptides generated by cleavage at the lysyl bonds. Nucleotide sequence analysis revealed that the isolated cDNA clone contained the 5'-untranslated sequence of 186 nucleotides, the translated sequence of 1176 nucleotides, and the 3'-untranslated sequence of 306 nucleotides with a poly(A) tail. A type AATAAA polyadenylation signal was located 17 nucleotides upstream of the start of a poly(A) tail. Comparison of the amino acid sequence predicted from the nucleotide sequence of the cDNA insert of the clone with the partial amino acid sequence of the mature BCKDH E1 beta subunit showed that the cDNA insert encodes for a 342 amino acid subunit with Mr 37,745 and that the subunit is synthesized as the precursor with a leader sequence of 50 amino acids and processed at the N-terminus. Northern blot analysis using the cDNA insert as a probe showed the presence of a 1.8-1.9-kb mRNA in bovine liver, suggesting that the insert covers nearly a full length of mRNA. Alignment of the deduced amino acid sequence of bovine BCKDH E1 beta with that of the human pyruvate dehydrogenase (PDH) complex E1 beta subunit revealed a high degree of sequence homology throughout the two enzymes.(ABSTRACT TRUNCATED AT 250 WORDS)     

Nucleic acid composition,codon usage,and the rate of synonymous substitution in protein-coding genes

Summary Based on the rates of synonymous substitution in 42 protein-codin gene pairs from rat and human, a correlation is shown to exist between the frequency of the nucleotides in all positions of the codon and the synonymous substitution rate. The correlation coefficients were positive for A and T and negative for C and G. This means that AT-rich genes accumulate more synonymous substitutions than GC-rich genes. Biased patterns of mutation could not account for this phenomenon. Thus, the variation in synonymous substitution rates and the resulting unequal codon usage must be the consequence of selection against A and T in synonymous positions. Most of the varition in rates of synonymous substitution can be explained by the nucleotide composition in synonymous positions. Codon-anticodon interactions, dinucleotide frequencies, and contextual factors influence neither the rates of synonymous substitution nor codon usage. Interestingly, the nucleotide in the second position of codons (always a nonsynonymous position) was found to affect the rate of synonymous substitution. This finding links the rate of nonsynonymous substitution with the synonymous rate. Consequently, highly conservative proteins are expected to be encoded by genes that evolve slowly in terms of synonymous substitutions, and are consequently highly biased in their codon usage.     

Sequence analysis and comparison of cDNAs of the zein multigene family .

The nucleotide sequence of two zein cDNAs in hybrid plasmids A20 and B49 have been determined. The insert in A20 is 921 bp long including a 5' non-coding region of 60 nucleotides, preceded by what is believed to be an artifactual sequence of 41 nucleotides, and a 3' non-coding region of 87 nucleotides. The B49 insert is 467 bp long and includes approximately one-half the protein coding sequence as well as a 3' non-coding region of 97 nucleotides. These sequences have been compared with the previously published sequence of another zein clone, A30 . A20 and A30 , both encoding 19 000 mol. wt. zeins , have approximately 85% homology at the nucleotide level. The B49 sequence, corresponding to a 22 000 mol. wt. zein, has approximately 65% homology to either A20 or A30 . All three zeins share common features including nearly identical amino acid compositions. In addition, the tandem repeats of 20 amino acids first seen in A30 are also present in A20 and B49 .     

Cloning and comparative mapping of a human class III (chi) alcohol dehydrogenase cDNA

A cDNA encoding human class III (chi ADH5) alcohol dehydrogenase was isolated, sequenced and used to comparatively map this unusual ADH. In their coding sequences, the three major ADH classes were approximately equisimilar, class II and III ADHs sharing the highest sequence identity (67%). A class III-like ADH was mapped to mouse chromosome 3, site of the ADH gene complex, and synteny of ADH5 with four other ADH loci on human chromosome 4 was confirmed. The nearly full-length 1613 nucleotide cDNA contained 433 nucleotides of 3' nontranslated sequence and two possible initiation sites for translation. A protein of 374 amino acid residues could be synthesized using the potential initiation codon at nucleotide 59. However, use of the likely initiation codon at nucleotide 5 would produce a protein of 392 residues with 19 additional N-terminal residues as compared to the known protein sequence. The derived protein sequence also differs at residue 166, where Tyr is found. This difference, due to a single base substitution, could result from cloning artifact, polymorphism, or two expressed class III ADH genes.     

Human glyceraldehyde-3-phosphate dehydrogenase pseudogenes: Molecular evolution and a possible mechanism for amplification

We screened two human genomic libraries and isolated 14 different clones, designated λG1 and EG1-EG13, homologous to human glyceraldehyde-3-phosphate dehydrogenase (GAPD) cDNA. Subcloning and sequencing these recombinant phages led us to classify them as five different pseudogenes (ψG1–ψG5). All these sequences show such features typical of processed pseudogenes as numerous mutations, insertions, and deletions. The identity of numerous mutated sites among these pseudogenes and the presence of two Alu sequences flanking both ends of ψG1 suggest that GAPD pseudogenes originated from a unique reverse transcribed mRNA followed by gene duplication. The rate of nucleotide substitutions per site per year for known GAPD functional genes is low both for the synonymous substitutions (1.87×10⁻⁹) and for the nonsynonymous substitutions (0.12￠10⁻⁹) and indicates that the GAPD cDNA sequence is well conserved not only at the amino acid level, but also at the nucleotide level. The rate of nucleotide substitutions per site per year for GAPD pseudogenes shows a higher value (5.9×10⁻⁹) and suggests that these pseudogenes do not have any functional role. This work was supported by grants from the Consiglio Nazionale delle Ricerche and the Ministero Pubblica Istruzione (Rome, Italy). Special acknowledgment is given to the “Progetto Finalizzato Ingegneria Genetica e Basi Molecolari delle Malattie Ereditarie.”     

Mammalian gene evolution: Nucleotide sequence divergence between mouse and rat

As a paradigm of mammalian gene evolution, the nature and extent of DNA sequence divergence between homologous protein-coding genes from mouse and rat have been investigated. The data set examined includes 363 genes totalling 411 kilobases, making this by far the largest comparison conducted between a single pair of species. Mouse and rat genes are on average 93.4% identical in nucleotide sequence and 93.9% identical in amino acid sequence. Individual genes vary substantially in the extent of nonsynonymous nucleotide substitution, as expected from protein evolution studies; here the variation is characterized. The extent of synonymous (or silent) substitution also varies considerably among genes, though the coefficient of variation is about four times smaller than for nonsynonymous substitutions. A small number of genes mapped to the X-chromosome have a slower rate of molecular evolution than average, as predicted if molecular evolution is male-driven. Base composition at silent sites varies from 33% to 95% G + C in different genes; mouse and rat homologues differ on average by only 1.7% in silent-site G + C, but it is shown that this is not necessarily due to any selective constraint on their base composition. Synonymous substitution rates and silent site base composition appear to be related (genes at intermediate G + C have on average higher rates), but the relationship is not as strong as in our earlier analyses. Rates of synonymous and nonsynonymous substitution are correlated, apparently because of an excess of substitutions involving adjacent pairs of nucleotides. Several factors suggest that synonymous codon usage in rodent genes is not subject to selection.     

Complete nucleotide sequence of ovine alpha-lactalbumin mRNA

The nucleotide sequence of ovine alpha-lactalbumin mRNA has been determined by chemical sequencing of two cDNA recombinant plasmids and a primer extension product. Ovine alpha-lactalbumin mRNA contains 723 nucleotides (excluding the poly(A) tail), with a 5' non-coding region of 26 nucleotides, followed by the 426 nucleotides of the coding region which determines a sequence signal of 19 amino acid residues and the 123 amino acid residues of mature alpha-lactalbumin. The coding region is followed by a 3' untranslated sequence of 271 nucleotides. The derived amino acid sequence of ovine pre-alpha-lactalbumin differs from that of its bovine counterpart by 8 amino acid substitutions, all but one originating from single mutations. Comparison of sequences of guinea pig, rat and human alpha-lactalbumin mRNAs with their ovine and bovine counterparts has revealed that these molecules have rapidly evolved. The highest degree of conservation was observed in the region coding for the mature protein and corresponds essentially to sequences which interact with UDP-galactosyltransferase and Ca2+ ions.     

Molecular evolution of hemagglutinin genes of H1N1 swine and human influenza A viruses

Summary The hemagglutinin (HA) genes of influenza type A (H1N1) viruses isolated from swine were cloned into plasmid vectors and their nucleotide sequences were determined. A phylogenetic tree for the HA genes of swine and human influenza viruses was constructed by the neighbor-joining method. It showed that the divergence between swine and human HA genes might have occurred around 1905. The estimated rates of synonymous (silent) substitutions for swine and human influenza viruses were almost the same. For both viruses, the rate of synonymous substitution was much higher than that of nonsynonymous (amino acid altering) substitution. It is the case even for only the antigenic sites of the HA. This feature is consistent with the neutral theory of molecular evolution. The rate of nonsynonymous substitution for human influenza viruses was three times the rate for swine influenza viruses. In particular, nonsynonymous substitutions at antigenic sites occurred less frequently in swine than in humans. The difference in the rate of nonsynonymous substitution between swine and human influenza viruses can be explained by the different degrees of functional constraint operating on the amino acid sequence of the HA in both hosts.     

Accelerated evolution in the protein-coding region of galectin cDNAs, congerin I and congerin II, from skin mucus of conger eel (Conger myriaster).

Two cDNAs encoding galectins named congerins I and II from the skin mucus of conger eel (Conger myriaster) were isolated and sequenced. Comparison of the nucleotide sequences of congerins I and II showed that the sequence similarities of the 5' and 3' untranslated regions (86 and 88%, respectively) were much higher than those of the protein-coding region (73%). The numbers of nucleotide substitutions per site (KN) for the untranslated regions are smaller than the numbers of nucleotide substitutions per synonymous site (KS) for the protein coding region. Furthermore, nonsynonymous nucleotide substitutions have accelerated more frequently than synonymous nucleotide substitutions in the protein coding region (KA/KS = 2.57). These results suggest that accelerated substitutions have occurred in the protein-coding regions of galectin genes to generate diverse galectins with different molecular properties. Northern blot analysis showed that both congerins were expressed not only in the skin tissues but also in the stomach of conger eel.     

Resolution of phylogenetic relationships among recently evolved species as a function of amount of DNA sequence: an empirical study based on woodpeckers (Aves: Picidae)

Synonymous substitutions in the 13 mitochondrial encoded protein genes form a large pool of characters that should approach the ideal for phylogenetic analysis of being independently and identically distributed. Pooling sequences from multiple mitochondrial protein-coding genes should result in statistically more powerful estimates of relationships among species that diverged sufficiently recently that most nucleotide substitutions are synonymous. Cytochrome oxidase I (COI) was sequenced for woodpecker species for which cytochrome b (cyt b) sequences were available. A pairing-design test based on the normal distribution indicated that cyt b evolves more rapidly than COI when all nucleotides are compared but their rates are equal for synonymous substitutions. Nearly all of the phylogenetically informative substitutions among woodpeckers are synonymous. Statistical support for relationships, as measured by bootstrap proportions, increased as the number of nucleotides increased from 1047 (cyt b) to 1512 (COI) to 2559 nucleotides (aggregate data set). Pseudo-bootstrap replicates showed the same trend and increasing the amount of sequence beyond the actual length of 2559 nucleotides to 5120 (2x) resulted in stronger bootstrap support, even though the amount of phylogenetic information was the same. However, the amount of sequence required to resolve an internode depends on the length of the internode and its depth in the phylogeny.     

Nucleotide sequence encoding the precursor of the small subunit of ribulose 1,5-bisphosphate carboxylase from Lemna gibba L.G-3

We have sequenced a cDNA clone, pLgSSU, which encodes the small subunit of ribulose 1,5-bisphosphate carboxylase of Lemna gibba L.G-3 a monocot plant. This clone contains a 832 basepair insert which encodes the entire 120 amino acids of the mature small subunit polypeptide (Mr = 14,127). In addition this clone encodes 53 amino acids of the amino terminal transit peptide of the precursor polypeptide and 242 nucleotides of the 3' non-coding region. Comparison of the nucleotide sequence of pLgSSU with Lemna gibba genomic sequences homologous to the 5' end of the cDNA clone suggests that nucleotides encoding four amino-terminal amino acids of the transit peptide are not included in the cDNA clone. The deduced amino acid sequence of the Lemna gibba mature small subunit polypeptide shows 70-75% homology to the reported sequences of other species. The transit peptide amino acid sequence shows less homology to other species. There is 50% homology to the reported soybean sequence and only 25% homology to the transit sequence of another monocot, wheat.