Comprehensive analysis of insertion sequences within rRNA genes of CPR bacteria and biochemical characterization of a homing endonuclease encoded by these sequences.

The Candidate Phyla Radiation (CPR) represents an extensive bacterial clade comprising primarily uncultured lineages and is distinguished from other bacteria by a significant prevalence of insertion sequences (ISs) within their rRNA genes. However, our understanding of the taxonomic distribution and characteristics of these ISs remains limited. In this study, we used a comprehensive approach to systematically determine the nature of the rRNA ISs in CPR bacteria. The analysis of hundreds of rRNA gene sequences across 65 CPR phyla revealed that ISs are present in 48% of 16S rRNA genes and 82% of 23S rRNA genes, indicating a broad distribution across the CPR clade, with exceptions in the 16S and 23S rRNA genes of Candidatus (Ca.) Saccharibacteria and the 16S rRNA genes of Ca. Peregrinibacteria. Over half the ISs display a group-I-intron-like structure, whereas specific 16S rRNA gene ISs display features reminiscent of group II introns. The ISs frequently encode proteins with homing endonuclease (HE) domains, centered around the LAGLIDADG motif. The LAGLIDADG HE (LHE) proteins encoded by the rRNA ISs of CPR bacteria predominantly have a single-domain structure, deviating from the usual single- or double-domain configuration observed in typical prokaryotic LHEs. Experimental analysis of one LHE protein, I-ShaI from Ca. Shapirobacteria, confirmed that its endonuclease activity targets the DNA sequence of its insertion site, and chemical cross-linking experiments demonstrated its capacity to form homodimers. These results provide robust evidence supporting the hypothesis that the explosive proliferation of rRNA ISs in CPR bacteria was facilitated by mechanisms involving LHEs. IMPORTANCE: Insertion sequences (ISs) in rRNA genes are relatively limited and infrequent in most bacterial phyla. With a comprehensive bioinformatic analysis, we show that in CPR bacteria, these ISs occur in 48% of 16S rRNA genes and 82% of 23S rRNA genes. We also report the systematic and biochemical characterization of the LAGLIDADG homing endonucleases (LHEs) encoded by these ISs in the first such analysis of the CPR bacteria. This study significantly extends our understanding of the phylogenetic positions of rRNA ISs within CPR bacteria and the biochemical features of their LHEs.

Microbiome analysis of the restricted bacteria in radioactive element-containing water at the Fukushima Daiichi Nuclear Power Station.

A major incident occurred at the Fukushima Daiichi Nuclear Power Station following the tsunami triggered by the Tohoku-Pacific Ocean Earthquake in March 2011, whereby seawater entered the torus room in the basement of the reactor building. Here, we identify and analyze the bacterial communities in the torus room water and several environmental samples. Samples of the torus room water (1 × 109 Bq137Cs/L) were collected by the Tokyo Electric Power Company Holdings from two sampling points between 30 cm and 1 m from the bottom of the room (TW1) and the bottom layer (TW2). A structural analysis of the bacterial communities based on 16S rRNA amplicon sequencing revealed that the predominant bacterial genera in TW1 and TW2 were similar. TW1 primarily contained the genus Limnobacter, a thiosulfate-oxidizing bacterium. γ-Irradiation tests on Limnobacter thiooxidans, the most closely related phylogenetically found in TW1, indicated that its radiation resistance was similar to ordinary bacteria. TW2 predominantly contained the genus Brevirhabdus, a manganese-oxidizing bacterium. Although bacterial diversity in the torus room water was lower than seawater near Fukushima, ~70% of identified genera were associated with metal corrosion. Latent environment allocation-an analytical technique that estimates habitat distributions and co-detection analyses-revealed that the microbial communities in the torus room water originated from a distinct blend of natural marine microbial and artificial bacterial communities typical of biofilms, sludge, and wastewater. Understanding the specific bacteria linked to metal corrosion in damaged plants is important for advancing decommissioning efforts. IMPORTANCE: In the context of nuclear power station decommissioning, the proliferation of microorganisms within the reactor and piping systems constitutes a formidable challenge. Therefore, the identification of microbial communities in such environments is of paramount importance. In the aftermath of the Fukushima Daiichi Nuclear Power Station accident, microbial community analysis was conducted on environmental samples collected mainly outside the site. However, analyses using samples from on-site areas, including adjacent soil and seawater, were not performed. This study represents the first comprehensive analysis of microbial communities, utilizing meta 16S amplicon sequencing, with a focus on environmental samples collected from the radioactive element-containing water in the torus room, including the surrounding environments. Some of the identified microbial genera are shared with those previously identified in spent nuclear fuel pools in countries such as France and Brazil. Moreover, our discussion in this paper elucidates the correlation of many of these bacteria with metal corrosion.

