tRNA Pro Research Articles

Efficient CRISPR/Cas9 system based on autonomously replicating plasmid with an AMA1 sequence and precisely targeted gene deletion in the edible fungus, Cordyceps militaris.

Cordyceps militaris is a popular edible fungus with important economic value worldwide. In this study, an efficient CRISPR/Cas9 genome‐editing system based on an autonomously replicating plasmid with an AMA1 sequence was constructed. Further, a precisely targeted gene deletion via homology‐directed repair was effectively introduced in C. militaris. Gene editing was successful, with efficiencies of 55.1% and 89% for Cmwc‐1 and Cmvvd, respectively. Precisely targeted gene deletion was achieved at an efficiency of 73.9% by a single guide RNA supplementation with donor DNAs. Double genes, Cmwc‐1 and Cmvvd, were edited simultaneously with an efficiency of 10%. Plasmid loss was observed under non‐selective culture conditions, which could permit recycling of the selectable marker and avoid the adverse effects of the CRISPR/Cas9 system on the fungus, which is beneficial for the generation of new cultivars. RNA Pol III promoters, endogenous tRNA Pro of C. militaris, and chimeric AfU6‐tRNA Gly can be used to improve the efficiency. Polyethylene glycol‐mediated protoplast transformation was markedly more efficient than Agrobacterium tumefaciens‐mediated transformation of C. militaris. To our knowledge, this is the first description of genome editing and precisely targeted gene deletion in mushrooms based on AMA1 plasmids. Our findings will enable the modification of multiple genes in both functional genomics research and strain breeding.

In silico model of mtDNA mutations effect on secondary and 3D structure of mitochondrial rRNA and tRNA in Leber’s hereditary optic neuropathy

The Leber’s hereditary optic neuropathy (LHON) is a rare disease caused by mitochondrial DNA (mtDNA) mutations. Beside primary mutations, the effect of secondary mtDNA mutations in still unclear. We examined the effect of secondary mtDNA mutations on secondary structure of different mitochondrial RNAs. Whole mitochondrial genome sequence of LHON patients has been obtained from in six non related pedigrees by Sanger sequencing method. The effect of mutations located in mitochondrial RNA genes was examined by creating in silico models of RNA secondary and regional 3D structure, accompanied by sequence conservation analysis. All three primary LHON mutations (m.3460G>A, m.11778G>A and m.14484 T>C) were revealed in study families. Four mutations in MT-RNR1 gene (m.750A>G, m.956delC, m.1438A>G and m.1555A>G) were identified and only an m.1555A>G causes significant changes of secondary structure of mitochondrial 12S ribosomal RNA (rRNA), while it is the only mutation which does not alter its 3D structure. Five mutations (m.1811A>G, m.2706A>G, m.2831G>A, m.3010G>A and m.3197T>C) were discovered in MT-RNR2 gene and all of them induced substantial alterations of mitochondrial 16S rRNA secondary structure. Significant changes of mitochondrial 16S rRNA 3D structure are caused by m.1811A>G, m.2706A>G, m.3010G>A and m.3197T>C. A single insertion variant (m.15986insG) has been found in the MT-TP gene which encodes mitochondrial transfer RNA for Proline (tRNA Pro). This mutation does not cause substantial changes of tRNA for Proline secondary structure, while the 3D geometry remains without major changes. Most of the mutation loci exhibited high level of sequence conservation. Presence of multiple mutations in a single family appears to cause more extensive changes in mitochondrial 12S and 16S rRNA, then their individual influence. The effect of discovered mutations on in silico modelled RNA structure is in a significant correlation with the present knowledge about the potential of these mutation to participate in the pathophysiology of LHON and other human diseases. The presence of certain multiple mitochondrial RNA mutations could be a possible explanation of LHON clinical presentation in some families. All revealed mutations have been evaluated for the first time in terms of in silico structural modelling. The application of bioinformatics tools such as secondary and 3D RNA structure prediction can have a great advantage in better understanding of the molecular standpoint of the LHON pathophysiology and clinical phenotype.

Genetic variation and population structure of spottybelly greenling (Hexagrammos agrammus) in Korean coasts analyzed by DNA markers emphasizing on microsatellites

