Human U1 small nuclear RNA pseudogenes do not map to the site of the U1 genes in 1p36 but are clustered in 1q12-q22.

Human U1 small nuclear RNA is encoded by approximately 30 gene copies. All of the U1 genes share several kilobases of essentially perfect flanking homology both upstream and downstream from the U1 coding region, but remarkably, for many U1 genes excellent flanking homology extends at least 24 kilobases upstream and 20 kilobases downstream. Class I U1 RNA pseudogenes are abundant in the human genome. These pseudogenes contain a complete but imperfect U1 coding region and possess extensive flanking homology to the true U1 genes. We mapped four class I pseudogenes by in situ hybridization to the long arm of chromosome 1, bands q12-q22, a region distinct from the site on the distal short arm of chromosome 1 to which the U1 genes have been previously mapped (Lund et al., Mol. Cell. Biol. 3:2211-2220, 1983; Naylor et al., Somat. Cell Mol. Genet. 10:307-313, 1984). We confirmed our in situ hybridization results by genomic blotting experiments with somatic cell hybrid lines with translocation products of human chromosome 1. These experiments provide further evidence that class I U1 pseudogenes and the true U1 genes are not interspersed. The results, along with those published elsewhere (Bernstein et al., Mol. Cell. Biol. 5:2159-2171, 1985), suggest that gene amplification may be responsible for the sequence homogeneity of the human U1 gene family.     

Human U1 RNA pseudogenes may be generated by both DNA- and RNA-mediated mechanisms.

Analysis of cloned human genomic loci homologous to the small nuclear RNA U1 established that such sequences are abundant and dispersed in the human genome and that only a fraction represent bona fide genes. The majority of genomic loci bear defective gene copies, or pseudogenes, which contain scattered base mismatches and in some cases lack the sequence corresponding to the 3' end of U1 RNA. Although all of the U1 genes examined to date are flanked by essentially identical sequences and therefore appear to comprise a single multigene family, we present evidence for the existence of at least three structurally distinct classes of U1 pseudogenes. Class I pseudogenes had considerable flanking sequence homology with the U1 gene family and were probably derived from it by a DNA-mediated event such as gene duplication. In contrast, the U1 sequence in class II and III U1 pseudogenes was flanked by single-copy genomic sequences completely unrelated to those flanking the U1 gene family; in addition, short direct repeats flanked the class III but not the class II pseudogenes. We therefore propose that both class II and III U1 pseudogenes were generated by an RNA-mediated mechanism involving the insertion of U1 sequence information into a new chromosomal locus. We also noted that two other types of repetitive DNA sequences in eucaryotes, the Alu family in vertebrates and the ribosomal DNA insertions in Drosophila, bore a striking structural resemblance to the classes of U1 pseudogenes described here and may have been created by an RNA-mediated insertion event.     

Examination of four HLA class I pseudogenes. Common events in the evolution of HLA genes and pseudogenes.

The HLA class I gene family in lymphoblastoid cell line 721 has been studied in detail and a number of sequences in addition to the classical genes have been identified. The cloning, characterization, and nucleotide sequences of four sequences, all full length HLA class I pseudogenes, are described in this report. These pseudogenes, contained within 5.4-, 5.9-, 7.0-, and 9.2-kb HindIII fragments, each have the class I exon-intron structure as well as class I homology in their 5' and 3' flanking regions. However, all four sequences have one or more substitutions that perturb the coding region, leaving little doubt that they are in fact pseudogenes. Comparisons among these sequences and the HLA class I genes revealed that their homology with the class I genes is patchwork. Thus, although some regions have diverged, other contiguous intron-exon sequences are highly conserved. Comparisons in the 5' regions indicate that the pseudogene promoters more closely resemble the classical HLA promoters than the nonclassical promoters as none of the unique structural features found in the HLA-E, -F, or -G regulatory regions are present in any of the pseudogene promoters. Further comparisons revealed that at least two putative gene conversion events, similar to those hypothesized to have occurred in the evolution of some HLA genes, may have occurred in the evolution of some of the pseudogenes. These and other hypothetical events in the evolution of the class I gene family are discussed.     

