Age and Detection of Retroprocessed Pseudogenes in Murine Rodents

Retroprocessed pseudogenes, calmodulin II (ψ₁, ψ₂, and ψ₃ CALMII), ψα-tubulin, π-glutathione S-transferase (ψπ-GST) from rat, lactic acid dehydrogenase (ψ LDH) from mouse, and heat shock protein 60 chaperonin (ψ HSP60) from Chinese hamster, were examined for their presence in these species by polymerase chain reaction (PCR). Pseudogenes of these murine rodents were detected by PCR only in those species in which the genes were originally identified, suggesting that the selected pseudogene of one species arose too recently to be detected in the genomes of the other rodent species. The calculated ages of the rodent pseudogenes ranged from 1.7 Myr (ψα-tubulin) to 7.5 Myr (ψ₃ CALMII) when employing a homologous functional gene of the taxon as a reference in the relative rate test with the mouse or rat as the outgroup. Given the high rate of divergence of the genes of rodents relative to other species, selection of an outgroup with similar mutation rates seems warranted. To justify further the conclusion that the selected pseudogenes were indeed retroprocessed after these three taxa diverged, the presence of the pseudogenes in the genome of different rat species was examined. The existence of ψ₃ CALMII and ψα-tubulin pseudogenes of Rattus norvegicus among species belonging to Rattus sensu stricto is evidence for the common ancestry of this group. Received: 23 May 2001 / Accepted: 31 May 2001     

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     

Characterization of α-gliadin genes from diploid wheats and the comparative analysis with those from polyploid wheats

To carry out comparative analysis of the α-gliadin genes on A genomes of diploid and polyploid wheats, 8 full-length α-gliadin genes, including 3 functional genes and 5 pseudogenes, were obtained from diploid wheats, among which 2, 2 and 4 α-gliadin genes were isolated from T. urartu, T. monococcum, and T. boeoticum, respectively. The results indicated that higher number of α-gliadin pseudogenes have been present in diploid wheats before the formation of polyploid wheats. Amino acid sequence comparative analysis among 26 α-gliadin genes, including 16 functional genes and 10 pseudogenes, from diploid and polyploid wheats was conducted. The results indicated that all α-gliadins contained four coeliac toxic peptide sequences (i.e., PSQQ, QQQP, QQPY, and QPYP). The polyglutamine domains are highly variable, and the second polyglutamine stretch is usually disrupted by the lysine or arginine residue at the fourth position. The unique domain I is the most conserved domain. There are four and two conserved cysteine residues in the unique domains I and II, respectively. Comparative analysis indicated that the functional α-gliadin genes from A genome are highly conserved, whereas the identity of pseudogenes in diploid wheats are higher than those in hexaploid wheats. Phylogenetic analysis indicated that all the analyzed functional α-gliadin genes could be clustered into two major groups, among which one group could be further divided into 5 subgroups. The origin of α-gliadin pseudogene and functional genes were also discussed. The text was submitted by the authors in English.     

Mutation Pattern Variation Among Regions of the Primate Genome