Two nuclear microsatellite DNA loci along with partial sequences of mitochondrial DNA tRNA PRO -D-loop region were analyzed to assess the genetic diversity and population genetic structure of spottybelly greenling ( Hexagrammos agrammus ). A total of 85 individuals were investigated from two different locations at the East and the South coasts of Korea which are Imwonjin (IM) and Tongyeong (TN), respectively. In the analysis of genetic diversity, nucleotide diversities were low showing 0.01 whereas haplotype diversities were as high as 0.92 and 0.97 in IM and TN populations, respectively indicating population bottleneck after rapid growth of these populations. No significant genealogical branches or clusters were recognized on the neighbor-joining phylogenetic tree. Pairwise population statistics FST and the exact test of population differentiation from the analysis of microsatellite DNA loci demonstrated no significant genetic difference between two populations investigated in the Korean coasts. In addition, significant deviations from Hardy-Weinberg equilibrium and heterozygote deficiencies were found from the microsatellite DNA loci. The results of the study will help to make a plan for fisheries management of the species.

Complete mitochondrial genome of the chocolate hind Cephalopholis boenak (Pisces: Perciformes)

The complete mitochondrial (mt) genome of the chocolate hind Cephalopholis boenak was obtained in this study, which is the first complete mitogenome in the genus Cephalopholis. The circular mtDNA molecule was 16 771 bp in size and the overall nucleotide composition of the H-strand was 29.58% A, 27.36% T, 15.84% G, and 27.22% C, with an A + T bias. The complete mitogenome encoded 13 protein-coding genes, 2 rRNAs, 22 tRNAs and 1 control region (D-loop), with the gene arrangement and translation direction basically identical to other typical vertebrate mitogenomes. The control region located between the tRNA Pro and tRNA Phe genes and the length was 1 064 bp, rich in A + T (70.11%); it was composed of 13 complete continuity tandem repeat units, 5′-TACATATCTATGTACTTAA-3′, and one imperfect tandem repeat.

Phylogeography of sable (Martes zibellina L. 1758) in the southeast portion of its range based on mitochondrial DNA variation: highlighting the evolutionary history of the sable

To reveal phylogeographic features of sable (Martes zibellina) in the southeast part of its range, we analyzed variability of the mitochondrial DNA (mtDNA) cytochrome b gene, tRNA (Pro), tRNA (Thr) and control region (D-loop) sequences from 78 specimens in populations of the Russian Far East, northeast China, and Mongolia. Our results revealed the presence of 49 different haplotypes split into two major phylogenetic groups—clades A and B, the latter separated into two clades, B1 and B2. Comparative analysis of D-loop haplotypes in populations originating from the southeast (Russian Far East, China and Mongolia) and the west (northern Urals) portions of sable range indicated that all three mtDNA clades were present in different regional groups. However, highest diversity of clade B1 in northeast China and its nearly complete absence from the Urals suggest that the southeast sable range, being a refuge during Pleistocene glacial periods, can be considered the center of genetic diversification and possibly origin of this species. All divergence estimates fall within the Pleistocene suggesting that Quarternary glaciations played an important role in phylogeographic differentiation of sable.

Analysis of the complete mitochondrial genome sequence of Pielomastax zhengi

The complete mitochondrial genome sequence of Pielomastax zhengi was determined by using long PCR and conserved primer walking approaches. The results showed that the entire mitochondrial genome of Pielomastax zhengi is 15 602 bp long with A+T content of 71.8%. All 37 genes are conserved in the positions observed in those of Locusta migratoria. All the genes are closely assembled by leaving 10 intergenic spacers in between. Those intergenic spacers are 47 bp in total (excluding the A+T rich region), with individual size ranges from 1 bp to 20 bp. In addition, there are a total of 52 bp overlapping sequences among 14 genes, ranging from 1-8 bp. All protein-coding genes start with a typical initiation codon in insects, ATN. Twelve protein-coding genes use the usual TAA and TAG termination codons, whereas, the ND5 genes have an incomplete termination codon (T). Except the tRNASer (AGN), whose DHU arm is absent; all the other 21 tRNA genes have typical clover-leaf secondary structures. But in Pielomastax zhengi, the secondary structures of five tRNA (tRNA(Cys), tRNA(Lys), tRNA(Phe), tRNA(Pro), tRNA(Arg)) genes can not be predicted by the conventional methods as they have only 3-4, rather than 5 base pairs in the TψC arm, while, the tRNA(Lys) and tRNA(Arg) have only 4 base pairs in the anticode arm. The predicted secondary structures of lrRNA and srRNA have 6 domains with 44 helices and 3 domains with 30 helices, respectively. The results of the rRNA secondary structure comparison showed that Pielomastax zhengi is more closely related with Locusta migratoria tibetensis than Ruspolia dubia. Like most insects, the mitochondrial genome of Pielomastax zhengi has a non-coding A+T rich region containing a polythymidine stretch, which may be involved in the replication and/or translation initiation of other genes.

