Conserved sequences in both coding and 5'' flanking regions of mammalian opal suppressor tRNA genes.

The rabbit genome encodes an opal suppressor tRNA gene. The coding region is strictly conserved between the rabbit gene and the corresponding gene in the human genome. The rabbit opal suppressor gene contains the consensus sequence in the 3' internal control region but like the human and chicken genes, the rabbit 5' internal control region contains two additional nucleotides. The 5' flanking sequences of the rabbit and the human opal suppressor genes contain extensive regions of homology. A subset of these homologies is also present 5' to the chicken opal suppressor gene. Both the rabbit and the human genomes also encode a pseudogene. That of the rabbit lacks the 3' half of the coding region. Neither pseudogene has homologous regions to the 5' flanking regions of the genes. The presence of 5' homologies flanking only the transcribed genes and not the pseudogenes suggests that these regions may be regulatory control elements specifically involved in the expression of the eukaryotic opal suppressor gene. Moreover the strict conservation of coding sequences indicates functional importance for the opal suppressor tRNA genes.     

Structural conservation and variation in the mitochondrial control region of fringilline finches (Fringilla spp.) and the greenfinch (Carduelis chloris)

We sequenced the entire control region and portions of flanking genes (tRNA(Phe), tRNA(Glu), and ND6) in the common chaffinch (Fringilla coelebs), blue chaffinch (F. teydea), brambling (F. montifringilla), and greenfinch (Carduelis chloris). In these finches the control region is similar in length (1,223-1,237 bp) and has the same flanking gene order as in other birds, and contains a putative TAS element and the highly conserved CSB-1 and F, D, and C boxes recognizable in most vertebrates. Cloverleaf-like structures associated with the TAS element at the 5' end and CSB-1 at the 3' end of the control region may be involved with the stop and start of D-loop synthesis, respectively. The pattern of nucleotide and substitution bias is similar to that in other vertebrates, and consequently the finch control region can be subdivided into a central, conserved G-rich domain (domain II) flanked by hypervariable 5'-C-rich (domain I) and 3'-AT-rich (domain III) segments. In pairwise comparisons among finch species, the central domain has unusually low transition/transversion ratios, which suggests that increased G + T content is a functional constraint, possibly for DNA primase efficiency. In finches the relative rates of evolution vary among domains according to a ratio of 4.2 (domain III) to 2.2 (domain I) to 1 (domain II), and extensively among sites within domains I and II. Domain I and III sequences are extremely useful in recovering intraspecific phylogeographic splits between populations in Africa and Europe, Madeira, and a basal lineage in Nefza, Tunisia. Domain II sequences are highly conserved, and are therefore only useful in conjunction with sequences from domains I and III in phylogenetic studies of closely related species.      

Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster

Mitochondrial DNA (mtDNA) regions corresponding to two major tRNA gene clusters were amplified and sequenced for the Japanese pit viper, himehabu. In one of these clusters, which in most vertebrates characterized to date contains three tightly connected genes for tRNA(Ile), and tRNA(Gln), and tRNA(Met), a sequence of approximately 1.3 kb was found to be inserted between the genes for tRNA(Ile) and tRNA(Gln). The insert consists of a control-region-like sequence possessing some conserved sequence blocks, and short flanking sequences which may be folded into tRNA(Pro), tRNA(Phe), and tRNA(Leu) genes. Several other snakes belonging to different families were also found to possess a control-region-like sequence and tRNA(Leu) gene between the tRNA(Ile)and tRNA(Gln) genes. We also sequenced a region surrounded by genes for cytochrome b and 12S rRNA, where the control region and genes for tRNA(Pro) and tRNA(Phe) are normally located in the mtDNAs of most vertebrates. In this region of three examined snakes, a control-region- like sequence exists that is almost completely identical to the one found between the tRNA(Ile) and tRNA(Gln) genes. The mtDNAs of these snakes thus possess two nearly identical control-region-like sequences which are otherwise divergent to a large extent between the species. These results suggest that the duplicate state of the control-region- like sequences has long persisted in snake mtDNAs, possibly since the original insertion of the control-region-like sequence and tRNA(Leu) gene into the tRNA gene cluster, which occurred in the early stage of the divergence of snakes. It is also suggested that the duplicated control-region-like sequences at two distant locations of mtDNA have evolved concertedly by a mechanism such as frequent gene conversion. The secondary structures of the determined tRNA genes point to the operation of simplification pressure on the T psi C arm of snake mitochondrial tRNAs.      

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.     

Cloning and sequence comparison of the mouse, human, and chicken engrailed genes reveal potential functional domains and regulatory regions.

We have isolated and characterized genomic DNA clones for the human and chicken homologues of the mouse En-1 and En-2 genes and determined the genomic structure and predicted protein sequences of both En genes in all three species. Comparison of these vertebrate En sequences with the Xenopus En-2 [Hemmati-Brivanlou et al., 1991) and invertebrate engrailed-like genes showed that the two previously identified highly conserved regions within the En protein ]reviewed in Joyner and Hanks, 1991] can be divided into five distinct subregions, designated EH1 to EH5. Sequences 5' and 3' to the predicted coding regions of the vertebrate En genes were also analyzed in an attempt to identify cis-acting DNA sequences important for the regulation of En gene expression. Considerable sequence similarity was found between the mouse and human homologues both within the putative 5' and 3' untranslated as well as 5' flanking regions. Between the mouse and Xenopus En-2 genes, shorter stretches of sequence similarity were found within the 3' untranslated region. The 5' untranslated regions of the mouse, chicken and Xenopus En-2 genes, however, showed no similarly conserved stretches. In a preliminary analysis of the expression pattern of the human En genes, En-2 protein and RNA were detected in the embryonic and adult cerebellum respectively and not in other tissues tested. These patterns are analogous to those seen in other vertebrates. Taken together these results further strengthen the suggestion that En gene function and regulation has been conserved throughout vertebrate evolution and, along with the five highly conserved regions within the En protein, raise an interesting question about the presence of conserved genetic pathways.     

