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.     

Nonallelic members of the cytochrome c multigene family of the rat may arise through different messenger RNAs

We determined the nucleotide sequences of three nonallelic cytochrome c genes (from recombinant clones Ch4A-RC5, 6 and 8) isolated from the rat cytochrome c gene family. In contrast with a fourth gene (from Ch4A-RC4), which has an intron and correctly encodes rat cytochrome c, these three appear to be pseudogenes and resemble mRNA molecules in two respects: they are all missing the intron of clone 4, and sequence homology with clone 4 in their 3' noncoding regions abruptly ends at two different A-rich tracts reminiscent of poly(A) tails. We also detect three cytochrome c mRNAs of sizes 1400, 1100 and 700 nucleotides in several tissues of the adult rat. The size differences among the mRNAs can be accounted for by length heterogeneity in their 3' noncoding regions. Two of the 3' ends map to the two points where the mRNA-like genes diverge from clone 4 at poly(A) tracts. Furthermore, short direct repeats flank the genes of clones 5, 6 and 8 at the positions where their sequences diverge. The observations suggest that these members of the cytochrome c multigene family may arise through insertion into the genome of DNA copies of cytochrome c mRNAs.     

人类基因组上的假基因

假基因是基因组上与编码基因序列非常相似的非功能性基因组DNA拷贝,一般情况都不被转录,且没有明确生理意义。假基因根据其来源可分为复制假基因和已加工假基因。迄今为止,明确鉴定的人类假基因多为已加工假基因,有8000个之多。在Swiss-Prot/TrEMBL收录的编码蛋白质的将近25500个基因序列中,约10％在基因组中有一个或多个近全长已加工假基因。其余的功能基因都没有已加工假基因。核糖体蛋白基因具有最多数量的已加工假基因,约有l700个(占已加工假基因数的22％),少数基因,如cyclophilinA、肌动蛋白(actin)、角蛋白(keratin)、GAPDH、细胞色素C(cytochromec)和nucleophosmin等则有很多份已加工假基因。总体上讲,假基因在人类染色体上的分布与染色体长度成比例,但已加工假基因在GC含量为41％～46％的染色体区域密度最高。已加工假基因的拷贝数和功能基因在生殖器官中的表达高度一致,说明许多假基因发生在胚胎阶段,另外也和基因中GC含量和基因大小密切相关。假基因的准确鉴定对基因组进化、分子医学研究和医学应用具有重要意义。     

Rate variation during molecular evolution: creationism and the cytochrome c molecular clock

Molecular clocks based upon amino acid sequences in proteins have played a major role in the clarification of evolutionary phylogenies. Creationist criticisms of these methods sometimes rely upon data that might initially seem to be paradoxical. For example, human cytochrome c differs from that of an alligator by 13 amino acids but differs by 14 amino acids from a much more closely related primate, Otolemur garnettii. The apparent anomaly is resolved by taking into consideration the variable substitution rate of cytochrome c, particularly among primates. This paper traces some of the history of extensive research into the topic of rate heterogeneity in cytochrome c including data from cytochrome c pseudogenes.     

Cytochrome c: gene structure, homology and ancestral relationships.

In this paper, the author notes the recommended definition of the word "homology" (i.e., indicating an ancestral relationship) and the recommended stipulation that "evidence for homology should be explicitly laid out". The postulated homology for somatic and testes-specific isozymes of cytochrome c is then examined, using recent data obtained from the study of cytochrome c genes. Consideration is also given to some newer findings of molecular biology and possibilities are considered for various types of change in the genome of an organism. Possible roles of introns, pseudogenes and multigene families are considered. The relationship of testes-specific cytochrome c to somatic cytochrome c is carefully considered from data obtained in experimental studies of genes of these two isozymes. If one assumes that these isozymes arose as a consequence of a gene duplication, data from rat and mouse genes indicate that the testes-specific isozyme has incorporated more amino acid changes than the somatic isozyme since the time of their divergence. However, when the 15 amino acid differences (testes-specific vs. somatic isozyme) are considered, there is virtually no similarity in these 15 positions of the testes-specific isozyme with any of the hypothetical ancestral sequences of the somatic isozyme. Nucleotide differences in cytochrome c genes have been evaluated by comparing genes for the two rodent cytochrome c isozymes to cytochrome c genes of fruit flies, chickens and humans. Comparisons of nucleotide substitution rates in genes for the two cytochrome c isozymes in rodents confirm the conclusions from amino acid sequence comparisons; namely, that more rapid nucleotide changes have occurred in the testes-specific cytochrome c gene, than in the somatic cytochrome c gene. Possible explanations for these findings are considered.     

Chromosomal localization of a cytochrome b5 gene to human chromosome 18 and a cytochrome b5 pseudogene to the X chromosome.

We have isolated cDNA clones that code for human cytochrome b5. Owing to the high degree of evolutionary conservation of cytochrome b5 sequences and the existence of human and rodent cytochrome b5 processed pseudogenes, we were unable to map unambiguously the chromosomal localization of the human gene(s) by Southern blot hybridization of DNA from human-rodent somatic cell hybrids. An alternative approach, based on restriction enzyme digestion of PCR-amplified DNA, enabled us to map the human cytochrome b5 gene(s) to chromosome 18 and one of its processed pseudogenes to the X chromosome. We propose the designations CYB5 and CYB5P1 for the gene and pseudogene loci, respectively.     

Frequency Matrix Approach Demonstrates High Sequence Quality in Avian BARCODEs and Highlights Cryptic Pseudogenes

