Evolution of human cytoplasmic actin gene sequences: chromosome mapping and structural characterizations of three cytoplasmic actin-like pseudogenes including one on the Y chromosome

Eight recombinant phage clones containing cytoplasmic actin-like gene sequences have been isolated from a human genomic library for structural characterization. Kpn I family repeat sequences flank six of these actin genes isolated, and Alu family repeats are scattered throughout the DNA inserts of all eight phage clones. Three of these genes are γ actin-like, and the other five are β actin-like. The complete nucleotide sequence analysis of one β and one γ actin-like genes and their flanking regions demonstrates that they both are processed pseudogenes. Using unique DNA sequences flanking these two pseudogenes as hybridization probes for human-mouse somatic cell hybrid DNAs, we have mapped the two actin pseudogenes on human chromosomes 8 and 3, respectively. We have also determined the DNA sequence of a human Y chromosome-linked, processed actin pseudogene. The different values of sequence divergence of these processed pseudogenes and their functional counterparts allow us to estimate the time of generation of the pseudogenes. The results suggest that the cDNA insertion events generating the human cytoplasmic actin-like pseudogenes have occurred at significantly different times during the evolution of primates, after their separation from other mammalian species.     

Juxta-centromeric region of human chromosome 21 is enriched for pseudogenes and gene fragments

A physical map including four pseudogenes and 10 gene fragments and spanning 500 kb in the juxta-centromeric region of the long arm of human chromosome 21 is presented. cDNA fragments isolated from a selected cDNA library were characterized and mapped to the 831B6 YAC and to two BAC contigs that cover 250 kb of the region. An 85 kb genomic sequence located in the proximal region of the map was analyzed for putative exons. Four pseudogenes were found, including psiIGSF3, psiEIF3, psiGCT-rel whose functional copies map to chromosome 1p13, chromosome 2 and chromosome 22q11, respectively. The TTLL1 pseudogene corresponds to a new gene whose functional copy maps to chromosome 22q13. Ten gene fragments represent novel sequences that have related sequences on different human chromosomes and show 97-100% nucleotide identity to chromosome 21. These may correspond to pseudogenes on chromosome 21 and to functional genes in other chromosomes. The 85 kb genomic sequence was analyzed also for GC content, CpG islands, and repetitive sequence distribution. A GC-poor L isochore spanning 40 kb from satellite 1 was observed in the most centromeric region, next to a GC-rich H isochore that is a candidate region for the presence of functional genes. The pericentric duplication of a 7.8 kb region that is derived from the 22q13 chromosome band is described. We showed that the juxta-centromeric region of human chromosome 21 is enriched for retrotransposed pseudogenes and gene fragments transferred by interchromosome duplications, but we do not rule out the possibility that the region harbors functional genes also.     

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.     

Characterization of the functional gene and several processed pseudogenes in the human triosephosphate isomerase gene family.

The functional gene and three intronless pseudogenes for human triosephosphate isomerase were isolated from a recombinant DNA library and characterized in detail. The functional gene spans 3.5 kilobase pairs and is split into seven exons. Its promoter contains putative TATA and CCAAT boxes and is extremely rich in G and C residues (76%). The pseudogenes share a high degree of homology with the functional gene but contain mutations that preclude the synthesis of an active triosephosphate isomerase enzyme. Sequence divergence calculations indicate that these pseudogenes arose approximately 18 million years ago. We present evidence that there is a single functional gene in the human triosephosphate isomerase gene family.     

Genome-Wide Survey and Comparative Analysis of LTR Retrotransposons and Their Captured Genes in Rice and Sorghum

