Comparative Genomics Search for Losses of Long-Established Genes on the Human Lineage

doi:10.1371/journal.pcbi.0030247

PLoS Comput Biol. 2007 December; 3(12): e247.

Published online 2007 December 14. Prepublished online 2007 October 30. doi: 10.1371/journal.pcbi.0030247

PMCID: PMC2134963

Copyright : © 2007 Zhu et al. This is an open-access article distributed under the terms of the Creative Commons Public Domain declaration which stipulates that, once placed in the public domain, this work may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose.

Comparative Genomics Search for Losses of Long-Established Genes on the Human Lineage

Jingchun Zhu,¹ J. Zachary Sanborn,¹ Mark Diekhans,¹ Craig B Lowe,¹ Tom H Pringle,¹ and David Haussler^1,²^*

¹ Center for Biomolecular Science and Engineering, University of California Santa Cruz, Santa Cruz, California, United States of America

² Howard Hughes Medical Institute, University of California Santa Cruz, Santa Cruz, California, United States of America

Gary Stormo, Editor

Washington University, United States of America

* To whom correspondence should be addressed. E-mail: haussler/at/soe.ucsc.edu

Received July 31, 2007; Accepted October 30, 2007.

Abstract

Taking advantage of the complete genome sequences of several mammals, we developed a novel method to detect losses of well-established genes in the human genome through syntenic mapping of gene structures between the human, mouse, and dog genomes. Unlike most previous genomic methods for pseudogene identification, this analysis is able to differentiate losses of well-established genes from pseudogenes formed shortly after segmental duplication or generated via retrotransposition. Therefore, it enables us to find genes that were inactivated long after their birth, which were likely to have evolved nonredundant biological functions before being inactivated. The method was used to look for gene losses along the human lineage during the approximately 75 million years (My) since the common ancestor of primates and rodents (the euarchontoglire crown group). We identified 26 losses of well-established genes in the human genome that were all lost at least 50 My after their birth. Many of them were previously characterized pseudogenes in the human genome, such as GULO and UOX. Our methodology is highly effective at identifying losses of single-copy genes of ancient origin, allowing us to find a few well-known pseudogenes in the human genome missed by previous high-throughput genome-wide studies. In addition to confirming previously known gene losses, we identified 16 previously uncharacterized human pseudogenes that are definitive losses of long-established genes. Among them is ACYL3, an ancient enzyme present in archaea, bacteria, and eukaryotes, but lost approximately 6 to 8 Mya in the ancestor of humans and chimps. Although losses of well-established genes do not equate to adaptive gene losses, they are a useful proxy to use when searching for such genetic changes. This is especially true for adaptive losses that occurred more than 250,000 years ago, since any genetic evidence of the selective sweep indicative of such an event has been erased.

Author Summary

One of the most important questions in biology is to identify the genetic changes underlying evolution, especially those along the lineage leading to the modern human. Although counterintuitive, losing a gene might actually bring a selective advantage to the organism. This type of gene loss is called adaptive gene loss. Although a few cases have been characterized in the literature, this is the first study to address adaptive gene losses on a scale of the whole human genome and a time period of up to 75 million years. The difficulty of identifying adaptive gene losses is in part the large number of pseudogenes in the human genome. To circumvent this problem, we used two methods to enrich the process for the adaptive candidates. The first is a novel approach for pseudogene detection that is highly sensitive in identifying single-copy pseudogenes that bear no apparent sequence homology to any functional human genes. Second, we used the length of time a gene is functional before loss as a proxy for biological importance, which allows us to differentiate losses of long-established genes from mere losses due to functional redundancy after gene duplication.

Introduction

It is intuitive to think that changes leading to increased complexity, adaptation, and intelligence are achieved by the gain and improvement of genetic components such as genes and regulatory elements. However, in certain scenarios, a loss of function can also bring a selective advantage. The best-known examples are losses of cell surface receptors to confer pathogenic resistance, such as the inactivation of the DUFFY gene contributing to malaria resistance [1] and homozygosity for a null allele of chemokine receptor CCR5 conveying resistance to infection by various pathogens, including HIV [2]. In addition, the loss of an existing biological component can open new developmental opportunities. The human-specific loss of a myosin heavy chain isoform expressed in the masticatory muscles has been linked to the weakening of human jaw muscles, possibly allowing the increase of cranial capacity in humans, although this is still quite speculative [3]. Adaptive gene loss is the type of genetic change that leads to better fitness for an organism by inactivating a functional gene. As argued by the “less-is-more” hypothesis, gene losses may be an important engine of evolutionary innovation [4].

In addition to adaptive evolution, gene losses can play an important role in human diseases where conditionally advantageous mutations improve fitness in a particular environment. For example, deleterious mutations affecting hemoglobin and other red blood cell proteins are common in many human populations due to a heterozygote advantage in malaria epidemic environments. This improved fitness comes at a cost for those born with deleterious mutations on both alleles, since the homozygous state causes anemia including sickle cell disease [5–7]. Other human diseases such as glucose-6-phosphate dehydrogenase deficiency [7] and cystic fibrosis [8,9] have also been associated with the heterozygote advantage.

Despite the apparent importance of adaptive gene loss, we know surprisingly little about its contribution and significance at the genomic level and over a broad time scale. Most research on adaptive evolution in mammals focuses on new genes or regulatory elements as well as on modifications to known genes, such as amino acid substitutions [10,11]. With the complete genomes of human and several other mammals including chimp, rhesus, mouse, rat, and dog [12–16], it is now feasible to systematically identify adaptive gene losses in the human lineage through the course of mammalian evolution.

A claim for adaptive genetic change typically requires evidence of DNA signatures indicating directional selection, and is accompanied by the identification of selective pressures acting on the organisms that are consistent with DNA, fossil, or historical evidence. Methods for detecting amino acid or DNA signatures left by natural selection are not generally applicable for identifying adaptive gene loss [17,18]. An inactivated gene is no longer maintained through the forces of natural selection, and secondary mutations begin to accumulate at the neutral rate. Therefore, methods based on sequence conservation or ratio of synonymous versus nonsynonymous mutations are not suitable to detect adaptive gene losses [11,19–21]. Recent adaptive losses can be detected by the distinct DNA signatures left by positive selection; however, those signatures only persist for a narrow evolutionary window of at most 250,000 years [22–24]. To detect adaptive gene losses further back into the evolutionary past, it is reasonable to assume that a nonredundant gene that was functional for a long time and then inactivated is a good candidate for adaptive gene loss. While not every loss of a well-established gene is adaptive, searching for those candidates can be used to enrich for adaptive gene losses.

