Multiple calmodulin mRNA species are derived from two distinct genes.

We have observed three calmodulin mRNA species in rat tissues. In order to know from how many expressed genes they are derived, we have investigated the genomic organization of calmodulin genes in the rat genome. From a rat brain cDNA library, we obtained two kinds of cDNAs (pRCM1 and pRCM3) encoding authentic calmodulin. DNA sequence analysis of these cDNA clones revealed substitutions of nucleotides at 73 positions of 450 nucleotides in the coding region, although the amino acid sequences of these calmodulins are exactly the same. DNA sequences in the 5' and 3' noncoding regions are quite different between these two cDNAs. From these results, we conclude that they are derived from two distinct bona fide calmodulin genes, CaMI (pRCM1) and CaMII (pRCM3). Total genomic Southern hybridization suggested four distinct calmodulin-related genes in the rat genome. By cloning and sequencing the calmodulin-related genes from rat genomic libraries, we demonstrated that the other two genes are processed pseudogenes generated from the CaMI (lambda SC9) and CaMII (lambda SC8) genes, respectively, through an mRNA-mediated process of insertions. Northern blotting showed that the CaMI gene is transcribed in liver, muscle, and brain in similar amounts, whereas the CaMII gene is transcribed mainly in brain. S1 nuclease mapping indicated that the CaMI gene produced two mRNA species (1.7 and 4 kilobases), whereas the CaMII gene expressed a single mRNA species (1.4 kilobases).     

Structure and evolution of a mouse tRNA gene cluster encoding tRNAAsp, tRNAGly and tRNAGlu and an unlinked, solitary gene encoding tRNAAsp

We have sequenced mouse tRNA genes from two recombinant lambda phage. An 1800 bp sequence from one phage contains 3 tRNA genes, potentially encoding tRNAAsp, tRNAGly, and tRNAGlu, separated by spacer sequences of 587 bp and 436 bp, respectively. The mouse tRNA gene cluster is homologous to a rat sequence (Sekiya et al., 1981, Nucleic Acids Res. 9, 2239-2250). The mouse and rat tRNAAsp and tRNAGly coding regions are identical. The tRNAGlu coding regions differ at two positions. The flanking sequences contain 3 non-homologous areas: a c. 100 bp insertion in the first mouse spacer, short tandemly repeated sequences in the second spacers and unrelated sequences at the 3' ends of the clusters. In contrast, most of the flanking regions are homologous, consisting of strings of consecutive, identical residues (5-17 bp) separated by single base differences and short insertions/deletions. The latter are often associated with short repeats. The homology of the flanking regions is c. 75%, similar to other murine genes. The second lambda clone contains a solitary mouse tRNAAsp gene. The coding region is identical to that of the clustered tRNAAsp gene. The 5' flanking regions of the two genes contain homologous areas (10-25 bp) separated by unrelated sequences. Overall, the flanking regions of the two mouse tRNAAsp genes are less homologous than those of the mouse and rat clusters.     

Structural organization of multiple rat calmodulin genes

Elsewhere, we have reported the structure of a rat calmodulin gene and two distinct rat calmodulin cDNAs, pRCM1 and pRCM3. Here, I report the cloning and sequencing of the third calmodulin cDNA (pRCM4) and two additional rat calmodulin genes. The original calmodulin gene is named CaM I (pRCM1) and the newly discovered calmodulin genes are named CaM II (pRCM3) and CaM III (pRCM4). CaM II spans about 10 x 10(3) base-pairs and consisted of five exons, while CaM III spans about 7.2 x 10(3) base-pairs and consisted of six exons. One of the introns (intron 3) observed in CaM I and CaM III is lost in CaM II. Otherwise, the intron/exon organization of these genes is exactly the same. In all calmodulin genes, the first intron separates the initiation codon (ATG) from the coding region of the protein. Northern blotting showed that CaM I is transcribed primarily into 1.7 x 10(3) base-pair mRNA in various tissues examined and 4.0 x 10(3) base-pair mRNA mainly in skeletal muscle, CaM II is transcribed into 1.4 x 10(3) base-pair mRNA almost exclusively in brain and CaM III is transcribed predominantly into 2.3 x 10(3) base-pair mRNA and faintly into 1.0 x 10(3) base-pair mRNA mainly in skeletal muscle and brain. DNA sequences in the promoter-regulator regions of these genes are partly homologous but essentially distinct and possess a number of direct repeats, palindromes and feasible stem-loop structures. Together with these, I report here the structures of the third and fourth calmodulin retropseudogenes.     

