Genomic organization of eukaryotic tRNAs

A genome resource for the marine annelid Platynereis dumerilii

The marine annelid Platynereis dumerilii is a model organism used in many research areas including evolution and development, neurobiology, ecology and regeneration. Here we present the genomes of P. dumerilii and of the closely related P. massiliensis and P. megalops, to facilitate comparative genomic approaches and help explore Platynereis biology. We used long-read sequencing technology and chromosomal-conformation capture along with extensive transcriptomic resources to obtain and annotate a draft genome assembly of ~1.47 Gbp for P. dumerilii, of which more than half represent repeat elements. We predict around 29,000 protein-coding genes, with relatively large intron sizes, over 38,000 non-coding genes, and 580 miRNA loci. We further explore the high genetic variation (~3% heterozygosity) within the Platynereis species complex. Gene ontology reveals the most variable loci to be associated with pigmentation, development and immunity. The current work sets the stage for further development of Platynereis genomic resources.

Toward a Categorization of Virus-ncRNA Interactions in the World of RNA to Disentangle the Tiny Secrets of Dengue Virus

In recent years, the function of noncoding RNAs (ncRNAs) as regulatory molecules of cell physiology has begun to be better understood. Advances in viral molecular biology have shown that host ncRNAs, cellular factors, and virus-derived ncRNAs and their interplay are strongly disturbed during viral infections. Nevertheless, the folding of RNA virus genomes has also been identified as a critical factor in regulating canonical and non-canonical functions. Due to the influence of host ncRNAs and the structure of RNA viral genomes, complex molecular and cellular processes in infections are modulated. We propose three main categories to organize the current information about RNA-RNA interactions in some well-known human viruses. The first category shows examples of host ncRNAs associated with the immune response triggered in viral infections. Even though miRNAs introduce a standpoint, they are briefly presented to keep researchers moving forward in uncovering other RNAs. The second category outlines interactions between virus-host ncRNAs, while the third describes how the structure of the RNA viral genome serves as a scaffold for processing virus-derived RNAs. Our grouping may provide a comprehensive framework to classify ncRNA-host-cell interactions for emerging viruses and diseases. In this sense, we introduced them to organize DENV-host-cell interactions.

The genome of the colonial hydroid Hydractinia reveals that their stem cells use a toolkit of evolutionarily shared genes with all animals

Hydractinia is a colonial marine hydroid that shows remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, Hydractinia symbiolongicarpus and Hydractinia echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell-type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from nonself.

Systematic computational hunting for small RNAs derived from ncRNAs during dengue virus infection in endothelial HMEC-1 cells

In recent years, a population of small RNA fragments derived from non-coding RNAs (sfd-RNAs) has gained significant interest due to its functional and structural resemblance to miRNAs, adding another level of complexity to our comprehension of small-RNA-mediated gene regulation. Despite this, scientists need more tools to test the differential expression of sfd-RNAs since the current methods to detect miRNAs may not be directly applied to them. The primary reasons are the lack of accurate small RNA and ncRNA annotation, the multi-mapping read (MMR) placement, and the multicopy nature of ncRNAs in the human genome. To solve these issues, a methodology that allows the detection of differentially expressed sfd-RNAs, including canonical miRNAs, by using an integrated copy-number-corrected ncRNA annotation was implemented. This approach was coupled with sixteen different computational strategies composed of combinations of four aligners and four normalization methods to provide a rank-order of prediction for each differentially expressed sfd-RNA. By systematically addressing the three main problems, we could detect differentially expressed miRNAs and sfd-RNAs in dengue virus-infected human dermal microvascular endothelial cells. Although more biological evaluations are required, two molecular targets of the hsa-mir-103a and hsa-mir-494 (CDK5 and PI3/AKT) appear relevant for dengue virus (DENV) infections. Here, we performed a comprehensive annotation and differential expression analysis, which can be applied in other studies addressing the role of small fragment RNA populations derived from ncRNAs in virus infection.

Selective gene expression maintains human tRNA anticodon pools during differentiation

Transfer RNAs are essential for translating genetic information into proteins. The human genome contains hundreds of predicted tRNA genes, many in multiple copies. How their expression is regulated to control tRNA repertoires is unknown. Here we combined quantitative tRNA profiling and chromatin immunoprecipitation with sequencing to measure tRNA expression following the differentiation of human induced pluripotent stem cells into neuronal and cardiac cells. We find that tRNA transcript levels vary substantially, whereas tRNA anticodon pools, which govern decoding rates, are more stable among cell types. Mechanistically, RNA polymerase III transcribes a wide range of tRNA genes in human induced pluripotent stem cells but on differentiation becomes constrained to a subset we define as housekeeping tRNAs. This shift is mediated by decreased mTORC1 signalling, which activates the RNA polymerase III repressor MAF1. Our data explain how tRNA anticodon pools are buffered to maintain decoding speed across cell types and reveal that mTORC1 drives selective tRNA expression during differentiation.

Plasmodium, the Apicomplexa Outlier When It Comes to Protein Synthesis

Plasmodium is an obligate intracellular parasite that has numerous interactions with different hosts during its elaborate life cycle. This is also the case for the other parasites belonging to the same phylum Apicomplexa. In this study, we bioinformatically identified the components of the multi-synthetase complexes (MSCs) of several Apicomplexa parasites and modelled their assembly using AlphaFold2. It appears that none of these MSCs resemble the two MSCs that we have identified and characterized in Plasmodium. Indeed, tRip, the central protein involved in the association of the two Plasmodium MSCs is different from its homologues, suggesting also that the tRip-dependent import of exogenous tRNAs is not conserved in other apicomplexan parasites. Based on this observation, we searched for obvious differences that could explain the singularity of Plasmodium protein synthesis by comparing tRNA genes and amino acid usage in the different genomes. We noted a contradiction between the large number of asparagine residues used in Plasmodium proteomes and the single gene encoding the tRNA that inserts them into proteins. This observation remains true for all the Plasmodia strains studied, even those that do not contain long asparagine homorepeats.

tRNA-Cys gene clusters exhibit high variability in Arabidopsis thaliana

Although most of the genes encoding tRNAs in plants are dispersed throughout the genome, a fraction of them form tRNA gene clusters. In Arabidopsis thaliana, the smallest of tRNA clusters on chromosome 5 consists of four tRNA-Cys-GCA genes placed within repeating units of 0.4 kbp. A systematic analysis of the genomic sequences of syntenic regions from various ecotypes of A. thaliana showed that the general structure of the cluster, consisting of a tRNA-Cys pseudogene followed by repeating units containing tRNA-Cys genes, is well conserved. However, there is significant heterogeneity in the number of repeating units between different ecotypes. A unique feature of this cluster is the presence of putative transposable elements (Helitron). In addition, two further tRNA-Cys gene mini-clusters (gene pairs) in A. thaliana were identified. RNA-seq-based evaluation of expression of tRNA-Cys-GCA genes showed a positive signal for 11 out of 13 unique transcripts. An analysis of the conservation of the tRNA-Cys clusters from A. thaliana with the corresponding regions from four other Arabidopsis species suggests a sequence of events that led to the divergence of these regions.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

The genome of the colonial hydroid Hydractinia reveals their stem cells utilize a toolkit of evolutionarily shared genes with all animals

Hydractinia is a colonial marine hydroid that exhibits remarkable biological properties, including the capacity to regenerate its entire body throughout its lifetime, a process made possible by its adult migratory stem cells, known as i-cells. Here, we provide an in-depth characterization of the genomic structure and gene content of two Hydractinia species, H. symbiolongicarpus and H. echinata, placing them in a comparative evolutionary framework with other cnidarian genomes. We also generated and annotated a single-cell transcriptomic atlas for adult male H. symbiolongicarpus and identified cell type markers for all major cell types, including key i-cell markers. Orthology analyses based on the markers revealed that Hydractinia's i-cells are highly enriched in genes that are widely shared amongst animals, a striking finding given that Hydractinia has a higher proportion of phylum-specific genes than any of the other 41 animals in our orthology analysis. These results indicate that Hydractinia's stem cells and early progenitor cells may use a toolkit shared with all animals, making it a promising model organism for future exploration of stem cell biology and regenerative medicine. The genomic and transcriptomic resources for Hydractinia presented here will enable further studies of their regenerative capacity, colonial morphology, and ability to distinguish self from non-self.

