Molecular cloning and sequencing of a murine pgk-1 pseudogene family

Seven genomic mouse DNA fragments carrying pgk-1-homologous regions have been cloned and sequenced. They have to be classified as processed genes because intervening sequences, present in their productive counterpart, are absent. Four pseudogenes (I-IV) represent nearly the complete sequence of pgk-1 cDNA. Two of these genes (I and II), although rather different from the published mouse pgk-1 cDNA in the 3'-untranslated region, represent the actual mouse pgk-1 cDNA sequence in the coding part except for substitutions in the third position of three codons. These genes can code for a functional PGK protein but, lacking as they do classical promoter structures, are probably not expressed. They show the typical characteristics of retroposons, being flanked by A-rich regions and direct repeats which are localized at the positions where the homology with the mouse pgk-1 cDNA is interrupted. Pseudogenes III and IV have numerous mutations. Gene III is also flanked by direct repeats, whereas gene IV is flanked by inverted repeats. The other three genes are flanked by direct repeats localized further inside the target sites. They are truncated and mutated extensively as usually observed with pseudogenes.     

Comparative analysis of the structural organization of two closely related vitellogenin genes in X. laevis

The structural organization of the two closely related vitellogenin genes A1 and A2 has been determined and compared by electron microscopy. In both genes the mRNA-coding sequence of 6 kb is interrupted 33 times, leading to a total gene length of 21 kb for gene A1 and 16 kb for gene A2. Thus both genes have a mean exon length of 0.175 kb, while the mean intron length is 0.45 kb in gene A1 and 0.31 kb in gene A2. Because the introns interrupt the structural sequence at homologous positions in genes A1 and A2, we suggest that these two genes are the products of a duplication of an ancestral gene which had an intron-exon arrangement similar to that of the extant genes. Since the duplication event, the sequence and length of the analogous introns have changed rapidly, whereas homologous exons have diverged to an extent of only 5% of their sequences. The results suggest different mechanisms of evolution for exons and introns. While the exons evolved primarily by point mutations, such mutations, as well as deletion, insertion and duplication events, were important in the evolution of the introns.     

Genomic organization of the human multidrug resistance (MDR1) gene and origin of P-glycoproteins

The MDR1 gene, responsible for multidrug resistance in human cells, encodes a broad specificity efflux pump (P-glycoprotein). P-glycoprotein consists of two similar halves, each half including a hydrophobic transmembrane region and a nucleotide-binding domain. On the basis of sequence homology between the N-terminal and C-terminal halves of P-glycoprotein, we have previously suggested that this gene arose by duplication of a primordial gene. We have now determined the complete intron/exon structure of the MDR1 gene by direct sequencing of cosmid clones and enzymatic amplification of genomic DNA segments. The MDR1 gene includes 28 introns, 26 of which interrupt the protein-coding sequence. Although both halves of the protein-coding sequence are composed of approximately the same number of exons, only two intron pairs, both within the nucleotide-binding domains, are located at conserved positions in the two halves of the protein. The other introns occur at different locations in the two halves of the protein and in most cases interrupt the coding sequence at different positions relative to the open reading frame. These results suggest that the P-glycoprotein arose by fusion of genes for two related but independently evolved proteins rather than by internal duplication.     

Vacuolar type H(+)-ATPase genes: presence of four genes including pseudogenes for the 16-kDa proteolipid subunit in the human genome.

Genes for the human vacuolar type H(+)-ATPase proteolipid (16-kDa) subunit were cloned and their nucleotide sequences were determined. Comparison of the deduced sequences indicated that at least four genes including pseudogenes are present in the human genome. One of them corresponded to that for the 16-kDa subunit expressed in HeLa cells. The coding sequence was separated by two introns. The second intron was located in the DNA segment giving a loop between the second and third transmembrane helices, supporting the idea that the 16-kDa subunit was evolved by gene duplication. The primary sequence determined from the second clone had a termination codon behind the third transmembrane helix. Possible translation products from the other two clones had no putative acidic residues essential for proton transport function of the 16-kDa subunit. Thus, it is interesting to know whether these genes are transcribed, since they may have unique cellular functions.     

