Identification, expression, and molecular evolution of microRNAs in the "living fossil" Triops cancriformis (tadpole shrimp)

Genome of tropical bed bug Cimex hemipterus (Cimicidae, Hemiptera) reveals tetraspanin expanded in bed bug ancestor

Cimex species are ectoparasites that exclusively feed on warm-blooded animals such as birds and mammals. Three cimicid species are known to be persistent pests for humans, including the tropical bed bug Cimex hemipterus, common bed bug Cimex lectularius, and Eastern bat bug Leptocimex boueti. To date, genomic information is restricted to the common bed bug C. lectularius, which limits understanding their biology and to provide controls of bed bug infestations. Here, a chromosomal-level genome assembly of C. hemipterus (495 Mb [megabase pairs]) contained on 16 pseudochromosomes (scaffold N50 = 34 Mb), together with 9 messenger RNA and small RNA transcriptomes were obtained. In comparison between hemipteran genomes, we found that the tetraspanin superfamily was expanded in the Cimex ancestor. This study provides the first genome assembly for the tropical bed bug C. hemipterus, and offers an unprecedented opportunity to address questions relating to bed bug infestations, as well as genomic evolution to hemipterans more widely.

Identification of novel microRNAs in the embryonic mouse brain using deep sequencing

Many advances in small RNA-seq technology and bioinformatics pipelines have been made recently, permitting the discovery of novel miRNAs in the embryonic day 15.5 (E15.5) mouse brain. We aimed to improve miRNA discovery in this tissue to expand our knowledge of the regulatory networks that underpin normal neurodevelopment, find new candidates for neurodevelopmental disorder aetiology, and deepen our understanding of non-coding RNA evolution. A high-quality small RNA-seq dataset of 458 M reads was generated. An unbiased miRNA discovery pipeline identified fifty putative novel miRNAs, six of which were selected for further validation. A combination of conservation analysis and target functional prediction was used to determine the authenticity of novel miRNA candidates. These findings demonstrate that miRNAs remain to be discovered, particularly if they have the features of other small RNA species.

Chromosome-scale genome assembly of the freshwater cladoceran crustacean Chydorus sphaericus: A resource for discovery of genes responsive to ecological challenges

The genus Chydorus Leach 1816 (family Chydoridae) is a diverse and ecologically important taxon within freshwater ecosystems. Despite having been widely used in ecological, evolutionary and eco-toxicological studies, no high-quality genomic resource is available for any member of the genus. Here, we present a high-quality chromosome-level assembly of the C. sphaericus genome by combining 7.40 Gb (∼ 50 × coverage) PacBio reads, 19.28 Gb (∼ 135 × coverage) Illumina paired-end reads, and 34.04 Gb Hi-C reads. Our genome assembly is approximately 151 Mb, with contig and scaffold N50 lengths of 1.09 Mb and 13.70 Mb, respectively. The assembly captured 94.9% of the complete eukaryotic BUSCO. Repetitive elements accounted for 17.6% of the genome, and 13,549 protein-coding genes were predicted (based on transcriptome sequencing data, ab-initio or homology-based prediction), of which 96.4% have been functionally annotated in the NCBI-NR database. We identified 303 gene families specific to C. sphaericus, mainly families enriched in functions related to immune response, visual senses and detoxification. Interestingly, we also found 53 significantly expanded gene families in C. sphaericus, mostly with functions related to detoxification. This high-quality assembly genome will act as a reference genome for C. sphaericus and benefit studies on functional and comparative genomics of Chydorus and other crustaceans.

Integrated microRNA-mRNA analysis reveals a possible molecular mechanism of enteritis susceptibility in Litopenaeus vannamei