Homo-trimeric structure of the ribonuclease for rRNA processing, FAU-1, from Pyrococcus furiosus.

Crystal structure of a ribonuclease for ribosomal RNA processing, FAU-1, from Pyrococcus furiosus was determined with the resolution of 2.57 Å in a homo-trimeric form. The monomer structure consists of two domains: N-terminal and C-terminal domains. C-terminal domain forms trimer and each N-terminal domain locates outside of the trimer core. In the obtained crystal, a dinucleotide, pApUp, was bound to the N-terminal domain, indicating that N-terminal domain has the RNA-binding ability. The affinities to RNA of FAU-1 and a fragment corresponding to the N-terminal domain, FAU-ΔC, were confirmed by polyacrylamide gel electrophoresis and nuclear magnetic resonance (NMR). Interestingly, well-dispersed NMR signals were observed at 318K, indicating that the FAU-ΔC-F18 complex form an ordered structure at higher temperature. As predicted in our previous works, FAU-1 and ribonuclease (RNase) E show a structural similarity in their RNA-binding regions. However, structural similarity between RNase E and FAU-1 could be found in the limited regions of the N-terminal domain. On the other hand, structural similarity between C-terminal domain and some proteins including a phosphatase was found. Thus, it is possible that the catalytic site is located in C-terminal domain.

Systematic Analysis of Diverse Polynucleotide Kinase Clp1 Family Proteins in Eukaryotes: Three Unique Clp1 Proteins of Trypanosoma brucei.

The Clp1 family proteins, consisting of the Clp1 and Nol9/Grc3 groups, have polynucleotide kinase (PNK) activity at the 5' end of RNA strands and are important enzymes in the processing of some precursor RNAs. However, it remains unclear how this enzyme family diversified in the eukaryotes. We performed a large-scale molecular evolutionary analysis of the full-length genomes of 358 eukaryotic species to classify the diverse Clp1 family proteins. The average number of Clp1 family proteins in eukaryotes was 2.3 ± 1.0, and most representative species had both Clp1 and Nol9/Grc3 proteins, suggesting that the Clp1 and Nol9/Grc3 groups were already formed in the eukaryotic ancestor by gene duplication. We also detected an average of 4.1 ± 0.4 Clp1 family proteins in members of the protist phylum Euglenozoa. For example, in Trypanosoma brucei, there are three genes of the Clp1 group and one gene of the Nol9/Grc3 group. In the Clp1 group proteins encoded by these three genes, the C-terminal domains have been replaced by unique characteristics domains, so we designated these proteins Tb-Clp1-t1, Tb-Clp1-t2, and Tb-Clp1-t3. Experimental validation showed that only Tb-Clp1-t2 has PNK activity against RNA strands. As in this example, N-terminal and C-terminal domain replacement also contributed to the diversification of the Clp1 family proteins in other eukaryotic species. Our analysis also revealed that the Clp1 family proteins in humans and plants diversified through isoforms created by alternative splicing.

Identification of Attenuators of Transcriptional Termination: Implications for RNA Regulation in Escherichia coli.