Molecular analysis of elements inserted into mouse gamma-actin processed pseudogenes.

DNA from ten mouse genomic clones, each containing distinct gamma-actin processed pseudogenes, was subjected to electron microscopic heteroduplex analysis, and in three cases (lambda mA36, lambda mA118 and lambda mA119) the heteroduplex formed with the DNA of a reference clone was found to be interrupted by a single-stranded loop. The genomic regions corresponding to these loops were subjected to structural analysis and they were found to represent different elements (IEs) inserted into the pseudogenes in a manner that gave rise to short target-site direct repeats. IE 36 (500 base-pairs in length) was found to be an intercisternal A-particle solo long terminal repeat (LTR), a 46 nucleotide region of which had undergone five-fold tandem amplification and subsequent mutation. IE 119 (501 base-pairs in length) was also a solo LTR, bearing similarity to the recently-described GLN-3 class of murine retroviral-like elements. IE 118 (865 base-pairs in length) is repeated 1000-2000 times in the mouse genome. It is not related to any known class of mobile elements, but does possess some sequence motifs that suggest it may be an LTR of a hitherto unrecognized family of retroviral-like elements. It also possesses a 26 out of 27 nucleotide identity to a region of the flanking pseudogene, suggesting that it may have suffered gene conversion.     

Two novel frameshift mutations associated with the presence of direct repeats of the LDL receptor gene in familial hypercholesterolemia

Two novel frameshift mutations were detected in the mutant LDL receptor genes responsible for familial hypercholesterolemia. One was a 5-bp insertion at codon 395 in exon 9, and the other was a one nucleotide deletion at codon 531 in exon 11. Both mutations alter the reading frame and consequently produce a premature stop codon in the region of the mature LDL receptor homologous to the epidermal growth factor (EGF) precursor. With regard to the mechanism responsible for the generation of these frameshift mutations, strand slipped mispairing mediated by short direct repeats is considered to be the most likely. The findings seem to support the hypothesis that a short direct repeat in DNA sequence can have a profound influence on the stability of a given gene and promote human gene mutations.     

Nucleotide sequences at the boundaries between gene and insertion regions in the rDNA of Drosophila melanogaster

Ribosomal RNA genes interrupted by type 1 insertions of 1 kb and 0.5 kb have been sequenced through the insertion region and compared with an uninterrupted gene. The 0.5 kb insertion is flanked by a duplication of a 14 bp segment that is present once in the uninterrupted gene; the 1 kb insertion is flanked by a duplication of 11 of these 14 bp. Short insertions are identical in their entire length to downstream regions of long insertions. No internal repeats occur in the insertion. The presence of target site duplications suggests that type 1 insertions arose by the introduction of transposable elements into rDNA. Short sequence homologies between the upstream ends of the insertions and the 28S' boundaries of the rRNA coding region suggest that short type 1 insertions may have arisen by recombination from longer insertions.We have sequenced both boundaries of two molecules containing type 2 insertions and the upstream boundary of a third; the points of interruption at the upstream boundary (28S' site) differ from each other in steps of 2 bp. Between the boundary in the 0.5 kb type 1 insertion and the type 2 boundaries there are distances of 74, 76, and 78 bp. At the downstream boundary (28S' site) the two sequenced type 2 insertions are identical. The rRNA coding region of one molecule extends across the insertion without deletion or duplication, but a 2 bp deletion in the RNA coding region is present in the second molecule. Stretches of 13 or 22 adenine residues occur at the downstream (28S') end of the two type 2 insertions.     

The complete nucleotide sequence of wild rice (Oryza nivara) chloroplast genome: first genome wide comparative sequence analysis of wild and cultivated rice

We determined the complete nucleotide sequence of the chloroplast genome of wild rice, Oryza nivara and compared it with the corresponding published sequence of relative cultivated rice, Oryza sativa. The genome was 134,494 bp long with a large single-copy region of 80,544 bp, a small single-copy region of 12,346 bp and two inverted repeats of 20,802 bp each. The overall A+T content was 61.0%. The O. nivara chloroplast genome encoded identical functional genes to O. sativa in the same order along the genome. On the other hand, detailed analysis revealed 57 insertion, 61 deletion and 159 base substitution events in the entire chloroplast genome of O. nivara. Among substitutions, transversions were much higher than transitions with the former even more frequent than the latter in the coding region. Most of the insertions/deletions were single-base but a few large length mutations were also detected. The frequency of insertion/deletion events was more in the coding region within inverted repeats. In contrast, a very few substitution events were identified in the coding region. Polymorphism was observed among rice cultivars at loci of large insertion/deletion events. This is the first report describing comparative and genome wide chloroplast analysis between a wild and cultivated crop.     