The accuracy of DNA barcode databases is critical for research and practical applications. Here we apply a frequency matrix to assess sequencing errors in a very large set of avian BARCODEs. Using 11,000 sequences from 2,700 bird species, we show most avian cytochrome c oxidase I (COI) nucleotide and amino acid sequences vary within a narrow range. Except for third codon positions, nearly all (96%) sites were highly conserved or limited to two nucleotides or two amino acids. A large number of positions had very low frequency variants present in single individuals of a species; these were strongly concentrated at the ends of the barcode segment, consistent with sequencing error. In addition, a small fraction (0.1%) of BARCODEs had multiple very low frequency variants shared among individuals of a species; these were found to represent overlooked cryptic pseudogenes lacking stop codons. The calculated upper limit of sequencing error was 8×10(-5) errors/nucleotide, which was relatively high for direct Sanger sequencing of amplified DNA, but unlikely to compromise species identification. Our results confirm the high quality of the avian BARCODE database and demonstrate significant quality improvement in avian COI records deposited in GenBank over the past decade. This approach has potential application for genetic database quality control, discovery of cryptic pseudogenes, and studies of low-level genetic variation.     

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.      

Processed pseudogenes for rat cytochrome c are preferentially derived from one of three alternate mRNAs.

Three cytochrome c mRNAs (1,400, 1,100 and 700 nucleotides) are colinear with RC4, a gene that has introns and correctly encodes cytochrome c. A comparison of RC4 to six nonallelic clones isolated from the rat cytochrome c multigene family demonstrates that all three mRNAs are represented in the genome as processed pseudogenes. Four of the six pseudogenes are derived from the 1,100-nucleotide mRNA, and genomic hybridizations further establish that nearly all of the 30 or so gene family members are also genomic copies of this mRNA despite the equimolar ratio of the three messages in rat tissues. Thus, the surprising multiplicity of cytochrome c sequences in the rat genome is mainly accounted for by the selective use of the 1,100-nucleotide mRNA for the formation of processed pseudogenes. In contrast to 700- and 1,400-nucleotide species which are polyadenylated downstream from AAGUAAA and AAUUAAA, respectively, the 1,100-nucleotide mRNA uses the ubiquitous AAUAAA and also displays a unique stem and loop structure (delta G = -59.4 kJ) centered 37 base pairs upstream from this sequence.     

Characterization of a mouse somatic cytochrome c gene and three cytochrome c pseudogenes.

Mouse contains two functional, but differentially expressed, cytochrome c genes. One of these genes is expressed in all somatic tissues so far examined. The other gene is expressed only in testis and is assumed to be spermatogenesis-specific. The nucleotide sequence of four mouse cytochrome c-like genes has been determined. One of these genes (MC1) contains an intron and encodes a polypeptide sequence identical to the published mouse somatic cytochrome c amino acid sequence. The other three genes can not properly encode a mouse cytochrome c protein and appear to be pseudogenes which have arisen via an insertion into the mouse genome of a cDNA copy of a cytochrome c mRNA molecule.     

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

Nucleotide sequences of immunoglobulin-epsilon pseudogenes in man and apes and their phylogenetic relationships

To understand the phylogenetic relationships between hominoids, the nucleotide sequences of immunoglobulin-epsilon processed pseudogenes from chimpanzee, gorilla and orangutan were determined. The basic structures of these processed pseudogenes agreed with their human counterpart. Although the degrees of nucleotide differences between man and the African apes had no statistical significance, all the analytical data examined supported the theory that chimpanzee is the closest relative of man. This result was consistent with that deduced by our recent qualitative study. Studies on the nucleotide sequences of globin genes have suggested that the molecular clock runs more slowly in hominoids than in non-hominoid primates. According to the present data, however, further retardation of the evolutionary rate was not observed in the human lineage. Assuming that orangutan diverged 14 million years ago and that the evolutionary rate between the orangutan lineage and the lineage leading to the other three species is constant, the divergence dates of chimpanzee and gorilla were estimated to be 4.9(+/- 0.9) and 5.9(+/- 0.9) million years ago, respectively.     

Three pseudogenes for human U13 snRNA belong to class III.

The nucleotide sequences of three pseudogenes for the small nucleolar RNA, U13, were determined from three human DNA clones. The sequences are reported 50 bp 5' and 3' to each gene. These pseudogenes belong to class III because they contain dispersed mismatches when compared to the previously determined U13 RNA sequence, an adenine-rich region at the 3' end, and short imperfect repeats flanking the 5' end of the coding sequence and the 3' end of the adenine-rich region.     

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.     

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.     

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.     

Molecular cloning and characterization of mutant and wild-type human beta-actin genes.

There are more than 20 beta-actin-specific sequences in the human genome, many of which are pseudogenes. To facilitate the isolation of potentially functional beta-actin genes, we used the new method of B. Seed (Nucleic Acids Res. 11:2427-2446, 1983) for selecting genomic clones by homologous recombination. A derivative of the pi VX miniplasmid, pi AN7 beta 1, was constructed by insertion of the 600-base-pair 3' untranslated region of the beta-actin mRNA expressed in human fibroblasts. Five clones containing beta-actin sequences were selected from an amplified human fetal gene library by homologous recombination between library phage and the miniplasmid. One of these clones contained a complete beta-actin gene with a coding sequence identical to that determined for the mRNA of human fibroblasts. A DNA fragment consisting of mostly intervening sequences from this gene was then used to identify 13 independent recombinant copies of the analogous gene from two specially constructed gene libraries, each containing one of the two types of mutant beta-actin genes found in a line of neoplastic human fibroblasts. The amino acid and nucleotide sequences encoded by the unmutated gene predict that a guanine-to-adenine transition is responsible for the glycine-to-aspartic acid mutation at codon 244 and would also result in the loss of a HaeIII site. Detection of this HaeIII polymorphism among the fibroblast-derived clones verified the identity of the beta-actin gene expressed in human fibroblasts.