Molecular structures of human argininosuccinate synthetase pseudogenes. Evolutionary and mechanistic implications

In the human genome there is one expressed gene for argininosuccinate synthetase and 14 pseudogenes. A cDNA coding for human argininosuccinate synthetase was used to screen a human genomic library. Twenty-five unique genomic clones were isolated and extensively characterized. At least seven clones represented processed argininosuccinate synthetase pseudogenes that lost the introns in the expressed gene. Restriction mapping demonstrated that these processed pseudogenes were located in distinct regions of the human genome. Complete nucleotide sequences of two processed pseudogenes, psi AS-1 and psi AS-3, and a partial sequence of psi AS-7 were determined. Both psi AS-1 and psi AS-3 had an adenine-rich region at their 3' end and were flanked by distinct imperfect direct repeats. A comparison of these pseudogene sequences to that of the cDNA demonstrated that psi AS-1 and psi AS-3 were 93% homologous to the cDNA, whereas psi AS-7 was 89% homologous to the cDNA. Therefore, it is estimated that psi AS-1 and psi AS-3 were created 10-11 million years ago, whereas psi AS-7 arose approximately 21 million years ago. We have estimated the evolutionary rate for the expressed argininosuccinate synthetase gene based on the sequences of psi AS-1 and psi AS-3. These data indicate that the expressed argininosuccinate synthetase gene is evolving at a rate similar to that of the beta-globin gene and much faster than the alpha-tubulin gene. Furthermore, a comparison of the sequences of psi AS-1 and psi AS-3 suggests the possibility that these pseudogenes arose from a common intermediate.     

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     

Identification of three human pseudogenes for subunit VIb of cytochrome c oxidase: a molecular record of gene evolution.

Three pseudogenes for the nuclear-encoded subunit VIb of cytochrome c oxidase (COX) were isolated by screening a human genomic library with cloned human cDNA coding for COX subunit VIb. The nucleotide sequences of the pseudogenes, designated psi COX6b-1, psi COX6b-2 and psi COX6b-3, were determined. Pseudogene psi COX6b-1 bears all the hallmarks of a processed pseudogene and diverged from the parental gene after the divergence of man and cow. Alu repetitive elements were integrated into the structural sequences of the other two pseudogenes. Comparison with the human and bovine cDNA sequences encoding COX subunit VIb suggests that psi COX6b-2 and psi COX6b-3 were formed earlier in evolution than psi COX6b-1. Genomic Southern analysis indicated that a few more pseudogenes for COX subunit VIb are likely to be present in the human genome. Identical nt differences with respect to the human cDNA sequence in the pseudogenes provide some clues on the evolution of the ancestral gene coding for COX subunit VIb.     

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.     

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.”     

A human dihydrofolate reductase intronless pseudogene with an Alu repetitive sequence: multiple DNA insertions at a single chromosomal site

A dihydrofolate reductase (DHFR) pseudogene, hDHFR-psi 3 has been isolated from a human genomic DNA fragment library. Sequence analysis of this gene revealed a lack of introns and the presence of a tract of nine adenines, 90 bp downstream from the end of the coding sequence. These features suggest that hDHFR-psi 3 was derived from a processed RNA molecule that has been converted into DNA and inserted into a chromosome, analogous to the origin of three intronless human DHFR genes previously described. An interesting feature of hDHFR-psi 3 is the presence of a member of the Alu moderately repetitive DNA sequence family within the DHFR coding region. This Alu element is flanked by a 16 bp directly repeated DNA segment derived from DHFR coding sequences. The Alu element apparently has been inserted into the intronless DHFR pseudogene and thus, there have been two insertions at a single chromosomal locus. The hDHFR-psi 3 contains only the 3' half of the DHFR coding sequence. Immediately upstream from the directly repeated sequence before the Alu element is an adenine-rich tract. The DNA farther upstream is moderately repetitive and is related to neither DHFR nor Alu DNA sequence. Therefore, it seems possible that a third insertion has occurred at the same site further disrupting the hDHFR coding sequences.     

Structure of the goat psi beta y beta-globin pseudogene. Analysis of goat pseudogene evolutionary patterns