The regulatory function of many bacterial small RNAs (sRNAs) requires the binding of the RNA chaperone Hfq to the 3' portion of the sRNA intrinsic terminator, and therefore sRNA signaling might be regulated by modulating its terminator. Here, using a multicopy screen developed with the terminator of sRNA SgrS, we identified an sRNA gene (cyaR) and three protein-coding genes (cspD, ygjH, and rof) that attenuate SgrS termination in Escherichia coli. Analyses of CyaR and YgjH, a putative tRNA binding protein, suggested that the CyaR activity was indirect and the effect of YgjH was moderate. Overproduction of the protein attenuators CspD and Rof resulted in more frequent readthrough at terminators of SgrS and two other sRNAs, and regulation by SgrS of target mRNAs was reduced. The effect of Rof, a known inhibitor of Rho, was mimicked by bicyclomycin or by a rho mutant, suggesting an unexpected role for Rho in sRNA termination. CspD, a member of the cold shock protein family, bound both terminated and readthrough transcripts, stabilizing them and attenuating termination. By RNA sequencing analysis of the CspD overexpression strain, we found global effects of CspD on gene expression across some termination sites. We further demonstrated effects of endogenous CspD under slow growth conditions where cspD is highly expressed. These findings provided evidence of changes in the efficiency of intrinsic termination, confirming this as an additional layer of the regulation of sRNA signaling. IMPORTANCE Growing evidence suggests that the modulation of intrinsic termination and readthrough of transcription is more widespread than previously appreciated. For small RNAs, proper termination plays a critical role in their regulatory function. Here, we present a multicopy screen approach to identify factors that attenuate small RNA termination and therefore abrogate signaling dependent on the small RNA. This study highlights a new aspect of regulation of small RNA signaling as well as the modulation of intrinsic termination.

Features of smaller ribosomes in candidate phyla radiation (CPR) bacteria revealed with a molecular evolutionary analysis.

The candidate phyla radiation (CPR) is a large bacterial group consisting mainly of uncultured lineages. They have small cells and small genomes, and they often lack ribosomal proteins uL1, bL9, and/or uL30, which are basically ubiquitous in non-CPR bacteria. Here, we comprehensively analyzed the genomic information on CPR bacteria and identified their unique properties. The distribution of protein lengths in CPR bacteria peaks at around 100-150 amino acids, whereas the position of the peak varies in the range of 100-300 amino acids in free-living non-CPR bacteria, and at around 100-200 amino acids in most symbiotic non-CPR bacteria. These results show that the proteins of CPR bacteria are smaller, on average, than those of free-living non-CPR bacteria, like those of symbiotic non-CPR bacteria. We found that ribosomal proteins bL28, uL29, bL32, and bL33 have been lost in CPR bacteria in a taxonomic lineage-specific manner. Moreover, the sequences of approximately half of all ribosomal proteins of CPR differ, in part, from those of non-CPR bacteria, with missing regions or specifically added regions. We also found that several regions in the 16S, 23S, and 5S rRNAs of CPR bacteria are lacking, which presumably caused the total predicted lengths of the three rRNAs of CPR bacteria to be smaller than those of non-CPR bacteria. The regions missing in the CPR ribosomal proteins and rRNAs are located near the surface of the ribosome, and some are close to one another. These observations suggest that ribosomes are smaller in CPR bacteria than those in free-living non-CPR bacteria, with simplified surface structures.

Distinct Expansion of Group II Introns During Evolution of Prokaryotes and Possible Factors Involved in Its Regulation.

Group II introns (G2Is) are ribozymes that have retroelement characteristics in prokaryotes. Although G2Is are suggested to have been an important evolutionary factor in the prokaryote-to-eukaryote transition, comprehensive analyses of these introns among the tens of thousands of prokaryotic genomes currently available are still limited. Here, we developed a bioinformatic pipeline that systematically collects G2Is and applied it to prokaryotic genomes. We found that in bacteria, 25% (447 of 1,790) of the total representative genomes had an average of 5.3 G2Is, and in archaea, 9% (28 of 296) of the total representative genomes had an average of 3.0 G2Is. The greatest number of G2Is per genome was 101 in Arthrospira platensis (phylum Cyanobacteriota). A comprehensive sequence analysis of the intron-encoded protein (IEP) in each G2I sequence was conducted and resulted in the addition of three new IEP classes (U1-U3) to the previous classification. This analysis suggested that about 30% of all IEPs are non-canonical IEPs. The number of G2Is per genome was defined almost at the phylum level, and at least in the following two phyla, Firmicutes, and Cyanobacteriota, the type of IEP was largely associated as a factor in the G2I increase, i.e., there was an explosive increase in G2Is with bacterial C-type IEPs, mainly in the phylum Firmicutes, and in G2Is with CL-type IEPs, mainly in the phylum Cyanobacteriota. We also systematically analyzed the relationship between genomic signatures and the mechanism of these increases in G2Is. This is the first study to systematically characterize G2Is in the prokaryotic phylogenies.

Proteomic and metabolomic analyses uncover sex-specific regulatory pathways in mouse fetal germline differentiation†.

