The Muridae glyceraldehyde-3-phosphate dehydrogenase family

Summary Although only one gene is known to be functional, numerous glyceraldehyde-3-phosphate dehydrogenase (GAPDH) related sequences are scattered throughoutMus musculus andRattus rattus genomes. In this report we show that: (1) GAPDH pseudogenes are repeated to comparable extents, at least 400 copies, in 12 other Muridae species; (2) the complete, or nearly so, sequence of GAPDH messenger RNA is amplified, and a high proportion, if not all of these copies, are intronless; (3) GAPDH pseudogenes are preferentially located in heavily methylated and DNAse I-insensitive regions of chromatin; and (4) the presence of atypical GAPDH-related mRNAs in different cellular contexts raises the possibility that more than one GAPDH gene is transcribed.     

How Many Processed Pseudogenes Are Accumulated in a Gene Family?

A simple kinetic model is developed that describes the accumulation of processed pseudogenes in a functional gene family. Insertion of new pseudogenes occurs at rate ν per gene and is countered by spontaneous deletion (at rate δ per DNA segment) of segments containing processed pseudogenes. If there are k functional genes in a gene family, the equilibrium number of processed pseudogenes is k(ν/δ), and the percentage of functional genes in the gene family at equilibrium is 1/[1 + (ν/δ)]. ν/δ values estimated for five gene families ranged from 1.7 to 15. This fairly narrow range suggests that the rates of formation and deletion of processed pseudogenes may be positively correlated for these families. If δ is sufficiently large relative to the per nucleotide mutation rate µ (δ > 20µ), processed pseudogenes will show high homology with each other, even in the absence of gene conversion between pseudogenes. We argue that formation of processed pseudogenes may share common pathways with transposable elements and retroviruses, creating the potential for correlated responses in the evolution of processed pseudogenes due to direct selection for control of transposable elements and/or retroviruses. Finally, we discuss the nature of the selective forces that may act directly or indirectly to influence the evolution of processed pseudogenes.

Anything produced by evolution is bound to be a bit of a mess—S. Brenner

     

Comparative Analysis of Evolution in a Rodent Histone H2a Pseudogene

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

The MGEA6 Multigene Family Has an Active Locus on 14q and at Least Nine Pseudogenes on Different Chromosomes

The meningioma expressed antigen-6 (MGEA6) was originally identified as an immunogenic antigen in meningioma patients. Somatic hybrid panel mapping and fluorescence in situ hybridization revealed MGEA6-related sequences on different human chromosomes. Here we carry out database analysis to investigate the complexity of the MGEA6-related sequences and demonstrate the existence of a multigene family. We localized the active gene (spanning over 83 kb) to chromosome 14q and elucidated its exon/intron structure. We identified and characterized 9 processed pseudogenes on 9 different chromosomes including chromosomes 2, 3, 6, 7, 9, 10, 12, 13, and 18. We performed phylogenetic analysis and concluded that the MGEA6 pseudogenes may result from more than one retrotransposition event; we calculated divergence times of the pseudogenes to be between 21.5 and 28.9 million years ago.     

Marmot phylogeny revisited: molecular evidence for a diphyletic origin of sociality

We established the phylogeny of 11 species of the genus Marmota based on the entire sequence of the mitochondrial cytochrome b ( cyt-b ) gene (1.1 kb) and a partial sequence of the NADH dehydrogenase subunit 4 ( ND4 ) gene (1.2 kb). In three species ( Marmota caligata , Marmota olympus , and Marmota bobac ) full-sized nuclear pseudogenes of the mitochondrial cyt-b were identified. The mitochondrial cyt-b genes and the three pseudogenes form separate clusters in the maximum parsimony dendrogram. This finding suggests that the pseudogenes originated from a single transfer to the nucleus that may have occurred prior to the radiation of the genus Marmota . Notably, compared with their functional mitochondrial equivalents the pseudogenes show a much lower substitution rate. In the dendrograms deduced from the mitochondrial sequences two distinct clusters become apparent: one cluster consists of the North-west American species, the other contains the Eurasian species together with the North American species Marmota monax . The position of M. monax as a member of the Eurasian clade is in accordance with the evolution of chromosome numbers. The results are of special interest with respect to the evolution of social systems in the genus that vary from solitary species ( M. monax ) to highly social species living in family groups (e.g. Marmota marmota ). The molecular phylogeny suggests a diphyletic origin of high sociality in the genus Marmota .     