The complete mitochondrial genome of Macrobrachium nipponense

The complete mitochondrial (mt) genome sequence plays an important role in the accurate determination of phylogenetic relationships among metazoans. Herein, we determined the complete mt genome sequence, structure and organization of Macrobrachium nipponense ( M. nipponense) (GenBank ID: NC_015073.1) and compared it to that of Macrobrachium lanchesteri ( M. lanchesteri) and Macrobrachium rosenbergii ( M. rosenbergii). The 15,806 base pair (bp) M. nipponense mt genome, which is comprised of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and 2 ribosomal RNAs (rRNAs), is slightly larger than that of M. lanchesteri (15,694 bp, GenBank ID: NC_012217.1) and M. rosenbergii (15,772 bp, GenBank ID: NC_006880.1). The M. nipponense genome contains a high AT content (66.0%), which is a common feature among metazoan mt genomes. Compared with M. lanchesteri and M. rosenbergii, we found a peculiar non-coding region of 950 bp with a microsatellite-like (TA) 6 element and many hairpin structures. The 13 PCGs are comprised of a total of 3707 codons, excluding incomplete termination codons, and the most frequently used amino acid is Leu (16.0%). The predicted start codons in the M. nipponense mt genome include ATG, ATC and ATA. Seven PCGs use TAA as a stop codon, whereas two use TAG, three use T and only one uses TA. Twenty-three of the genes are encoded on the L strand, and ND1, ND4, ND5, ND4L, 12S rRNA, 16S rRNA, tRNA His , tRNA Pro , tRNA Phe , tRNA Val , tRNA Gln , tRNA Cys , tRNA Tyr and a tRNA Leu are encoded on the H strand. The two rRNAs of M. nipponense and M. rosenbergii are encoded on the H strand, whereas the M. lanchesteri rRNAs are encoded on the L stand.

The complete mitochondrial genome of the marbled rockfish Sebastiscus marmoratus (Scorpaeniformes, Scorpaenidae): Genome characterization and phylogenetic considerations

The complete mitochondrial genome sequence of the marbled rockfish Sebastiscus marmoratus (Scorpaeniformes, Scorpaenidae) was determined and phylogenetic analysis was conducted to elucidate the evolutionary relationship of the marbled rockfish with other Sebastinae species. This mitochondrial genome, consisting of 17301 bp, is highly similar to that of most other vertebrates, containing the same gene order and an identical number of genes or regions, including 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and one putative control region. Most of the genes are encoded on the H-strand, while the ND6 and seven tRNA genes (for Gln, Ala, Asn, Tyr, Ser (UCA), Glu, and Pro) are encoded on the L-strand. The reading frame of two pairs of genes overlapped on the same strand (the ATPase 8 and 6 genes overlapped by ten nucleotides; ND4L and ND4 genes overlapped by seven nucleotides). The possibly nonfunctional light-strand replication origin folded into a typical stem-loop secondary structure and a conserved motif (5'-GCCGG-3') was found at the base of the stem within the tRNA(Cys) gene. An extent termination-associated sequence (ETAS) and conserved sequence blocks (CSB) were identified in the control region, except for CSB-1; unusual long tandem repeats were found at the 3' end of the control region. Phylogenetic analyses supported the view that Sebastinae comprises four genera (Sebates, Hozukius, Helicolenus, and Sebasticus).

Two groups of phenylalanine biosynthetic operon leader peptides genes: a high level of apparently incidental frameshifting in decoding Escherichia coli pheL