Regulatory mechanisms of germline differentiation have generally been explained via the function of signaling pathways, transcription factors, and epigenetic regulation; however, little is known regarding proteomic and metabolomic regulation and their contribution to germ cell development. Here, we conducted integrated proteomic and metabolomic analyses of fetal germ cells in mice on embryonic day (E)13.5 and E18.5 and demonstrate sex- and developmental stage-dependent changes in these processes. In male germ cells, RNA processing, translation, oxidative phosphorylation, and nucleotide synthesis are dominant in E13.5 and then decline until E18.5, which corresponds to the prolonged cell division and more enhanced hyper-transcription/translation in male primordial germ cells and their subsequent repression. Tricarboxylic acid cycle and one-carbon pathway are consistently upregulated in fetal male germ cells, suggesting their involvement in epigenetic changes preceding in males. Increased protein stability and oxidative phosphorylation during female germ cell differentiation suggests an upregulation of aerobic energy metabolism, which likely contributes to the proteostasis required for oocyte maturation in subsequent stages. The features elucidated in this study shed light on the unrevealed mechanisms of germ cell development.

Behavioral and brain- transcriptomic synchronization between the two opponents of a fighting pair of the fish Betta splendens.

Conspecific male animals fight for resources such as food and mating opportunities but typically stop fighting after assessing their relative fighting abilities to avoid serious injuries. Physiologically, how the fighting behavior is controlled remains unknown. Using the fighting fish Betta splendens, we studied behavioral and brain-transcriptomic changes during the fight between the two opponents. At the behavioral level, surface-breathing, and biting/striking occurred only during intervals between mouth-locking. Eventually, the behaviors of the two opponents became synchronized, with each pair showing a unique behavioral pattern. At the physiological level, we examined the expression patterns of 23,306 brain transcripts using RNA-sequencing data from brains of fighting pairs after a 20-min (D20) and a 60-min (D60) fight. The two opponents in each D60 fighting pair showed a strong gene expression correlation, whereas those in D20 fighting pairs showed a weak correlation. Moreover, each fighting pair in the D60 group showed pair-specific gene expression patterns in a grade of membership analysis (GoM) and were grouped as a pair in the heatmap clustering. The observed pair-specific individualization in brain-transcriptomic synchronization (PIBS) suggested that this synchronization provides a physiological basis for the behavioral synchronization. An analysis using the synchronized genes in fighting pairs of the D60 group found genes enriched for ion transport, synaptic function, and learning and memory. Brain-transcriptomic synchronization could be a general phenomenon and may provide a new cornerstone with which to investigate coordinating and sustaining social interactions between two interacting partners of vertebrates.

Large-Scale Molecular Evolutionary Analysis Uncovers a Variety of Polynucleotide Kinase Clp1 Family Proteins in the Three Domains of Life.

Clp1, a polyribonucleotide 5'-hydroxyl kinase in eukaryotes, is involved in pretRNA splicing and mRNA 3'-end formation. Enzymes similar in amino acid sequence to Clp1, Nol9, and Grc3, are present in some eukaryotes and are involved in prerRNA processing. However, our knowledge of how these Clp1 family proteins evolved and diversified is limited. We conducted a large-scale molecular evolutionary analysis of the Clp1 family proteins in all living organisms for which protein sequences are available in public databases. The phylogenetic distribution and frequencies of the Clp1 family proteins were investigated in complete genomes of Bacteria, Archaea and Eukarya. In total, 3,557 Clp1 family proteins were detected in the three domains of life, Bacteria, Archaea, and Eukarya. Many were from Archaea and Eukarya, but a few were found in restricted, phylogenetically diverse bacterial species. The domain structures of the Clp1 family proteins also differed among the three domains of life. Although the proteins were, on average, 555 amino acids long (range, 196-2,728), 122 large proteins with >1,000 amino acids were detected in eukaryotes. These novel proteins contain the conserved Clp1 polynucleotide kinase domain and various other functional domains. Of these proteins, >80% were from Fungi or Protostomia. The polyribonucleotide kinase activity of Thermus scotoductus Clp1 (Ts-Clp1) was characterized experimentally. Ts-Clp1 preferentially phosphorylates single-stranded RNA oligonucleotides (Km value for ATP, 2.5 µM), or single-stranded DNA at higher enzyme concentrations. We propose a comprehensive assessment of the diversification of the Clp1 family proteins and the molecular evolution of their functional domains.