Copy number variation in tRNA isodecoder genes impairs mammalian development and balanced translation

The number of tRNA isodecoders has increased dramatically in mammals, but the specific molecular and physiological reasons for this expansion remain elusive. To address this fundamental question we used CRISPR editing to knockout the seven-membered phenylalanine tRNA gene family in mice, both individually and combinatorially. Using ATAC-Seq, RNA-seq, ribo-profiling and proteomics we observed distinct molecular consequences of single tRNA deletions. We show that tRNA-Phe-1-1 is required for neuronal function and its loss is partially compensated by increased expression of other tRNAs but results in mistranslation. In contrast, the other tRNA-Phe isodecoder genes buffer the loss of each of the remaining six tRNA-Phe genes. In the tRNA-Phe gene family, the expression of at least six tRNA-Phe alleles is required for embryonic viability and tRNA-Phe-1-1 is most important for development and survival. Our results reveal that the multi-copy configuration of tRNA genes is required to buffer translation and viability in mammals.

High production of transfer RNAs identifies the presence of developing oocytes in ovaries and intersex testes of teleost fish

5S rRNA is highly transcribed in fish oocytes and this transcription levels can be used to identify the presence of oocytes in the intersex testes of fish exposed to xenoestrogens. Similar to 5S rRNA, tRNAs are transcribed by RNA polymerase III (Pol-III) in eukaryotes, so this study focuses in the analysis of the levels of expression of tRNAs in the gonads (ovaries and testes) of eight teleost species as a possible new oocyte molecular marker. Total RNA extracted from gonads of six commercial teleost species in the Biscay Bay, from the pollution sentinel species thicklip grey mullet (Chelon labrosus) known present intersex testes in response to xenoestrogens in Gernika estuary and from the laboratory model species Danio rerio were analysed through capillary electrophoresis. Bioanalyzer electropherograms were used to quantify the concentrations of tRNAs, 5S and 5.8S rRNA. All studied ovaries expressed significantly higher levels of tRNAs and 5S rRNA than testes. A tRNA to 5.8S rRNA index was calculated which differentiates ovaries from testes, and identifies some intersex testes in between testes and ovaries in mullets. The tRNA/5.8S ratio was highest in ovaries in previtellogenic stage, decreasing towards maturity. Thus, strong oocyte expression of tRNAs is an additional proof of high activity levels of Pol-III during early stages of oocyte development in teleost ovaries. Incidentally, we observed that miRNA concentrations were always higher in testes than ovaries. The indexing approach developed in the present study could have multiple applications in teleost reproduction research and in the development of early molecular markers of intersex condition.

The Evolution of tRNA Copy Number and Repertoire in Cellular Life

tRNAs are universal decoders that bridge the gap between transcriptome and proteome. They can also be processed into small RNA fragments with regulatory functions. In this work, we show that tRNA copy number is largely controlled by genome size in all cellular organisms, in contrast to what is observed for protein-coding genes that stop expanding between ~20,000 and ~35,000 loci per haploid genome in eukaryotes, regardless of genome size. Our analyses indicate that after the bacteria/archaea split, the tRNA gene pool experienced the evolution of increased anticodon diversity in the archaeal lineage, along with a tRNA gene size increase and mature tRNA size decrease. The evolution and diversification of eukaryotes from archaeal ancestors involved further expansion of the tRNA anticodon repertoire, additional increase in tRNA gene size and decrease in mature tRNA length, along with an explosion of the tRNA gene copy number that emerged coupled with accelerated genome size expansion. Our findings support the notion that macroscopic eukaryotes with a high diversity of cell types, such as land plants and vertebrates, independently evolved a high diversity of tRNA anticodons along with high gene redundancy caused by the expansion of the tRNA copy number. The results presented here suggest that the evolution of tRNA genes played important roles in the early split between bacteria and archaea, and in eukaryogenesis and the later emergence of complex eukaryotes, with potential implications in protein translation and gene regulation through tRNA-derived RNA fragments.

PlantRNA 2.0: an updated database dedicated to tRNAs of photosynthetic eukaryotes