The bacterial pheL gene encodes the leader peptide for the phenylalanine biosynthetic operon. Translation of pheL mRNA controls transcription attenuation and, consequently, expression of the downstream pheA gene. Fifty-three unique pheL genes have been identified in sequenced genomes of the gamma subdivision. There are two groups of pheL genes, both of which are short and contain a run(s) of phenylalanine codons at an internal position. One group is somewhat diverse and features different termination and 5′-flanking codons. The other group, mostly restricted to Enterobacteria and including Escherichia coli pheL, has a conserved nucleotide sequence that ends with UUC_CCC_UGA. When these three codons in E. coli pheL mRNA are in the ribosomal E-, P- and A-sites, there is an unusually high level, 15%, of +1 ribosomal frameshifting due to features of the nascent peptide sequence that include the penultimate phenylalanine. This level increases to 60% with a natural, heterologous, nascent peptide stimulator. Nevertheless, studies with different tRNAPro mutants in Salmonella enterica suggest that frameshifting at the end of pheL does not influence expression of the downstream pheA. This finding of incidental, rather than utilized, frameshifting is cautionary for other studies of programmed frameshifting.

Pathogenic mutations in antisense mitochondrial tRNAs

Pathogenic mutations in mitochondrial tRNAs are 6.5 times more frequent than in other mitochondrial genes. This suggests that tRNA mutations perturb more than one function. A potential additional tRNA gene function is that of templating for antisense tRNAs. Pathogenic mutations weaken cloverleaf secondary structures of sense tRNAs. Analyses here show similar effects for most antisense tRNAs, especially after adjusting for associations between sense and antisense cloverleaf stabilities. These results imply translational activity by antisense tRNAs. For sense tRNAs Ala and Ser UCN, pathogenicity associates as much with sense as with antisense cloverleaf formation. For tRNA Pro, pathogenicity seems associated only with antisense, not sense tRNA cloverleaf formation. Translational activity by antisense tRNAs is expected for the 11 antisense tRNAs processed by regular sense RNA maturation, those recognized by their cognate amino acid’s tRNA synthetase, and those forming relatively stable cloverleaves as compared to their sense counterpart. Most antisense tRNAs probably function routinely in translation and extend the tRNA pool (extension hypothesis); others do not (avoidance hypothesis). The greater the expected translational activity of an antisense tRNA, the more pathogenic mutations weaken its cloverleaf secondary structure. Some evidence for RNA interference, a more classical role for antisense tRNAs, exists only for tRNA Ser UCN. Mutation pathogenicity probably frequently results from a mixture of effects due to sense and antisense tRNA translational activity for many mitochondrial tRNAs. Genomic studies should routinely explore for translational activity by antisense tRNAs.

The complete mitochondrial genome sequence of the cutlassfish Trichiurus japonicus (Perciformes: Trichiuridae): Genome characterization and phylogenetic considerations

Mitochondrial genome sequence and structure analysis has become a powerful tool for studying molecular evolution and phylogenetic relationships. To understand the systematic status of Trichiurus japonicus in suborder Scombroidei, we determined the complete mitochondrial genome (mitogenome) sequence using the long-polymerase chain reaction (long-PCR) and shotgun sequencing method. The entire mitogenome is 16,796 bp in length and has three unusual features, including (1) the absence of tRNA Pro gene, (2) the possibly nonfunctional light-strand replication origin (O L) showing a shorter loop in secondary structure and no conserved motif (5'-GCCGG-3'), (3) two sets of the tandem repeats at the 5' and 3' ends of the control region. The three features seem common for Trichiurus mitogenomes, as we have confirmed them in other three T. japonicus individuals and in T. nanhaiensis. Phylogenetic analysis does not support the monophyly of Trichiuridae, which is against the morphological result. T. japonicus is most closely related to those species of family Scombridae; they in turn have a sister relationship with Perciformes members including suborders Acanthuroidei, Caproidei, Notothenioidei, Zoarcoidei, Trachinoidei, and some species of Labroidei, based on the current dataset of complete mitogenome. T. japonicus together with T. brevis, T. lepturus and Aphanopus carbo form a clade distinct from Lepidopus caudatus in terms of the complete Cyt b sequences. T. japonicus mitogenome, as the first discovered complete mitogenome of Trichiuridae, should provide important information on both genomics and phylogenetics of Trichiuridae.

