Molecular Nature of Spontaneous Mutations in Mouse Lactate Dehydrogenase-a Processed Pseudogenes

The presence of at least ten mouse LDH-A pseudogenes was demonstrated in the genomic blot analysis, and four different processed pseudogenes have thus far been isolated and characterized. In this report, the nucleotide sequences to two different mouse lactate dehydrogenase-A processed pseudogenes, M11 and M14, were determined and compared with the protein-coding sequences of the mouse and rat LDH-A functional genes. In the pseudogene M11, the sequence of 64 nucleotides from codon no. 257 to 278 was tandemly duplicated. In the pseudogene M14, the sequence of 22 nucleotides from codon no. 68 to 75 was replaced by an inserted repetitive sequence of 242 nucleotides homologous to a mouse truncated R element. The pattern of nucleotide substitutions accumulated in mouse LDH-A pseudogenes M11 and M14, as well as that of pseudogene M10 identified previously, was analyzed, and the substitution frequencies of the C or G at the CG dinucleotide were found to be high.     

Comparative Analysis of Evolution in a Rodent Histone H2a Pseudogene

Sequences were obtained from five species of rodents that are orthologous to an H2a histone pseudogene from Mus musculus. The pseudogene is part of the cluster of replication-dependent histone genes found on Mus musculus chromosome 13. Comparative analysis of these five sequences together with the previously published sequence from M. musculus shows that this gene has likely been a pseudogene throughout the evolution of the genus Mus, while the gene from Rattus norvegicus is likely functional. Three large (>20 bp) deletions were found among the Mus pseudogenes, a feature that is very unusual compared to surveys of processed pseudogenes. In addition, there are two single-base deletions and one 4-bp insertion among the Mus pseudogenes. The species distributions of one of the large deletions and the 4-bp insertion require either independent insertions of an identical sequence, independent deletions with identical boundaries, or a deletion followed by precise reintegration of the original sequence. The evidence favors the hypothesis of multiple deletions with identical boundaries. The ``coding' regions of the Mus pseudogenes show a much reduced level of among-species variability in the 3′ half of the pseudogene, compared both to the 5′ half and to flanking sequences. This supports a hypothesis that the 3′ end of the pseudogene is the target of frequent gene conversion by functional H2a genes. Received: 1 April 1997 / Accepted: 12 June 1997     

Evolution of human cytoplasmic actin gene sequences: chromosome mapping and structural characterizations of three cytoplasmic actin-like pseudogenes including one on the Y chromosome

Eight recombinant phage clones containing cytoplasmic actin-like gene sequences have been isolated from a human genomic library for structural characterization. Kpn I family repeat sequences flank six of these actin genes isolated, and Alu family repeats are scattered throughout the DNA inserts of all eight phage clones. Three of these genes are γ actin-like, and the other five are β actin-like. The complete nucleotide sequence analysis of one β and one γ actin-like genes and their flanking regions demonstrates that they both are processed pseudogenes. Using unique DNA sequences flanking these two pseudogenes as hybridization probes for human-mouse somatic cell hybrid DNAs, we have mapped the two actin pseudogenes on human chromosomes 8 and 3, respectively. We have also determined the DNA sequence of a human Y chromosome-linked, processed actin pseudogene. The different values of sequence divergence of these processed pseudogenes and their functional counterparts allow us to estimate the time of generation of the pseudogenes. The results suggest that the cDNA insertion events generating the human cytoplasmic actin-like pseudogenes have occurred at significantly different times during the evolution of primates, after their separation from other mammalian species.     

Loci for human U1 RNA: structural and evolutionary implications