Gene loss normally leaves behind a pseudogene. However, the vast majority of pseudogenes in a genome did not bring a selective advantage to the organism. Most pseudogenes arise through a gene copying operation of either retrotransposition (reverse-transcribing a processed mRNA back to DNA, which is reinserted in the genome at a different location) [25], or by segmental or tandem duplication of a genomic region [26]. These are called processed or unprocessed pseudogenes, respectively. While processed pseudogenes in general have a single exon and a polyadenine tail, unprocessed pseudogenes typically have multiple exons and preserve the intron–exon structures of the parental gene. The vast majority of processed pseudogenes are “dead on arrival,” due to the lack of complete coding regions or necessary transcription and translation signals in the new genomic location. Even when a functional gene is formed by segmental duplication, one copy often becomes silenced by degenerative mutations due to functional redundancy [27]. In contrast, adaptive gene losses arise from degradation of genes with well-established functions, which often do not have close homologs in the genome. Taking advantage of the genomic signatures left behind by retrotransposition or gene duplication, several genomic surveys identified tens of thousands of pseudogenes in the human genome using sequence homology to a functional parental gene [28–32]. However, because they lack close homologs, many losses of well-established genes were missed by these studies. More importantly, these analyses focused on cataloging pseudogenes in the human genome, but not on addressing whether the pseudogenizations played a role in evolution.

This study identified losses of well-established protein-coding genes in the human lineage since the common ancestor of euarchontoglires (primates, lemurs, tree shrews, rodents, and lagomorphs such as rabbits). We applied a novel comparative genomic method to identify pseudogenes by syntenic mapping of gene structures between the human–mouse–dog trio of genomes. This approach is able to systematically detect the sequence signature left by losses of well-established genes, distinguishing true losses from mere loss of redundant genes following duplication or retrotransposition. Our analysis was able to differentiate the losses of well-established genes from the large background of human pseudogenes. Twenty six losses of well-established genes were identified in the human lineage since the common ancestor of euarchontoglires, approximately 75 million years ago (Mya). Sixteen of those were previously uncharacterized gene losses in the human genome, such as the loss of acyltransferase 3 during great ape evolution.

Results

Detecting Gene Losses Using Gene Structure Conservation

After a mutation inactivates a functional gene, the signature of the intron–exon structure can still be detected for some time before neutral decay erases it from the genome. Based on the observation that mammalian gene structures are typically conserved between species, a gene prediction program called TransMap was developed that exploits the large-scale conservation of gene order and orientation on mammalian chromosomes to map gene structures between genomes. TransMap is essentially a cross-species mRNA alignment grogram (Text S1) that relies upon the “syntenic” alignments produced by the BLASTZ program [33]. TransMap is highly sensitive in detecting gene structures for both genes and pseudogenes (Table S2). Unlike most existing pseudogene detection methods [30,31,34–36], TransMap does not rely on sequence homology to a parental gene from the same genome; therefore, it is well-suited for detecting losses of well-established genes whose functional precursor has not been recently duplicated. To identify gene losses, the mapped coding region is conceptually translated and scanned for ORF-disrupting mutations. A TransMap prediction is a candidate for gene loss if any of the following ORF-disrupting mutations are detected in the mapped coding regions: stop codons, frameshifts, a splice junction that is not GT-AG, or when less than 50% of the coding region can be mapped to the target (Figure S1). To decrease the number of false positives, a valid conceptual translation is required in an outgroup genome. For example, if a mouse mRNA TransMaps from mouse to the dog genome with a valid conceptual translation, but fails to do so from mouse to the human genome, the gene is a candidate for a human gene loss (Figure 1).

Figure 1

Algorithm for Identifying Losses of Well-Established Genes in the Human Lineage Since the Common Ancestor of Euarchontoglires