The complete mitochondrial genome of the ridgetail white prawn Exopalaemon carinicauda Holthuis, 1950 (Crustacean: Decapoda: Palaemonidae) revealed a novel rearrangement of tRNA genes

The complete mitochondrial (mt) DNA sequence was determined for a ridgetail white prawn, Exopalaemon carinicauda Holthuis, 1950 (Crustacea: Decopoda: Palaemonidae). The mt genome is 15,730 bp in length, encoding a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes, which is typical for metazoans. The majority-strand consists of 33.6% A, 23.0% C, 13.4% G, and 30.0% T bases (AT skew = 0.057; GC skew = − 0.264). A total of 1045 bp of non-coding nucleotides were observed in 16 intergenic regions, including a major A + T rich (79.7%) noncoding region (886 bp). A novel translocation of tRNA Pro and tRNA Thr was found when comparing this genome with the pancrustacean ground pattern indicating that gene order is not conserved among caridean mitochondria. Furthermore, the rate of Ka/ Ks in 13 protein-coding genes between three caridean species is much less than 1, which indicates a strong purifying selection within this group. To investigate the phylogenetic relationship within Malacostraca, phylogenetic trees based on currently available malacostracan complete mitochondrial sequences were built with the maximum likelihood and Bayesian models. All analyses based on nucleotide and amino acid data strongly support the monophyly of Decapoda. The Penaeidae, Reptantia, Caridea, and Meiura clades were also recovered as monophyletic groups with strong statistical support. However, the phylogenetic relationships within Pleocyemata are unstable, as represented by the inclusion or exclusion of Caridea.

Structural Determinants Characteristic to AARS Subclasses and tRNA-splicing Endonuclease in Eukaryotes

We compare and analyse the whole set of cytoplasmic, nonorganellar transfer RNA genes from 22 eukaryal genomes. Grouping this whole set into sets of isoacceptors we have elucidated structural elements that are characteristic to individual isoacceptor tRNAs within the subclasses of aminoacyl tRNA synthetases. Further, we have observed structural motifs straddling the exon-intron boundaries, which includes selective occurrence of both symmetric- 3-4-3 and asymmetric-3-4-(4, 5) bulge-helix-bulge-like structural motifs. Among all the tRNA isoacceptors, Ile, Leu, Ser, Pro, Met, Arg, and Tyr harbor BHB-like secondary structures at exon-intron boundaries. The structural signatories at the exon-intron boundaries appear to contribute to the specificity of splice site recognition.

Functional guanine–arginine interaction between tRNA Pro and prolyl-tRNA synthetase that couples binding and catalysis

Aminoacyl-tRNA synthetases catalyze the attachment of specific amino acids to their cognate tRNAs. Specific aminoacylation is dictated by a set of recognition elements that mark tRNA molecules as substrates for particular synthetases. Escherichia coli prolyl-tRNA synthetase (ProRS) has previously been shown to recognize specific bases of tRNA Pro in both the anticodon domain, which mediate initial complex formation, and in the acceptor stem, which is proximal to the site of catalysis. In this work, we unambiguously define the molecular interaction between E. coli ProRS and the acceptor stem of cognate tRNA Pro. Oxidative cross-linking studies using 2'-deoxy-8-oxo-7,8-dihydroguanosine-containing proline tRNAs identify a direct interaction between a critical arginine residue (R144) in the active site of E. coli ProRS and the G72 residue in the acceptor stem of tRNA Pro. Assays conducted with motif 2 loop variants and tRNA mutants wherein specific atomic groups of G72 were deleted, are consistent with a functionally important hydrogen-bonding network between R144 and the major groove of G72. These results taken together with previous studies suggest that breaking this key contact uncouples the allosteric interaction between the anticodon domain and the aminoacylation active site, providing new insights into the communication network that governs the synthetase–tRNA interaction.

Evolution of acceptor stem tRNA recognition by class II prolyl-tRNA synthetase