Three clones U1-1, U1-6, and U1-8 containing sequences related to human U1 RNA have been studied by sequence analysis. The results show that each of the three clones represents a distinct locus. The U1-6 locus is closely related to the HU1-1 locus, which is believed to represent a functional U1 gene. The U1-1 and U1-8 loci are pseudogenes by definition, since they contain sequences that are closely related to but not identical with the human U1 RNA sequence. The U1-6 locus contains the sequence T-A-T-A-T close to the 5'-end of the U1 sequence but it is unclear if this represents the promoter. When the U1-8 locus was compared to the U1-6 locus, it was observed that the 5'-flanking sequences, except in the immediate vicinity of the pseudogene, are as well-conserved as the U1-related sequence itself, at least up to position -220. The high degree of homology in the 5'-flanking region suggests that U1 genes have a much more strict sequence requirement with regard to 5'-flanking sequences than most other eukaryotic genes. The U1-6 and U1-8 loci contain the sequence T-A-T-G-T-A-G-A-T-G-A between positions -211 and -221. An identical sequence is present in the equivalent position in the HU1-1 locus, and may represent the promoter. The high degree of conservation in the postulated promoter region indicates that pseudogenes like U1-8 possibly could be expressed. A truncated U1-related sequence is present between 106 to 150 nucleotides upstream from the U1 gene/pseudogene in the U1-6, the U1-8 and the HU1-1 loci, suggesting that the U1 genes may have been clustered early in evolution. The U1-1 locus has a strikingly different structure from the U1-8 locus; the pseudogene itself is as closely related to the U1 RNA sequence as is the U1-8 pseudogene but the flanking sequences, both on the 5' and the 3' side, share no detectable homology with the corresponding regions in the U1-6 or U1-8 loci. It may therefore be postulated that small nuclear RNA pseudogenes are created by several different mechanisms.     

Nucleotide sequences of mouse genomic loci including a gene or pseudogene for U6 (4.8S) nuclear RNA.

We have isolated four clones which hybridize with U6 (4.8S) nuclear RNA, a mammalian small nuclear RNA(nRNA), from DNA of BALB/C mouse liver. Their restriction maps are totally different from each other, indicating that they derived from different loci in the mouse genome. The nucleotide sequences around the hybridizing region in the three clones have been determined. One clone gives a gene that is co-linear with the U6 RNA. There is a sequence TATAAAT beginning 31 nucleotides upstream of the gene, which may suggest that the U6 RNA is transcribed by RNA polymerase II. The other two clones contain a pseudogene for the U6 RNA which has 7 or 9 nucleotide changes from the RNA. The pseudogenes are surrounded by radically different sequences from those surrounding the gene, and they are closely linked to a pseudogene for another snRNA, 4.5S-I RNA, or a part of highly repetitive an interspersed sequence B1.     

Nonrandom integration of human U4 RNA pseudogenes.

Four loci for human U4 RNA have been characterized by DNA sequence analysis. The results show that all four loci represent pseudogenes, which are flanked by direct repeats. Three of the pseudogenes, designated U4/5, U4/6, and U4/8, have very similar structures; they are all truncated and contain the first 67 to 68 nucleotides of the U4 RNA sequence. Their properties suggest that they were created by integration of truncated cDNA copies of the U4 RNA into new chromosomal sites. An interesting observation was that their flanking regions exhibit sequence homology. A purine-rich 5'-flanking sequence 12 to 13 nucleotides long is almost perfectly conserved in all three loci. Boxes of homology were also found on the 3' side when the U4/6 and U4/8 loci were compared. The U4/4 locus has a slightly different structure; the pseudogene matches the first 79 nucleotides of U4 RNA, but contains a greater number of mutations than the other pseudogenes. Taken together, the results suggest that a frequently occurring type of pseudogene for human U4 was created by a RNA-mediated mechanism and that the integration sites have features in common.     

Conserved sequences in both coding and 5'' flanking regions of mammalian opal suppressor tRNA genes.

The rabbit genome encodes an opal suppressor tRNA gene. The coding region is strictly conserved between the rabbit gene and the corresponding gene in the human genome. The rabbit opal suppressor gene contains the consensus sequence in the 3' internal control region but like the human and chicken genes, the rabbit 5' internal control region contains two additional nucleotides. The 5' flanking sequences of the rabbit and the human opal suppressor genes contain extensive regions of homology. A subset of these homologies is also present 5' to the chicken opal suppressor gene. Both the rabbit and the human genomes also encode a pseudogene. That of the rabbit lacks the 3' half of the coding region. Neither pseudogene has homologous regions to the 5' flanking regions of the genes. The presence of 5' homologies flanking only the transcribed genes and not the pseudogenes suggests that these regions may be regulatory control elements specifically involved in the expression of the eukaryotic opal suppressor gene. Moreover the strict conservation of coding sequences indicates functional importance for the opal suppressor tRNA genes.     

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.     