We sequenced three argininosuccinate-synthetase-processed pseudogenes (ΨAS-A1, ΨAS-A3, ΨAS-3) and their noncoding flanking sequences in human, orangutan, baboon, and colobus. Our data showed that these pseudogenes were incorporated into the genome of the Old World monkeys after the divergence of the Old World and New World monkey lineages. These pseudogene flanking regions show variable mutation rates and patterns. The variation in the G/C to A/T mutation rate (u) can account for the unequal GC contents at equilibrium: 34.9, 36.9, and 41.7% in the pseudogene ΨAS-A1, ΨAS-A3, and ΨAS-3 flanking regions, respectively. The A/T to G/C mutation rate (v) seems stable and the u/v ratios equal 1.9, 1.7, and 1.4 in the flanking regions of ΨAS-A1, ΨAS-A3, and ΨAS-3, respectively. These ``regional' variations of the mutation rate affect the evolution of the pseudogenes, too. The ratio u/v being greater than 1.0 in each case, the overall mutation rate in the GC-rich pseudogenes is, as expected, higher than in their GC-poor flanking regions. Moreover, a ``sequence effect' has been found. In the three cases examined u and v are higher (at least 20%) in the pseudogene than in its flanking region—i.e., the pseudogene appears as mutation ``hot' spots embedded in ``cold' regions. This observation could be partly linked to the fact that the pseudogene flanking regions are long-standing unconstrained DNA sequences, whereas the pseudogenes were relieved of selection on their coding functions only around 30–40 million years ago. We suspect that relatively more mutable sites maintained unchanged during the evolution of the argininosuccinate gene are able to change in the pseudogenes, such sites being eliminated or rare in the flanking regions which have been void of strong selective constraints over a much longer period. Our results shed light on (1) the multiplicity of factors that tune the spontaneous mutation rate and (2) the impact of the genomic position of a sequence on its evolution. Received: 10 February 1997 / Accepted: 21 April 1997     

Loci for human U1 RNA: structural and evolutionary implications

Three clones U1-1, U1-6, and U1-8 containing sequences related to human U1 RNA have been studied by sequence analysis. The results show that each of the three clones represents a distinct locus. The U1-6 locus is closely related to the HU1-1 locus, which is believed to represent a functional U1 gene. The U1-1 and U1-8 loci are pseudogenes by definition, since they contain sequences that are closely related to but not identical with the human U1 RNA sequence. The U1-6 locus contains the sequence T-A-T-A-T close to the 5'-end of the U1 sequence but it is unclear if this represents the promoter. When the U1-8 locus was compared to the U1-6 locus, it was observed that the 5'-flanking sequences, except in the immediate vicinity of the pseudogene, are as well-conserved as the U1-related sequence itself, at least up to position -220. The high degree of homology in the 5'-flanking region suggests that U1 genes have a much more strict sequence requirement with regard to 5'-flanking sequences than most other eukaryotic genes. The U1-6 and U1-8 loci contain the sequence T-A-T-G-T-A-G-A-T-G-A between positions -211 and -221. An identical sequence is present in the equivalent position in the HU1-1 locus, and may represent the promoter. The high degree of conservation in the postulated promoter region indicates that pseudogenes like U1-8 possibly could be expressed. A truncated U1-related sequence is present between 106 to 150 nucleotides upstream from the U1 gene/pseudogene in the U1-6, the U1-8 and the HU1-1 loci, suggesting that the U1 genes may have been clustered early in evolution. The U1-1 locus has a strikingly different structure from the U1-8 locus; the pseudogene itself is as closely related to the U1 RNA sequence as is the U1-8 pseudogene but the flanking sequences, both on the 5' and the 3' side, share no detectable homology with the corresponding regions in the U1-6 or U1-8 loci. It may therefore be postulated that small nuclear RNA pseudogenes are created by several different mechanisms.     

A Translocated Mitochondrial Cytochrome b Pseudogene in Voles (Rodentia: Microtus)

A full-length cytochrome b pseudogene was found in rodents; it has apparently been translocated from a mitochondrion to the nuclear genome in the subfamily Arvicolinae. The pseudogene (ψcytb) differed from its mitochondrial counterpart at 201 of 1143 sites (17.6%) and by four indels. Cumulative evidence suggests that the pseudogene has been translocated to the nucleus. Phylogenetic reconstruction indicates that the pseudogene arose before the diversification of M. arvalis/M. rossiaemeridionalis from M. oeconomus, but after the divergence of the peromyscine/sigmodontine/arvicoline clades some ∼10 MYA. Published rates of divergence between mitochondrial genes and their nuclear pseudogenes suggest that the translocation of this mitochondrial gene to the nuclear genome occurred some 6 MYA, in agreement with the phylogenetic evidence. Received: 16 January 1998 / Accepted: 18 July 1998     