人类基因组上的假基因

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

Isolation and characterization of rat and human glyceraldehyde-3-phosphate dehydrogenase cDNAs: genomic complexity and molecular evolution of the gene.

Full length cDNAs encoding the glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from rat and man have been isolated and sequenced. Many GAPDH gene-related sequences have been found in both genomes based on genomic blot hybridization analysis. Only one functional gene product is known. Results from genomic library screenings suggest that there are 300-400 copies of these sequences in the rat genome and approximately 100 in the human genome. Some of these related sequences have been shown to be processed pseudogenes. We have isolated several rat cDNA clones corresponding to these pseudogenes indicating that some pseudogenes are transcribed. Rat and human cDNAs are 89% homologous in the coding region, and 76% homologous in the first 100 base pairs of the 3'-noncoding region. Comparison of these two cDNA sequences with those of the chicken, Drosophila and yeast genes allows the analysis of the evolution of the GAPDH genes in detail.     

Pseudogenization in pathogenic fungi with different host plants and lifestyles might reflect their evolutionary past

Pseudogenes are genes with significant homology to functional genes, but contain disruptive mutations (DMs) leading to the production of non‐ or partially functional proteins. Little is known about pseudogenization in pathogenic fungi with different lifestyles. Here, we report the identification of DMs causing pseudogenes in the genomes of the fungal plant pathogens Botrytis cinerea, Cladosporium fulvum, Dothistroma septosporum, Mycosphaerella fijiensis, Verticillium dahliae and Zymoseptoria tritici. In these fungi, we identified 1740 gene models containing 2795 DMs obtained by an alignment‐based gene prediction method. The contribution of sequencing errors to DMs was minimized by analyses of resequenced genomes to obtain a refined dataset of 924 gene models containing 1666 true DMs. The frequency of pseudogenes varied from 1% to 5% in the gene catalogues of these fungi, being the highest in the asexually reproducing fungus C. fulvum (4.9%), followed by D. septosporum (2.4%) and V. dahliae (2.1%). The majority of pseudogenes do not represent recent gene duplications, but members of multi‐gene families and unitary genes. In general, there was no bias for pseudogenization of specific genes in the six fungi. Single exceptions were those encoding secreted proteins, including proteases, which appeared more frequently pseudogenized in C. fulvum than in D. septosporum. Most pseudogenes present in these two phylogenetically closely related fungi are not shared, suggesting that they are related to adaptation to a different host (tomato versus pine) and lifestyle (biotroph versus hemibiotroph).     

Concerted and birth-and-death evolution of 26S ribosomal DNA in Camellia L.

Background and AimsThe ribosomal DNA (rDNA) gene family, encoding ribosomal RNA (rRNA), has long been regarded as an archetypal example illustrating the model of concerted evolution. However, controversy is arising, as rDNA in many eukaryotic species has been proved to be polymorphic. Here, a metagenomic strategy was applied to detect the intragenomic polymorphism as well as the evolutionary patterns of 26S rDNA across the genus Camellia.MethodsDegenerate primer pairs were designed to amplify the 26S rDNA fragments from different Camellia species. The amplicons were then paired-end sequenced on the Illumina MiSeq platform.Key ResultsAn extremely high level of rDNA polymorphism existed universally in Camellia. However, functional rDNA was still the major component of the family, and was relatively conserved among different Camellia species. Sequence variations mainly came from rRNA pseudogenes and favoured regions that are rich in GC. Specifically, some rRNA pseudogenes have existed in the genome for a long time, and have even experienced several expansion events, which has greatly enriched the abundance of rDNA polymorphism.Conclusions Camellia represents a group in which rDNA is subjected to a mixture of concerted and birth-and-death evolution. Some rRNA pseudogenes may still have potential functions. Conversely, when released from selection constraint, they can evolve in the direction of decreasing GC content and structural stability through a methylation-induced process, and finally be eliminated from the genome.     