PlantRNA (http://plantrna.ibmp.cnrs.fr/) is a comprehensive database of transfer RNA (tRNA) gene sequences retrieved from fully annotated nuclear, plastidial and mitochondrial genomes of photosynthetic organisms. In the first release (PlantRNA 1.0), tRNA genes from 11 organisms were annotated. In this second version, the annotation was implemented to 51 photosynthetic species covering the whole phylogenetic tree of photosynthetic organisms, from the most basal group of Archeplastida, the glaucophyte Cyanophora paradoxa, to various land plants. tRNA genes from lower photosynthetic organisms such as streptophyte algae or lycophytes as well as extremophile photosynthetic species such as Eutrema parvulum were incorporated in the database. As a whole, about 37 000 tRNA genes were accurately annotated. In the frame of the tRNA genes annotation from the genome of the Rhodophyte Chondrus crispus, non-canonical splicing sites in the D- or T-regions of tRNA molecules were identified and experimentally validated. As for PlantRNA 1.0, comprehensive biological information including 5'- and 3'-flanking sequences, A and B box sequences, region of transcription initiation and poly(T) transcription termination stretches, tRNA intron sequences and tRNA mitochondrial import are included.

A Cocktail of Three Virulent Phages Controls Multidrug-Resistant Salmonella Enteritidis Infection in Poultry

Salmonella enterica is not only the most common pathogen of poultry and poultry-derived products but is also a significant foodborne pathogen. In recent years, many S. enterica isolates have exhibited multi-drug resistance, which places huge pressure on global economy and health. Since phages are an attractive alternative to biocontrol pathogens, we isolated a total of 15 Salmonella phages from sewage effluent, sediment, and chicken manure. The GRNsp1, GRNsp3, GRNsp6, GRNsp21, GRNsp27, GRNsp30, GRNsp50, and GRNsp51 phages exhibited a wide host range against S. enterica serovars Enteritidis and Typhimurium in vitro. In particular, GRNsp51 exerted highly efficient lytic effects against a large proportion of S. Enteritidis and S. Typhimurium strains isolated from different regions of China. Meanwhile, GRNsp8 expanded the host range of GRNsp6 and GRNsp51. Based on their host ranges and lytic capacities, GRNsp6, GRNssp8, and GRNsp51 were selected for further investigation. Morphology, one-step growth curves, and stability assays revealed that GRNsp6, GRNsp8, and GRNsp51 all belong to the Caudovirales order and display relatively short latency periods with broad pH and thermal stability. Genomic analysis indicated that the genomes of these three phages contained no genes related to virulence, antibiotic resistance, or lysogeny. In addition, we tested the effectiveness of a cocktail composed of these three phages against S. Enteritidis in a chicken model. Treatment with the oral phage cocktail 24 h before or alongside Salmonella challenge significantly reduced colonization of the intestinal tract and decreased the mRNA expression of IL-6, IFN-γ, and IL-1β in the duodenum. Together, these findings indicate that a cocktail of the GRNsp6, GRNsp8, and GRNsp51 phages could serve as an effective antimicrobial therapeutic agent against multidrug-resistant Salmonella in animal production to mitigate infections by multiple zoonotic Salmonella species.

Highly contiguous genome assemblies of the Guinea paper wasp (Polistes exclamans) and Mischocyttarus mexicanus

Paper wasps are a model system for the study of social evolution due to a high degree of inter- and intraspecific variation in cooperation, aggression, and visual signals of social status. Increasing the taxonomic coverage of genomic resources for this diverse clade will aid comparative genomic approaches for testing predictions about the molecular basis of social evolution. Here, we provide draft genome assemblies for two well-studied species of paper wasps, Polistes exclamans and Mischocyttarus mexicanus. The P. exclamans genome assembly is 221.5 Mb in length with a scaffold N50 of 4.11 Mb. The M. mexicanus genome assembly is 227 Mb in length with a scaffold N50 of 1.1 Mb. Genomes have low repeat content (9.54–10.75%) and low GC content (32.06–32.4%), typical of other social hymenopteran genomes. The DNA methyltransferase gene, Dnmt3 , was lost early in the evolution of Polistinae. We identified a second independent loss of Dnmt3 within hornets (genus: Vespa).

Neural networks enable efficient and accurate simulation-based inference of evolutionary parameters from adaptation dynamics

The rate of adaptive evolution depends on the rate at which beneficial mutations are introduced into a population and the fitness effects of those mutations. The rate of beneficial mutations and their expected fitness effects is often difficult to empirically quantify. As these 2 parameters determine the pace of evolutionary change in a population, the dynamics of adaptive evolution may enable inference of their values. Copy number variants (CNVs) are a pervasive source of heritable variation that can facilitate rapid adaptive evolution. Previously, we developed a locus-specific fluorescent CNV reporter to quantify CNV dynamics in evolving populations maintained in nutrient-limiting conditions using chemostats. Here, we use CNV adaptation dynamics to estimate the rate at which beneficial CNVs are introduced through de novo mutation and their fitness effects using simulation-based likelihood-free inference approaches. We tested the suitability of 2 evolutionary models: a standard Wright-Fisher model and a chemostat model. We evaluated 2 likelihood-free inference algorithms: the well-established Approximate Bayesian Computation with Sequential Monte Carlo (ABC-SMC) algorithm, and the recently developed Neural Posterior Estimation (NPE) algorithm, which applies an artificial neural network to directly estimate the posterior distribution. By systematically evaluating the suitability of different inference methods and models, we show that NPE has several advantages over ABC-SMC and that a Wright-Fisher evolutionary model suffices in most cases. Using our validated inference framework, we estimate the CNV formation rate at the GAP1 locus in the yeast Saccharomyces cerevisiae to be 10-4.7 to 10-4 CNVs per cell division and a fitness coefficient of 0.04 to 0.1 per generation for GAP1 CNVs in glutamine-limited chemostats. We experimentally validated our inference-based estimates using 2 distinct experimental methods-barcode lineage tracking and pairwise fitness assays-which provide independent confirmation of the accuracy of our approach. Our results are consistent with a beneficial CNV supply rate that is 10-fold greater than the estimated rates of beneficial single-nucleotide mutations, explaining the outsized importance of CNVs in rapid adaptive evolution. More generally, our study demonstrates the utility of novel neural network-based likelihood-free inference methods for inferring the rates and effects of evolutionary processes from empirical data with possible applications ranging from tumor to viral evolution.

Characterization of tRNA expression profiles in large offspring syndrome

Background

Assisted Reproductive Technologies (ART) use can increase the risk of congenital overgrowth syndromes, such as large offspring syndrome (LOS) in ruminants. Epigenetic variations are known to influence gene expression and differentially methylated regions (DMRs) were previously determined to be associated with LOS in cattle. We observed DMRs overlapping tRNA clusters which could affect tRNA abundance and be associated with tissue specificity or overgrowth. Variations in tRNA expression have been identified in several disease pathways suggesting an important role in the regulation of biological processes. Understanding the role of tRNA expression in cattle offers an opportunity to reveal mechanisms of regulation at the translational level. We analyzed tRNA expression in the skeletal muscle and liver tissues of day 105 artificial insemination-conceived, ART-conceived with a normal body weight, and ART-conceived bovine fetuses with a body weight above the 97^th percentile compared to Control-AI.

Results

Despite the centrality of tRNAs to translation, in silico predictions have revealed dramatic differences in the number of tRNA genes between humans and cattle (597 vs 1,659). Consistent with reports in human, only a fraction of predicted tRNA genes are expressed. We detected the expression of 474 and 487 bovine tRNA genes in the muscle and liver with the remainder being unexpressed. 193 and 198 unique tRNA sequences were expressed in all treatment groups within muscle and liver respectively. In addition, an average of 193 tRNA sequences were expressed within the same treatment group in different tissues. Some tRNA isodecoders were differentially expressed between treatment groups. In the skeletal muscle and liver, we categorized 11 tRNA isoacceptors with undetected expression as well as an isodecoder that was unexpressed in the liver (Ser^GGA). Our results identified variation in the proportion of tRNA gene copies expressed between tissues and differences in the highest contributing tRNA anticodon within an amino acid family due to treatment and tissue type. Out of all amino acid families, roughly half of the most highly expressed tRNA isoacceptors correlated to their most frequent codon in the bovine genome.

Conclusion

Although the number of bovine tRNA genes is nearly triple of that of the tRNA genes in human, there is a shared occurrence of transcriptionally inactive tRNA genes in both species. We detected differential expression of tRNA genes as well as tissue- and treatment- specific tRNA transcripts with unique sequence variations that could modulate translation during protein homeostasis or cellular stress, and give rise to regulatory products targeting genes related to overgrowth in the skeletal muscle and/or tumor development in the liver of LOS individuals. While the absence of certain isodecoders may be relieved by wobble base pairing, missing tRNA species could increase the likelihood of mistranslation or mRNA degradation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08496-7.

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson–Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

BOPAL 2.0 and a study of tRNA and rRNA gene evolution in Clostridium

We present BOPAL 2.0, an improved version of the BOPAL algorithm for the evolutionary history inference of tRNA and rRNA genes in bacterial genomes. Our approach can infer complete evolutionary scenarios and ancestral gene orders on a phylogeny and considers a wide range of events such as duplications, deletions, substitutions, inversions and transpositions. It is based on the fact that tRNA and rRNA genes are often organized in operons/clusters in bacteria, and this information is used to help identify orthologous genes for each genome comparison. BOPAL 2.0 introduces new features, such as a triple-wise alignment step, context-aware singleton matching and a second pass of the algorithm. Evaluation on simulated datasets shows that BOPAL 2.0 outperforms the original BOPAL in terms of the accuracy of inferred events and ancestral genomes. We also present a study of the tRNA/rRNA gene evolution in the Clostridium genus, in which the organization of these genes is very divergent. Our results indicate that tRNA and rRNA genes in Clostridium have evolved through numerous duplications, losses, transpositions and substitutions, but very few inversions were inferred.

Role of a cryptic tRNA gene operon in survival under translational stress

As compared to eukaryotes, bacteria have a reduced tRNA gene set encoding between 30 and 220 tRNAs. Although in most bacterial phyla tRNA genes are dispersed in the genome, many species from distinct phyla also show genes forming arrays. Here, we show that two types of arrays with distinct evolutionary origins exist. This work focuses on long tRNA gene arrays (L-arrays) that encompass up to 43 genes, which disseminate by horizontal gene transfer and contribute supernumerary tRNA genes to the host. Although in the few cases previously studied these arrays were reported to be poorly transcribed, here we show that the L-array of the model cyanobacterium Anabaena sp. PCC 7120, encoding 23 functional tRNAs, is largely induced upon impairment of the translation machinery. The cellular response to this challenge involves a global reprogramming of the transcriptome in two phases. tRNAs encoded in the array are induced in the second phase of the response, directly contributing to cell survival. Results presented here show that in some bacteria the tRNA gene set may be partitioned between a housekeeping subset, which constantly sustains translation, and an inducible subset that is generally silent but can provide functionality under particular conditions.

Transfer RNA-derived fragments in aging Caenorhabditis elegans originate from abundant homologous gene copies

Small RNAs that originate from transfer RNA (tRNA) species, tRNA-derived fragments (tRFs), play diverse biological functions but little is known for their association with aging. Moreover, biochemical aspects of tRNAs limit discovery of functional tRFs by high throughput sequencing. In particular, genes encoding tRNAs exist as multiple copies throughout genome, and mature tRNAs have various modified bases, contributing to ambiguities for RNA sequencing-based analysis of tRFs. Here, we report age-dependent changes of tRFs in Caenorhabditis elegans. We first analyzed published RNA sequencing data by using a new strategy for tRNA-associated sequencing reads. Our current method used unique mature tRNAs as a reference for the sequence alignment, and properly filtered out false positive enrichment for tRFs. Our analysis successfully distinguished de novo mutation sites from differences among homologous copies, and identified potential RNA modification sites. Overall, the majority of tRFs were upregulated during aging and originated from 5′-ends, which we validated by using Northern blot analysis. Importantly, we revealed that the major source of tRFs upregulated during aging was the tRNAs with abundant gene copy numbers. Our analysis suggests that tRFs are useful biomarkers of aging particularly when they originate from abundant homologous gene copies.

Molecular features similarities between SARS-CoV-2, SARS, MERS and key human genes could favour the viral infections and trigger collateral effects

In December 2019, rising pneumonia cases caused by a novel β-coronavirus (SARS-CoV-2) occurred in Wuhan, China, which has rapidly spread worldwide, causing thousands of deaths. The WHO declared the SARS-CoV-2 outbreak as a public health emergency of international concern, since then several scientists are dedicated to its study. It has been observed that many human viruses have codon usage biases that match highly expressed proteins in the tissues they infect and depend on the host cell machinery for the replication and co-evolution. In this work, we analysed 91 molecular features and codon usage patterns for 339 viral genes and 463 human genes that consisted of 677,873 codon positions. Hereby, we selected the highly expressed genes from human lung tissue to perform computational studies that permit to compare their molecular features with those of SARS, SARS-CoV-2 and MERS genes. The integrated analysis of all the features revealed that certain viral genes and overexpressed human genes have similar codon usage patterns. The main pattern was the A/T bias that together with other features could propitiate the viral infection, enhanced by a host dependant specialization of the translation machinery of only some of the overexpressed genes. The envelope protein E, the membrane glycoprotein M and ORF7 could be further benefited. This could be the key for a facilitated translation and viral replication conducting to different comorbidities depending on the genetic variability of population due to the host translation machinery. This is the first codon usage approach that reveals which human genes could be potentially deregulated due to the codon usage similarities between the host and the viral genes when the virus is already inside the human cells of the lung tissues. Our work leaded to the identification of additional highly expressed human genes which are not the usual suspects but might play a role in the viral infection and settle the basis for further research in the field of human genetics associated with new viral infections. To identify the genes that could be deregulated under a viral infection is important to predict the collateral effects and determine which individuals would be more susceptible based on their genetic features and comorbidities associated.

Genome Size Reduction and Transposon Activity Impact tRNA Gene Diversity While Ensuring Translational Stability in Birds

As a highly diverse vertebrate class, bird species have adapted to various ecological systems. How this phenotypic diversity can be explained genetically is intensively debated and is likely grounded in differences in the genome content. Larger and more complex genomes could allow for greater genetic regulation that results in more phenotypic variety. Surprisingly, avian genomes are much smaller compared to other vertebrates but contain as many protein-coding genes as other vertebrates. This supports the notion that the phenotypic diversity is largely determined by selection on non-coding gene sequences. Transfer RNAs (tRNAs) represent a group of non-coding genes. However, the characteristics of tRNA genes across bird genomes have remained largely unexplored. Here, we exhaustively investigated the evolution and functional consequences of these crucial translational regulators within bird species and across vertebrates. Our dense sampling of 55 avian genomes representing each bird order revealed an average of 169 tRNA genes with at least 31% being actively used. Unlike other vertebrates, avian tRNA genes are reduced in number and complexity but are still in line with vertebrate wobble pairing strategies and mutation-driven codon usage. Our detailed phylogenetic analyses further uncovered that new tRNA genes can emerge through multiplication by transposable elements. Together, this study provides the first comprehensive avian and cross-vertebrate tRNA gene analyses and demonstrates that tRNA gene evolution is flexible albeit constrained within functional boundaries of general mechanisms in protein translation.

Structural and Genetic Determinants of Convergence in the Drosophila tRNA Structure-Function Map

The evolution of tRNA multigene families remains poorly understood, exhibiting unusual phenomena such as functional conversions of tRNA genes through anticodon shift substitutions. We improved FlyBase tRNA gene annotations from twelve Drosophila species, incorporating previously identified ortholog sets to compare substitution rates across tRNA bodies at single-site and base-pair resolution. All rapidly evolving sites fell within the same metal ion-binding pocket that lies at the interface of the two major stacked helical domains. We applied our tRNA Structure-Function Mapper (tSFM) method independently to each Drosophila species and one outgroup species Musca domestica and found that, although predicted tRNA structure-function maps are generally highly conserved in flies, one tRNA Class-Informative Feature (CIF) within the rapidly evolving ion-binding pocket-Cytosine 17 (C17), ancestrally informative for lysylation identity-independently gained asparaginylation identity and substituted in parallel across tRNAAsn paralogs at least once, possibly multiple times, during evolution of the genus. In D. melanogaster, most tRNALys and tRNAAsn genes are co-arrayed in one large heterologous gene cluster, suggesting that heterologous gene conversion as well as structural similarities of tRNA-binding interfaces in the closely related asparaginyl-tRNA synthetase (AsnRS) and lysyl-tRNA synthetase (LysRS) proteins may have played a role in these changes. A previously identified Asn-to-Lys anticodon shift substitution in D. ananassae may have arisen to compensate for the convergent and parallel gains of C17 in tRNAAsn paralogs in that lineage. Our results underscore the functional and evolutionary relevance of our tRNA structure-function map predictions and illuminate multiple genomic and structural factors contributing to rapid, parallel and compensatory evolution of tRNA multigene families.

Combining tRNA sequencing methods to characterize plant tRNA expression and post-transcriptional modification

Differences in tRNA expression have been implicated in a remarkable number of biological processes. There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of high-throughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana, which exhibited many anciently conserved similarities with humans but also positions that appear to be 'hot spots' for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

Diversity of tRNA Clusters in the Chloroviruses

Viruses rely on their host’s translation machinery for the synthesis of their own proteins. Problems belie viral translation when the host has a codon usage bias (CUB) that is different from an infecting virus due to differences in the GC content between the host and virus genomes. Here, we examine the hypothesis that chloroviruses adapted to host CUB by acquisition and selection of tRNAs that at least partially favor their own CUB. The genomes of 41 chloroviruses comprising three clades, each infecting a different algal host, have been sequenced, assembled and annotated. All 41 viruses not only encode tRNAs, but their tRNA genes are located in clusters. While differences were observed between clades and even within clades, seven tRNA genes were common to all three clades of chloroviruses, including the tRNA^Arg gene, which was found in all 41 chloroviruses. By comparing the codon usage of one chlorovirus algal host, in which the genome has been sequenced and annotated (67% GC content), to that of two of its viruses (40% GC content), we found that the viruses were able to at least partially overcome the host’s CUB by encoding tRNAs that recognize AU-rich codons. Evidence presented herein supports the hypothesis that a chlorovirus tRNA cluster was present in the most recent common ancestor (MRCA) prior to divergence into three clades. In addition, the MRCA encoded a putative isoleucine lysidine synthase (TilS) that remains in 39/41 chloroviruses examined herein, suggesting a strong evolutionary pressure to retain the gene. TilS alters the anticodon of tRNA^Met that normally recognizes AUG to then recognize AUA, a codon for isoleucine. This is advantageous to the chloroviruses because the AUA codon is 12–13 times more common in the chloroviruses than their host, further helping the chloroviruses to overcome CUB. Among large DNA viruses infecting eukaryotes, the presence of tRNA genes and tRNA clusters appear to be most common in the Phycodnaviridae and, to a lesser extent, in the Mimiviridae.

Epigenetic silencing of clustered tRNA genes in Arabidopsis

Beyond their key role in translation, cytosolic transfer RNAs (tRNAs) are involved in a wide range of other biological processes. Nuclear tRNA genes (tDNAs) are transcribed by the RNA polymerase III (RNAP III) and cis-elements, trans-factors as well as genomic features are known to influence their expression. In Arabidopsis, besides a predominant population of dispersed tDNAs spread along the 5 chromosomes, some clustered tDNAs have been identified. Here, we demonstrate that these tDNA clusters are transcriptionally silent and that pathways involved in the maintenance of DNA methylation play a predominant role in their repression. Moreover, we show that clustered tDNAs exhibit repressive chromatin features whilst their dispersed counterparts contain permissive euchromatic marks. This work demonstrates that both genomic and epigenomic contexts are key players in the regulation of tDNAs transcription. The conservation of most of these regulatory processes suggests that this pioneering work in Arabidopsis can provide new insights into the regulation of RNA Pol III transcription in other organisms, including vertebrates.

The multi-faceted regulation of nuclear tRNA gene transcription

Transfer RNAs are among the most ancient molecules of life on earth. Beyond their crucial role in protein synthesis as carriers of amino acids, they are also important players in a plethora of other biological processes. Many debates in term of biogenesis, regulation and function persist around these fascinating non-coding RNAs. Our review focuses on the first step of their biogenesis in eukaryotes, i.e. their transcription from nuclear genes. Numerous and complementary ways have emerged during evolution to regulate transfer RNA gene transcription. Here, we will summarize the different actors implicated in this process: cis-elements, trans-factors, genomic contexts, epigenetic environments and finally three-dimensional organization of nuclear genomes. © 2019 IUBMB Life, 2019 © 2019 IUBMB Life, 71(8):1099-1108, 2019.

Mycobacterium genus and tRNA arrays

The presence of tRNA array, a region with high tRNA gene number and density, has been demonstrated in Mycobacterium genus. However, a recent phylogenomic study revealed the existence of five distinct monophyletic groups (genera) within this genus. Considering this new scenario, and based on in-silico analyses, we have identified and characterised the abundance and diversity of tRNA array units within Mycobacterium, Mycolicibacterium gen. nov., Mycolicibacillus gen. nov., and Mycobacteroides gen. nov. The occurrence and prevalence of tRNA arrays among the genera belonging to Actinobacteria indicate their possible role in the organismal fitness.

The cellular landscape of mid-size noncoding RNA

Noncoding RNA plays an important role in all aspects of the cellular life cycle, from the very basic process of protein synthesis to specialized roles in cell development and differentiation. However, many noncoding RNAs remain uncharacterized and the function of most of them remains unknown. Mid-size noncoding RNAs (mncRNAs), which range in length from 50 to 400 nucleotides, have diverse regulatory functions but share many fundamental characteristics. Most mncRNAs are produced from independent promoters although others are produced from the introns of other genes. Many are found in multiple copies in genomes. mncRNAs are highly structured and carry many posttranscriptional modifications. Both of these facets dictate their RNA-binding protein partners and ultimately their function. mncRNAs have already been implicated in translation, catalysis, as guides for RNA modification, as spliceosome components and regulatory RNA. However, recent studies are adding new mncRNA functions including regulation of gene expression and alternative splicing. In this review, we describe the different classes, characteristics and emerging functions of mncRNAs and their relative expression patterns. Finally, we provide a portrait of the challenges facing their detection and annotation in databases. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Global In-Silico Scenario of tRNA Genes and Their Organization in Virus Genomes

Viruses are known to be highly dependent on the host translation machinery for their protein synthesis. However, tRNA genes are occasionally identified in such organisms, and in addition, few of them harbor tRNA gene clusters comprising dozens of genes. Recently, tRNA gene clusters have been shown to occur among the three domains of life. In such a scenario, the viruses could play a role in the dispersion of such structures among these organisms. Thus, in order to reveal the prevalence of tRNA genes as well as tRNA gene clusters in viruses, we performed an unbiased large-scale genome survey. Interestingly, tRNA genes were predicted in ssDNA (single-stranded DNA) and ssRNA (single-stranded RNA) viruses as well in many other dsDNA viruses of families from Caudovirales order. In the latter group, tRNA gene clusters composed of 15 to 37 tRNA genes were characterized, mainly in bacteriophages, enlarging the occurrence of such structures within viruses. These bacteriophages were from hosts that encompass five phyla and 34 genera. This in-silico study presents the current global scenario of tRNA genes and their organization in virus genomes, contributing and opening questions to be explored in further studies concerning the role of the translation apparatus in these organisms.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation

BACKGROUND

Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes.

FINDINGS

Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated.

CONCLUSION

The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

Exploring tRNA gene cluster in archaea

BACKGROUND

Shared traits between prokaryotes and eukaryotes are helpful in the understanding of the tree of life evolution. In bacteria and eukaryotes, it has been shown a particular organisation of tRNA genes as clusters, but this trait has not been explored in the archaea domain.

OBJECTIVE

Explore the occurrence of tRNA gene clusters in archaea.

METHODS

In-silico analyses of complete and draft archaeal genomes based on tRNA gene isotype and synteny, tRNA gene cluster content and mobilome elements.

FINDINGS

We demonstrated the prevalence of tRNA gene clusters in archaea. tRNA gene clusters, composed of archaeal-type tRNAs, were identified in two Archaea class, Halobacteria and Methanobacteria from Euryarchaeota supergroup. Genomic analyses also revealed evidence of the association between tRNA gene clusters to mobile genetic elements and intra-domain horizontal gene transfer.

MAIN CONCLUSIONS

tRNA gene cluster occurs in the three domains of life, suggesting a role of this type of tRNA gene organisation in the biology of the living organisms.

Accurate mapping of tRNA reads

Motivation

Many repetitive DNA elements are transcribed at appreciable expression levels. Mapping the corresponding RNA sequencing reads back to a reference genome is notoriously difficult and error-prone task, however. This is in particular true if chemical modifications introduce systematic mismatches, while at the same time the genomic loci are only approximately identical, as in the case of tRNAs.

Results

We therefore developed a dedicated mapping strategy to handle RNA-seq reads that map to tRNAs relying on a modified target genome in which known tRNA loci are masked and instead intronless tRNA precursor sequences are appended as artificial 'chromosomes'. In a first pass, reads that overlap the boundaries of mature tRNAs are extracted. In the second pass, the remaining reads are mapped to a tRNA-masked target that is augmented by representative mature tRNA sequences. Using both simulated and real life data we show that our best-practice workflow removes most of the mapping artefacts introduced by simpler mapping schemes and makes it possible to reliably identify many of chemical tRNA modifications in generic small RNA-seq data. Using simulated data the FDR is only 2%. We find compelling evidence for tissue specific differences of tRNA modification patterns.

Availability and implementation

The workflow is available both as a bash script and as a Galaxy workflow from https://github.com/AnneHoffmann/tRNA-read-mapping.

Contact

fabian@tbi.univie.ac.at.

Supplementary information

Supplementary data are available at Bioinformatics online.

Repertoires of tRNAs: The Couplers of Genomics and Proteomics

The pool of transfer RNA (tRNA) molecules in cells allows the ribosome to decode genetic information. This repertoire of molecular decoders is positioned in the crossroad of the genome, the transcriptome, and the proteome. Omics and systems biology now allow scientists to explore the entire repertoire of tRNAs of many organisms, revealing basic exciting biology. The tRNA gene set of hundreds of species is now characterized, in addition to the tRNA genes of organelles and viruses. Genes encoding tRNAs for certain anticodon types appear in dozens of copies in a genome, while others are universally absent from any genome. Transcriptome measurement of tRNAs is challenging, but in recent years new technologies have allowed researchers to determine the dynamic expression patterns of tRNAs. These advances reveal that availability of ready-to-translate tRNA molecules is highly controlled by several transcriptional and posttranscriptional regulatory processes. This regulation shapes the proteome according to the cellular state. The tRNA pool profoundly impacts many aspects of cellular and organismal life, including protein expression level, translation accuracy, adequacy of folding, and even mRNA stability. As a result, the shape of the tRNA pool affects organismal health and may participate in causing conditions such as cancer and neurological conditions.

Association of Transfer RNA Fragments in White Blood Cells With Antibody Response to Bovine Leukemia Virus in Holstein Cattle

Bovine leukemia virus (BLV) affects cattle health and productivity worldwide, causing abnormal immune function and immunosuppression. Transfer RNA fragments (tRFs) are known to be involved in inhibition of gene expression and have been associated with stress and immune response, tumor growth, and viral infection. The objective of this study was to identify tRFs associated with antibody response to BLV in Holstein cattle. Sera from 14 animals were collected to establish IgG reactivity to BLV by ELISA. Seven animals were seropositive (positive group) and seven were seronegative (negative group) for BLV exposure. Leukocytes from each animal were collected and tRFs were extracted for sequencing. tRF5^GlnCTG, tRF5^GlnTTG, and tRF5^HisGTG, were significantly different between seropositive and seronegative groups (P < 0.0067). In all cases the positive group had a lower number of normalized sequences for tRFs when compared to the negative group. Result suggests that tRF5s could potentially be used as biomarkers to establish exposure of cattle to BLV.

Genome-Wide Analysis Reveals Ancestral Lack of Seventeen Different tRNAs and Clade-Specific Loss of tRNA-CNNs in Archaea

Transfer RNA (tRNA) is a category of RNAs that specifically decode messenger RNAs (mRNAs) into proteins by recognizing a set of 61 codons commonly adopted by different life domains. The composition and abundance of tRNAs play critical roles in shaping codon usage and pairing bias, which subsequently modulate mRNA translation efficiency and accuracy. Over the past few decades, effort has been concentrated on evaluating the specificity and redundancy of different tRNA families. However, the mechanism and processes underlying tRNA evolution have only rarely been investigated. In this study, by surveying tRNA genes in 167 completely sequenced genomes, we systematically investigated the composition and evolution of tRNAs in Archaea from a phylogenetic perspective. Our data revealed that archaeal genomes are compact in both tRNA types and copy number. Generally, no more than 44 different types of tRNA are present in archaeal genomes to decode the 61 canonical codons, and most of them have only one gene copy per genome. Among them, tRNA-Met was significantly overrepresented, with an average of three copies per genome. In contrast, the tRNA-UAU and 16 tRNAs with A-starting anticodons (tRNA-ANNs) were rarely detected in all archaeal genomes. The conspicuous absence of these tRNAs across the archaeal phylogeny suggests they might have not been evolved in the common ancestor of Archaea, rather than have lost independently from different clades. Furthermore, widespread absence of tRNA-CNNs in the Methanococcales and Methanobacteriales genomes indicates convergent loss of these tRNAs in the two clades. This clade-specific tRNA loss may be attributing to the reductive evolution of their genomes. Our data suggest that the current tRNA profiles in Archaea are contributed not only by the ancestral tRNA composition, but also by differential maintenance and loss of redundant tRNAs.

Novel Genomic and Evolutionary Perspective of Cyanobacterial tRNAs

Transfer RNA (tRNA) plays a central role in protein synthesis and acts as an adaptor molecule between an mRNA and an amino acid. A tRNA has an L-shaped clover leaf-like structure and contains an acceptor arm, D-arm, D-loop, anti-codon arm, anti-codon loop, variable loop, Ψ-arm and Ψ-loop. All of these arms and loops are important in protein translation. Here, we aimed to delineate the genomic architecture of these arms and loops in cyanobacterial tRNA. Studies from tRNA sequences from 61 cyanobacterial species showed that, except for few tRNAs (tRNA^Asn, tRNA^Leu, tRNA^Gln, and tRNA^Met), all contained a G nucleotide at the 1st position in the acceptor arm. tRNA^Leu and tRNA^Met did not contain any conserved nucleotides at the 1st position whereas tRNA^Asn and tRNA^Gln contained a conserved U¹ nucleotide. In several tRNA families, the variable region also contained conserved nucleotides. Except for tRNA^Met and tRNA^Glu, all other tRNAs contained a conserved A nucleotide at the 1st position in the D-loop. The Ψ-loop contained a conserved U¹-U²-C³-x-A⁵-x-U⁷ sequence, except for tRNA^Gly, tRNA^Ala, tRNA^Val, tRNA^Phe, tRNA^Thr, and tRNA^Gln in which the U⁷ nucleotide was not conserved. However, in tRNA^Asp, the U⁷ nucleotide was substituted with a C⁷ nucleotide. Additionally, tRNA^Arg, tRNA^Gly, and tRNA^Lys of cyanobacteria contained a group I intron within the anti-codon loop region. Maximum composite likelihood study on the transition/transversion of cyanobacterial tRNA revealed that the rate of transition was higher than the rate of transversion. An evolutionary tree was constructed to understand the evolution of cyanobacterial tRNA and analyses revealed that cyanobacterial tRNA may have evolved polyphyletically with high rate of gene loss.

SMORE: Synteny Modulator of Repetitive Elements

Several families of multicopy genes, such as transfer ribonucleic acids (tRNAs) and ribosomal RNAs (rRNAs), are subject to concerted evolution, an effect that keeps sequences of paralogous genes effectively identical. Under these circumstances, it is impossible to distinguish orthologs from paralogs on the basis of sequence similarity alone. Synteny, the preservation of relative genomic locations, however, also remains informative for the disambiguation of evolutionary relationships in this situation. In this contribution, we describe an automatic pipeline for the evolutionary analysis of such cases that use genome-wide alignments as a starting point to assign orthology relationships determined by synteny. The evolution of tRNAs in primates as well as the history of the Y RNA family in vertebrates and nematodes are used to showcase the method. The pipeline is freely available.

An artificial intelligence approach fit for tRNA gene studies in the era of big sequence data

Unsupervised data mining capable of extracting a wide range of knowledge from big data without prior knowledge or particular models is a timely application in the era of big sequence data accumulation in genome research. By handling oligonucleotide compositions as high-dimensional data, we have previously modified the conventional self-organizing map (SOM) for genome informatics and established BLSOM, which can analyze more than ten million sequences simultaneously. Here, we develop BLSOM specialized for tRNA genes (tDNAs) that can cluster (self-organize) more than one million microbial tDNAs according to their cognate amino acid solely depending on tetra- and pentanucleotide compositions. This unsupervised clustering can reveal combinatorial oligonucleotide motifs that are responsible for the amino acid-dependent clustering, as well as other functionally and structurally important consensus motifs, which have been evolutionarily conserved. BLSOM is also useful for identifying tDNAs as phylogenetic markers for special phylotypes. When we constructed BLSOM with 'species-unknown' tDNAs from metagenomic sequences plus 'species-known' microbial tDNAs, a large portion of metagenomic tDNAs self-organized with species-known tDNAs, yielding information on microbial communities in environmental samples. BLSOM can also enhance accuracy in the tDNA database obtained from big sequence data. This unsupervised data mining should become important for studying numerous functionally unclear RNAs obtained from a wide range of organisms.

Functionality of tRNAs encoded in a mobile genetic element from an acidophilic bacterium

The genome of the acidophilic, bioleaching bacterium Acidithiobacillus ferrooxidans, strain ATCC 23270, contains 95 predicted tRNA genes. Thirty-six of these genes (all 20 species) are clustered within an actively excising integrative-conjugative element (ICEAfe1). We speculated that these tRNA genes might have a role in adapting the bacterial tRNA pool to the codon usage of ICEAfe1 genes. To answer this question, we performed theoretical calculations of the global tRNA adaptation index to the entire A. ferrooxidans genome with and without the ICEAfe1 encoded tRNA genes. Based on these calculations, we observed that tRNAs encoded in ICEAfe1 negatively contribute to adapt the tRNA pool to the codon use in A. ferrooxidans. Although some of the tRNAs encoded in ICEAfe1 are functional in aminoacylation or protein synthesis, we found that they are expressed at low levels. These findings, along with the identification of a tRNA-like RNA encoded in the same cluster, led us to speculate that tRNA genes encoded in the mobile genetic element ICEAfe1 might have acquired mutations that would result in either inactivation or the acquisition of new functions.

Analyses of Genomic tRNA Reveal Presence of Novel tRNAs in Oryza sativa

Transfer rRNAs are important molecules responsible for the translation event during protein synthesis. tRNAs are widespread found in unicellular to multi-cellular organisms. Analysis of tRNA gene family members in Oryza sativa revealed the presence of 750 tRNA genes distributed unevenly in different chromosomes. The length of O. sativa tRNAs genes were ranged from 66 to 91 nucleotides encoding 52 isoacceptor in total. tRNA^Ser found in chromosome 8 of O. sativa encoded only 66 nucleotides which is the smallest tRNA of O. sativa and to our knowledge, this is the smallest gene of eukaryotic lineage reported so far. Analyses revealed the presence of several novel/pseudo tRNA genes in O. sativa which are reported for the first time. Multiple sequence alignment of tRNAs revealed the presence of family specific conserved consensus sequences. Functional study of these novel tRNA and family specific conserved consensus sequences will be crucial to decipher their importance in biological events. The rate of transition of O. sativa tRNA was found to be higher than the rate of transversion. Evolutionary study revealed, O. sativa tRNAs were evolved from the lineages of multiple common ancestors. Duplication and loss study of tRNAs genes revealed, majority of the O. sativa tRNA were duplicated and 17 of them were found to be undergone loss during the evolution. Orthology and paralogy study showed, the majority of O. sativa tRNA were paralogous and only a few of tRNA^Ser were found to contain orthologous tRNAs.

The mitochondrial genomes of the acoelomorph worms Paratomella rubra, Isodiametra pulchra and Archaphanostoma ylvae

Acoels are small, ubiquitous - but understudied - marine worms with a very simple body plan. Their internal phylogeny is still not fully resolved, and the position of their proposed phylum Xenacoelomorpha remains debated. Here we describe mitochondrial genome sequences from the acoels Paratomella rubra and Isodiametra pulchra, and the complete mitochondrial genome of the acoel Archaphanostoma ylvae. The P. rubra and A. ylvae sequences are typical for metazoans in size and gene content. The larger I. pulchra  mitochondrial genome contains both ribosomal genes, 21 tRNAs, but only 11 protein-coding genes. We find evidence suggesting a duplicated sequence in the I. pulchra mitochondrial genome. The P. rubra, I. pulchra and A. ylvae mitochondria have a unique genome organisation in comparison to other metazoan mitochondrial genomes. We found a large degree of protein-coding gene and tRNA overlap with little non-coding sequence in the compact P. rubra genome. Conversely, the A. ylvae and I. pulchra genomes have many long non-coding sequences between genes, likely driving genome size expansion in the latter. Phylogenetic trees inferred from mitochondrial genes retrieve Xenacoelomorpha as an early branching taxon in the deuterostomes. Sequence divergence analysis between P. rubra sampled in England and Spain indicates cryptic diversity.

The Echinococcus canadensis (G7) genome: a key knowledge of parasitic platyhelminth human diseases

BACKGROUND

The parasite Echinococcus canadensis (G7) (phylum Platyhelminthes, class Cestoda) is one of the causative agents of echinococcosis. Echinococcosis is a worldwide chronic zoonosis affecting humans as well as domestic and wild mammals, which has been reported as a prioritized neglected disease by the World Health Organisation. No genomic data, comparative genomic analyses or efficient therapeutic and diagnostic tools are available for this severe disease. The information presented in this study will help to understand the peculiar biological characters and to design species-specific control tools.

RESULTS

We sequenced, assembled and annotated the 115-Mb genome of E. canadensis (G7). Comparative genomic analyses using whole genome data of three Echinococcus species not only confirmed the status of E. canadensis (G7) as a separate species but also demonstrated a high nucleotide sequences divergence in relation to E. granulosus (G1). The E. canadensis (G7) genome contains 11,449 genes with a core set of 881 orthologs shared among five cestode species. Comparative genomics revealed that there are more single nucleotide polymorphisms (SNPs) between E. canadensis (G7) and E. granulosus (G1) than between E. canadensis (G7) and E. multilocularis. This result was unexpected since E. canadensis (G7) and E. granulosus (G1) were considered to belong to the species complex E. granulosus sensu lato. We described SNPs in known drug targets and metabolism genes in the E. canadensis (G7) genome. Regarding gene regulation, we analysed three particular features: CpG island distribution along the three Echinococcus genomes, DNA methylation system and small RNA pathway. The results suggest the occurrence of yet unknown gene regulation mechanisms in Echinococcus.

CONCLUSIONS

This is the first work that addresses Echinococcus comparative genomics. The resources presented here will promote the study of mechanisms of parasite development as well as new tools for drug discovery. The availability of a high-quality genome assembly is critical for fully exploring the biology of a pathogenic organism. The E. canadensis (G7) genome presented in this study provides a unique opportunity to address the genetic diversity among the genus Echinococcus and its particular developmental features. At present, there is no unequivocal taxonomic classification of Echinococcus species; however, the genome-wide SNPs analysis performed here revealed the phylogenetic distance among these three Echinococcus species. Additional cestode genomes need to be sequenced to be able to resolve their phylogeny.

Automated detection of ncRNAs in the draft genome sequence of a colonial tunicate: the carpet sea squirt Didemnum vexillum

BACKGROUND

The colonial ascidian Didemnum vexillum, sea carpet squirt, is not only a key marine organism to study morphological ancestral patterns of chordates evolution but it is also of great ecological importance due to its status as a major invasive species. Non-coding RNAs, in particular microRNAs (miRNAs), are important regulatory genes that impact development and environmental adaptation. Beyond miRNAs, not much in known about tunicate ncRNAs.

RESULTS

We provide here a comprehensive homology-based annotation of non-coding RNAs in the recently sequenced genome of D. vexillum. To this end we employed a combination of several computational approaches, including blast searches with a wide range of parameters, and secondary structured centered survey with infernal. The resulting candidate set was curated extensively to produce a high-quality ncRNA annotation of the first draft of the D. vexillum genome. It comprises 57 miRNA families, 4 families of ribosomal RNAs, 22 isoacceptor classes of tRNAs (of which more than 72 % of loci are pseudogenes), 13 snRNAs, 12 snoRNAs, and 1 other RNA family. Additionally, 21 families of mitochondrial tRNAs and 2 of mitochondrial ribosomal RNAs and 1 long non-coding RNA.

CONCLUSIONS

The comprehensive annotation of the D. vexillum non-coding RNAs provides a starting point towards a better understanding of the restructuring of the small RNA system in ascidians. Furthermore it provides a valuable research for efforts to establish detailed non-coding RNA annotations for other recently published and recently sequences in tunicate genomes.

Orthologs, turn-over, and remolding of tRNAs in primates and fruit flies

Background

Transfer RNAs (tRNAs) are ubiquitous in all living organism. They implement the genetic code so that most genomes contain distinct tRNAs for almost all 61 codons. They behave similar to mobile elements and proliferate in genomes spawning both local and non-local copies. Most tRNA families are therefore typically present as multicopy genes. The members of the individual tRNA families evolve under concerted or rapid birth-death evolution, so that paralogous copies maintain almost identical sequences over long evolutionary time-scales. To a good approximation these are functionally equivalent. Individual tRNA copies thus are evolutionary unstable and easily turn into pseudogenes and disappear. This leads to a rapid turnover of tRNAs and often large differences in the tRNA complements of closely related species. Since tRNA paralogs are not distinguished by sequence, common methods cannot not be used to establish orthology between tRNA genes.

Results

In this contribution we introduce a general framework to distinguish orthologs and paralogs in gene families that are subject to concerted evolution. It is based on the use of uniquely aligned adjacent sequence elements as anchors to establish syntenic conservation of sequence intervals. In practice, anchors and intervals can be extracted from genome-wide multiple sequence alignments. Syntenic clusters of concertedly evolving genes of different families can then be subdivided by list alignments, leading to usually small clusters of candidate co-orthologs. On the basis of recent advances in phylogenetic combinatorics, these candidate clusters can be further processed by cograph editing to recover their duplication histories. We developed a workflow that can be conceptualized as stepwise refinement of a graph of homologous genes. We apply this analysis strategy with different types of synteny anchors to investigate the evolution of tRNAs in primates and fruit flies. We identified a large number of tRNA remolding events concentrated at the tips of the phylogeny. With one notable exception all phylogenetically old tRNA remoldings do not change the isoacceptor class.

Conclusions

Gene families evolving under concerted evolution are not amenable to classical phylogenetic analyses since paralogs maintain identical, species-specific sequences, precluding the estimation of correct gene trees from sequence differences. This leaves conservation of syntenic arrangements with respect to “anchor elements” that are not subject to concerted evolution as the only viable source of phylogenetic information. We have demonstrated here that a purely synteny-based analysis of tRNA gene histories is indeed feasible. Although the choice of synteny anchors influences the resolution in particular when tight gene clusters are present, and the quality of sequence alignments, genome assemblies, and genome rearrangements limits the scope of the analysis, largely coherent results can be obtained for tRNAs. In particular, we conclude that a large fraction of the tRNAs are recent copies. This proliferation is compensated by rapid pseudogenization as exemplified by many very recent alloacceptor remoldings.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2927-4) contains supplementary material, which is available to authorized users.