Simple DNA sequences in homologous flanking regions near immunoglobulin VH genes: a role in gene interaction?

Five closely related immunoglobulin VH genes (subgroup II) were compared by sequencing of several kb of DNA. In three of the genes homology greater than 75% was found along an area of 4 kb that includes the coding region. The homology in flanking regions is only slightly lower than that in the coding sequences. Two other genes, which are located on the same EcoRI fragment, show high homology to the first three genes in the coding and immediately flanking regions. In more distant flanking regions no homology is found with the first three genes. This indicates that their evolutionary history differs from that of the other three genes. A region of simple DNA sequence composed of repetitive TCC and TCA elements was found at a distance of approximately 380 bp upstream from the initiator ATG of these VH genes. This region is the site where the two sets of genes abruptly start to diverge. The structure of the simple DNA sequence in the various VH genes suggests that it may be involved in gene interaction. We propose that both simple DNA sequences and homology in flanking regions serve a function in the correction of VH genes, which seem to be rather free to diverge and drift into pseudogenes. A correction mechanism may help this gene family to maintain its two major features, multiplicity and diversity.     

Human DNA sequences homologous to a protein coding region conserved between homeotic genes of Drosophila

Several human DNA sequences were isolated by virtue of homology to a highly conserved region that has been identified in a number of homeotic genes in Drosophila. Structural analysis of the human DNAs indicate that two separate and distinct regions sharing a high degree of homology with the homeo box sequences of Drosophila are separated by only 5 kb in the human genome. Sequence determination of these regions reveals that both human DNA sequences contain a region capable of coding 61 amino acids, which shares greater than 90% homology with the peptide sequences specified by the homeo box domain of Drosophila homeotic genes, Antennapedia, fushi tarazu, and Ultrabithorax. By contrast, the human DNA sequences lying outside of the 190 nucleotide homeo box region share virtually no sequence homology, either with the flanking sequences of the other human clones or with flanking regions of the known Drosophila homeotic genes.     

Four genes in two diverged subfamilies encode the ribulose-1,5-bisphosphate carboxylase small subunit polypeptides of Arabidopsis thaliana

The multigene family encoding the small subunit polypeptides of ribulose-1,5-bisphosphate carboxylase/oxygenase in the crucifer Arabidopsis thaliana has been isolated and the organization and structure of the individual members determined. The family consists of four genes which have been divided into two subfamilies on the basis of linkage and DNA and amino acid sequence similarities. Three of the genes, designated ats1B, ats2B, and ats3B, reside in tandem on an 8 kb stretch of the chromosome. These genes share greater than 95% similarity in DNA sequence and encode polypeptides identical in length and 96.7% similar in amino acid sequence. The fourth gene, ats1A, is at least 10 kb removed from, or completely unlinked to the B subfamily. The B subfamily genes are more similar to each other than to ats1A in nucleotide and amino acid sequence. All four genes are interupted by two introns whose placement within the coding region of the genes is conserved. The introns of the B subfamily genes are similar in length and nucleotide sequence, but show no similarity to the introns of ats1A. Comparison of the DNA sequences within the immediate 5 and 3 flanking sequences among the genes revealed only limited regions of homology. S1 analysis shows that all four genes are expressed.     

Analysis of a drosophila tRNA gene cluster: two tRNALeu genes contain intervening sequences