Identification of a heat-shock pseudogene from Caenorhabditis elegans

While characterizing the hsp70 gene family from Caenorhabditis elegans we encountered an unusual member of this family. Sequence data reveal that the hsp-2ps gene is a pseudogene of the constitutively expressed, heat-inducible hsp-1 gene. Two stop codons generated near the 5' end of the sequence as well as several frameshift mutations and a large internal deletion confirm the identification of hsp-2ps as a pseudogene. The nucleotide substitution rate of the third codon position was twice that of the first and second codon positions, suggesting that the hsp-2ps gene was nonfunctional since the time of the duplication event. The hsp-2ps gene duplicates a region of the hsp-1 gene that lies exclusively within the transcribed region and retains the introns. We feel that the hsp-2ps gene was produced by a transpositional duplication event, which occurred approximately 8.5 million years ago.     

Genomic structure and promoter activity of the testis haploid germ cell-specific intronless genes,Tact1 and Tact2

The Tact1 and Tact2 genes, each of which encodes an actin-like protein, are exclusively expressed and translated in haploid germ cells in testis. To characterize the haploid germ cell-specific gene structure, a mouse genomic library was screened with a Tact1 cDNA as a probe, and four independent phage clones containing the Tact1 gene were isolated. Southern hybridization and sequencing analyses revealed that Tact1 and Tact2 were single copy genes contained on a common fragment in a head-to-head orientation, and that the distance between these genes was less than 2 kb. Comparison of the nucleotide sequences of genomic DNA and cDNA demonstrated that Tact1 and Tact2 lack introns, although all known actin or actin-related genes in mammals contain introns. Human Tact orthologues also lack introns and are located within 6.4 kb in a head-to-head orientation. These findings indicate that Tact1 and Tact2 or one of these genes arose by retroposition of a spliced mRNA transcribed from an actin progenitor gene prior to the divergence of rodents and primates. The Tact1 and Tact2 genes are unusual retroposons in that they have retained an open reading frame and are expressed in testicular germ cells, because almost all retroposons become pseudogenes. It was revealed that a 2kb sequence between the two genes bidirectionally controls haploid germ-cell specific expression by analyzing transgenic mice. Comparison of the murine Tact genes with their human orthologues showed a high level of identity between the two species in the 5'-upstream and non-coding sequences as well as in the coding region, indicating that conserved elements in these regions may be involved in the regulation of haploid germ cell-specific expression. The promoter region contains no TATA-, CCAAT- or GC-boxes, although there are potential cAMP response element (CRE)-like motifs in the 5'-upstream region and the 5'-untranslated region in Tact1 and Tact2, respectively. Transient promoter analyses indicate that CREMtau may activate Tact1 and Tact2 expression in germ cells.     

Molecular evolution of serpins: homologous structure of the human alpha 1-antichymotrypsin and alpha 1-antitrypsin genes

alpha 1-Antichymotrypsin belongs to a supergene family that includes alpha 1-antitrypsin, antithrombin III, ovalbumin, and angiotensinogen. The human chromosomal alpha 1-antichymotrypsin gene has been cloned and its molecular structure established. The gene is approximately 12 kb in length and contains five exons and four introns. The locations of the introns within the alpha 1-antichymotrypsin gene are identical with those of the human alpha 1-antitrypsin and angiotensinogen genes. Other members of this supergene family contain introns located at nonhomologous positions of the genes. The homologous organization of the alpha 1-antichymotrypsin and alpha 1-antitrypsin genes corresponds with the high degree of homology between their protein sequences and suggests that these loci arose by recent gene duplication. A model is presented for the evolution of both the genomic structure and the protein sequences of the serine protease inhibitor superfamily.     