Tissue- and Time-Specific Expression of Otherwise Identical tRNA Genes

Codon usage bias affects protein translation because tRNAs that recognize synonymous codons differ in their abundance. Although the current dogma states that tRNA expression is exclusively regulated by intrinsic control elements (A- and B-box sequences), we revealed, using a reporter that monitors the levels of individual tRNA genes in Caenorhabditis elegans, that eight tryptophan tRNA genes, 100% identical in sequence, are expressed in different tissues and change their expression dynamically. Furthermore, the expression levels of the sup-7 tRNA gene at day 6 were found to predict the animal's lifespan. We discovered that the expression of tRNAs that reside within introns of protein-coding genes is affected by the host gene's promoter. Pairing between specific Pol II genes and the tRNAs that are contained in their introns is most likely adaptive, since a genome-wide analysis revealed that the presence of specific intronic tRNAs within specific orthologous genes is conserved across Caenorhabditis species.

A Comprehensive tRNA Genomic Survey Unravels the Evolutionary History of tRNA Arrays in Prokaryotes

Considering the importance of tRNAs in the translation machinery, scant attention has been paid to tRNA array units defined as genomic regions containing at least 20 tRNA genes with a minimal tRNA gene density of two tRNA genes per kilobase. Our analysis of Acidithiobacillus ferrivorans CF27 and Acidithiobacillus ferrooxidans ATCC 23270(T) genomes showed that both display a tRNA array unit with syntenic conservation which mainly contributed to the tRNA gene redundancy in these two organisms. Our investigations into the occurrence and distribution of tRNA array units revealed that 1) this tRNA organization is limited to few phyla and mainly found in Gram-positive bacteria; and 2) the presence of tRNA arrays favors the redundancy of tRNA genes, in particular those encoding the core tRNA isoacceptors. Finally, comparative array organization revealed that tRNA arrays were acquired through horizontal gene transfer (from Firmicutes or unknown donor), before being subjected to tRNA rearrangements, deletions, and duplications. In Bacilli, the most parsimonious evolutionary history involved two common ancestors and the acquisition of their arrays arose late in evolution, in the genera branches. Functional roles of the array units in organism lifestyle, selective genetic advantage and translation efficiency, as well as the evolutionary advantages of organisms harboring them were proposed. Our study offers new insight into the structural organization and evolution of tRNA arrays in prokaryotic organisms.