Molecular analysis of glyceraldehyde-3-phosphate dehydrogenase in Trypanoplasma borelli: An evolutionary scenario of subcellular compartmentation in Kinetoplastida

In Trypanoplasma borelli, a representative of the Bodonina within the Kinetoplastida, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) activity was detected in both the cytosol and glycosomes. This situation is similar to that previously found in Trypanosomatidae, belonging to a different Kinetoplastida suborder. In Trypanosomatidae different isoenzymes, only distantly related, are responsible for the activity in the two cell compartments. In contrast, immunoblot analysis indicated that the GAPDH activity in cytosol and glycosomes of T. borelli should be attributed to identical or at least very similar proteins related to the glycosomal GAPDH of Trypanosomatidae. Moreover, only genes related to the glycosomal GAPDH genes of Trypanosomatidae could be detected. All attempts to identify a gene related to the one coding for the trypanosomatid cytosolic GAPDH remained unsuccessful. Two tandemly arranged genes were found which are 95% identical. The two encoded polypeptides differ in 17 residues. Their sequences are 72–77% identical to the glycosomal GAPDH of the other Kinetoplastida and share with them some characteristic features: an excess of positively charged residues, specific insertions, and a small carboxy-terminal extension containing the sequence -AKL. This tripeptide conforms to the consensus signal for targeting of proteins to glycosomes. One of the two gene copies has undergone some mutations at positions coding for highly conserved residues of the active site and the NAD⁺-binding domain of GAPDH. Modeling of the protein's three-dimensional structure suggested that several of the substitutions compensate each other, retaining the functional coenzyme-binding capacity, although this binding may be less tight. The presented analysis of GAPDH in T. borelli gives further support to the assertion that one isoenzyme, the cytosolic one, was acquired by horizontal gene transfer during the evolution of the Kinetoplastida, in the lineage leading to the suborder Trypanosomatina (Trypanosome, Leishmania), after the divergence from the Bodonina (Trypanoplasma). Furthermore, the data clearly suggest that the original GAPDH of the Kinetoplastida has been compartmentalized during evolution.Abbreviations GAPDH glyceraldehyde-3-phosphate dehydrogenase (EC 1.2.1.12) - HK hexokinase (EC 2.7.1.1) - PGI glucosephosphate isomerase (EC 5.3.1.9) - PGK phosphoglycerate kinase (EC 2.7.2.3) - PYK pyruvate kinase (EC 2.7.1.40) - TIM triosephosphate isomerase (EC 5.3.1.1) - SDS sodium dodecyl sulfate - SSC saline sodium citrate (0.15 M NaCl, 15 mM sodium citrate, pH 7.0) - MYR millions of years Nucleotide sequence data reported in this paper have been submitted to the EMBL/Genbank/DDBJ nucleotide sequence databases under accession number X74535 Correspondence to: P.A.M. Michels     

Evolution of the NANOG pseudogene family in the human and chimpanzee genomes

Background  

The NANOG gene is expressed in mammalian embryonic stem cells where it maintains cellular pluripotency. An unusually large family of pseudogenes arose from it with one unprocessed and ten processed pseudogenes in the human genome. This article compares the NANOG gene and its pseudogenes in the human and chimpanzee genomes and derives an evolutionary history of this pseudogene family.     

The Old World sparrows (genus Passer) phylogeography and their relative abundance of nuclear mtDNA pseudogenes