The Human Actin-Regulatory Protein Cap G: Gene Structure and Chromosome Location

Cap G (formerly called macrophage capping protein or gCap39) is a member of the gelsolin/villin fanlily of actin-regulatory proteins. Unlike all other members of this family, Cap G caps the barbed ends of actin filaments, but does not sever them. This protein is half the molecular weight and contains half the number of repeat subunits (3 vs 6) of gelsolin and villin, suggesting that these two proteins may have arisen by gene duplication of the Cap G gene. To investigate this possibility we have cloned and sequenced the human Cap G gene (gene symbol CAPG). The gene is 16.6 kb in size, contains 10 exons and 9 introns and is located on the proximal short arm of chromosome 2. The open reading frame is 6.9 kb, having 9 exons and 8 introns. This region contains 3 splice sites that are nearly identical to the human gelsolin gene, but shares only one with villin, indicating that CAPG is more closely related to gelsolin. Further comparisons of these three genes, however, indicate that the evolutionary steps resulting in human gelsolin and villin are likely to have been more complex than a simple tandem duplication of the Cap G gene.     

Duplicated cytoglobin genes in teleost fishes

Cytoglobin is a recently discovered myoglobin-related O2-binding protein of vertebrates with uncertain function. It occurs as single-copy gene in mammals. Here, we demonstrate the presence of two paralogous cytoglobin genes (Cygb-1 and Cygb-2) in the teleost fishes Danio rerio, Oryzias latipes, Tetraodon nigroviridis, and Takifugu rubripes. The globin-typical introns at positions B12.2 and G7.0 are conserved in both genes, whereas the C-terminal exon found in mammalian cytoglobin is absent in the fish genes. Phylogenetic analyses show that the two cytoglobin genes diverged early in teleost evolution. This is confirmed by gene synteny analyses, which suggest a large-scale duplication event. Although both cytoglobin genes are highly conserved and have evolved under purifying selection, substitution rates are significantly higher in Cygb-1 than in Cygb-2. Similar to their mammalian ortholog, both fish cytoglobins are expressed in a broad range of tissues. However, Cygb-2 is more than 250-fold stronger expressed in neuronal tissues, suggesting a subfunctionalization of the two cytoglobin paralogs after gene duplication.     

Intervening sequences in paralogous genes: a comparative genomic approach to study the evolution of X chromosome introns

The enlargement of the genome size and the decrease in genome compactness with increase in the number and size of introns is a general pattern during the evolution of eukaryotes. Among the possible mechanisms for modifying intron size, it has been suggested that the insertion of transposable elements might have an important role in driving intron evolution. The analysis of large portions of the human genome demonstrated that a relatively recent (50 to 100 MYA) accumulation of transposable elements appears to be biased, favoring a preferential insertion of LINE1 transposons into sex chromosomes rather than into autosomes. In the present work, the effect of chromosomal location on the increase in size of introns was evaluated with a comparative analysis performed on pairs of human paralogous genes, one located on the X chromosome and the second on an autosome. A phylogenetic analysis was also performed on the X-encoded proteins and their paralogs to confirm orthology-paralogy and to approximately estimate the time of gene duplication. Statistical analysis of total intron length for each pair of paralogous genes provided no evidence for a larger size of introns in the gene copies located on the X chromosome. On the opposite, introns of autosomal genes were found to be significantly longer than introns of their X-linked paralogs. Likewise, LINE1 elements were not significantly more frequent in X-chromosome introns, whereas the frequency of SINE elements showed a marginally significant bias toward autosomal introns.     

Human placental protein 14 gene: sequence and characterization of a short duplication

