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.     

A maximum likelihood method for analyzing pseudogene evolution: implications for silent site evolution in humans and rodents.

We present a new likelihood method for detecting constrained evolution at synonymous sites and other forms of nonneutral evolution in putative pseudogenes. The model is applicable whenever the DNA sequence is available from a protein-coding functional gene, a pseudogene derived from the protein-coding gene, and an orthologous functional copy of the gene. Two nested likelihood ratio tests are developed to test the hypotheses that (1) the putative pseudogene has equal rates of silent and replacement substitutions; and (2) the rate of synonymous substitution in the functional gene equals the rate of substitution in the pseudogene. The method is applied to a data set containing 74 human processed-pseudogene loci, 25 mouse processed-pseudogene loci, and 22 rat processed-pseudogene loci. Using the informatics resources of the Human Genome Project, we localized 67 of the human-pseudogene pairs in the genome and estimated the GC content of a large surrounding genomic region for each. We find that, for pseudogenes deposited in GC regions similar to those of their paralogs, the assumption of equal rates of silent and replacement site evolution in the pseudogene is upheld; in these cases, the rate of silent site evolution in the functional genes is approximately 70% the rate of evolution in the pseudogene. On the other hand, for pseudogenes located in genomic regions of much lower GC than their functional gene, we see a sharp increase in the rate of silent site substitutions, leading to a large rate of rejection for the pseudogene equality likelihood ratio test.     

Eleven daughters of NANOG

Nanog is a recently discovered ANTP class homeobox gene. Mouse Nanog is expressed in the inner cell mass and in embryonic stem cells and has roles in self-renewal and maintenance of pluripotency. Here we describe the location, genomic organization, and relative ages of all human NANOG pseudogenes, comprising ten processed pseudogenes and one tandem duplicate. These are compared to the original, intact human NANOG gene. Eleven is an unusually high number of pseudogenes for a homeobox gene and must reflect expression in the human germ line. A pseudogene orthologous to NANOGP4 was found in chimpanzee and an expressed pseudogene in macaque. Examining pseudogenes of differing ages gives insight into pseudogene decay, which involves an excess of deletion mutations over insertions. The mouse genome has two processed pseudogenes, which are not clear orthologues of the primate pseudogenes.     

A systematic investigation identifies a significant number of probable pseudogenes in the Escherichia coli genome

Pseudogenes are open reading frames (ORFs) encoding dysfunctional proteins with high homology to known protein-coding genes. Although pseudogenes were reported to exist in the genomes of many eukaryotes and bacteria, no systematic search for pseudogenes in the Escherichia coli genome has been carried out. Genome comparisons of E. coli strains K-12 and O157 revealed that many protein-coding sequences have prematurely terminated orthologs encoding unstable proteins. To systematically screen for pseudogenes, we selected ORFs generated by premature termination of the orthologous protein-coding genes and subsequently excluded those possibly arising from sequence errors. Lastly we eliminated those with close homologs in this and other species, as these shortened ORFs may actually have functions. The process produced 95 and 101 pseudogene candidates in K-12 and O157, respectively. The assigned three-dimensional structures suggest that most of the encoded proteins cannot fold properly and thus are dysfunctional, indicating that they are probably pseudogenes. Therefore, the existence of a significant number of probable pseudogenes in E. coli is predicted, awaiting experimental verification. Most of them were found to be genes with paralogs or horizontally transferred genes or both. We suggest that pseudogenes constitute a small fraction of the genomes of free-living bacteria in general, reflecting the faster elimination than production of pseudogenes.     

More reliable estimates of divergence times in Pan using complete mtDNA sequences and accounting for population structure

Here, we report the sequencing and analysis of eight complete mitochondrial genomes of chimpanzees (Pan troglodytes) from each of the three established subspecies (P. t. troglodytes, P. t. schweinfurthii and P. t. verus) and the proposed fourth subspecies (P. t. ellioti). Our population genetic analyses are consistent with neutral patterns of evolution that have been shaped by demography. The high levels of mtDNA diversity in western chimpanzees are unlike those seen at nuclear loci, which may reflect a demographic history of greater female to male effective population sizes possibly owing to the characteristics of the founding population. By using relaxed-clock methods, we have inferred a timetree of chimpanzee species and subspecies. The absolute divergence times vary based on the methods and calibration used, but relative divergence times show extensive uniformity. Overall, mtDNA produces consistently older times than those known from nuclear markers, a discrepancy that is reduced significantly by explicitly accounting for chimpanzee population structures in time estimation. Assuming the human–chimpanzee split to be between 7 and 5 Ma, chimpanzee time estimates are 2.1–1.5, 1.1–0.76 and 0.25–0.18 Ma for the chimpanzee/bonobo, western/(eastern + central) and eastern/central chimpanzee divergences, respectively.     

Pseudogene evolution and natural selection for a compact genome