Aminoacyl-tRNA synthetases (AARS) are an essential family of enzymes that catalyze the attachment of amino acids to specific tRNAs during translation. Previously, we showed that base-specific recognition of the tRNAPro acceptor stem is critical for recognition by Escherichia coli prolyl-tRNA synthetase (ProRS), but not for human ProRS. To further delineate species-specific differences in acceptor stem recognition, atomic group mutagenesis was used to probe the role of sugar–phosphate backbone interactions in recognition of human tRNAPro. Incorporation of site-specific 2′-deoxynucleotides, as well as phosphorothioate and methylphosphonate modifications within the tRNA acceptor stem revealed an extensive network of interactions with specific functional groups proximal to the first base pair and the discriminator base. Backbone functional groups located at the base of the acceptor stem, especially the 2′-hydroxyl of A66, are also critical for aminoacylation catalytic efficiency by human ProRS. Therefore, in contrast to the bacterial system, backbone-specific interactions contribute significantly more to tRNA recognition by the human enzyme than base-specific interactions. Taken together with previous studies, these data show that ProRS-tRNA acceptor stem interactions have co-adapted through evolution from a mechanism involving ‘direct readout’ of nucleotide bases to one relying primarily on backbone-specific ‘indirect readout’.

The many applications of acid urea polyacrylamide gel electrophoresis to studies of tRNAs and aminoacyl-tRNA synthetases

Here we describe the many applications of acid urea polyacrylamide gel electrophoresis (acid urea PAGE) followed by Northern blot analysis to studies of tRNAs and aminoacyl-tRNA synthetases. Acid urea PAGE allows the electrophoretic separation of different forms of a tRNA, discriminated by changes in bulk, charge, and/or conformation that are brought about by aminoacylation, formylation, or modification of a tRNA. Among the examples described are (i) analysis of the effect of mutations in the Escherichia coli initiator tRNA on its aminoacylation and formylation; (ii) evidence of orthogonality of suppressor tRNAs in mammalian cells and yeast; (iii) analysis of aminoacylation specificity of an archaeal prolyl-tRNA synthetase that can aminoacylate archaeal tRNA Pro with cysteine, but does not aminoacylate archaeal tRNA Cys with cysteine; (iv) identification and characterization of the AUA-decoding minor tRNA Ile in archaea; and (v) evidence that the archaeal minor tRNA Ile contains a modified base in the wobble position different from lysidine found in the corresponding eubacterial tRNA.

Identification and characterization of a new conjugation/type IVA secretion system ( trb/tra) of Legionella pneumophila Corby localized on two mobile genomic islands

Horizontal gene transfer probably contributes to evolution of Legionella pneumophila and its adaptation to different environments. Although horizontal gene transfer was observed in Legionella, the mechanism is still not specified. In this study we identified and analysed a new type of conjugation/type IVA secretion system ( trb/tra) of L. pneumophila Corby, a virulent human isolate. Two similar versions of this conjugation system were identified, localized on two different genomic islands (Trb-1, 42,710 bp and Trb-2, 34,434 bp). Trb-1 and Trb-2 are integrated within the tRNA Pro gene ( lpc2778) and the tmRNA gene ( lpc0164), respectively. Both islands exhibit an oriT region and both can be excised from the chromosome forming episomal circles. Trb-1 was analysed in more detail. It is active and can be horizontally transferred to other Legionella strains by conjugation and then integrated into the genome in a site-specific manner within the tRNA Pro gene. We characterized the sequence of the excision and integration sites of Trb-1 in three different L. pneumophila strains. Here we demonstrate that L. pneumophila exhibits a functional oriT region and that genomic islands in Legionella can be mobilized and conjugated to other species of Legionella. Thus, we describe for the first time a mechanism that may explain the observed horizontal transfer of chromosomal DNA in Legionella.

The complete mitochondrial genome of the whiting, Merlangius merlangus and the haddock, Melanogrammus aeglefinus: A detailed genomic comparison among closely related species of the Gadidae family