Evolutionary Analysis of HIV-1 Pol Proteins Reveals Representative Residues for Viral Subtype Differentiation.

RNA viruses have been used as model systems to understand the patterns and processes of molecular evolution because they have high mutation rates and are genetically diverse. Human immunodeficiency virus 1 (HIV-1), the etiological agent of acquired immune deficiency syndrome, is highly genetically diverse, and is classified into several groups and subtypes. However, it has been difficult to use its diverse sequences to establish the overall phylogenetic relationships of different strains or the trends in sequence conservation with the construction of phylogenetic trees. Our aims were to systematically characterize HIV-1 subtype evolution and to identify the regions responsible for HIV-1 subtype differentiation at the amino acid level in the Pol protein, which is often used to classify the HIV-1 subtypes. In this study, we systematically characterized the mutation sites in 2,052 Pol proteins from HIV-1 group M (144 subtype A; 1,528 subtype B; 380 subtype C), using sequence similarity networks. We also used spectral clustering to group the sequences based on the network graph structures. A stepwise analysis of the cluster hierarchies allowed us to estimate a possible evolutionary pathway for the Pol proteins. The subtype A sequences also clustered according to when and where the viruses were isolated, whereas both the subtype B and C sequences remained as single clusters. Because the Pol protein has several functional domains, we identified the regions that are discriminative by comparing the structures of the domain-based networks. Our results suggest that sequence changes in the RNase H domain and the reverse transcriptase (RT) connection domain are responsible for the subtype classification. By analyzing the different amino acid compositions at each site in both domain sequences, we found that a few specific amino acid residues (i.e., M357 in the RT connection domain and Q480, Y483, and L491 in the RNase H domain) represent the differences among the subtypes. These residues were located on the surface of the RT structure and in the vicinity of the amino acid sites responsible for RT enzymatic activity or function.

An archaeal RNA binding protein, FAU-1, is a novel ribonuclease related to rRNA stability in Pyrococcus and Thermococcus.

Ribosome biogenesis and turnover are processes necessary for cell viability and proliferation, and many kinds of proteins are known to regulate these processes. However, many questions still remain, especially in the Archaea. Generally, several ribonucleases are required to process precursor rRNAs to their mature forms, and to degrade rRNAs for quality control. Here, we found that FAU-1, which is known to be an RNA binding protein, possesses an RNase activity against precursor 5S rRNA derived from P. furiosus and T. kodakarensis in the order Thermococcales in vitro. An in vitro analysis revealed that UA sequences in the upstream of 5S rRNA were preferentially degraded by addition of FAU-1. Moreover, a fau-1 gene deletion mutant of T. kodakarensis showed a delay of exponential phase, reduction of maximum cell number and significant changes in the nucleotide sequence lengths of its 5S, 16S, and 23S rRNAs in early exponential phase. Our results suggest that FAU-1 is a potential RNase involved in rRNA stability through maturation and/or degradation processes.

Distinct requirements for energy metabolism in mouse primordial germ cells and their reprogramming to embryonic germ cells.

Primordial germ cells (PGCs), undifferentiated embryonic germ cells, are the only cells that have the ability to become gametes and to reacquire totipotency upon fertilization. It is generally understood that the development of PGCs proceeds through the expression of germ cell-specific transcription factors and characteristic epigenomic changes. However, little is known about the properties of PGCs at the metabolite and protein levels, which are directly responsible for the control of cell function. Here, we report the distinct energy metabolism of PGCs compared with that of embryonic stem cells. Specifically, we observed remarkably enhanced oxidative phosphorylation (OXPHOS) and decreased glycolysis in embryonic day 13.5 (E13.5) PGCs, a pattern that was gradually established during PGC differentiation. We also demonstrate that glycolysis and OXPHOS are important for the control of PGC reprogramming and specification of pluripotent stem cells (PSCs) into PGCs in culture. Our findings about the unique metabolic property of PGCs provide insights into our understanding of the importance of distinct facets of energy metabolism for switching PGC and PSC status.

Systematic Analysis of the Binding Surfaces between tRNAs and Their Respective Aminoacyl tRNA Synthetase Based on Structural and Evolutionary Data.