Pseudogenes are nonfunctional copies of protein-coding genes that are presumed to evolve without selective constraints on their coding function. They are of considerable utility in evolutionary genetics because, in the absence of selection, different types of mutations in pseudogenes should have equal probabilities of fixation. This theoretical inference justifies the estimation of patterns of spontaneous mutation from the analysis of patterns of substitutions in pseudogenes. Although it is possible to test whether pseudogene sequences evolve without constraints for their protein-coding function, it is much more difficult to ascertain whether pseudogenes may affect fitness in ways unrelated to their nucleotide sequence. Consider the possibility that a pseudogene affects fitness merely by increasing genome size. If a larger genome is deleterious--for example, because of increased energetic costs associated with genome replication and maintenance--then deletions, which decrease the length of a pseudogene, should be selectively advantageous relative to insertions or nucleotide substitutions. In this article we examine the implications of selection for genome size relative to small (1-400 bp) deletions, in light of empirical evidence pertaining to the size distribution of deletions observed in Drosophila and mammalian pseudogenes. There is a large difference in the deletion spectra between these organisms. We argue that this difference cannot easily be attributed to selection for overall genome size, since the magnitude of selection is unlikely to be strong enough to significantly affect the probability of fixation of small deletions in Drosophila.     

Comparative analysis of human masculinity

To study rapidly evolving male specific Y (MSY) genes we retrieved and analyzed nine such genes. VCY, HSFY and RBMY were found to have functional X gametologs, but the rest did not. Using chimpanzee orthologs for XKRY, CDY, HSFY, PRY, and TSPY, the average silent substitution is estimated as 0.017 +/- 0.006/site and the substitution rate is 1.42 x 10(-9)/site/year. Except for VCY, all other loci possess two or more pseudogenes on the Y chromosome. Sequence differences from functional genes show that BPY2, DAZ, XKRY, and RBMY each have one pseudogene for each one that is human specific, while others were generated well before the human-chimpanzee split, by means of duplication, retro-transposition or translocation. Some functional MSY gene duplication of VCY, CDY and HSFY, as well as X-linked VCX and HSFX duplication, occurred in the lineage leading to humans; these duplicates have accumulated nucleotide substitutions that permit their identification.     

Structure and evolution of human and African ape rDNA pseudogenes

We discuss the evolutionary significance of four aberrant 18S rDNA clones that were obtained from human, chimpanzee, and gorilla DNA libraries. We show that these clones carry representatives of a small 18S rDNA pseudogene family that arose in a common ancestor of these species. Aspects of their structure and phylogenetic distribution suggest that the 18S pseudogenes no longer interact genetically with normal ribosomal genes and therefore may not be linked to nucleolus organizer regions.      

Phylogenetic analysis and identification of pseudogenes reveal a progressive loss of zona pellucida genes during evolution of vertebrates

Vertebrate eggs are surrounded by an extracellular matrix with similar functions and conserved individual components: the zona pellucida (ZP) glycoproteins. In mammals, chickens, frogs, and some fish species, we established an updated list of the ZP genes, studied the relationships within the ZP gene family using phylogenetic analysis, and identified ZP pseudogenes. Our study confirmed the classification of ZP genes in six subfamilies: ZPA/ZP2, ZPB/ZP4, ZPC/ZP3, ZP1, ZPAX, and ZPD. The identification of a Zpb pseudogene in the mouse genome, Zp1 pseudogenes in the dog and bovine genomes, and Zpax pseudogenes in the human, chimpanzee, macaque, and bovine genomes showed that the evolution of ZP genes mainly occurs by death of genes. Our study revealed that the extracellular matrix surrounding vertebrate eggs contains three to at least six ZP glycoproteins. Mammals can be classified in three categories. In the mouse, the ZP is composed of three ZP proteins (ZPA/ZP2, ZPC/ZP3, and ZP1). In dog, cattle and, putatively, pig, cat, and rabbit, the zona is composed of three ZP proteins (ZPA/ZP2, ZPB/ZP4, and ZPC/ZP3). In human, chimpanzee, macaque, and rat, the ZP is composed of four ZP proteins (ZPA/ZP2, ZPB/ZP4, ZPC/ZP3, and ZP1). Our review provides new directions to investigate the molecular basis of sperm-egg recognition, a mechanism which is not yet elucidated.     

Genomic relationships and speciation times of human, chimpanzee, and gorilla inferred from a coalescent hidden Markov model