A recombinant DNA phage containing a cluster of Drosophila melanogaster tRNA genes has been isolated and analyzed. The insert of this phage has been mapped by in situ hybridization to chromosomal region 50AB, a known tRNA site. Nucleotide sequencing of the entire Drosophila tRNA coding region reveals seven tRNA genes spanning 2.5 kb of chromosomal DNA. This cluster is separated from other tRNA regions on the chromosome by at least 2.7 kb on one side, and 9.6 kb on the other. Two tRNA genes are nearly identical and contain intervening sequences of length 38 and 45 bases, respectively, in the anticodon loop. These two genes are assigned to be tRNALeu genes because of significant sequence homology with yeast tRNA3Leu, and secondary structure homology with yeast tRNA3Leu intervening sequence. In addition, an 8 base sequence (AAAAUCUU) is conserved in the same location in the intervening sequences of Drosophila tRNALeu genes and a yeast tRNA3Leu gene. Similar sequenes occur in all other tRNAs containing intervening sequences. The remaining five genes are identical tRNAIle genes, which are also identical to a tRNAIle gene from chromosomal region 42A. The 5' flanking regions are only weakly homologous, but each set of isoacceptors contains short regions of strong homology approximately 20 nucleotides preceding the tRNA coding sequences: GCNTTTTG preceding tRNAIle genes; and GANTTTGG preceding tRNALeu genes. The genes are irregularly distributed on both DNA strands; spacing regions are divergent in sequence and length.     

The keratin BIIIB gene family: isolation of cDNA clones and structure of a gene and a related pseudogene

The nucleotide sequence of cDNA clones encoding the three major BIIIB high-sulfur wool keratin proteins (BIIIB2, 3, and 4) and the structure of a BIIIB4 gene and a BIIIB3 pseudogene are reported. Although Southern blot analysis indicates that the BIIIB genes comprise a multigene family in the sheep genome, they are poorly represented in genomic DNA libraries. The family sequence homology of the coding region extends into the 5' and 3' untranslated regions and the near 5' flanking region of the BIIIB3 and 4 genes. These homologies suggest that the BIIIB3 and 4 genes represent the latest gene duplication event in the evolution of the BIIIB multigene family. Like the genes coding for other wool keratin matrix protein components, the BIIIB genes have the conserved 18-bp sequence immediately 5' to the initiation codon and also appear to lack introns.     

Structural analysis of gene loci for rat U1 small nuclear RNA.

Four phage clones which hybridize with U1 small nuclear RNA were obtained from a rat gene library. Two clones contain a presumed pseudogene. A third clone includes two gene candidates that are co-linear with the rat U1-RNA, 3.6kb apart and in the opposite orientation. The two genes are surrounded by identical sequences of 491bp upstream and 178bp downstream. The upstream sequences do not contain a TATA box, but share many block homologies with those for the human U1-RNA gene(1-3). A 101bp "identifier (ID) sequence", which was reported to be specifically expressed in rat brain (4), is inserted immediately after the shared sequence downstream of one of the genes. In the fourth clone, there are two putative pseudogenes, which have one or three nucleotide changes, 3kb apart and in the same orientation. Southern blot analysis of total rat DNA reveals about 50 U1-RNA genes/pseudogenes in the genome.     

Evolution of the functional human beta-actin gene and its multi-pseudogene family: conservation of noncoding regions and chromosomal dispersion of pseudogenes.

We have assigned six members of the human beta-actin multigene family to specific human chromosomes. The functional gene, ACTB, is located on human chromosome 7, and the other assigned beta-actin-related sequences are dispersed over at least four different chromosomes including one locus assigned to the X chromosome. Using intervening sequence probes, we showed that the functional gene is single copy and that all of the other beta-actin related sequences are recently generated in evolution and are probably processed pseudogenes. The entire nucleotide sequence of the functional gene has been determined and is identical to cDNA clones in the coding and 5' untranslated regions. We have previously reported that the 3' untranslated region is well conserved between humans and rats (Ponte et al., Nucleic Acids Res. 12:1687-1696, 1984). Now we report that four additional noncoding regions are evolutionarily conserved, including segments of the 5' flanking region, 5' untranslated region, and, surprisingly, intervening sequences I and III. These conserved sequences, especially those found in the introns, suggest a role for internal sequences in the regulation of beta-actin gene expression.     

Drosophila melanogaster U1 and U2 small nuclear RNA genes contain common flanking sequences

The flanking sequences of three U2 genes (or pseudogenes) and one U1 gene of Drosophila melanogaster have been determined. Comparison of the sequences reveals a remarkable homology between position ?30 and ?65 upstream from the structural genes, starting with a TATA box-like sequence. The 3′ flanking regions are also conserved in all genes and contain a canonical A-A-T-A-A-A polyadenylation signal.     

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.     