Enteritis is one of the main diseases affecting Pacific whiteleg shrimp (Litopenaeus vannamei) in recent years, and it has resulted in huge losses to the aquaculture industry. Prior to this study, the molecular mechanism underlying enteritis in L. vannamei was unclear, and comprehensive multi-omics analysis had not been conducted. In this study, 1209 differentially expressed genes (DEGs) were identified from the hepatopancreas of L. vannamei with and without enteritis. Kyoto Encyclopedia of Genes and Genomes analysis showed that genes were significantly enriched in immune, metabolic, and endocrine regulatory pathways. Forty-eight significantly different microRNAs (miRNAs) were identified in the miRNA-Seq analysis. Further functional annotation analysis showed that the regulatory pathway of target gene enrichment of differentially expressed miRNAs was consistent with DEGs. Through miRNA-mRNA integration analysis, 47 meaningful miRNA-mRNA pairs were obtained, of which melanogenesis and pancreatic secretion were considered key pathways. Subsequent miRNA-mRNA interaction network analysis revealed that mja-miR-6493-3p, Mja-miR-6494, novel-198, novel-272, novel-261, novel-200, novel-183, novel-184, novel-237, and novel-192 may be key miRNAs involved in the regulation of these two signaling pathways. Finally, the RAS signaling pathway was found to inhibit the translation level of proteins in the hepatopancreas. These results suggest that target gene integration analysis of mRNA-miRNA can reveal the molecular mechanism underlying enteritis in L. vannamei and also provide valuable new insights for resisting enteritis.

Comparative genomics of Hox and ParaHox genes among major lineages of Branchiopoda with emphasis on tadpole shrimps

Hox and ParaHox genes (HPHGs) are key developmental genes that pattern regional identity along the anterior–posterior body axis of most animals. Here, we identified HPHGs in tadpole shrimps (Pancrustacea, Branchiopoda, Notostraca), an iconic example of the so-called “living fossils” and performed a comparative genomics analysis of HPHGs and the Hox cluster among major branchiopod lineages. Notostraca possess the entire Hox complement, and the Hox cluster seems to be split into two different subclusters, although we were not able to support this finding with chromosome-level assemblies. However, the genomic structure of Hox genes in Notostraca appears more derived than that of Daphnia spp., which instead retains the plesiomorphic condition of a single compact cluster. Spinicaudata and Artemia franciscana show instead a Hox cluster subdivided across two or more genomic scaffolds with some orthologs either duplicated or missing. Yet, branchiopod HPHGs are similar among the various clades in terms of both intron length and number, as well as in their pattern of molecular evolution. Sequence substitution rates are in fact generally similar for most of the branchiopod Hox genes and the few differences we found cannot be traced back to natural selection, as they are not associated with any signals of diversifying selection or substantial switches in selective modes. Altogether, these findings do not support a significant stasis in the Notostraca Hox cluster and further confirm how morphological evolution is not tightly associated with genome dynamics.

Daphnia magna egg piRNA cluster expression profiles change as mothers age

BACKGROUND

PiRNAs prevent transposable elements wreaking havoc on the germline genome. Changes in piRNA expression over the lifetime of an individual may impact on ageing through continued suppression, or release, of transposable element expression. We identified piRNA producing clusters in the genome of Daphnia magna by a combination of bioinformatic methods, and then contrasted their expression between parthenogenetically produced eggs representing maternally-deposited germline piRNAs of young (having their 1st clutch) and old (having their 5th clutch) mothers. Results from eggs were compared to cluster expression in three generations of adults.

RESULTS

As for other arthropods, D. magna encodes long uni-directionally transcribed non-coding RNAs consisting of fragmented transposable elements which account for most piRNAs expressed. Egg tissues showed extensive differences between clutches from young mothers and those from old mothers, with 578 and 686 piRNA clusters upregulated, respectively. Most log fold-change differences for significant clusters were modest, however. When considering only highly expressed clusters, there was a bias towards 1st clutch eggs at 41 upregulated versus eight clusters in the eggs from older mothers. F0 generation differences between young and old mothers were fewer than eggs, as 179 clusters were up-regulated in young versus 170 old mothers. This dropped to 31 versus 22 piRNA clusters when comparing adults in the F1 generation, and no differences were detected in the F3 generation. Inter-generational losses of differential piRNA cluster were similar to that observed for D. magna micro-RNA expression.

CONCLUSIONS