Complete mitochondrial genome of Trypauchen vagina (Perciformes, Gobioidei)

The complete mitochondrial genome of Trypauchen vagina was determined first. The genome is 16,686 bp in length and consists of 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, and 2 main non-coding regions [the control region (CR) and the origin of the light strand replication], the gene composition and order of which was similar to most other vertebrates. The overall base composition of T. vagina is T 27.6%, C 27.6%, A 29.5%, and G 15.3%, with a slight A+T bias of 57.1%. In addition to the discrete and conserved sequence blocks, an incomplete tandem repeat unit is detected within the CR. This mitogenome sequence data would play an important role in population genetics and phylogenetic analysis of the Gobioidei.     

Complete mitochondrial genome of Oxuderces dentatus (Perciformes, Gobioidei)

The complete mitochondrial genome (mitogenome) of Oxuderces dentatus was determined first. The genome was 17,116?bp in length and consisted of 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, and 2 main non-coding regions [the control region (CR) and the origin of the light strand replication], the gene composition and order of which was similar to most other vertebrates. The overall base composition of the heavy strand was T 27.9%, C 26.8%, A 30.2%, and G 15.1%, with a slight A+T bias of 58.1%. In addition to the discrete and conserved sequence blocks, unusual long tandem repeat unit (three 150-bp tandem repeat units and an incomplete copy of 146?bp) was also detected within CR. This mitogenome sequence data would play an important role in population genetics and phylogenetic analysis of the Gobioidei.     

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.     

Complete mitochondrial genome of the nerippe fritillary butterfly, Argynnis nerippe (Lepidoptera: Nymphalidae)

The complete mitochondrial genome sequence of the nerippe fritillary butterfly, Argynnis nerippe, which is listed as an endangered species in Korea, is described with an emphasis on the A+T-rich region. The 15,140-bp long circular molecule consisted of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and 1 control region, known in insect as the A+T-rich region, as found in typical metazoans. The 329-bp long A+T-rich region located between srRNA and tRNA(Met) possessed the highest A/T content (95.7%) than any other region of the genome. Along with the several conserved sequences found typically in the lepidopteran insects the genome contained one tRNA(Met)-like and tRNA(Leu)(UUR)-like sequence in the A+T-rich region.     

Mutations in mitochondrial tRNA genes: a frequent cause of neuromuscular diseases.

We have sequenced the tRNA genes of mtDNA from patients with chronic progressive external ophthalmoplegia (CPEO) without detectable mtDNA deletions. Four point mutations were identified, located within highly conserved regions of mitochondrial tRNA genes, namely tRNA(Leu)(UAG), tRNA(Ser)(GCU), tRNA(Gly) and tRNA(Lys). One of these mutations (tRNA(Leu)(UAG)) was found in four patients with different forms of mitochondrial myopathy. An accumulation of three different tRNA point mutations (tRNA(Leu)(UAG)), tRNA(Ser)(GCU) and tRNA(Gly) was observed in a single patient, suggesting that mitochondrial tRNA genes represent hotspots for point mutations causing neuromuscular diseases.     

Nested Evolution of a tRNALeu(UAA) Group I Intron by both Horizontal Intron Transfer and Recombination of the Entire tRNA Locus

The origin and evolution of bacterial introns are still controversial issues. Here we present data on the distribution and evolution of a recently discovered divergent tRNA(Leu)(UAA) intron. The intron shows a higher sequence affiliation with introns in tRNA(Ile)(CAU) and tRNA(Arg)(CCU) genes in alpha- and beta-proteobacteria, respectively, than with other cyanobacterial tRNA(Leu)(UAA) group I introns. The divergent tRNA(Leu)(UAA) intron is sporadically distributed both within the Nostoc and the Microcystis radiations. The complete tRNA gene, including flanking regions and intron from Microcystis aeruginosa strain NIVA-CYA 57, was sequenced in order to elucidate the evolutionary pattern of this intron. Phylogenetic reconstruction gave statistical evidence for different phylogenies for the intron and exon sequences, supporting an evolutionary model involving horizontal intron transfer. The distribution of the tRNA gene, its flanking regions, and the introns were addressed by Southern hybridization and PCR amplification. The tRNA gene, including the flanking regions, were absent in the intronless stains but present in the intron-containing strains. This suggests that the sporadic distribution of this intron within the Microcystis genus cannot be attributed to intron mobility but rather to an instability of the entire tRNA(Leu)(UAA) intron-containing genome region. Taken together, the complete data set for the evolution of this intron can best be explained by a model involving a nested evolution of the intron, i.e., wherein the intron has been transferred horizontally (probably through a single or a few events) to a tRNA(Leu)(UAA) gene which is located within a unstable genome region.