The wheat ω-gliadin genes: structure and EST analysis

A survey and analysis is made of all available ω-gliadin DNA sequences including ω-gliadin genes within a large genomic clone, previously reported gene sequences, and ESTs identified from the large wheat EST collection. A contiguous portion of the Gli-B3 locus is shown to contain two apparently active ω-gliadin genes, two pseudogenes, and four fragments of the 3′ portion of ω-gliadin sequences. Comparison of ω-gliadin sequences allows a phylogenetic picture of their relationships and genomes of origin. Results show three groupings of ω-gliadin active gene sequences assigned to each of the three hexaploid wheat genomes, and a fourth group thus far consisting of pseudogenes assigned to the A-genome. Analysis of ω-gliadin ESTs allows reconstruction of two full-length model sequences encoding the AREL- and ARQL-type proteins from the Gli-A3 and Gli-D3 loci, respectively. There is no DNA evidence of multiple active genes from these two loci. In contrast, ESTs allow identification of at least three to four distinct active genes at the Gli-B3 locus of some cultivars. Additional results include more information on the position of cysteines in some ω-gliadin genes and discussion of problems in studying the ω-gliadin gene family. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evolution of Growth Hormone in Primates: The GH Gene Clusters of the New World Monkeys Marmoset (Callithrix jacchus) and White-Fronted Capuchin (Cebus albifrons)

The GH gene cluster in marmoset, Callithrix jacchus, comprises eight GH-like genes and pseudogenes and appears to have arisen as a consequence of gene duplications occurring independently of those leading to the human GH gene cluster. We report here the complete sequence of the marmoset GH gene locus, including the intergenic regions and 5′ and 3′ flanking sequence, and a study of the multiple GH-like genes of an additional New World monkey (NWM), the white-fronted capuchin, Cebus albifrons. The marmoset sequence includes 945 nucleotides (nt) of 5′ flanking sequence and 1596 nt of 3′ flanking sequence that are “unique”; between these are eight repeat units, including the eight GH genes/pseudogenes. The breakpoints between these repeats are very similar, indicating a regular pattern of gene duplication. These breakpoints do not correspond to those found in the much less regular human GH gene cluster. This and phylogenetic analysis of the repeat units within the marmoset gene cluster strongly support the independent origin of these gene clusters, and the idea that the episode of rapid evolution that occurred during GH evolution in primates preceded the gene duplications. The marmoset GH gene cluster also differs from that of human in having fewer and more evenly distributed Alu sequences (a single pair in each repeat unit) and a “P-element” upstream of every gene/pseudogene. In human there is no P-element upstream of the gene encoding pituitary GH, and these elements have been implicated in placental expression of the other genes of the cluster. The GH gene clusters in marmoset and capuchin appear to have arisen as the consequence of a single-gene duplication event, but in capuchin there was then a remarkable expansion of the GH locus, giving at least 40 GH-like genes and pseudogenes. Thus even among NWMs the GH gene cluster is very variable. [Reviewing Editor: Nicolas Galtier]     

Genes and pseudogenes for human U2 RNA: Implications for the mechanism of pseudogene formation

Three loci, designated U2/4, U2/6 and U2/7, which contain sequences related to human U2 RNA, have been studied. The U2/6 locus contains a tandem array of bona fide U2 genes. U2/4 and U2/7, in contrast, contain pseudogenes whose sequences deviate significantly from that of mammalian U2 RNA. The two pseudogenes appear to have been created by different mechanisms. The sequences that flank the pseudogene in the U2/4 locus lack homology to the corresponding sequences in functional human U2 genes, except for 10 base-pairs immediately following the 3′ end. The conserved 3′-flanking segment is homologous to those nucleotides that are present in a U2 RNA precursor. No direct repeats flank the pseudogene in the U2/4 locus. The observations thus suggest that a complementary DNA copy of the U2 RNA precursor was inserted into a blunt-ended chromosomal break to generate the U2/4 locus.The U2/7 locus, in contrast, revealed flanking sequence homology when compared to functional U2 genes, both on the 5′ and 3′ sides of the pseudogene. The homology was interrupted on both sides by repetitive sequences belonging to the Alu family. On the 5′ side the homology continues beyond the Alu repeats whereas on the 3′ side it ends precisely at the Alu repeat. This Alu repeat is inserted in a region where a homocopolymeric region of alternating C and T residues is located in functional U2 loci. The observed organization of the U2/7 locus suggests that a previously functional U2 locus was invaded by Alu repeats and subsequently accumulated base substitutions to become a pseudogene.     