Little overlap in differentially expressed clusters was found between adults containing mixed somatic and germline (ovary) tissues and germ-line representing eggs. A cluster encompassing a Tudor domain containing gene important in the piRNA pathway was upregulated in the eggs from old mothers. We hypothesise that regulation of this gene could form part of a feedback loop that reduces piRNA pathway activity explaining the reduced number of highly-expressed clusters in eggs from old mothers.

Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis

This analysis presents five genome assemblies of four Notostraca taxa. Notostraca origin dates to the Permian/Upper Devonian and the extant forms show a striking morphological similarity to fossil taxa. The comparison of sequenced genomes with other Branchiopoda genomes shows that, despite the morphological stasis, Notostraca share a dynamic genome evolution with high turnover for gene families' expansion/contraction and a transposable elements content comparable to other branchiopods. While Notostraca substitutions rate appears similar or lower in comparison to other branchiopods, a subset of genes shows a faster evolutionary pace, highlighting the difficulty of generalizing about genomic stasis versus dynamism. Moreover, we found that the variation of Triops cancriformis transposable elements content appeared linked to reproductive strategies, in line with theoretical expectations. Overall, besides providing new genomic resources for the study of these organisms, which appear relevant for their ecology and evolution, we also confirmed the decoupling of morphological and molecular evolution.

Immunohistochemical Localization of a GnRH-Like Peptide in the Nerve Ganglion of Three Classes of Crustaceans, the Tadpole Shrimp Triops longicaudatus (Branchiopoda), the Barnacle Balanus crenatus (Hexanauplia), and the Hermit Crab Pagurus filholi (Malacostraca)

In vertebrates, gonadotropin-releasing hormone (GnRH) regulates gonadal maturation by stimulating the synthesis and release of pituitary gonadotropins. GnRH has also been identified in invertebrates. Crustacea consists of several classes including Cephalocarida, Remipedia, Branchiopoda (e.g., tadpole shrimp), Hexanauplia (e.g., barnacle) and Malacostraca (e.g., shrimp, crab). In the malacostracan crustaceans, the presence of GnRH has been detected in several species, mainly by immunohistochemistry. In the present study, we examined whether a GnRH-like peptide exists in the brain and/or nerve ganglion of three classes of crustaceans, the tadpole shrimp Triops longicaudatus (Branchiopoda), the barnacle Balanus crenatus (Hexanauplia), and the hermit crab Pagurus filholi (Malacostraca), by immunohistochemistry using a rabbit polyclonal antibody raised against chicken GnRH-II (GnRH2). This antibody was found to recognize the giant freshwater prawn Macrobrachium rosenbergii GnRH (MroGnRH). In the tadpole shrimp, GnRH-like-immunoreactive (ir) cell bodies were located in the circumesophageal connective of the deuterocerebrum, and GnRH-like-ir fibers were detected also in the ventral nerve cord. In the barnacle, GnRH-like-ir cell bodies and fibers were located in the supraesophageal ganglion (brain), the subesophageal ganglion, and the circumesophageal connective. In the hermit crab, GnRH-like-ir cell bodies were detected in the anterior-most part of the supraesophageal ganglion and the subesophageal ganglion. GnRH-like-ir fibers were observed also in the thoracic ganglion and the eyestalk. These results suggest that a GnRH-like peptide exists widely in crustacean species.

MicroRNAs in Daphnia magna identified and characterized by deep sequencing, genome mapping and manual curation

MicroRNAs (miRNAs) are small non-coding RNAs that function in RNA silencing and post-transcriptional regulation of gene expression in most organisms. The water flea, Daphnia magna is a key model to study phenotypic, physiological and genomic responses to environmental cues and miRNAs can potentially mediate these responses. By using deep sequencing, genome mapping and manual curations, we have characterised the miRNAome of D. magna. We identified 66 conserved miRNAs and 13 novel miRNAs; all of these were found in the three studied life stages of D. magna (juveniles, subadults, adults), but with variation in expression levels between stages. Forty-one of the miRNAs were clustered into 13 genome clusters also present in the D. pulex genome. Most miRNAs contained sequence variants (isomiRs). The highest expressed isomiRs were 3′ template variants with one nucleotide deletion or 3′ non-template variants with addition of A or U at the 3′ end. We also identified offset RNAs (moRs) and loop RNAs (loRs). Our work extends the base for further work on all species (miRNA, isomiRs, moRNAs, loRNAs) of the miRNAome of Daphnia as biomarkers in response to chemical substances and environment cues, and underline age dependency.