Using mouse mRNAs as queries and the dog genome as an outgroup, we identified gene losses in the human lineage since the common ancestor of euarchontoglires. We chose the RefSeq database because it is one of the most comprehensive collections of mouse transcripts [37]. Based on the differential prediction status in the dog and human genomes, 1,008 of the 19,541 mouse RefSeq transcripts were identified as potential human losses (Table S3). The candidates were reduced to 90 after filtering out overlaps with GenBank human mRNAs in order to identify pseudogenes where no transcriptional activity in the form of mature mRNA has been observed in humans, and to remove the false positive pseudogenes predicted by TransMap ((Text S1). A visual examination of an alignment of the mouse mRNA sequence, mouse, human, and dog genomic sequences and their three-frame translations further reduced the number of candidates to 72, eliminating 18 that we are not confident represent true pseudogenes in the human genome. The pipeline is illustrated in Figure 1.

Gene Losses in the Human Lineage

Of the 72 candidates, 27 are predicted to be olfactory receptors (ORs) and ten are members of large gene families, such as keratins. These gene families are organized as tandem gene clusters that have experienced copy number changes and/or complex local rearrangements since the common ancestor of euarchontoglires. The dynamics of gene clusters make it difficult to unambiguously discern ortholog/paralog relationships among species and analyze the evolutionary history of the lost genes. Therefore, we focused on the remaining 35 (non-OR, non-cluster) candidates that were confirmed by visual inspection, which we referred to as the “definite” losses in the human lineage. Among them, 21 gene losses are annotated with some biological functions in mouse and 14 have not been characterized functionally (Table 1).

Table 1

Gene Losses in the Human Lineage since the Common Ancestor of Euarchontoglires

The large number of ORs found by this method is consistent with observations that human OR genes experienced a rapid acceleration of pseudogene formation [38,39]. Previous studies have shown evidence that human genes involved in olfaction have a significant tendency to be under positive selection, indicating ORs have undergone directional selection in humans, including by pseudogenization [18,39]. Many previously identified non-olfactory gene losses were confirmed as well. These include gulonolactone (L-) oxidase (GULO), a vitamin C biosynthesis enzyme that is the genetic basis for scurvy [40], and urate oxidase (UOX), an enzyme converting uric acid to allantoin [41]. In addition, a human-specific nonsense mutation was confirmed in an orphan chemoattractant G protein–coupled receptor 33 (GPR33). Additionally, confirmed gene losses included the human-specific loss of cardiotrophin-2 (CTF2) due to a 8 bp deletion [42], cytochrome c oxidase subunit VIIIb (COX8B) with only 40% mapping to the human genome yielding an ORF of eight amino acids [43], and others listed in Table 1, note a.

Our analysis also identified 23 previously uncharacterized losses, of which 21 belong to the definite loss group and two are members of complex gene clusters. One example of a previously unknown definite loss is acyltransferase 3 (ACYL3; NM_177028), identified by the Riken mouse cDNA project [37,44] and whose function in mammals has not been characterized experimentally. We conducted protein profile analysis and determined that this gene has a highly conserved acyltransferase 3 domain (hence we annotated the gene ACYL3) [45]. Further structural modeling (SAM, TMHMM, SignalP) revealed mouse Acyl3 is a multipass transmembrane protein with its C-terminal domain forming a helix bundle. The N-terminal is extracellular and hydrophilic with conserved cysteine residues able to form disulfide bonds [46–49]. The extracellular domain might be involved in cellular response to external signals. Thus, Acyl3 might be a membrane protein with acyltransferase activity or a multipass transporter to pass molecules across the membrane upon external signals [50,51]. Phylogenetically, ACYL3 is ancient and conserved in archaea, bacteria, fungi, worms, flies, and mammals. While numerous copies of ACYL3 are encoded in fly and worm genomes, mammalian genomes have only one copy. A nonsense mutation (TGG to TGA) located in one of the transmembrane helices is shared by the human and chimpanzee genomes. However, the ancestral TGG codon is present in two orangutan trace sequences, and a valid conceptual translation is present in the rhesus genome. To narrow down the timing of the inactivation, we sequenced a PCR product amplified from the corresponding region in a gorilla DNA sample. The sequencing result showed the TGG (W) codon is present in the gorilla genome. Based on this evidence, it appears that the nonsense mutation inactivated ACYL3 after the divergence of gorillas from the human lineage, and before the divergence of humans and chimpanzees (Figure 2). It is intriguing that the last copy of such an ancient enzyme as ACYL3 was lost during the evolution of great apes. Although we do not know the precise evolutionary impact of the loss, its expression pattern in the mouse pituitary gland and developmental abnormalities observed in Drosophila null mutants suggests the loss might be related to development or hormonal regulation [52–54].

Figure 2

A Nonsense Mutation Inactivated ACYL3 During Great Ape Evolution

As shown in Table 1, other previously unknown gene losses include CETN4, a mammalian centrin expressed in ciliated cells including those present in the cerebellum [55], NEPN (nephrocan), an inhibitor of TGF-β signaling pathway [56], and NRADD, a death domain containing membrane protein involved in mediating apoptosis in response to ER stress [57]. It is worth noting that ER-mediated apoptosis triggers a cascade leading to the activation of Caspase 12 (CASP12), a gene that is also lost in humans. The loss of CASP12 is still polymorphic in humans and has been shown to have experienced a recent selective sweep [58,59].

Timing of the Gene Losses

The timing of the gene losses was determined by finding the branch interval that encloses the earliest shared ORF-disrupting mutations between humans and other mammals on a phylogenetic tree. The branch intervals on the human lineage for the 35 definite losses are illustrated on a mammalian phylogeny (Figure 3).

Figure 3

Timing of the Gene Losses in the Human Lineage Since the Common Ancestor of Euarchontoglires, Estimated Based on Shared Mutation Analysis

Using complete genome sequences of human, chimp, rhesus, mouse, and dog, six genes were determined to be human-specific losses, i.e., lost after the divergence of humans and chimpanzees. Ten genes were found to be lost during the period between the human–chimp split and the divergence of old world monkeys from the human lineage. Among the ten genes, seven were observed to have independent ORF-disrupting mutations in the rhesus lineage that are not shared with humans or chimps. Seventeen genes were determined to be lost prior to the divergence of old world monkeys from the human lineage and after the common ancestor of euarchontoglires. Due to insufficient sequence information in the rhesus genome, two genes could only be determined to be lost at some point during the 70 million years (My) prior to the human–chimp split and after the common ancestor of euarchontoglires. To refine the timing of the gene losses, we extracted trace sequences from several additional primates (orangutan, marmoset, tarsier, galago) and tree shrew. Using these trace sequences, we were able to narrow down 50% of gene losses to a much more precise branch on a phylogenetic tree. For example, the timing of gene losses for GUCY2D, NEPN, and others (number 19 to 22) was narrowed down to the period of approximately 25 to 40 Mya, between the dates when old world and new world monkeys split off from the human lineage (Figure 3).

In addition to identifying previously unknown gene losses in the human genome, our analysis refined the timing of several previously known gene losses. For example, SULT1D1, a sulfotransferase, and GSTA4, a glutathione S-transferase alpha 4 isozyme are known to be pseudogenes in human while their mouse orthologs remain functional [60,61]; however, it was unclear when the inactivation occurred. Our analysis discovered that the pseudogenization of GSTA4 and SULT1D1 occurred approximately 14 to 25 Mya, between the dates when orangutans and old world monkeys split off from the human lineage. GULO is known to be inactive in primates with the inactivation dating to some time prior to the separation of apes and old world monkeys (>25 Mya) [62]. Frameshift indels, nonsense mutations, and genomic deletions are observed in the human GULO sequence, indicating an older pseudogene that has experienced numerous secondary mutations. The shared mutation analysis has shown that human, chimp, rhesus, and marmoset share a nonsense mutation, while galago, mouse, and dog share the TGG tryptophan codon. Therefore, the inactivation of GULO occurred before the separation of new world and old world monkeys (>40 Mya).

Losses of Well-Established Genes versus Duplication-Induced Pseudogenes

A significant contribution of this analysis is to differentiate losses of well-established genes from the large background of pseudogenes caused by retrotransposition or formed shortly after segmental or tandem duplication. The method of syntenically mapping gene structures to both a target and outgroup genome is likely to filter out almost all processed pseudogenes. However, TransMap does not fully eliminate those genes that were silenced soon after duplication, which we referred to as duplication-induced pseudogenes. To identify duplication-induced gene losses, we need to determine when the duplication occurred.

Recent segmental duplications can be detected by within-genome sequence homology. If there is a self-alignment chain in the UCSC human genome browser [63] enclosing the gene loss region, the region is determined to have been recently duplicated. We determined when the duplications had occurred by tracing along the human lineage through a seven-species syntenic alignment (human, chimp, rhesus, mouse, rat, dog, opossum) to determine the origin of each duplicate in the best self-alignment (the one with the highest alignment score recorded in the UCSC genome browser). If a gene and its duplicate trace back to a single region in an outgroup genome, the duplication was determined to have occurred on the branch immediately after the outgroup split off from the human lineage. If the gene and its duplicate consistently traced back to different regions through the series of outgroups all the way back to opossum, the duplication was determined to occur prior to the common ancestor of human and opossum. Figure 4A is a schematic illustration of this procedure. The branch of gene duplication is the branch of gene birth (by duplication).

Figure 4

Timing of Gene Birth Is Estimated by Determining Duplication Histories of Genomic Regions Surrounding the Gene Losses

In many cases, there are no detectable self-alignments, indicating an ancient duplication had formed the functional precursor to the pseudogene. We presume the functional precursor of the pseudogene existed prior to the earliest common ancestor of human and the species whose genomic sequence can be aligned to the human exons, therefore providing a lower bound timing of the gene birth event. To narrow down the timing of gene birth on the long branch between dog and opossum, we included scaffold assemblies of the elephant, tenrec, and armadillo genomes. To infer gene birth events that occurred further back in the evolutionary past, we included the chicken genome in the analysis (Figure 4B). Using the above method, the gene birth by duplication branch for the 35 non-OR, non-cluster definite gene losses was determined and is shown in Table 2 (Gene Birth Branch).

Table 2

Losses of Well-Established Genes in the Human Genome, Classified Using Estimated Functional Time Length

Subsequently, we estimated the length of time (in My) a gene remained functional before its pseudogenization using the separation of the gene birth and death branches. Since the timings of both events are estimated as branch intervals, an upper and lower bound estimation of the separation was obtained. Using a 50 My threshold, we classified the candidates based on their functional time lengths as losses of well-established genes, duplication-induced pseudogenes, or undetermined. If the lower estimation of functional time length is greater than 50 My, the candidate is classified as a loss of a well-established gene. If the upper estimation of functional time length is smaller than 50 My, the gene loss is classified as a duplication-induced pseudogene. The gene loss is classified as undetermined if its functional time length overlaps the 50 My threshold.

Table 2 gives the estimated functional time length for the 35 definite losses. Among them, 26 are classified to be losses of well-established genes, which accounted for the majority (74%) of the definite losses. Five are classified as duplication-induced losses—CYP2G1, SORD, S100a15, CXCL7, and UNC93a— (labeled “**” in Tables 1 and and2).2). The remaining four are undetermined. Of these 26 losses of well-established genes, 16 have not been previously characterized as human pseudogenes in the literature. Among these 16, four have been functionally characterized in mouse, which are NRADD, NEPN, CETN4, and GUCY2D. Table S5 describes various subsets constructed using the 35 “definite” losses. All four candidates do not have detectable homologs in the human genome. Most strikingly, NRADD, NEPN, and CETN4 remained functional for more than 300 My before being inactivated.