Eukaryotic tRNA paradox

tRNAs are widely believed to segregate into two classes, I and II. Computational analysis of eukaryotic tRNA entries in Genomic tRNA Database, however, leads to new, albeit paradoxical, presence of more than a thousand class-I tRNAs with uncharacteristic long variable arms (V-arms), like in class-II. Out of 62,202 tRNAs from 69 eukaryotes, as many as 1431 class-I tRNAs have these novel extended V-arms, and we refer to them as paradoxical tRNAs (pxtRNAs). A great majority of these 1431 pxtRNA genes are located in intergenic regions, about 18% embedded in introns of genes or ESTs, and just one in 3'UTR. A check on the conservations of 2D and 3D base pairs for each position of these pxtRNAs reveals a few variations, but they seem to have almost all the known features (already known identity and conserved elements of tRNA). Analyses of the A-Box and B-Box of these pxtRNA genes in eukaryotes display salient deviations from the previously annotated conserved features of the standard promoters, whereas the transcription termination signals are just canonical and non-canonical runs of thymidine, similar to the ones in standard tRNA genes. There is just one such pxtRNA(ProAGG) gene in the entire human genome, and the availability of data allows epigenetic analysis of this human pxtRNA(ProAGG) in three different cell lines, H1 hESC, K562, and NHEK, to assess the level of its expression. Histone acetylation and methylation of this lone pxtRNA(ProAGG) gene in human differ from that of the nine standard human tRNA(ProAGG) genes. The V-arm nucleotide sequences and their secondary structures in pxtRNA differ from that of class-II tRNA. Considering these differences, hypotheses of alternative splicing, non-canonical intron and gene transfer are examined to partially improve the Cove scores of these pxtRNAs and to critically question their antecedence and novelty.