Evolutionary emergence of Hairless as a novel component of the Notch signaling pathway

Suppressor of Hairless [Su(H)], the transcription factor at the end of the Notch pathway in Drosophila, utilizes the Hairless protein to recruit two co-repressors, Groucho (Gro) and C-terminal Binding Protein (CtBP), indirectly. Hairless is present only in the Pancrustacea, raising the question of how Su(H) in other protostomes gains repressive function. We show that Su(H) from a wide array of arthropods, molluscs, and annelids includes motifs that directly bind Gro and CtBP; thus, direct co-repressor recruitment is ancestral in the protostomes. How did Hairless come to replace this ancestral paradigm? Our discovery of a protein (S-CAP) in Myriapods and Chelicerates that contains a motif similar to the Su(H)-binding domain in Hairless has revealed a likely evolutionary connection between Hairless and Metastasis-associated (MTA) protein, a component of the NuRD complex. Sequence comparison and widely conserved microsynteny suggest that S-CAP and Hairless arose from a tandem duplication of an ancestral MTA gene.

Integrated analysis of differentially expressed microRNAs and mRNAs to screen miRNAs and genes related to reproduction in Macrobrachium nipponense

Female Macrobrachium nipponense has the characteristic of short sexual maturity during the breeding season, which can increase breeding risk and lead to prevalent female individual miniaturization. In this study, we characterized micro (mi)RNA-seq data of the eyestalk (E) and cerebral ganglia (B) of female M. nipponense during breeding and non-breeding seasons. A total of 393 and 189 differentially expressed miRNAs (DE miRNAs) were identified in BSE vs. NBSE and BSB vs. NBSB, respectively. The most abundant up- and down-regulated DE miRNAs were miR-124, miR-14, and miR-7. Enrichment analysis showed that DE miRNA target genes were mainly involved in 'metabolic process' and 'binding', and were associated with 'neurohormonal regulation' and 'photoreceptor activity' signaling pathways. Integrated analysis of miRNA-mRNA expression showed that the most abundant DE miRNAs were miR-14 and miR-278 in BSE vs. NBSE and BSB vs. NBSB, respectively. Four pairs of DE miRNAs and their corresponding target annotated genes were selected from the DE miRNA-mRNA interaction network (bmo-miR-316-5p/opsin protein, ame-miR-125/skeletal muscle actin 8, dmo-miR-278/sugar transporter, and tca-miR-3885-5p/5-HT1 receptor). Gene expression analysis of these four pairs in different ovary development stages showed their potential regulatory roles in ovary maturation.

Characterization of the salmon louse Lepeophtheirus salmonis miRNome: Sex-biased differences related to the coding and non-coding RNA interplay

The salmon louse Lepeophtheirus salmonis is a marine ectoparasite that has a detrimental impact on salmon farms. Genomic knowledge of adult stages is critical to understand the reproductive success and lifecycle completion of this species. Here, we report a comprehensive characterization of the L. salmonis miRNome with emphasis on the sex-differences of the parasite. Small-RNA sequencing was conducted on males and females, and mRNA-sequencing was also conducted to identify miRNA-targets at these stages. Based on bioinformatics analyses, 3101 putative miRNAs were found in L. salmonis, including precursors and variants. The most abundant and over-expressed miRNAs belonged to the bantam, mir-100, mir-1, mir-263a and mir-276 families, while the most differentially expressed mRNAs corresponded to genes related to reproduction and other biological processes involved in cell-differentiation. Target analyses revealed that the most up-regulated miRNAs in males can act by inhibiting the expression of genes related to female differentiation such as vitellogenin genes. Target prediction and expression patterns suggested a pivotal role of miRNAs in the reproductive development of L. salmonis.

microRNAs in Macrobrachium olfersii embryos: Identification, their biogenesis components and potential targets