Discussion

This study presents the first attempt to systematically identify adaptive gene losses in the human genome since the common ancestor of euarchontoglires, approximately 75 Mya. Using losses of well-established genes as the proxy for adaptive gene losses, we focused on identifying a class of pseudogenes that were once functional and retained this function through tens of millions of years of evolution. We confidently identified 26 losses of well-established genes, including 16 that were not previously known in the literature. The highlight of this analysis is the ability to automatically detect losses of genes bearing no significant homology to any functional paralog in the human genome. Their functional precursors had an ancient origin, but enough evolutionary time has elapsed to erase any significant homology with other genes in the human genome. These genes were functioning for hundreds of millions of years and silenced recently within the past 75 My.

It has been proposed that the majority of pseudogenes are either dead-on-arrival [58] or inactivated quickly after duplication [27]. Therefore, it is not surprising that we have identified a much smaller number of pseudogenes as compared to the thousands identified by previous whole genome analysis that aimed to catalog the human genome for unprocessed pseudogenes [30,31,36,64]. We overlapped our results with two well-known pseudogene databases, Yale pseudogene database, composed of mostly various computational predictions [64], and VEGA pseudogene collection, compiled by manual curation [65]. We found limited overlap between the losses identified in Table 1 with both pseudogene sets (Table S4). Only two out of 31 annotated, zero out of 14 hypothetical, and five out of 27 ORs were found by all three analyses. Neither database has GULO, Cardiotropin 2, or many others listed in Table 1 (see note b). A recent genome scan identified 67 human-specific gene losses, including 36 ORs [58]. Excluding ORs, only one out of the six human-specific gene losses identified in Table 1, Gpr33, was also discovered in that study [58]. Another possible overlap is Ugt2b1, which belongs to a tandem cluster of Ugt2B genes on Chromosome 4. The limited overlap in part reflects the difference in methodology used to identify the pseudogenes, but also makes apparent that none of these methods in their present state are able to form the complete set of losses of genes with ancient origins. It also confirms that we have identified some unprocessed pseudogenes derived from functional precursors of ancient origin, where evolution has erased any significant homology to their current functional paralogs.

The gene loss candidates shown in Table 1 are by no means a complete list of losses of well-established genes in the human lineage during the past 75 My. TransMap gene model prediction methodology is not perfect, many factors can introduce prediction errors including uncertainties in sequence alignments, errors generated by the gene model prediction and evaluation procedures, and evolutionary changes of the gene structures across mammalian species (Text S1). For example, the well-known human specific loss of CMAH (a CMP-sialic acid hydroxylase) [66] was not found by this analysis due to the strictness of TransMap gene model predictions, causing a valid CMAH gene model in the dog genome to be excluded because it featured a noncanonical GC-AG splice junction. However, the use of an outgroup genome and the mRNA filter makes the analysis far more likely to produce false negatives than false positives. Several other factors also contribute to this incompleteness. First, our method using human–mouse–dog comparison relied upon well-defined mouse genes to seed the search and valid dog predictions for outgroup confirmation. Problems in either one will return a false negative result. Our analysis missed MYH16 because it is not in mouse RefSeq, which could be due to an independent loss or a misannotation. We further investigated its absence and found that the MYH16 syntenic region is not present in the mouse genome, indicating an independent loss in mouse via genomic deletion. Our analysis required a valid conceptual translation in the dog genome, which may fail to occur due to TransMap prediction errors, sequencing gaps, or an independent loss in dog. However, the chance of producing a valid mapping increases if multiple outgroups, such as the opossum genome [67] or a computationally reconstructed ancestral genome [68], were used and the resultant gene loss predictions were combined. For example, the previously documented human specific loss of Htr5b [69] can be identified using a reconstructed boreoeutherian genome as the outgroup (Haussler lab, unpublished data). Our analysis can also be improved by extending our seed mRNAs to include those from other species and by using multiple outgroup genomes. For example, using chimpanzee MYH16 mRNA as a seed could have found this pseudogene in human.

Our analysis may not identify human polymorphic gene losses. For example, the human-specific loss of CASP12 [58,59] was not identified by our analysis because the latest human genome assembly (NCBI release 36) has the functional allele. Several other human polymorphic losses were also missed by our analysis for the same reason [70,71]. These polymorphic null alleles are potentially crucial to human diseases, e.g., CASP12 in sepsis and CCR5 in HIV infection. Incorporating human EST and mRNA information, as was done by Hahn et al. [70,71], or the human SNP dataset [72], could help our method identify human polymorphic gene losses. Overlapping those alleles with human disease loci, such as those documented in OMIM database [73] or identified by genetic association studies, might lead to the identification of new human disease associated genes. Another factor that may cause the method to overlook gene losses is related to segmental duplication. After a gene is duplicated, both the ancestral copy (the copy in the original genomic context) and the daughter copy (the copy duplicated in the new genomic context) are equally subject to degenerative mutations. Since our analysis evaluates based on the status of the ancestral copy, if evolution silences the daughter copy, it will not be identified by our method. However, this type of false negative is quite limited in our results because it only applies when a segmental duplication occurred after the boreoeutherian common ancestor. Treating the daughter copy in the same way as the ancestral copy will solve this problem, except in the case of a tandem segmental duplication, where it is difficult to distinguish the ancestral copy from the daughter copy.