GENCODE pseudogenes

Historically pseudogenes were believed to represent nonfunctional genomic fossils; however, there is emerging evidence that many of them could be biologically active. This possibility has ignited interest in pseudogene loci and made the need for their high-quality annotation more pressing as an accurate knowledge of all pseudogenes in the human reference genome sequence facilitates confident functional analysis. GENCODE have undertaken the first genome-wide pseudogene assignment for protein-coding genes combining both large-scale manual annotation and computational pseudogene prediction pipelines. Multiple computational predictions provide an unbiased set of hints for manual annotators to investigate, both during first-pass annotation and as part of QC to identify any potential missing pseudogene loci. Where a pseudogene is identified, the extent of its homology to the parent locus is fully investigated by a manual annotator; a pseudogene model is built and assigned to one of eight pseudogene biotypes depending on the mechanism of creation and on the presence of locus-specific transcriptional or proteomic data. The high-quality, information-rich set of pseudogenes created has been integrated with ENCODE functional genomics data, specifically expression level, transcription factor and RNA polymerase II binding, and chromatin marks. In this way we have been able to identify some pseudogenes that possess conventional characteristics of functionality as well as others with interesting patterns of partial activity, which might suggest that putatively inactive loci could be gaining a novel function, for example as long noncoding RNAs. The activity data associated with every pseudogene is stored in the psiDR resource.