Long terminal repeat (LTR) retrotransposons are the major class I mobile elements in plants. They play crucial roles in gene expansion, diversification and evolution. However, their captured genes are yet to be genome-widely identified and characterized in most of plants although many genomes have been completely sequenced. In this study, we have identified 7,043 and 23,915 full-length LTR retrotransposons in the rice and sorghum genomes, respectively. High percentages of rice full-length LTR retrotransposons were distributed near centromeric region in each of the chromosomes. In contrast, sorghum full-length LTR retrotransposons were not enriched in centromere regions. This dissimilarity could be due to the discrepant retrotransposition during and after divergence from their common ancestor thus might be contributing to species divergence. A total of 672 and 1,343 genes have been captured by these elements in rice and sorghum, respectively. Gene Ontology (GO) and gene set enrichment analysis (GSEA) showed that no over-represented GO term was identified in LTR captured rice genes. For LTR captured sorghum genes, GO terms with functions in DNA/RNA metabolism and chromatin organization were over-represented. Only 36% of LTR captured rice genes were expressed and expression divergence was estimated as 11.9%. Higher percentage of LTR captured rice genes have evolved into pseudogenes under neutral selection. On the contrary, higher percentage of LTR captured sorghum genes were under purifying selection and 72.4% of them were expressed. Thus, higher percentage of LTR captured sorghum genes was functional. Small RNA analysis suggested that some of LTR captured genes in rice and sorghum might have been involved in negative regulation. On the other hand, positive selection has been observed in both rice and sorghum LTR captured genes and some of them were still expressed and functional. The data suggest that some of these LTR captured genes might have evolved into new gene functions.     

The protein kinase ERK3 is encoded by a single functional gene: genomic analysis of the ERK3 gene family

Extracellular signal-regulated kinase 3 (ERK3) is a distantly related member of the mitogen-activated protein (MAP) kinase family of serine/threonine kinases. Here, we report the characterization of the genomic loci encoding ERK3 in mice and humans. The mouse ERK3 gene (Mapk6) spans more than 20 kb and is split into six exons. Its structure is similar to that of the human MAPK6 gene, which extends over 40 kb. We also identified and characterized a mouse Mapk6 processed pseudogene. In humans, database analysis has revealed the presence of six MAPK6 processed pseudogenes localized on four different chromosomes. We further show that the structure of MAPK6 is closely related to that of the gene encoding the homologous protein kinase p63(MAPK) (MAPK4), suggesting that the two genes arose by duplication. Our analysis demonstrates that the ERK3 subfamily of MAP kinase genes is composed of two functional genes, MAPK6 and MAPK4, and several pseudogenes.     

Characterization of the bovine BCL2L1 gene and related pseudogenes

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

The ornithine aminotransferase-encoding gene family of rat: cloning, characterization, and evolutionary relationships between a single expressed gene and three pseudogenes

As a first step towards understanding the molecular mechanisms through which the expression of the gene (OAT) encoding ornithine aminotransferase (OAT) is regulated in a tissue-specific manner, we have used a near full length OAT cDNA to isolate related sequences from a rat genomic DNA library. Twenty-one unique clones representing five contigs and spanning approximately 140 kb of genomic DNA were isolated and characterized. From these clones we have identified a single expressed OAT gene and three processed pseudogenes. The comparison of the EcoRI, BamHI, and HindIII fragments contained within these genomic clones with those detected in total genomic DNA by the cDNA probe suggests that essentially all of the OAT-related sequences in the rat genome have been isolated. Thus, the tissue-specific regulation of OAT gene expression appears to be effected through a single expressed gene. Data are presented which suggest that the OAT-1, OAT-2, and OAT-3 pseudogenes arose approximately 28.5, 7.3, and 25.1 Myr ago, respectively. Mutation rates are presented for each codon position of the expressed rat and human OAT genes. The region of the rat genome flanking the boundary of the OAT-3 pseudogene is of additional interest as it shares considerable identity to sequences contained within expressed genes and flanking other processed pseudogenes.     

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

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

A neutral model of transcriptome evolution

Microarray technologies allow the identification of large numbers of expression differences within and between species. Although environmental and physiological stimuli are clearly responsible for changes in the expression levels of many genes, it is not known whether the majority of changes of gene expression fixed during evolution between species and between various tissues within a species are caused by Darwinian selection or by stochastic processes. We find the following: (1) expression differences between species accumulate approximately linearly with time; (2) gene expression variation among individuals within a species correlates positively with expression divergence between species; (3) rates of expression divergence between species do not differ significantly between intact genes and expressed pseudogenes; (4) expression differences between brain regions within a species have accumulated approximately linearly with time since these regions emerged during evolution. These results suggest that the majority of expression differences observed between species are selectively neutral or nearly neutral and likely to be of little or no functional significance. Therefore, the identification of gene expression differences between species fixed by selection should be based on null hypotheses assuming functional neutrality. Furthermore, it may be possible to apply a molecular clock based on expression differences to infer the evolutionary history of tissues.     

人类基因组上的假基因

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