Among the 26 losses of well-established genes, six were identified to be lost independently in the human and old world monkey lineages (numbers 8, 11, 12, 13, 15, 25 in Table 1). This can be interpreted as a confirmation for adaptive evolution, if we believe that a common selection pressure forced these genes to be lost in separate clades. Other known independent losses such as Caspase15 and Gpr33 seem to confirm this hypothesis [74,75]. An alternative interpretation is that the gene function is no longer needed, such as the loss of GULO in guinea pigs and humans [40]. However, it is also quite probable that the original loss did not occur independently on different lineages, but rather a common mutation that was missed by the analysis might have occurred earlier on a shared ancestor to inactivate the gene. This might have been a mutation in a noncoding region, or a mutation that was erased by secondary mutations such as genomic deletions. For example, a prior, noncoding mutation in any of the six cases we found could have disrupted the transcription, translation, or regulatory signals of the gene in the common ancestor of old world monkeys and apes, rendering the gene effectively inactive at the time that these lineages split. Since the gene is no longer under selective pressure to maintain its integrity, secondary ORF-disrupting mutations could follow, occurring independently in the separate lineages, as observed by our analysis.

To identify genes that are truly lost, we have focused on regions lacking any reported mRNA evidence, including in cell lines derived from cancer cells. A large number of candidates with differential mutational status in the human and dog gene predictions (918 out of 1,008) were filtered out because they overlap with some mRNA evidence in humans. The majority of these are likely to be TransMap prediction errors (Text S1, Table S1). However, some pseudogenes still generate transcripts if the transcription signal is intact, and these would be overlooked by our method. An example of a transcribed pseudogene in the human genome that appears on this list is CATSPER2 (chr15: 41815434–41825788), represented by GenBank mRNA BC066967, and BC047442. The mammalian gene collection annotates it as a transcribed pseudogene. If a pseudogene is transcribed and spliced, its mRNA transcript with ORF-disrupting mutations (i.e., premature stop codon) is targeted and degraded by the cell's RNA surveillance pathway of nonsense mediated decay [76], although this process may not be complete. Only with time will these pseudogenes will be completely silenced at the level of transcription. In addition, studies have shown that occasionally a pseudogene, like Makorin1, not only transcribes but also plays a vital biological role in stabilizing the mRNA of its homologous coding gene [77]. Thus it is difficult to prove that a transcribed pseudogene is completely nonfunctional.

Theories of molecular evolution suggest three outcomes for new genes arising from gene duplication: degeneration due to functional redundancy, evolution into a new function, or function sharing by both copies [27]. The expected time that elapses before a gene is inactivated is thought to be relatively short [27]. Lynch and Conery estimated the half-life of a new duplicate to be around 16 My in the human lineage [27,78,79]. Using this estimate, after our cutoff of 50 My, 11% of redundant genes caused by duplications are expected to be intact by chance. After 60 My (the shortest estimation that passes the cutoff in Table 2), only 7.5% will be left. Twenty-six candidates in Table 2 are classified as losses of established genes using the 50 My cutoff, and many have an estimated functional period after duplication that is much longer than 50 My. This suggests that they are likely to have evolved independent functions before pseudogenization and thus likely to be true losses of well-established genes. In addition, our method used the lower-bound estimation for the functional time length for this classification. Although the higher-bound estimations for four candidates (PFPL, ABCA14, LOC344492, BC018465 in Table 2) satisfy the 50 My cutoff, their low estimations do not. As complete genome sequences for additional mammals become available in the future, the timing of duplication and pseudogenization can be greatly refined, potentially classifying some of these four candidates as losses of established genes as well.

It is nontrivial to determine whether these losses we have found were truly adaptive. It is very likely that neutral losses at dispensable loci account for a subset of our results. For example, GULO, a vitamin C biosynthesis gene, is thought to have been lost in primates because primates have ample dietary supply of ascorbic acid, reducing or removing the selective pressure that maintains this gene. In general, it is difficult to differentiate between neutral loss due to removal of selective pressure, as proposed by the “use it or lose it hypothesis,” and positively selected adaptive loss, as by the “less is more” hypothesis, without knowing the gene's precise biological functions. Given our current knowledge of human genes, identifying the losses of established genes seem to be the best strategy in the search for more ancient (before 250 Mya) adaptive gene losses on a genomic scale. The resulting list is a much more enriched set of candidates.

In summary, our analysis identified a set of losses that are highly enriched for well-established genes in the human genome against a large background of pseudogenes. Expanding these results to include genes and genomes from the entire mammalian clade will generate a more accurate and comprehensive picture of adaptive gene losses in human evolution. From a theoretical standpoint, it will provide insight into the role that loss of functional genes plays in evolutionary adaptation [4]. The method presented here can also be generalized to discover gene losses in other organisms on a genomic scale.

Materials and Methods

Mapping mouse mRNAs to the human and dog genomes.

BLAT [80] alignments of mouse mRNAs sequences from Genbank [81] and RefSeq [37] to the cognate genome were obtained from the UCSC Genome Browser Database [82]. These alignments, along with the coding sequence annotations associated with the mRNAs, provide annotations of gene structure in the genome. BLASTZ [33] syntenic chained alignments [63] of cognate genome (mouse) to the target genome (human or dog) are used to project the mouse mRNA alignments to the target genome. This algorithm, known as TransMap and illustrated in Figure 1, results in predictions of orthologous gene structures in the target organism. The TransMap prediction methodology is described in detail in Text S1.

Sources of sequences, mRNAs, and pseudogene databases.

Genomic sequences used in this study were obtained from the UCSC genome browser [82]. Sequence release of the following species are human (NCBI release 36.1; UCSC hg18 March 2006), chimp (Chimpanzee Sequencing and Analysis Consortium Build 2 version 1; UCSC March 2006), rhesus macaque (BCM HGSC version 1.0, Mmul_051212; UCSC Jan 2006), mouse (NCBI release 36; UCSC Feb 2006), dog (Broad Institute assembly version 2.0; UCSC May 2005), rat (Baylor Human Genome Sequencing Center HGSC version 3.4; UCSC Nov 2004), opossum (draft assembly produced by the Broad Institute; UCSC Jan 2006), chicken (version 2.1 draft assembly produced by the Genome Sequencing Center at the Washington University; UCSC May 2006), elephant (Broad Institute version 1.0; UCSC May 2005), tenrec (Broad Institute echTel 1.0; UCSC Jul 2005), and armadillo (Broad Institute version 1.0; UCSC May 2005). Trace sequences of orangutan, marmoset, tarsier, galago, and tree shrew were downloaded from NCBI Trace Archive. Mouse RefSeq genes were obtained from the UCSC mouse genome browser, which is consistent with the NCBI mouse genome build 36.1. Yale and VEGA pseudogene datasets were also obtained from the UCSC human genome browser track “Yale Pseudo” and track “Vega Pseudogenes.” Human mRNAs filter is the GenBank human mRNAs, obtained from the UCSC human genome browser track “Human mRNAs”. The identifiers of the two orangutan trace sequences used to confirm the TGG (W) codon is present in the orangutan genome are ti865941905 and ti1012155976 in the NCBI trace archive.