In embryonic development, microRNAs (miRNAs) regulate the complex gene expression associated with the complexity of embryogenesis. Today, few studies have been conducted on the identification of miRNAs and components of miRNA biogenesis on embryonic development in crustaceans, especially in prawns. In this context, the aim of this study was to identify in silico components of miRNA biogenesis, and miRNAs and potential target genes during embryonic development in the prawn Macrobrachium olfersii through small RNAs and transcriptome analyses. Using the miRDeep2 program, we identified 17 miRNA precursors in M. olfersii, which seven (miR-9, miR-10, miR-92, miR-125, miR-305, miR-1175, and miR-2788) were reported in the miRBase database, indicating high evolutionary conservation of these sequences among animals. The other 10 miRNAs of M. olfersii were novel miRNAs and only similar to Macrobrachium niponnense miRNAs, indicating genus-specific miRNAs. In addition, eight key components of miRNA biogenesis (DROSHA, PASHA/DGCR8, XPO5, RAN, DICER, TRBP2, AGO, and PIWI) were identified in M. olfersii embryos unigenes. In the annotation of miRNA targets, 516 genes were similar to known sequences in the GenBank database. Regarding the conserved miRNAs, we verified that they were differentially expressed during embryonic development in M. olfersii. In conclusion, this is the first study that identifies conserved and novel miRNAs in the prawn M. olfersii with some miRNA target genes involved in embryonic development. Our results will allow further studies on the function of these miRNAs and miRNA biogenesis components during embryonic development in M. olfersii and other prawns of commercial interest.

Novel microRNAs are associated with population divergence in transcriptional response to thermal stress in an intertidal copepod

The role of gene expression in adaptation to differing thermal environments has been assayed extensively. Yet, in most natural systems, analyses of gene expression reveal only one level of the complexity of regulatory machineries. MicroRNAs (miRNAs) are small noncoding RNAs which are key components of many gene regulatory networks, and they play important roles in a variety of cellular pathways by modulating post-transcriptional quantities of mRNA available for protein synthesis. The characterization of miRNA loci and their regulatory dynamics in nonmodel systems are still largely understudied. In this study, we examine the role of miRNAs in response to high thermal stress in the intertidal copepod Tigriopus californicus. Allopatric populations of this species show varying levels of local adaptation with respect to thermal regimes, and previous studies showed divergence in gene expression between populations from very different thermal environments. We examined the transcriptional response to temperature stress in two populations separated by only 8 km by utilizing RNA-seq to quantify both mRNA and miRNA levels. Using the currently available genome sequence, we first describe the repertoire of miRNAs in T. californicus and assess the degree to which transcriptional response to temperature stress is governed by miRNA activity. The two populations showed large differences in the number of genes involved, the magnitude of change in commonly used genes and in the number of miRNAs involved in transcriptional modulation during stress. Our results suggest that an increased level of regulatory network complexity may underlie improved survivorship under thermal stress in one of the populations.

The identification of microRNAs involved in the response of Chinese shrimp Fenneropenaeus chinensis to white spot syndrome virus infection

MicroRNA (miRNA) is a class of small noncoding RNA, which is involved in the post-transcriptional regulation in all metazoan eukaryotes. MiRNAs might play an important role in the host response to virus infection. However, miRNAs in the aquatic crustacean species were not extensively investigated. To obtain a better understanding of the response of Chinese shrimp Fenneropenaeus chinensis to white spot syndrome virus (WSSV) infection, the sequence and expression profile of miRNAs in the hepatopancreas of WSSV-infected F. chinensis were obtained by the high-throughput Illumina HiSeq 2500 deep sequencing technique. A total number of 129 known miRNAs and 44 putative novel miRNAs were identified from the deep sequencing data. The peak size of miRNAs was 22 nt (37.0%). 25 miRNAs were significantly (P < 0.05) differentially expressed post WSSV infection. Six of the differentially expressed miRNAs were randomly selected for further verification by the real-time RT-PCR technique. The results showed that there was a consistency between the deep sequencing and real-time RT-PCR assay. The target genes of differentially expressed miRNAs were predicted. Each miRNA had 4 target genes on average. The results suggested that some specific miRNAs might be involved in the response of F. chinensis to WSSV infection, and further provided basic information for the investigation of specific miRNAs in F. chinensis.