Cytochrome oxidase subunit II gene of rice has an insertion sequence within the intron.

We have isolated and sequenced the cytochrome oxidase subunit II gene from rice (Oryza sativa L. var Labelle). The overall structural organization of this gene is very similar to that of the maize gene. This gene contains an intron in a position identical to the intron in the maize gene. However, the intron in the rice gene is longer than that of the maize gene largely due to a 461 bp insertion sequence, which has inverted repeats at its termini and is flanked by direct repeats, characteristic of transposable elements. Apart from this insertion sequence, the remainder of the intron sequence is strikingly homologous to that of maize (98.6% homology), suggesting a possible functional or structural role. The coding regions of the two genes exhibit 99.5% nucleotide sequence homology and their deduced amino acid sequences are identical. Similarly, the 3'-noncoding regions, except for several small insertions and deletions, show complete sequence homology. On the contrary, no sequence homology is detected in the 5'-noncoding regions.     

Complete nucleotide and encoded amino acid sequence of a mammalian myosin heavy chain gene. Evidence against intron-dependent evolution of the rod

The complete nucleotide sequence and exon/intron structure of the rat embryonic skeletal muscle myosin heavy chain (MHC) gene has been determined. This gene comprises 24 X 10(3) bases of DNA and is split into 41 exons. The exons encode a 6035 nucleotide (nt) long mRNA consisting of 90 nt of 5' untranslated, 5820 nt of protein coding and 125 nt of 3' untranslated sequence. The rat embryonic MHC polypeptide is encoded by exons 3 to 41 and contains 1939 amino acid residues with a calculated Mr of 223,900. Its amino acid sequence displays the structural features typical for all sarcomeric MHCs, i.e. an amino-terminal "globular" head region and a carboxy-terminal alpha-helical rod portion that shows the characteristics of a coiled coil with a superimposed 28-residue repeat pattern interrupted at only four positions by "skip" residues. The complex structure of the rat embryonic MHC gene and the conservation of intron locations in this and other MHC genes are indicative of a highly split ancestral sarcomeric MHC gene. Introns in the rat embryonic gene interrupt the coding sequence at the boundaries separating the proteolytic subfragments of the head, but not at the head/rod junction or between the 28-residue repeats present within the rod. Therefore, there is little evidence for exon shuffling and intron-dependent evolution by gene duplication as a mechanism for the generation of the ancestral MHC gene. Rather, intron insertion into a previously non-split ancestral MHC rod gene consisting of multiple tandemly arranged 28-residue-encoding repeats, or convergent evolution of an originally non-repetitive ancestral MHC rod gene must account for the observed structure of the rod-encoding portion of present-day MHC genes.     

Sequence simplicity and evolution of the 3′ untranslated region of the histone H1° Gene

The H1° gene has a long 3′ untranslated region (3′UTR) of 1,125 nucleotides in the rat and 1,310 in humans. Analysis of the sequences shows that they have features of simple DNA that suggest involvement of replication slippage in their evolution. These features include the length imbalance between the rat and human sequences; the abundance of single-base repeats, two-base runs and other simple motifs clustered along the sequence; and the presence of single-base repeat length polymorphisms in the rat and mouse sequences. Pairwise comparisons show numerous short insertions/deletions, often flanked by direct repeats. In addition, a proportion of short insertions/deletions results from length differences in conserved single-base repeats. Quantification of the sequence simplicity shows that simple sequences have been more actively incorporated in the human lineage than in the rodent lineage. The combination of insertions/deletions and nucleotide substitutions along the sequence gives rise to three main regions of homology: a highly variable central region flanked by more conserved regions nearest the coding region and the polyA addition site. Correspondence to: P. Suau     

Genetic analysis of the cell binding domain region of the chicken fibronectin gene