Evolution of a multigene family that encodes the Kunitz chymotrypsin inhibitor in winged bean: a possible intermediate in the generation of a new gene with a distinct pattern of expression

Winged bean Kunitz chymotrypsin inhibitor (WCI) accumulates in an organ-specific and temporally regulated manner. The protein is encoded by a multigene family that includes at least four putative inhibitor-coding genes and three pseudogenes. The structure of the WCI genes indicates that an insertion at a 5′ proximal site occurred after duplication of the ancestral WCI gene and that several gene conversion events subsequently contributed to the evolution of this gene family. Analysis of the promoter activity of the 5′ regions of the WCI genes in transgenic tobacco showed that only the 5′ regions of the WCI-3a and WCI-3b genes, which encode the major WCI protein in winged bean, promoted the organ-specific and temporally regulated expression of a reporter gene. The 5′ region of a pseudogene, the WCI-P1 gene which contains frameshift mutations, exhibited constitutive promoter activity in tobacco, an indication that the 5′ region of the WCI-P1 gene might spontaneously have acquired new regulatory sequences during evolution. Since gene conversion is a relatively frequent event and since the homology between the WCI-P1 and WCI-3a/b genes is disrupted at a 5′ proximal site by remnants of an inserted sequence, the WCI-P1 gene appears to be a possible intermediate that could be converted into a new functional gene with a distinct pattern of expression by a single gene-conversion event. Received: 17 April 1996 / Accepted: 23 October 1996     

Structure of the catfish IGH locus: analysis of the region including the single functional IGHM gene

The catfish IGH locus is large (∼1 Mb) and complex, having undergone multiple internal duplications and transpositions. To define the structure of the locus that contains the single expressed IGHM gene, two overlapping bacterial-artificial-chromosome (BAC) clones spanning the most 3′ end of the channel catfish immunoglobulin heavy (IGH) chain locus have been completely sequenced. The analyses created a contig of 257,153 bp containing 55 VH, 6 D, 12 JH genes and the IGH constant region genes encoding the functional secreted and membrane forms of IgM and the membrane form of IgD. This analysis revealed three major features. First, no C-region genes were found aside from the previously described IGHM1 and IGHD1, with the latter gene being the most 3′ C-region gene of the catfish IGH locus. There was no evidence in the region sequenced for genes that could encode an Ig class similar to the IgZ/IgT described in zebrafish, trout and pufferfish. Second, there are a high number of VH pseudogenes, 28 out of 55 (51%). In contrast, the entire zebrafish IGH locus has 40 functional VH genes and eight pseudogenes (17%). Third, an internal duplication of a 52.4-kb block of VH genes has occurred. These observations suggest that the IGH locus of teleost fish varies significantly from species to species in the diversity of C-region genes as well as the numbers of genes encoding V regions.     

Organization, Structure, and Evolution of the Nonadult Rat β-Globin Gene Cluster