To determine the mechanism underlying the flow of genetic information, it is important to understand the relationship between a tRNA and its binding enzyme, a member of the aminoacyl-tRNA synthetase (aaRS) family. We have developed a novel method to project the interacting regions of tRNA-aaRS complexes, obtained from their three-dimensional structures, onto two-dimensional space. The interacting surface between each tRNA and its aaRS was successfully identified by determining these interactions with an atomic distance threshold of 3.3 Å. We analyzed their interactions, using 60 mainly bacterial and eukaryotic tRNA-aaRS complexes, and showed that the tRNA sequence regions that interacted most strongly with each aaRS are the anticodon loop and the CCA terminal region, followed by the D-stem. A sequence conservation analysis of the canonical tRNAs was conducted in 83 bacterial, 182 archaeal, and 150 eukaryotic species. Our results show that the three tRNA regions that interact with the aaRS and two additional loop regions (D-loop and TΨC-loop) known to be important for formation of the tRNA L-shaped structure are broadly conserved. We also found sequence conservations near the tRNA discriminator in the Bacteria and Archaea, and an enormous number of noncanonical tRNAs in the Eukaryotes. This is the first global view of tRNA evolution based on its structure and an unprecedented number of sequence data.

Systematic characterization of artificial small RNA-mediated inhibition of Escherichia coli growth.

A new screening system for artificial small RNAs (sRNAs) that inhibit the growth of Escherichia coli was constructed. In this system, we used a plasmid library to express RNAs of ∼120 nucleotides, each with a random 30-nucleotide sequence that can recognize its target mRNA(s). After approximately 60,000 independent colonies were screened, several plasmids that inhibited bacterial growth were isolated. To understand the inhibitory mechanism, we focused on one sRNA, S-20, that exerted a strong inhibitory effect. A time-course analysis of the proteome of S-20-expressing E. coli and a bioinformatic analysis were used to identify potential S-20 target mRNAs, and suggested that S-20 binds the translation initiation sites of several mRNAs encoding enzymes such as peroxiredoxin (osmC), glycyl-tRNA synthetase α subunit (glyQ), uncharacterized protein ygiM, and tryptophan synthase ß chain (trpB). An in vitro translation analysis of chimeric luciferase-encoding mRNAs, each containing a potential S-20 target sequence, indicated that the translation of these mRNAs was inhibited in the presence of S-20. A gel shift analysis combined with the analysis of a series of S-20 mutants suggested that S-20 targets multiple mRNAs that are responsible for inhibiting E. coli growth. These data also suggest that S-20 acts like an endogenous sRNA and that E. coli can utilize artificial sRNAs.

Factors Associated with Changes in Brain Atrophy during a Three-Year Observation in Elderly Diabetic Patients: Effect of Renal Impairment on Hippocampal Atrophy.

BACKGROUND/AIMS: We conducted a 3-year longitudinal study concerning factors associated with changes in brain atrophy in elderly diabetic patients. METHODS: We evaluated hippocampal and global brain atrophy using automatic voxel-based morphometry of structural magnetic resonance images, 4 cognitive function tests, and cerebral small vessel disease (SVD) in 66 diabetic patients. RESULTS: During the 3-year follow-up, hippocampal and global brain atrophy advanced, and cognitive functions worsened. For changes in hippocampal atrophy, changes in estimated glomerular filtration rate (eGFR), albuminuria, and being an ApoE ε4 carrier were independent factors; change in the number of silent brain infarctions was an independent factor for changes in global brain atrophy. A significant association of changes in eGFR and albuminuria with hippocampal atrophy remained after adjusting for confounders including SVD. Both types of brain atrophy at baseline were significantly correlated with cognitive impairment at baseline and especially associated with changes in delayed word recall during the follow-up after adjusting for confounders. CONCLUSION: Changes in eGFR and albuminuria during follow-up were independent risk factors for hippocampal atrophy, which was associated with decline in delayed word recall, suggesting that management of chronic kidney disease may prevent the progression of hippocampal atrophy.

Involvement of nucleotide diphosphate kinase 2 in the reopening of the sensitive period of filial imprinting of domestic chicks (Gallus gallus domesticus).