The genealogical relationship of human, chimpanzee, and gorilla varies along the genome. We develop a hidden Markov model (HMM) that incorporates this variation and relate the model parameters to population genetics quantities such as speciation times and ancestral population sizes. Our HMM is an analytically tractable approximation to the coalescent process with recombination, and in simulations we see no apparent bias in the HMM estimates. We apply the HMM to four autosomal contiguous human–chimp–gorilla–orangutan alignments comprising a total of 1.9 million base pairs. We find a very recent speciation time of human–chimp (4.1 ± 0.4 million years), and fairly large ancestral effective population sizes (65,000 ± 30,000 for the human–chimp ancestor and 45,000 ± 10,000 for the human–chimp–gorilla ancestor). Furthermore, around 50% of the human genome coalesces with chimpanzee after speciation with gorilla. We also consider 250,000 base pairs of X-chromosome alignments and find an effective population size much smaller than 75% of the autosomal effective population sizes. Finally, we find that the rate of transitions between different genealogies correlates well with the region-wide present-day human recombination rate, but does not correlate with the fine-scale recombination rates and recombination hot spots, suggesting that the latter are evolutionarily transient.     

Identification of protein-coding sequences using the hybridization of 18S rRNA and mRNA during translation

We introduce a new approach in this article to distinguish protein-coding sequences from non-coding sequences utilizing a period-3, free energy signal that arises from the interactions of the 3′-terminal nucleotides of the 18S rRNA with mRNA. We extracted the special features of the amplitude and the phase of the period-3 signal in protein-coding regions, which is not found in non-coding regions, and used them to distinguish protein-coding sequences from non-coding sequences. We tested on all the experimental genes from Saccharomyces cerevisiae and Schizosaccharomyces pombe. The identification was consistent with the corresponding information from GenBank, and produced better performance compared to existing methods that use a period-3 signal. The primary tests on some fly, mouse and human genes suggests that our method is applicable to higher eukaryotic genes. The tests on pseudogenes indicated that most pseudogenes have no period-3 signal. Some exploration of the 3′-tail of 18S rRNA and pattern analysis of protein-coding sequences supported further our assumption that the 3′-tail of 18S rRNA has a role of synchronization throughout translation elongation process. This, in turn, can be utilized for the identification of protein-coding sequences.     

Mouse cytoskeletal gamma-actin: analysis and implications of the structure of cloned cDNA and processed pseudogenes

The nucleotide sequence corresponding to almost the whole of a mouse gamma-cytoskeletal actin mRNA was determined from overlapping cloned DNA copies derived from brain mRNA. Several gamma-actin processed pseudogenes were isolated from a library of cloned DBA mouse genomic DNA, and the nucleotide sequences of these were determined and compared with that of the cDNA. This showed that two of these pseudogenes had arisen from a gene duplication or amplification event, and indicated that they had subsequently undergone partial correction against one another. The relative ages of the pseudogenes were estimated on the basis of their percentage divergence from the cDNA sequence and these were compared with an estimation based on the number of presumed silent mutations in the cDNA since each pseudogene had arisen. Consistent results were obtained, except in the case of one pseudogene which also showed an anomalous regional distribution of differences from the cDNA sequence. One way of accounting for the features of this anomalous pseudogene is by postulating that it is derived from a second functional gene for gamma-actin, different from that represented by the cDNA described here.     

Tubulin pseudogenes as markers for hominoid divergence

Processed pseudogenes arise via unimolecular events that result in the integration of nonfunctional (and therefore non-selected) regions of DNA into the germ line. The sequence of such pseudogenes can be used as a novel form of evolutionary clock: the older a particular pseudogene, the more mutations it has acquired relative to the selectively constrained functional gene from which it was originally derived. We have used specific beta-tubulin gene probes to assay for the presence of fully sequenced processed pseudogenes in genomic DNA from various hominoid species. The data suggest that orangutan is more closely related to human, chimpanzee and gorilla than is generally believed.     

The chimpanzee GH locus: composition, organization, and evolution

In most mammals the growth hormone (GH) locus comprises a single gene expressed primarily in the anterior pituitary gland. However, in higher primates multiple duplications of the GH gene gave rise to a complex locus containing several genes. In man this locus comprises five genes, including GH-N (expressed in pituitary) and four genes expressed in the placenta, but in other species the number and organization of these genes vary. The situation in chimpanzee has been unclear, with suggestions of up to seven GH-like genes. We have re-examined the GH locus in chimpanzee and have deduced the complete sequence. The locus includes five genes apparently organized in a fashion similar to that in human, with two of these genes encoding GH-like proteins, and three encoding chorionic somatomammotropins/placental lactogens (CSHs/PLs). There are notable differences between the human and chimpanzee loci with regard to the expressed proteins, gene regulation, and gene conversion events. In particular, one human gene (hCSH-L) has changed substantially since the chimpanzee/human split, potentially becoming a pseudogene, while the corresponding chimpanzee gene (CSH-A1) has been conserved, giving a product almost identical to the adjacent CSH-A2. Chimpanzee appears to produce two CSHs, with potentially differing biological properties, whereas human produces a single CSH. The pattern of gene conversion in human has been quite different from that in chimpanzee. The region around the GH-N gene in chimpanzee is remarkably polymorphic, unlike the corresponding region in human. The results shed new light on the complex evolution of the GH locus in higher primates.