The 12-member beta-globin gene locus of the goat contains three beta(adult)-type pseudogenes, one in each of three four-gene subsets of the locus. We have determined the complete nucleotide sequence of psi beta y, the pseudogene present in the most downstream four-gene subset, which also contains the functional fetal gene, beta F. psi beta y contains, throughout its length, numerous incapacitating mutations in common with the previously sequenced goat psi beta x and psi beta z pseudogenes consistent with the model that all were descended from a common pseudogene ancestor which became defective prior to the expansion of the beta-globin locus in the goat lineage. Evolutionary analysis of the psi beta y sequence in comparison to psi beta x and psi beta z provides evidence that nucleotide substitutions were fixed in a random manner within these pseudogenes with respect to polarity, coding versus non-coding regions, and replacement sites versus silent sites. However, substitutions appear to have accumulated asymmetrically between different pseudogenes in a manner that provides evidence for partial gene conversion. Moreover, the presence of deletions in goat psi beta y, which are also observed in the cow pseudogene psi 2, but not in the cow psi 1 pseudogene, indicate that goat psi beta y and cow psi 2 are orthologous but cow psi 1 actually arose prior to the goat/cow divergence. The authentic goat orthologue to cow psi 1 temporarily existed in the goat lineage but was deleted, probably prior to the divergence of goats and sheep.     

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.     

Characterization of the bovine BCL2L1 gene and related pseudogenes

We have amplified and characterized partial regions of exons 2 and 3 of the bovine BCL2L1 gene, one of the anti-apoptotic members of the B-cell lymphoma 2 gene family. Cloning and sequencing of the amplified products revealed the existence of several BCL2L1-related sequences, including the bovine BCL2L1 gene and various processed pseudogenes. The bovine BCL2L1 gene revealed two polymorphic nucleotide sequences that resulted in two protein variants, with amino acid replacements at positions 60 and 69. In addition, we report three bovine BCL2L1-related sequences (BCL2L1psi) that probably correspond to intronless processed pseudogenes. These BCL2L1psi pseudogene sequences have accumulated multiple substitutions, deletions and insertions that translated into stop codons or changed the open reading frame of the functional gene. We provide evidence suggesting that the retro-transposition event that originated these processed pseudogenes took place before the divergence of the Cervidae and Bovidae families.     

Germ line maintenance of the pseudogene donor pool for somatic immunoglobulin gene conversion in chickens.

Somatic immunoglobulin diversity is generated in avian species by sequential gene conversion of variable (V) gene segments of the immunoglobulin heavy- and light-chain loci during B-cell development. The germ line pools of donor sequence information for somatic V-region gene conversion are found in families of V pseudogenes, located 5' of the single functional V gene of each locus. The sequence relationships among the pseudogenes (psi VL) and functional VL1 gene of the chicken light-chain alleles in three inbred strains were compared to determine the extent of diversity within the germ line pseudogene cluster. Numerous differences were observed. For example, compared with the previously reported CB allele and the G4 allele, the S3 allele contains two intact pseudogenes between psi VL16 and psi VL18. These two adjacent psi VL gene segments (psi VL17a and psi VL17b) could have given rise to the psi VL17 segment of the G4 and CB alleles by homologous recombination. The majority of other sequence polymorphisms among the psi VL alleles appear to be the result of meiotic gene conversion. The incidence of untemplated mutations within psi VL segments is significantly lower than the incidence of mutation within the pseudogene flanking regions. Together with the observations that most psi VL segments have open reading frames and lack stop codons, these data support the hypothesis that the psi VL cluster resembles a functional multigene family maintained by evolutionary selection for its functional role in generating somatic antibody diversity. Meiotic gene conversion events within the psi VL cluster serve both to introduce diversity by the exchange of short segments between family members and to prevent the accumulation of random mutations.     

Repetitive nucleotide sequence insertions into a novel calmodulin-related gene and its processed pseudogene

A gene containing a transposon-like human repeat element, called THE 1, has been isolated and characterized. The gene, termed T+, encodes a polypeptide resembling known calcium-binding proteins. The THE 1 element is present in the 3'-untranslated region of its message. The cDNA clone corresponding to the gene's mRNA product led to the identification of this gene. A processed RNA pseudogene related to the authentic gene has also been isolated. In addition to intron processing, this pseudogene differs from the gene in that it contains an interspersed Alu repeat instead of a THE 1 element in the 3'-untranslated region. Thus, we compare a site containing a THE 1 element to an ancestrally related transposon-less target site. The comparison suggests a retroviral-related mechanism of THE 1 insertion. This system is unusual in that the parent gene is associated with three distinct retrotransposition events: the parent gene was converted to a processed RNA pseudogene, an Alu repeat inserted into the pseudogene, and a THE 1 element inserted into the parent gene.     

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.