The β-globin gene cluster of Wistar rat was extensively cloned and the embryonic genes were mapped and sequenced. Four overlapping λ Dash recombinant clones cover about 31 kb and contain four nonadult β-globin genes, 5′–ε1–γ1–γ2–ψγ3–3′. The ε1 and γ2 are active genes, since their protein products were detected in the fetal stage of the rat (Iwahara et al., J Biochem 119:360–366, 1996). The γ1 locus might be a pseudogene, since the ATA box in the promoter region is mutated to GTA; however, no other defect is observed. The ψγ3 locus is a truncated pseudogene because a 19-base deletion, which causes a shift of the reading frame, is observed between the second nucleotide of the putative codon 68 and codon 76. A sequence comparison suggests that the ψγ3 might be produced by a gene conversion event of the proto-γ-globin gene set. Possible histories of the evolution of rat nonadult β-globin genes are discussed. Received: 6 August 1998 / Accepted: 12 February 1999     

Gene Conversions Can Generate Sequence Variants in the Late Chorion Multigene Families of Bombyx Mori

The 140-kbp late chorion locus of Bombyx mori strain 703 contains 15 divergently oriented gene pairs encoding the high cysteine (Hc) eggshell proteins. Sequence homology is approximately 91% for the 2-kb region of each gene pair, including the 5' flanking region, intron and exons. The homology rapidly disappears within a few hundred basepairs of the 3' end of most genes. Here we present the results of the nucleotide sequence and genomic blot comparison of Hc genes from different races of B. mori. Comparison of the nucleotide sequences of the same gene pair in two different races reveals that most of the nucleotide differences occur in clusters or patches and correspond to sequences present in other Hc genes in the locus. The number of nucleotide differences that have accumulated in the highly conserved regions of the gene pair (2.3/100 bp), most of which are attributable to patchwork exchanges, is significantly higher than the number of differences in the poorly conserved 3' flanking regions (0.6/100 bp), due primarily to new mutations. These data are consistent with a gene conversion process, which in the short-term generates new combinations of sequence variants, but in the long-term results in concerted evolution. Genomic blot analyses of different geographical races of B. mori reveal that there is variation in the number of Hc gene pairs (14-19 gene pairs), indicating that unequal crossovers also occur in the locus.     

Identification and characterization of plastid trnF(GAA) pseudogenes in four species of Solanum (Solanaceae)

Non-functional trnF pseudogenes that rarely occur in embryophytes have been found in Solanaceae. We have sequenced the trnL-F intergenic spacer of four species of Solanum, and found duplicated regions of the original trnF gene. These repeats were 94–260 bp long causing large length variation in the trnL-F intergenic spacer resulting from differences in pseudogene copy number (2–4). The duplicated trnF regions are comprised of several highly structured motifs, which were partial residues, or entire parts of the Anticodon, T- and D-domains of the original gene, but all lacked the acceptor stems at the 5′- or 3′-end. Pseudogenes included several transitions and transversions in their sequences compared to the original trnF gene. Among pseudogene copies, T-domains were more frequent and fragmented than D-domain elements. Our results demonstrate that although chloroplast evolution is uniform such structural duplications in the sequences used for phylogenetic reconstructions should be treated with great caution.     

Genomic organization of the sheep immunoglobulin JH segments and their contribution to heavy chain variable region diversity

 The sheep immunoglobulin heavy chain Igh-J locus has been characterized in order to determine the genomic organization of JH segments and their contribution to heavy chain diversity. The locus contains six segments, of which two are functional and four are apparently pseudogenes. These segments span a 1.8 kilobase (kb) region. The distance between JH-ps4 (the 3′-most segment) and the first domain of the μ-chain encoding constant gene is about 5 kb. The two functional JH segments have a standard upstream recombination signal sequence, including heptamer and nonamer sequences separated by a 22–23 nucleotide spacer, and end with a RNA donor splice site. These two segments possess all the characteristic JH invariant residues and are found in expressed μ heavy chain variable regions. The 5′ functional JH1 segment is used in more than 90% of the cDNAs sequenced to date. The contribution of JH segment germline multiplicity to variable regions diversity appears therefore to be minimal. Comparison with other mammalian JH segments shows that all loci are very closely related and probably have evolved from a common ancestral locus. Received: 19 November 1996 / Revised: 17 March 1997