Filial imprinting is a behavior characterized by the sensitive or critical period restricted to the first few days after hatching. Once the sensitive period is closed, it is widely believed that chicks can never be imprinted under natural conditions. Previously, we showed that the exogenous injection of T3 reopened the sensitive period which was already closed. That study suggested that T3 functioned by way of a rapid non-genomic action; however, the molecular mechanism of how T3 reopens the sensitive period remains unknown. Here, we show that the phosphorylation level of nucleotide diphosphate kinase 2 (NDPK2) was upregulated following T3 injection. Pharmacological deprivation of the kinase activity of NDPK hampered the molecular process prerequisite for the reopening of the sensitive period of filial imprinting. Moreover, it is shown that the kinase activity of NDPK2 participates in the priming process by T3 signaling which endows the potential for learning. Our data indicate that NDPK2 plays a crucial role downstream of T3 action and that its phosphorylation is involved in the non-genomic signaling during imprinting.

Expansion of Noncanonical V-Arm-Containing tRNAs in Eukaryotes.

Transfer RNA (tRNA) is essential for the translation of genetic information into proteins, and understanding its molecular evolution is important if we are to understand the genetic code. In general, long variable-arm (V-arm) structures form in tRNA(Leu), tRNA(Ser), and bacterial and organellar tRNA(Tyr). However, as we have previously reported, noncanonical V-arms occur in nematode tRNA(Gly) and tRNA(Ile), and potentially affect translational fidelity. Here, we comprehensively analyzed 69 eukaryotic genome sequences and examined the evolutionary divergence of the V-arm-containing tRNAs. In total, 253 V-arm-containing tRNAs, with neither leucine nor serine anticodons, were identified in organisms ranging from nematodes to fungi, plants, and vertebrates. We defined them as "noncanonical V-arm-containing tRNAs" (nov-tRNAs). Moreover, 2,415 nov-tRNA-like sequences lacking some of the conserved features of tRNAs were also identified, largely in vertebrate genomes. These nov-tRNA/nov-tRNA-like sequences can be categorized into three types, based on differences in their possible evolutionary origins. The type A nov-tRNAs in nematodes probably evolved not only from tRNA(Leu) but also from tRNA(Ser) and other isotypes on several independent occasions. The type B nov-tRNAs are dispersed abundantly throughout vertebrate genomes, and seem to have originated from retrotransposable elements. The type C nov-tRNAs may have been acquired from plant chloroplasts or from bacteria through horizontal transfer. Our findings provide unexpected insight into the evolution of the tRNA molecule, which was diverse and occurred independently in nematodes, vertebrates, and plants.

Precise mapping and dynamics of tRNA-derived fragments (tRFs) in the development of Triops cancriformis (tadpole shrimp).

BACKGROUND: In a deep sequencing analysis of small RNAs prepared from a living fossil, the tadpole shrimp Triops cancriformis, a 32-nt small RNA was specifically detected in the adult stage. A nucleotide sequence comparison between the 32-nt small RNA and predicted tRNA sequences in the draft nuclear genomic DNA showed that the small RNA was derived from tRNA(Gly)(GCC). To determine the overall features of the tRNA-derived fragments (tRFs) of T. cancriformis, the small RNA sequences in each of the six developmental stages (egg, 1st-4th instar larvae, and adult) were compared with the mitochondrial and nuclear tRNA sequences. RESULTS: We found that the tRFs were derived from mitochondrial and nuclear tRNAs corresponding to 16 and 39 anticodons, respectively. The total read number of nuclear tRFs was approximately 400 times larger than the number of mitochondrial tRFs. Interestingly, the main regions in each parental tRNA from which these tRFs were derived differed, depending on the parental anticodon. Mitochondrial tRF(Ser)(GCU)s were abundantly produced from the 5' half regions of the parental tRNA, whereas mitochondrial tRF(Val)(UAC)s were mainly produced from the 3' end regions. Highly abundant nuclear tRFs, tRF(Gly)(GCC)s, tRF(Gly)(CCC)s, tRF(Glu)(CUC)s, and tRF(Lys)(CUU)s were derived from the 5' half regions of the parental tRNAs. Further analysis of the tRF read counts in the individual developmental stages suggested that the expression of mitochondrial and nuclear tRFs differed during the six stages. Based on these data, we precisely summarized the positions of the tRFs in their parental tRNAs and their expression changes during development. CONCLUSIONS: Our results reveal the entire dynamics of the tRFs from both the nuclear and mitochondrial genomes of T. cancriformis and indicate that the majority of tRFs in the cell are derived from nuclear tRNAs. This study provides the first examples of developmentally expressed mitochondrial tRFs.