We determined the first complete mitochondrial DNA (mtDNA) sequences for the whiting ( Merlangius merlangus, family Gadidae, order Gadiformes) and the haddock ( Melanogrammus aeglefinus, family Gadidae, order Gadiformes). The entire mitogenomes were amplified and sequenced by primer walking using newly designed specific internal primers. Lengths were 16,569 and 16,585 bases for whiting and haddock respectively, lengths which lie within the range of previously reported gadiform sequences from Atlantic cod ( Gadus morhua, 16,696 bases) and walleye pollock ( Theragra chalcogramma, 16,570 bases). Gene arrangement in both species conformed to the order seen in most vertebrate mitochondrial genomes. We identified a long intergenic spacer located between the tRNA Thr and tRNA Pro genes (of 100 and 70 bp long for whiting and haddock, respectively), as previously described for other species of the order Gadiformes. Using nucleotide and amino acid divergence data of four complete gadoid mitogenomes ( M. merlangius, M. aeglefinus, G. morhua and T. chalcogramma), we examined in detail the relative mtDNA mutation patterns across genes and among Gadidae species and tested for the performance of each protein-coding, transfer RNA and ribosomal RNA gene in depicting the expected phylogeny among the four species, as compared with the whole genome dataset. This comparison may be particularly useful in phylogenetic analyses of such a diverse fish family, as well as for the understanding of the patterns of nucleotide substitution of the mtDNA at low levels of divergence.

Bacillus subtilis tRNA Pro with the anticodon mo 5UGG can recognize the codon CCC

In Bacillus subtilis, four codons, CCU, CCC, CCA, and CCG, are used for proline. There exists, however, only one proline-specific tRNA having the anticodon mo 5UGG. Here, we found that this tRNA Pro(mo 5UGG) can read not only the codons CCA, CCG and CCU but also CCC, using an in vitro assay system. This means that the first nucleoside of its anticodon, 5-methoxyuridine (mo 5U), recognizes A, G, U and C. On the other hand, it was reported that mo 5U at the first position of the anticodon of tRNA Val(mo 5UAC) can recognize A, G, and U but not C. A comparison of the structure of the anticodon stem and loop of tRNA Pro(mo 5UGG) with those of other tRNAs containing mo 5U at the first positions of the anticodons suggests that a modification of nucleoside 32 to pseudouridine (Ψ) enables tRNA Pro(mo 5UGG) to read the CCC codon.

Murine leukemia virus with a primer-binding site complementary to tRNA Lys,3 adapts to select new tRNAs for replication following extended in vitro culture

The preference of MuLV for the selection of tRNA Pro as a replication primer was investigated by altering the primer-binding site (PBS) to be complementary to tRNA Lys,3. MuLV-based vectors with a PBS complementary to tRNA Lys,3 were found to be approximately 2-fold less infectious than vectors with the wild-type PBS complementary to tRNA Pro. MuLV with a PBS complementary to tRNA Lys,3 was replication competent and maintained the PBS during early stages of in vitro culture. Upon extended culture, PBS were isolated which were complementary to tRNA Arg. A second MuLV was generated in which the region upstream of the PBS which is predicted to form an RNA stem loop structure was altered so that the nucleotide sequence within the loop would be complementary to the anticodon of tRNA Lys,3. The virus with both the U5 and PBS complementary to tRNA Lys,3 was also replication competent. Upon extended in vitro culture though, this virus reverted to utilize tRNA Lys1,2. Analysis of the infectivity and replication of the wild-type and mutant viruses revealed that tRNA Pro was the preferred tRNA for high-level replication. Viruses with a PBS complementary to tRNA Arg or tRNA Ly1,2 replicated at levels approximately 30% and 10% as effective as the wild-type virus, while virus with a PBS complementary to tRNA Lys,3 had the slowest replication kinetics and least infectivity. Comparison of the virion tRNA content of the wild-type and mutant viruses revealed similar ratios with respect to levels of tRNA Pro, tRNA Arg and tRNA Lys. Modeling of the U5-PBS region revealed that the predicted RNA structure for the virus that selected tRNA Arg was more similar to the wild type virus that uses tRNA Pro than the virus which use tRNA Lys1,2 or tRNA Lys,3; the virus that uses tRNA Lys,3 had the most profound disruption in the predicted RNA structure. The results of these studies demonstrate that MuLV has evolved to preferentially select tRNA Pro for high-level replication and are discussed with respect to common features of the primer selection process between MuLV and other retroviruses.