We have determined the nucleotide sequence of the cell binding domain region of the chicken fibronectin gene and analyzed it evolutionaly. We present here the complete nucleotide sequence of 4.3 kb HindIII/EcoRI segment from the clone lambda FC23 of the chicken fibronectin gene. There were five exons in this segment. When we lined up the amino acid of exons 28, 29 and 31, three alignments, known as the Type III repeat, appeared. Tetrapeptide, -RGDS-, called the cell binding domain, existed in the second repeat, coding exon 30. It was presumed that the Type III repeats were composed of two exons in the chicken gene, the same as in the rat and humans. We found repeatedly appearing amino-acid sequences such as -TIT- (three arrays in these Type III repeats) but also found one of the amino acids substituted in the tripeptide in these Type III repeats (seven arrays). We analyzed these repeats from the point of view of evolution. We used three of the nucleotide sequences (12-18 bp) coding such -TIT- repeats as a unit length for comparing the various homologies after dividing the coding region into 56 segments. The mutual homology of the divided segments to each one of three showed 53% on average. On the other hand, the mutual nucleotide homology of the Type III repeat was 44%. This suggested that the Type III repeat may have been developed by frequent duplication of small gene units.     

Random insertion mutagenesis of the intracellular pathogen Rhodococcus equi using transposomes

The identification of virulence factors in Rhodococcus equi has been severely hampered by the lack of a method for in vivo random insertion mutagenesis. This study reports the use of transposomes to generate random insertions of a gene conferring kanamycin resistance into the genome of R. equi ATCC 33701. Southern hybridisation using the kanamycin resistance gene as probe showed that insertion of transposome is random. This was confirmed following nucleotide sequence analysis of the junction between the transposome and chromosomal DNA. The presence of a 9 bp duplication of the target sequence showed that random integration of the transposome was due to a bona fide Tn5 transposition event.     

Sequence and comparative analysis of the rabbit alpha-like globin gene cluster reveals a rapid mode of evolution in a G + C-rich region of mammalian genomes

A sequence of 10,621 base-pairs from the alpha-like globin gene cluster of rabbit has been determined. It includes the sequence of gene zeta 1 (a pseudogene for the rabbit embryonic zeta-globin), the functional rabbit alpha-globin gene, and the theta 1 pseudogene, along with the sequences of eight C repeats (short interspersed repeats in rabbit) and a J sequence implicated in recombination. The region is quite G + C-rich (62%) and contains two CpG islands. As expected for a very G + C-rich region, it has an abundance of open reading frames, but few of the long open reading frames are associated with the coding regions of genes. Alignments between the sequences of the rabbit and human alpha-like globin gene clusters reveal matches primarily in the immediate vicinity of genes and CpG islands, while the intergenic regions of these gene clusters have many fewer matches than are seen between the beta-like globin gene clusters of these two species. Furthermore, the non-coding sequences in this portion of the rabbit alpha-like globin gene cluster are shorter than in human, indicating a strong tendency either for sequence contraction in the rabbit gene cluster or for expansion in the human gene cluster. Thus, the intergenic regions of the alpha-like globin gene clusters have evolved in a relatively fast mode since the mammalian radiation, but not exclusively by nucleotide substitution. Despite this rapid mode of evolution, some strong matches are found 5' to the start sites of the human and rabbit alpha genes, perhaps indicating conservation of a regulatory element. The rabbit J sequence is over 1000 base-pairs long; it contains a C repeat at its 5' end and an internal region of homology to the 3'-untranslated region of the alpha-globin gene. Part of the rabbit J sequence matches with sequences within the X homology block in human. Both of these regions have been implicated as hot-spots for recombination, hence the matching sequences are good candidates for such a function. All the interspersed repeats within both gene clusters are retroposon SINEs that appear to have inserted independently in the rabbit and human lineages.     

Coding region deletions associated with the major form of rDNA interruption in Drosophila

The nucleotide sequences at and around the termini of 5 kb type 1 interruptions in three separate clones of D. melanogaster rDNA repeats have been determined, and have been compared with the sequence of the corresponding region of an insertion-free rDNA repeat. All three interrupted rDNA repeats contain a small deletion of 28S rRNA coding material at the left coding/insertion sequence junction. A second deletion was found in one of the three clones, ad other aberrations were suggested by the results of restriction enzyme digestions of unfractionated rDNA. The termini of 5 kb type 1 rDNA insertions in D. melanogaster were also compared with the corresponding regions of 28S rDNA interruptions in D. virilis: the insertion site is identical in the two species, but the termini of the two species' interruptions show no homology. I sequenced a 1.1 kb region of the 5 kb type 1 D. melanogaster rDNA interruption that covers the sequences of the 1 kb and 0.5 kb insertions. There is 98% homology between the rightmost 1 kb of the 5 kb interruption and the sequences of the shorter insertions. Data suggest that Drosophila rDNA interruptions arose as a transposable element, and that divergence had included length alterations generated by unequal crossing over.