ACYL3 domain analysis and structural modeling.

The mouse gene ACYL3 (NM_177028) is predicted to contain an acyltransferase 3 (acyl3, IPR002656) and a nose resistant to fluoxetine-4 (NRF, IPR006621) domain in the InterPro database [45]. We predicted ACYL3 to have eight or nine transmembrane helixes in its C-terminal sequence using TMHMM v2.0 [48] and to have a 20 amino acid signal peptide in its N-terminal sequence by SignalP 3.0 [49]. We performed protein structure modeling using SAM-T05 [46]. SAM (dssp-eh12 model) predicted ACYL3 to be a multi-pass membrane protein with four conserved cysteine residues in its N-terminal sequence and 12 helixes in its C-terminal sequence. The nonsense stop codon shared by humans and chimpanzees is located in the tenth helix of the structural prediction.

DNA amplification and sequencing of ACYL3 region in gorilla.

Genomic DNA surrounding the ACYL3 TGG (W) codon was PCR-amplified from a gorilla DNA sample. Degenerate PCR primers were designed based on the conservation in human, chimp, and rhesus genomic sequences. The forward primer is 5′-GGTCACCCTATTTGCGGTGGCCGCTCTGGCATACA-3′ and the reverse primer is 5′-TGGGCTGGGTCCTCTTTGCGTGCCACNGAGGATATGGAGGTATGGA-3′. The PCR product was sequenced in both forward and reverse directions, and results were combined to generate a 162 bp gorilla sequence.

Determining earliest shared mutations.

The earliest shared mutations were determined by examining an alignment of genomic sequences of human, chimp, rhesus, mouse, and dog, plus 200 bp trace sequences surrounding the mutation site, when available in the NCBI trace archive, from orangutan, marmoset, tarsier, galago, and tree shrew. Trace sequence analysis was limited to candidates with only point mutations (i.e., stop codons, frameshift indels, or noncanonical splice sites). Candidates with coverage problems (>50% missing) were only analyzed to check whether the limited coverage is also shared by the chimp or rhesus genomes, and their point mutations were not analyzed except for a stop codon mutation in GULO.

Supporting Information

Figure S1: A Histogram of the Number of TransMap Gene Models Predicted in the Human Genome against the Fraction of the Length of the Mouse mRNA Coding Sequence That Can Be Aligned

(21 KB PDF)

Click here for additional data file.^{(22K, pdf)}

Table S1: Evaluation of the TransMap Gene Model Classifications against ENCODE Gene Annotations

(52 KB DOC)

Click here for additional data file.^{(53K, doc)}

Table S2: Comparison of Sequence Coverage by TransMap and Other Sequence Aligners

(22 KB DOC)

Click here for additional data file.^{(23K, doc)}

Table S3: The Percentage of Alignments or Gene Models with a “Valid” Code Assignment

(20 KB DOC)

Click here for additional data file.^{(20K, doc)}

Table S4: Genomic Coordinates of the Lost Genes and Their Overlap with YALE and VEGA Pseudogenes

(58 KB DOC)

Click here for additional data file.^{(58K, doc)}

Table S5: Description of the Various Subsets Constructed Using the 35 “Definite Losses”

(27 KB DOC)

Click here for additional data file.^{(27K, doc)}

Text S1: Description of the TransMap Gene Prediction Methodology

(49 KB DOC)

Click here for additional data file.^{(50K, doc)}

Acknowledgments

We thank Ting Wang for critical reading of the manuscript and scientific discussions. We thank Sofie Salama for helping to carry out the Acyl3 sequencing experiment. We thank Webb Miller for the BLASTZ alignments.

Abbreviations


My	million years

Mya	million years ago

OR	olfactory receptor

Footnotes

A previous version of this article appeared as an Early Online Release on October 30, 2007 (doi:10.1371/journal.pcbi.0030247.eor).

Author contributions. JZ and DH conceived and designed the experiments. JZ, JZS, MD, and CBL performed the experiments. JZ, JZS, and THP analyzed the data. JZ, JZS, and MD contributed reagents/materials/analysis tools. JZ, JZS, and DH wrote the paper.

Funding. This work is supported by the National Human Genome Research Institute (JZ, CBL), the NIH Training Grant T32 GM070386 (JZS), the National Cancer Institute NO1-CO-12400 (MD), and the Howard Hughes Medical Institute (DH).

Competing interests. The authors have declared that no competing interests exist.

References