The phylogenetic relationships of genus Passer (Old World sparrows) have been studied with species covering their complete world living range. Mitochondrial (mt) cyt b genes and pseudogenes have been analyzed, the latter being strikingly abundant in genus Passer compared with other studied songbirds. The significance of these Passer pseudogenes is presently unclear. The mechanisms by which mt cyt b genes become pseudogenes after nuclear translocation are discussed together with their mode of evolution, i.e., transition/transversion mitochondrial ratio is decreased in the nucleus, as is the constraint for variability at the three codon positions. However, the skewed base composition according to codon position (in 1^st position the percentage is very similar for the four bases, in 2^nd position there are fewer percentage of A and G and more percentage of T, and in 3^rd codon position fewer percentage of G and T and is very rich in A and C) is maintained in the translocated nuclear pseudogenes. Different nuclear internal mechanisms and/or selective pressures must exist for explaining this nuclear/mitochondrial differential DNA base evolutive variability. Also, the phylogenetic usefulness of pseudogenes for defining relationships between closely related lineages is stressed. The analyses suggest that the primitive genus Passer species comes from Africa, the Cape sparrow being the oldest: P. hispaniolensis italiae is more likely conspecific to P. domesticus than to P. hispaniolensis. Also, Passer species are not included within weavers or Estrildinae or Emberizinae, as previously suggested. European and American Emberizinae sparrows are closely related to each other and seem to be the earliest species that radiated among the studied songbirds (all in the Miocene Epoch). Received: 29 November 2000 / Accepted: 22 March 2001     

Molecular characterization of a novel,nuclear-encoded,NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase in plastids of the gymnosperm Pinus sylvestris L.

Angiosperms and algae possess two distinct glyceraldehyde-3-phosphate dehydrogenase (GAPDH) enzymes, an NAD⁺-dependent tetramer involved in cytosolic glycolysis and an NADP⁺-dependent enzyme of the Calvin cycle in chloroplasts. We have found that the gymnosperm Pinus sylvestris possesses, in addition to these, a nuclear-encoded, plastid-specific, NAD⁺-dependent GAPDH, designated GapCp, which has not previously been described from any plant. Several independent full-size cDNAs for this enzyme were isolated which encode a functional transit peptide and mature subunit very similar to that of cytosolic GAPDH of angiosperms and algae. A molecular phylogeny reveals that chloroplast GapCp and cytosolic GapC arose through gene duplication early in chlorophyte evolution. The GapCp gene is expressed as highly as that for GapC in light-grown pine seedlings. These findings suggest that aspects of compartmentalized sugar phosphate metabolism may differ in angiosperms and gymnosperms and furthermore underscore the contributions of endosymbiotic gene transfer and gene duplication to the nuclear complement of genes for enzymes of plant primary metabolism.     

Identification and genetic mapping of the murine gene and 20 related sequences encoding chromosomal protein HMG-17

HMG-17 is an abundant, nonhistone chromosomal protein that binds preferentially to nucleosomal core particles of mammalian chromatin. The human gene for HMG-17 has been localized to Chromosome (Chr) 1p, but the murine gene has not been previously mapped. Here we identify the murine functional gene, Hmg17, from among more than 25 related sequences (probably processed pseudogenes) and show that it is located on mouse Chr 4, in a region known to have conserved linkage relationships with human Chr 1p. We also report the map locations of 20 additional Hmg17-related sequences on mouse Chrs 1, 2, 3, 5, 7, 8, 9, 13, 15, 16, 17, 18, and X. The multiple, dispersed members of the Hmg17 multigene family can be detected efficiently with a single cDNA probe and provide useful markers for genetic mapping studies in mice.     

Evolutionary and functional diversity of green fluorescent proteins in cephalochordates

Green fluorescent protein (GFP) has been widely used as a molecular marker in modern biological research. Before the recent report of one GFP gene in Branchiostoma floridae, GFP family members were cloned only from other two groups of species: Cnidaria and Copepoda. Here we describe the complete GFP gene repertoire of B. floridae which includes 13 functional genes and 2 pseudogenes, representing the largest GFP family found so far. Coupling with nine other GFP sequences from another two species of genus Branchiostoma and the sequences from Cnidaria and Copepoda, we made a deep-level phylogenetic analysis for GFP genes in cephalochordates and found: 1) GFP genes have experienced a divergent evolution in cephalochordates; 2) all amphioxus GFP genes form four main clades on the tree which had diverged before the radiation of the last common ancestor of all extant cephalochordates; 3) GFP genes in amphioxus shared a common ancestor with that in Copepoda rather than being derived from horizontal gene transfer, which indicates that our ancestor was derived from a fluorescent organism and lost this ability after its separation from Cephalochordata, and also makes GFP a rare gene which has a rather unusual evolutionary path. In addition, we also provided evidence indicating that GFP genes have evolved divergent functions by specializing their expression profile, and different fluorescent spectra by changing their emission peaks. These findings spark two interesting issues: what are GFP in vivo functions in cephalochordates and why they are lost in other examined deuterostomes?     