Three types of rat U1 small nuclear RNA genes with different flanking sequences are induced to express in vivo

There are about 50 copies of U1 RNA genes/pseudogenes in the rat genome. To date, we have isolated so far 25 phage clones carrying a U1 RNA gene/pseudogene from two rat genomic libraries. The 12 clones were selected by hybridization with the U1 RNA coding region under a stringent condition, and were mapped and sequenced. Here, we report three types of U1 RNA genes with different flanking sequences, all of which were shown to be induced to express in vivo by transfection with their polylinker-inserted maxi U1 RNA genes into cultured rat cells. Although these three classes of U1 RNA genes have few homologous flanking sequences, they provide both upstream and downstream of the genes two conserved blocks, which may possibly play an important role in U1 RNA expression.     

DNase I sensitive domain of the gene coding for the glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase

We have cloned a 36-kilobase segment of chicken DNA containing the gene coding for glyceraldehyde-3-phosphate dehydrogenase [GAPDH (EC 1.2.1.12)], a glycolytic enzyme which is expressed constitutively in all cell types. Using defined segments of this cloned DNA as probes, we have determined the DNase I sensitive domain of the GAPDH natural gene in the hen oviduct. When nuclei isolated from hen oviduct were treated with DNase I under conditions known to preferentially degrade actively transcribed genes (i.e., 15-20% of the DNA rendered perchloric acid soluble), a region of approximately 12 kilobase(s) (kb) containing the GAPDH coding sequences and flanking DNA was found to be highly susceptible to digestion by DNase I. Approximately 4 kb downstream from the end of the coding sequences, there was an abrupt transition from the DNase I sensitive or "open" configuration to the resistant or "closed" configuration. The chromatin then remained in a closed conformation for at least 10 kb further downstream. On the 5' side of the gene, the transition from a sensitive to a resistant configuration was located about 4 kb upstream from the gene. In addition, we have localized two repeated sequences in the area of DNA that was cloned. One of these is of the CR1 family of middle repetitive elements. It is located about 18 kb 3' to the gene and as such lies well outside of the DNase I sensitive region which encompasses GAPDH. The other repetitive element is of an uncharacterized family. It is located upstream from the gene and appears to be within a region of transition from the DNase I sensitive to resistant states.     

Cloning and characterization of two processed pseudogenes and the cDNA for the murine U1 snRNP-specific protein C

Genes for the snRNP proteins U1-70K, U1-A, Sm-B′/B, Sm-D1 and Sm-E have been isolated from various metazoan species. The genes for Sm-D1 and Sm-E, which were isolated from a murine and human source respectively, appear to belong to a multigene family. It has been suggested that also for the mammalian U1-C protein such a multigene family exists. With the human U1-C cDNA as a probe, two genes containing sequences homologous to the probe sequence were isolated from a mouse genomic library. Simultaneously, a murine U1-C cDNA was isolated from a mouse cDNA library. This 0.74 kb cDNA contains an open reading frame (ORF) of 477 bp encoding a polypeptide of 159 amino acids (aa) which differs at only one position (position 65) from the human U1-C protein. One of the isolated U1-C genes contains an ORF as well and shares 92% nucleotide sequence identity with the mouse U1-C cDNA. The features of this gene, in particular the absence of introns, the acquisition of a 3′ poly(A) tail and flanking direct repeats, indicate that it represents a processed pseudogene. At the predicted aa sequence level, substitutions of conserved residues at functionally important positions are observed, strongly suggesting that expression of this gene would not lead to a functional polypeptide. The second U1-C gene appeared to be a pseudogene as well because it is also intronless and contains a frameshift mutation compared to the ORF in the mouse U1-C cDNA. The characterization of these two pseudogenes points to the existence of a U1-C multigene family in mice. Furthermore, comparison of aa sequences of the murine, human and Xenopus U1-C shows that the protein is highly conserved through evolution. Since the Xenopus U1-C differs from the two mammalian counterparts solely at a number of positions in the C-terminal region, it can be concluded that aa changes are less well tolerated in the N-terminal region of U1-C than in the rest of the protein.