Composition and genomic organization of arthropod Hox clusters

Background

The ancestral arthropod is believed to have had a clustered arrangement of ten Hox genes. Within arthropods, Hox gene mutations result in transformation of segment identities. Despite the fact that variation in segment number/character was common in the diversification of arthropods, few examples of Hox gene gains/losses have been correlated with morphological evolution. Furthermore, a full appreciation of the variation in the genomic arrangement of Hox genes in extant arthropods has not been recognized, as genome sequences from each major arthropod clade have not been reported until recently. Initial genomic analysis of the chelicerate Tetranychusurticae suggested that loss of Hox genes and Hox gene clustering might be more common than previously assumed. To further characterize the genomic evolution of arthropod Hox genes, we compared the genomic arrangement and general characteristics of Hox genes from representative taxa from each arthropod subphylum.

Results

In agreement with others, we find arthropods generally contain ten Hox genes arranged in a common orientation in the genome, with an increasing number of sampled species missing either Hox3 or abdominal-A orthologs. The genomic clustering of Hox genes in species we surveyed varies significantly, ranging from 0.3 to 13.6 Mb. In all species sampled, arthropod Hox genes are dispersed in the genome relative to the vertebrate Mus musculus. Differences in Hox cluster size arise from variation in the number of intervening genes, intergenic spacing, and the size of introns and UTRs. In the arthropods surveyed, Hox gene duplications are rare and four microRNAs are, in general, conserved in similar genomic positions relative to the Hox genes.

Conclusions

The tightly clustered Hox complexes found in the vertebrates are not evident within arthropods, and differential patterns of Hox gene dispersion are found throughout the arthropods. The comparative genomic data continue to support an ancestral arthropod Hox cluster of ten genes with a shared orientation, with four Hox gene-associated miRNAs, although the degree of dispersion between genes in an ancestral cluster remains uncertain. Hox3 and abdominal-A orthologs have been lost in multiple, independent lineages, and current data support a model in which inversions of the Abdominal-B locus that result in the loss of abdominal-A correlate with reduced trunk segmentation.

Electronic supplementary material

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

The identification and characteristics of salinity-related microRNAs in gills of Portunus trituberculatus

MicroRNAs (miRNAs) are small noncoding RNAs that post-transcriptionally regulate gene expression in organisms. To understand the underlying mechanisms behind the molecular response of the crab to low salt-stress, high-throughput Illumina/Solexa deep sequencing technology was used to investigate the expression profiles of miRNAs under low salinity challenged. Two mixed RNA pool libraries of gill tissues from low salinity challenged (LC) and the control groups (NC) were sequenced on the Illumina platform. A total of 6,166,057 and 7,032,973 high-quality reads were obtained from the NC and LC libraries, respectively. Sixty-seven miRNAs consisting of 16 known and 51 novel ones were identified, among which, 12 miRNAs were differentially expressed in LC compared to NC. Thirty-four of the target genes predicted were differentially expressed in the opposite direction to the miRNAs, which were involved in crucial processes related to osmoregulation by gene ontology (GO) functional enrichment analysis, such as anion transport processes (GO:0006820) and chitin metabolic process (GO:0006030). These results provide a basis for further investigation of the miRNA-modulating networks in osmoregulation of Portunus trituberculatus.

Roles of small RNAs in the immune defense mechanisms of crustaceans

Small RNAs, 21-24 nucleotides in length, are non-coding RNAs found in most multicellular organisms, as well as in some viruses. There are three main types of small RNAs including microRNA (miRNA), small-interfering RNA (siRNA), and piwi-interacting RNA (piRNA). Small RNAs play key roles in the genetic regulation of eukaryotes; at least 50% of all eukaryote genes are the targets of small RNAs. In recent years, studies have shown that some unique small RNAs are involved in the immune response of crustaceans, leading to lower or higher immune responses to infections and diseases. SiRNAs could be used as therapy for virus infection. In this review, we provide an overview of the diverse roles of small RNAs in the immune defense mechanisms of crustaceans.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0245-5) contains supplementary material, which is available to authorized users.