Relaxed purifying selection of rhodopsin gene within a Chinese endemic cavefish genus <Emphasis Type="Italic">Sinocyclocheilus</Emphasis> (Pisces: Cypriniformes)

The cyprinid fish genus Sinocyclocheilus, as the most cavefish rich genus, includes many species showing striking adaptation to caves and convergent reduction or even loss of eyes and pigmentation. RH1 is responsible for dim vision. In order to explore the evolution of RH1 gene in this genus, we sequenced the complete gene from 28 individuals of 16 representative species of Sinocyclocheilus, with cave and surface species included. Phylogenetic analyses supported the monophyly of Sinocyclocheilus and polyphyly of the cave species. Codon models implemented in PAML were used to infer the evolution of RH1. We found that Sinocyclocheilus had a significantly higher evolutionary rate for amino acids than other cyprinid fishes compared, which might be the result of relaxation of purifying selection and could be ascribed to cave habit of this genus. In contrast to previous hypotheses, both cave and surface lineages exhibited a similar rate of molecular evolution, so the RH1 of cave species may still be functional, although these species were highly adapted to cave environment. Two amino acid substitutions (D83G and E122V) that were not reported before were found, which may be useful for site-directed mutagenesis in the future. Handling editor: Christian Sturmbauer     

Olfactory receptors in aquatic and terrestrial vertebrates

In species representing different levels of vertebrate evolution, olfactory receptor genes have been identified by molecular cloning techniques. Comparing the deduced amino-acid sequences revealed that the olfactory receptor gene family of Rana esculenta resembles that of Xenopus laevis, indicating that amphibians in general may comprise two classes of olfactory receptors. Whereas teleost fish, including the goldfish Carassius auratus, possess only class I receptors, the `living fossil' Latimeria chalumnae is endowed with both receptor classes; interestingly, most of the class II genes turned out to be pseudogenes. Exploring receptor genes in aquatic mammals led to the discovery of a large array of only class II receptor genes in the dolphin Stenella Coeruleoalba; however, all of these genes were found to be non-functional pseudogenes. These results support the notion that class I receptors may be specialized for detecting water-soluble odorants and class II receptors for recognizing volatile odorants. Comparing the structural features of both receptor classes from various species revealed that they differ mainly in their extracellular loop 3, which may contribute to ligand specificity. Comparing the number and diversity of olfactory receptor genes in different species provides insight into the origin and the evolution of this unique gene family. Accepted: 29 July 1998     

Evolution of the bovine lysozyme gene family: Changes in gene expression and reversion of function

Recruitment of lysozyme to a digestive function in ruminant artiodactyls is associated with amplification of the gene. At least four of the approximately ten genes are expressed in the stomach, and several are expressed in nonstomach tissues. Characterization of additional lysozymelike sequences in the bovine genome has identified most, if not all, of the members of this gene family. There are at least six stomachlike lysozyme genes, two of which are pseudogenes. The stomach lysozyme pseudogenes show a pattern of concerted evolution similar to that of the functional stomach genes. At least four nonstomach lysozyme genes exist. The nonstomach lysozyme genes are not monophyletic. A gene encoding a tracheal lysozyme was isolated, and the stomach lysozyme of advanced ruminants was found to be more closely related to the tracheal lysozyme than to the stomach lysozyme of the camel or other nonstomach lysozyme genes of ruminants. The tracheal lysozyme shares with stomach lysozymes of advanced ruminants the deletion of amino acid 103, and several other adaptive sequence characteristics of stomach lysozymes. I suggest here that tracheal lysozyme has reverted from a functional stomach lysozyme. Tracheal lysozyme then represents a second instance of a change in lysozyme gene expression and function within ruminants. Correspondence to: D.M. Irwin