Tournamille C, Colin Y, Cartron JP, Le Van Kim C. Disruption of a GATA motif in the Duffy gene promoter abolishes erythroid gene expression in Duffy-negative individuals. Nat Genet. 1995;10:224–228. [PubMed]
Dean M, Carrington M, Winkler C, Huttley GA, Smith MW, et al. Genetic restriction of HIV-1 infection and progression to AIDS by a deletion allele of the CKR5 structural gene. Hemophilia Growth and Development Study, Multicenter AIDS Cohort Study, Multicenter Hemophilia Cohort Study, San Francisco City Cohort, ALIVE Study. Science. 1996;273:1856–1862. [PubMed]
Stedman HH, Kozyak BW, Nelson A, Thesier DM, Su LT, et al. Myosin gene mutation correlates with anatomical changes in the human lineage. Nature. 2004;428:415–418. [PubMed]
Olson MV. When less is more: gene loss as an engine of evolutionary change. Am J Hum Genet. 1999;64:18–23. [PubMed]
Ringelhann B, Hathorn MK, Jilly P, Grant F, Parniczky G. A new look at the protection of hemoglobin AS and AC genotypes against plasmodium falciparum infection: a census tract approach. Am J Hum Genet. 1976;28:270–279. [PubMed]
Tishkoff SA, Varkonyi R, Cahinhinan N, Abbes S, Argyropoulos G, et al. Haplotype diversity and linkage disequilibrium at human G6PD: recent origin of alleles that confer malarial resistance. Science. 2001;293:455–462. [PubMed]
Ruwende C, Khoo SC, Snow RW, Yates SN, Kwiatkowski D, et al. Natural selection of hemi- and heterozygotes for G6PD deficiency in Africa by resistance to severe malaria. Nature. 1995;376:246–249. [PubMed]
Schroeder SA, Gaughan DM, Swift M. Protection against bronchial asthma by CFTR delta F508 mutation: a heterozygote advantage in cystic fibrosis. Nat Med. 1995;1:703–705. [PubMed]
Gabriel SE, Brigman KN, Koller BH, Boucher RC, Stutts MJ. Cystic fibrosis heterozygote resistance to cholera toxin in the cystic fibrosis mouse model. Science. 1994;266:107–109. [PubMed]
Lai CS, Fisher SE, Hurst JA, Vargha-Khadem F, Monaco AP. A forkhead-domain gene is mutated in a severe speech and language disorder. Nature. 2001;413:519–523. [PubMed]
Pollard KS, Salama SR, Lambert N, Lambot MA, Coppens S, et al. An RNA gene expressed during cortical development evolved rapidly in humans. Nature. 2006;443:167–172. [PubMed]
Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005;437:69–87. [PubMed]
Gibbs RA, Rogers J, Katze MG, Bumgarner R, Weinstock GM, et al. Evolutionary and biomedical insights from the rhesus macaque genome. Science. 2007;316:222–234. [PubMed]
Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, et al. Initial sequencing and comparative analysis of the mouse genome. Nature. 2002;420:520–562. [PubMed]
Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005;438:803–819. [PubMed]
Gibbs RA, Weinstock GM, Metzker ML, Muzny DM, Sodergren EJ, et al. Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature. 2004;428:493–521. [PubMed]
Nielsen R. Molecular signatures of natural selection. Annu Rev Genet. 2005;39:197–218. [PubMed]
Clark AG, Glanowski S, Nielsen R, Thomas PD, Kejariwal A, et al. Inferring nonneutral evolution from human–chimp–mouse orthologous gene trios. Science. 2003;302:1960–1963. [PubMed]
Bejerano G, Pheasant M, Makunin I, Stephen S, Kent WJ, et al. Ultraconserved elements in the human genome. Science. 2004;304:1321–1325. [PubMed]
Felsenstein J. Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981;17:368–376. [PubMed]
Yang Z, Nielsen R. Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages. Mol Biol Evol. 2002;19:908–917. [PubMed]
Sabeti PC, Schaffner SF, Fry B, Lohmueller J, Varilly P, et al. Positive natural selection in the human lineage. Science. 2006;312:1614–1620. [PubMed]
Kreitman M. Methods to detect selection in populations with applications to the human. Annu Rev Genomics Hum Genet. 2000;1:539–559. [PubMed]
Bustamante CD, Wakeley J, Sawyer S, Hartl DL. Directional selection and the site-frequency spectrum. Genetics. 2001;159:1779–1788. [PubMed]
Vanin EF. Processed pseudogenes: characteristics and evolution. Annu Rev Genet. 1985;19:253–272. [PubMed]
Mighell AJ, Smith NR, Robinson PA, Markham AF. Vertebrate pseudogenes. FEBS Lett. 2000;468:109–114. [PubMed]
Lynch M, Conery JS. The evolutionary fate and consequences of duplicate genes. Science. 2000;290:1151–1155. [PubMed]
Zhang Z, Harrison PM, Liu Y, Gerstein M. Millions of years of evolution preserved: a comprehensive catalog of the processed pseudogenes in the human genome. Genome Res. 2003;13:2541–2558. [PubMed]
Zhang Z, Carriero N, Gerstein M. Comparative analysis of processed pseudogenes in the mouse and human genomes. Trends Genet. 2004;20:62–67. [PubMed]
Suyama M, Harrington E, Bork P, Torrents D. Identification and analysis of genes and pseudogenes within duplicated regions in the human and mouse genomes. PLoS Comput Biol. 2006. e76 doi: 10.1371/journal.pcbi.0020076. [PMC free article] [PubMed]
Torrents D, Suyama M, Zdobnov E, Bork P. A genome-wide survey of human pseudogenes. Genome Res. 2003;13:2559–2567. [PubMed]
Khelifi A, Duret L, Mouchiroud D. HOPPSIGEN: a database of human and mouse processed pseudogenes. Nucleic Acids Res. 2005;33:D59–D66. [PMC free article] [PubMed]
Schwartz S, Kent WJ, Smit A, Zhang Z, Baertsch R, et al. Human–mouse alignments with BLASTZ. Genome Res. 2003;13:103–107. [PubMed]
Zhang Z, Carriero N, Zheng D, Karro J, Harrison PM, et al. PseudoPipe: an automated pseudogene identification pipeline. Bioinformatics. 2006;22:1437–1439. [PubMed]
Zheng D, Gerstein MB. A computational approach for identifying pseudogenes in the ENCODE regions. Genome Biol. 2006;7(Supplement 1) [PMC free article] [PubMed]
Harrison PM, Hegyi H, Balasubramanian S, Luscombe NM, Bertone P, et al. Molecular fossils in the human genome: identification and analysis of the pseudogenes in chromosomes 21 and 22. Genome Res. 2002;12:272–280. [PubMed]
Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35:D61–D65. [PubMed]
Gilad Y, Man O, Paabo S, Lancet D. Human specific loss of olfactory receptor genes. Proc Natl Acad Sci U S A. 2003;100:3324–3327. [PubMed]
Gilad Y, Bustamante CD, Lancet D, Paabo S. Natural selection on the olfactory receptor gene family in humans and chimpanzees. Am J Hum Genet. 2003;73:489–501. [PubMed]
Nishikimi M, Fukuyama R, Minoshima S, Shimizu N, Yagi K. Cloning and chromosomal mapping of the human nonfunctional gene for L-gulono-gamma-lactone oxidase, the enzyme for L-ascorbic acid biosynthesis missing in man. J Biol Chem. 1994;269:13685–13688. [PubMed]
Wu XW, Lee CC, Muzny DM, Caskey CT. Urate oxidase: primary structure and evolutionary implications. Proc Natl Acad Sci U S A. 1989;86:9412–9416. [PubMed]
Derouet D, Rousseau F, Alfonsi F, Froger J, Hermann J, et al. Neuropoietin, a new IL-6-related cytokine signaling through the ciliary neurotrophic factor receptor. Proc Natl Acad Sci U S A. 2004;101:4827–4832. [PubMed]
Goldberg A, Wildman DE, Schmidt TR, Huttemann M, Goodman M, et al. Adaptive evolution of cytochrome c oxidase subunit VIII in anthropoid primates. Proc Natl Acad Sci U S A. 2003;100:5873–5878. [PubMed]
Okazaki Y, Furuno M, Kasukawa T, Adachi J, Bono H, et al. Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs. Nature. 2002;420:563–573. [PubMed]
Apweiler R, Attwood TK, Bairoch A, Bateman A, Birney E, et al. InterPro—an integrated documentation resource for protein families, domains and functional sites. Bioinformatics. 2000;16:1145–1150. [PubMed]
Krogh A, Brown M, Mian IS, Sjolander K, Haussler D. Hidden Markov models in computational biology. Applications to protein modeling. J Mol Biol. 1994;235:1501–1531. [PubMed]
Karchin R, Hughey R. Weighting hidden Markov models for maximum discrimination. Bioinformatics. 1998;14:772–782. [PubMed]
Krogh A, Larsson B, von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–580. [PubMed]
Bendtsen JD, Nielsen H, von Heijne G, Brunak S. Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004;340:783–795. [PubMed]
Choy RK, Thomas JH. Fluoxetine-resistant mutants in C. elegans define a novel family of transmembrane proteins. Mol Cell. 1999;4:143–152. [PubMed]
Choy RK, Kemner JM, Thomas JH. Fluoxetine-resistance genes in Caenorhabditis elegans function in the intestine and may act in drug transport. Genetics. 2006;172:885–892. [PubMed]
Lein ES, Hawrylycz MJ, Ao N, Ayres M, Bensinger A, et al. Genome-wide atlas of gene expression in the adult mouse brain. Nature. 2007;445:168–176. [PubMed]
Dzitoyeva S, Dimitrijevic N, Manev H. Identification of a novel Drosophila gene, beltless, using injectable embryonic and adult RNA interference (RNAi) BMC Genomics. 2003;4:33. [PMC free article] [PubMed]
Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, et al. A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci U S A. 2004;101:6062–6067. [PubMed]
Gavet O, Alvarez C, Gaspar P, Bornens M. Centrin4p, a novel mammalian centrin specifically expressed in ciliated cells. Mol Biol Cell. 2003;14:1818–1834. [PMC free article] [PubMed]
Mochida Y, Parisuthiman D, Kaku M, Hanai J, Sukhatme VP, et al. Nephrocan, a novel member of the small leucine-rich repeat protein family, is an inhibitor of transforming growth factor-beta signaling. J Biol Chem. 2006;281:36044–36051. [PubMed]
Wang X, Shao Z, Zetoune FS, Zeidler MG, Gowrishankar K, et al. NRADD, a novel membrane protein with a death domain involved in mediating apoptosis in response to ER stress. Cell Death Differ. 2003;10:580–591. [PubMed]
Wang X, Grus WE, Zhang J. Gene losses during human origins. PLoS Biol. 2006. e52 doi: 10.1371/journal.pbio.0040052. [PMC free article] [PubMed]
Xue Y, Daly A, Yngvadottir B, Liu M, Coop G, et al. Spread of an inactive form of caspase-12 in humans is due to recent positive selection. Am J Hum Genet. 2006;78:659–670. [PubMed]
Meinl W, Glatt H. Structure and localization of the human SULT1B1 gene: neighborhood to SULT1E1 and a SULT1D pseudogene. Biochem Biophys Res Commun. 2001;288:855–862. [PubMed]
Morel F, Rauch C, Coles B, Le Ferrec E, Guillouzo A. The human glutathione transferase alpha locus: genomic organization of the gene cluster and functional characterization of the genetic polymorphism in the hGSTA1 promoter. Pharmacogenetics. 2002;12:277–286. [PubMed]
Ohta Y, Nishikimi M. Random nucleotide substitutions in primate nonfunctional gene for L-gulono-gamma-lactone oxidase, the missing enzyme in L-ascorbic acid biosynthesis. Biochim Biophys Acta. 1999;1472:408–411. [PubMed]
Kent WJ, Baertsch R, Hinrichs A, Miller W, Haussler D. Evolution's cauldron: duplication, deletion, and rearrangement in the mouse and human genomes. Proc Natl Acad Sci U S A. 2003;100:11484–11489. [PubMed]
Karro JE, Yan Y, Zheng D, Zhang Z, Carriero N, et al. Pseudogene.org: a comprehensive database and comparison platform for pseudogene annotation. Nucleic Acids Res. 2007;35:D55–D60. [PubMed]
Ashurst JL, Chen CK, Gilbert JG, Jekosch K, Keenan S, et al. The Vertebrate Genome Annotation (Vega) database. Nucleic Acids Res. 2005;33:D459–D465. [PMC free article] [PubMed]
Chou HH, Hayakawa T, Diaz S, Krings M, Indriati E, et al. Inactivation of CMP-N-acetylneuraminic acid hydroxylase occurred prior to brain expansion during human evolution. Proc Natl Acad Sci U S A. 2002;99:11736–11741. [PubMed]
Mikkelsen TS, Wakefield MJ, Aken B, Amemiya CT, Chang JL, et al. Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences. Nature. 2007;447:167–177. [PubMed]
Blanchette M, Green ED, Miller W, Haussler D. Reconstructing large regions of an ancestral mammalian genome in silico. Genome Res. 2004;14:2412–2423. [PubMed]
Grailhe R, Grabtree GW, Hen R. Human 5-HT(5) receptors: the 5-HT(5A) receptor is functional but the 5-HT(5B) receptor was lost during mammalian evolution. Eur J Pharmacol. 2001;418:157–167. [PubMed]
Hahn Y, Lee B. Human-specific nonsense mutations identified by genome sequence comparisons. Hum Genet. 2006;119:169–178. [PubMed]
Hahn Y, Lee B. Identification of nine human-specific frameshift mutations by comparative analysis of the human and the chimpanzee genome sequences. Bioinformatics. 2005;21(Supplement 1):i186–194. [PubMed]
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29:308–311. [PMC free article] [PubMed]
Hamosh A, Scott AF, Amberger JS, Bocchini CA, McKusick VA. Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders. Nucleic Acids Res. 2005;33:D514–D517. [PMC free article] [PubMed]
Rompler H, Schulz A, Pitra C, Coop G, Przeworski M, et al. The rise and fall of the chemoattractant receptor GPR33. J Biol Chem. 2005;280:31068–31075. [PubMed]
Eckhart L, Uthman A, Sipos W, Tschachler E. Genome sequence comparison reveals independent inactivation of the caspase-15 gene in different evolutionary lineages of mammals. Mol Biol Evol. 2006;23:2081–2089. [PubMed]
Lewis BP, Green RE, Brenner SE. Evidence for the widespread coupling of alternative splicing and nonsense-mediated mRNA decay in humans. Proc Natl Acad Sci U S A. 2003;100:189–192. [PubMed]
Hirotsune S, Yoshida N, Chen A, Garrett L, Sugiyama F, et al. An expressed pseudogene regulates the messenger-RNA stability of its homologous coding gene. Nature. 2003;423:91–96. [PubMed]
Zhang L, Gaut BS, Vision TJ. Gene duplication and evolution. Science. 2001;293:1551. [PubMed]
Long M, Thornton K. Gene duplication and evolution. Science. 2001;293:1551. [PubMed]
Kent WJ. BLAT—the BLAST-like alignment tool. Genome Res. 2002;12:656–664. [PubMed]
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL. GenBank: update. Nucleic Acids Res. 2004;32:D23–D26. [PMC free article] [PubMed]
Karolchik D, Baertsch R, Diekhans M, Furey TS, Hinrichs A, et al. The UCSC Genome Browser Database. Nucleic Acids Res. 2003;31:51–54. [PMC free article] [PubMed]

Articles from PLoS Computational Biology are provided here courtesy of
Public Library of Science