Differential hybridization of cDNAs corresponding to mRNAs expressed in the human endometrium during the secretory phase or during the first trimester of pregnancy, but not during the proliferative phase, allowed us to isolate and characterize cDNAs encoding human placental protein 14 (PP14). The cDNA was used to isolate the PP14 gene from a human genomic library. The entire gene encompasses 5.05 kb divided into seven exons by six introns. The human PP14 gene shows identical organization with the ovine beta-lactoglobulin gene, as expected from protein homology. Sequencing of 3 kb of the 5'-flanking region of the gene allowed us to characterize a 400-bp duplication of the PP14 gene lying at position -2,660. This duplication was homologous to 100 bp of exon 4 and 300 bp of intron 4, including 180 bp corresponding exactly to the right arm of an Alu element lying on the complementary strand. This homology suggests that this duplication may have arisen through a retroposition event.     

Intron loss and gain during evolution of the catalase gene family in angiosperms.

Angiosperms (flowering plants), including both monocots and dicots, contain small catalase gene families. In the dicot, Arabidopsis thaliana, two catalase (CAT) genes, CAT1 and CAT3, are tightly linked on chromosome 1 and a third, CAT2, which is more similar to CAT1 than to CAT3, is unlinked on chromosome 4. Comparison of positions and numbers of introns among 13 angiosperm catalase genomic sequences indicates that intron positions are conserved, and suggests that an ancestral catalase gene common to monocots and dicots contained seven introns. Arabidopsis CAT2 has seven introns; both CAT1 and CAT3 have six introns in positions conserved with CAT2, but each has lost a different intron. We suggest the following sequence of events during the evolution of the Arabidopsis catalase gene family. An initial duplication of an ancestral catalase gene gave rise to CAT3 and CAT1. CAT1 then served as the template for a second duplication, yielding CAT2. Intron losses from CAT1 and CAT3 followed these duplications. One subclade of monocot catalases has lost all but the 5''-most and 3''-most introns, which is consistent with a mechanism of intron loss by replacement of an ancestral intron-containing gene with a reverse-transcribed DNA copy of a fully spliced mRNA. Following this event of concerted intron loss, the Oryza sativa (rice, a monocot) CAT1 lineage acquired an intron in a novel position, consistent with a mechanism of intron gain at proto-splice sites.     

The role of introns in repeat protein gene formation

Genes composed of tandem repetitive sequence motifs are abundant in nature and are enriched in eukaryotes. To investigate repeat protein gene formation mechanisms, we have conducted a large-scale analysis of their introns and exons. We find that a wide variety of repeat motifs exhibit a striking conservation of intron position and phase, and are composed of exons that encode one or two complete repeats. These results suggest a simple model of repeat protein gene formation from local duplications. This model is corroborated by amino acid sequence similarity patterns among neighboring repeats from various repeat protein genes. The distribution of one- and two-repeat exons indicates that intron-facilitated repeat motif duplication, in which the start and end points of duplication are located in consecutive intronic regions, significantly exceeds intron-independent duplication. These results suggest that introns have contributed to the greater abundance of repeat protein genes in eukaryotic versus prokaryotic organisms, a conclusion that is supported by taxonomic analysis.     

Gene structure and evolution of testicular haploid germ cell-specific genes, Oxct2a and Oxct2b

OXCT/SCOT is the rate-determining enzyme in ketolysis in mitochondria of many extrahepatic organs. Two testicular isoforms, Oxct2a and Oxct2b, are highly homologous and specifically expressed in haploid spermatids of the mouse. In this report, we analyzed the structure and evolution of Oxct2a and Oxct2b. Both Oxct2's are single-copy intronless genes, of which nucleotide sequences are conserved with Oxct, indicating that these genes are transposons generated from Oxct. A CpG island was found within both Oxct2's. Oxct2a and Oxct2b are located in the third introns of Bmp8a and Bmp8b, and they are positioned within a 240-kb region in a tail-to-tail orientation on chromosome 4. This structural feature was also conserved in a syntenic region of human 1p34.3. Structural similarity between mice and humans indicated that these two sets of genes were generated by a segmental gene duplication, which occurred before the primate-rodent split. Dot matrix and phylogenetic tree analyses demonstrated that multiple rounds of intrachromosomal gene conversion between the two loci occurred in each species independently.