Nucleotide sequences of the cDNA and an intronless pseudogene for human lactate dehydrogenase-A isozyme

Eight cDNA clones for lactate dehydrogenase-A isozyme (LDH-A) were isolated from a human fibroblast cDNA library, characterized, and no sequence heterogeneity was found. Four cDNA clones appear to contain nearly full-length cDNA inserts and the complete nucleotide sequence of 1710 base pairs consists of the protein-coding sequence (999 base pairs), the 5' (97 base pairs) and 3' (565 base pairs) untranslated regions and poly(dA) tail (49 base pairs). The predicted amino acid sequence of the human LDH-A polypeptide shows 92% homology (27 differences out of 331 amino acids compared) with that of the pig LDH-A subunit determined by direct protein sequencing [Kiltz et al. (1977) Hoppe-Seyler's Z. Physiol. Chem. 358, 123-127]. Human genomic clones containing an LDH-A pseudogene were isolated and the nucleotide sequence of 1635 base pairs from an intronless pseudogene was determined. The presence of two termination codons, two deletions of three nucleotides each and the replacement of three arginine residues at the active site (nos 98, 105 and 168) by other amino acids renders its coding region incapable of producing a functional LDH-A protein. A comparison between human LDH-A cDNA and the pseudogene sequences reveals 12.9% differences (114 transitions, 65 transversions and 36 deletions/insertions). Further, only four out of the 25 dCpdG dinucleotides present in the cDNA sequence remain unchanged, although the sequences possess 87.1% homology.     

Three types of rat U1 small nuclear RNA genes with different flanking sequences are induced to express in vivo

There are about 50 copies of U1 RNA genes/pseudogenes in the rat genome. To date, we have isolated so far 25 phage clones carrying a U1 RNA gene/pseudogene from two rat genomic libraries. The 12 clones were selected by hybridization with the U1 RNA coding region under a stringent condition, and were mapped and sequenced. Here, we report three types of U1 RNA genes with different flanking sequences, all of which were shown to be induced to express in vivo by transfection with their polylinker-inserted maxi U1 RNA genes into cultured rat cells. Although these three classes of U1 RNA genes have few homologous flanking sequences, they provide both upstream and downstream of the genes two conserved blocks, which may possibly play an important role in U1 RNA expression.     

Genes and pseudogenes for human U2 RNA: Implications for the mechanism of pseudogene formation

Three loci, designated U2/4, U2/6 and U2/7, which contain sequences related to human U2 RNA, have been studied. The U2/6 locus contains a tandem array of bona fide U2 genes. U2/4 and U2/7, in contrast, contain pseudogenes whose sequences deviate significantly from that of mammalian U2 RNA. The two pseudogenes appear to have been created by different mechanisms. The sequences that flank the pseudogene in the U2/4 locus lack homology to the corresponding sequences in functional human U2 genes, except for 10 base-pairs immediately following the 3′ end. The conserved 3′-flanking segment is homologous to those nucleotides that are present in a U2 RNA precursor. No direct repeats flank the pseudogene in the U2/4 locus. The observations thus suggest that a complementary DNA copy of the U2 RNA precursor was inserted into a blunt-ended chromosomal break to generate the U2/4 locus.The U2/7 locus, in contrast, revealed flanking sequence homology when compared to functional U2 genes, both on the 5′ and 3′ sides of the pseudogene. The homology was interrupted on both sides by repetitive sequences belonging to the Alu family. On the 5′ side the homology continues beyond the Alu repeats whereas on the 3′ side it ends precisely at the Alu repeat. This Alu repeat is inserted in a region where a homocopolymeric region of alternating C and T residues is located in functional U2 loci. The observed organization of the U2/7 locus suggests that a previously functional U2 locus was invaded by Alu repeats and subsequently accumulated base substitutions to become a pseudogene.     

Characterization of two novel human retropseudogenes related to the calmodulin-encoding gene, CaMII

Two human genomic clones (lambda hg22 and lambda hg29), containing two novel calmodulin (CaM) retropseudogenes, were isolated and characterized. The two pseudogenes show high similarity with the human CaMII cDNA, hCE1 [SenGupta et al., J. Biol. Chem. 262 (1987) 16663-16670] and the CaMII-type retropseudogene, hCE2 (CaMII-psi 1) [SenGupta et al., Nucleic Acids Res. 17 (1989) 2868]. One of them, in clone lambda hg22 (CaMII-psi 2), shows all the characteristics of a processed pseudogene. In clone lambda hg29 (CaMII-psi 3), however, an Alu repetitive sequence was detected immediately upstream from the ancestral 5'-untranslated region. Downstream from the truncated 3'-untranslated region, three additional copies of Alu repetitive sequences flanking about 750 nucleotides of unknown origin were found. Such a processed retropseudogene flanked by multiple Alu repeats may be a target for further recombination events. The three retropseudogenes CaMII-psi 1, CaMII-psi 2 and CaMII-psi 3 are estimated to be about 49, 21 and 25 million years old, respectively.     

Identification of ABCC6 pseudogenes on human chromosome 16p: implications for mutation detection in pseudoxanthoma elasticum

Pseudoxanthoma elasticum (PXE), a heritable disorder affecting the skin, eyes, and the cardiovascular system, has recently been linked to mutations in the ABCC6 gene on chromosome 16p13.1. The original mutation detection strategy employed by us consisted of the amplification of each exon of the ABCC6 gene with primer pairs placed on the flanking introns, followed by heteroduplex scanning and direct nucleotide sequencing. However, this approach suggested the presence of multiple copies of the 5'-region of the gene when total genomic DNA was used as a template. In this study, we have identified two pseudogenes containing sequences highly homologous to the 5'-end of ABCC6. First, by the use of allele-specific polymerase chain reaction (PCR), two bacterial artificial chromosome (BAC) clones containing a putative pseudogene of ABCC6, designated as ABCC6-psi 1, were isolated from the human BAC library. Sequence analysis of ABCC6-psi 1 revealed it to be a truncated copy of ABCC6, which contains the upstream region and exon 1 through intron 9 of the gene. Secondly, a homology search of a high-throughput sequence database revealed the presence of another truncated copy of ABCC6, which was designated as ABCC6-psi 2, and which was shown to harbor upstream sequences and a segment spanning exon 1 through intron 4 of ABCC6. In addition to several nucleotide differences in the flanking introns and the upstream region, both pseudogenes contain several nucleotide changes in the exonic sequences, including stop codon mutations, which complicate mutation analysis in patients with PXE. Nucleotide differences in flanking introns between these two pseudogenes and ABCC6 allowed us to design allele-specific primers that eliminated the amplification of both pseudogene sequences by PCR and provided reliable amplification of ABCC6-specific sequences only. The use of allele-specific PCR has revealed, thus far, two novel 5'-end PXE mutations, 179del9 and T364R in exons 2 and 9, respectively, and several polymorphisms within the upstream region and exons 1-9 of ABCC6. These strategies facilitate comprehensive analysis of ABCC6 for mutations in PXE.     

Mutation Pattern Variation Among Regions of the Primate Genome

We sequenced three argininosuccinate-synthetase-processed pseudogenes (ΨAS-A1, ΨAS-A3, ΨAS-3) and their noncoding flanking sequences in human, orangutan, baboon, and colobus. Our data showed that these pseudogenes were incorporated into the genome of the Old World monkeys after the divergence of the Old World and New World monkey lineages. These pseudogene flanking regions show variable mutation rates and patterns. The variation in the G/C to A/T mutation rate (u) can account for the unequal GC contents at equilibrium: 34.9, 36.9, and 41.7% in the pseudogene ΨAS-A1, ΨAS-A3, and ΨAS-3 flanking regions, respectively. The A/T to G/C mutation rate (v) seems stable and the u/v ratios equal 1.9, 1.7, and 1.4 in the flanking regions of ΨAS-A1, ΨAS-A3, and ΨAS-3, respectively. These ``regional' variations of the mutation rate affect the evolution of the pseudogenes, too. The ratio u/v being greater than 1.0 in each case, the overall mutation rate in the GC-rich pseudogenes is, as expected, higher than in their GC-poor flanking regions. Moreover, a ``sequence effect' has been found. In the three cases examined u and v are higher (at least 20%) in the pseudogene than in its flanking region—i.e., the pseudogene appears as mutation ``hot' spots embedded in ``cold' regions. This observation could be partly linked to the fact that the pseudogene flanking regions are long-standing unconstrained DNA sequences, whereas the pseudogenes were relieved of selection on their coding functions only around 30–40 million years ago. We suspect that relatively more mutable sites maintained unchanged during the evolution of the argininosuccinate gene are able to change in the pseudogenes, such sites being eliminated or rare in the flanking regions which have been void of strong selective constraints over a much longer period. Our results shed light on (1) the multiplicity of factors that tune the spontaneous mutation rate and (2) the impact of the genomic position of a sequence on its evolution. Received: 10 February 1997 / Accepted: 21 April 1997     

DNA sequences complementary to human 7 SK RNA show structural similarities to the short mobile elements of the mammalian genome

A complementary DNA clone of 7 SK RNA from HeLa cells was used to study the genomic organization of 7 SK sequences in the human genome. Genomic hybridizations and genomic clones show that 7 SK is homologous to a family of disperse repeated sequences most of which lack the 3' end of the 7 SK RNA sequence. Only few of the genomic K sequences are homologous to both 3' and 5' 7 SK probes and presumably include the gene(s) for 7 SK RNA. The sequence of four genomic 7 SK clones confirms that they are in most cases pseudogenes. Although Alu sequences are frequently found near the 3' and 5' end of K DNA, the sequences immediately flanking the pseudogenes are different in all clones studied. However, direct repeats were found flanking directly the K DNA or the K-Alu unit, suggesting that the K sequences alone or in conjunction with Alu DNA might constitute a mobile element.     

The nuclear autoantigen La/SS-B: mapping and sequencing of the gene and the three retropseudogenes

One target of autoantibodies in sera of patients with systemic lupus erythematosus or primary Sjögren's syndrome is the nuclear autoantigen La/SS-B. Lambda clones and cosmids were isolated, which contained the sequences of the La gene and the three La pseudogenes. They were used for preparation of a physical map. Finally, the La gene and pseudogenes were sequenced. The pseudogenes were characterized as retropseudogenes. Their evolutionary ages were estimated to be approx. 4, 4.5 and 5 million years. Inserts of 4, 16 and 24 nucleotides, which were mostly A-residues, were found in exon 7 of the respective pseudogene. The oldest pseudogene contained the longest insert, the youngest pseudogene contained the smallest insert. The oligonucleotides seem to be the result of repeated inserts of A-residues in a hot spot region of the La genes. Two La cDNAs were isolated which contained either a deletion or an insert of an A-residue at the same position.     

Association of a truncated cytochrome c processed pseudogene with a similarly truncated member from a long interspersed repeat family of rat.

The cytochrome c multigene family of rat contains approximately 30 processed pseudogenes that represent genomic DNA copies of three alternate mRNAs. Here, the DNA sequence of an unusual processed pseudogene reveals that it has a complete 3' noncoding region including a short poly A tail but unlike the others is abruptly truncated at its 5' end, 19 amino acid codons from the translation terminator. At this position the pseudogene is fused through 17 consecutive adenylic acid residues to a 1.3 kb repetitive sequence. This repetitive element is flanked by direct repeats and represents a truncated member from a major long interspersed repeat family. The rat element is a composite of sequences observed in long interspersed repeats from both rodents and primates. Comparison to the equivalent mouse sequences shows that the 5' half of the repeat distal to the pseudogene has an open reading frame and is highly conserved whereas the half adjacent to the pseudogene is evolutionarily unstable. The proportion of cytochrome c pseudogene recombinant clones containing this repetitive DNA is 3 fold greater than observed in random isolates and may reflect a general tendency of processed pseudogenes to associate with other repetitive sequences in the genome.