Differential expression patterns of conserved miRNAs and isomiRs during Atlantic halibut development

Benchmarking of bioinformatics tools for NGS-based microRNA profiling with RT-qPCR method

MicroRNAs are vital gene expression regulators, extensively studied worldwide. The large-scale characterization of miRNAomes is possible using next-generation sequencing (NGS). This technology offers great opportunities, but these cannot be fully exploited without proper and comprehensive bioinformatics analysis. This may be achieved by the use of reliable dedicated software; however, different programs may generate divergent results, leading to additional discrepancies. Thus, the aim of this study was to compare three bioinformatic algorithms dedicated to NGS-based microRNA profiling and validate them using an alternative method, namely RT-qPCR. The comparison analysis revealed differences in the number and sets of identified miRNAs. The qPCR confirmed the expression of the investigated microRNAs. The correlation analysis of NGS and qPCR measurements showed strong and significant coefficients for a subset of the tested miRNAs, including those detected by all three algorithms. Single miRNA variants (isomiRs) showed different levels of correlation with the qPCR data. The obtained results revealed the good performance of all tested programs, despite the observed differences. Moreover, they implied that some specific miRNAs may be differentially estimated using NGS technology and the qPCR method, regardless of the used bioinformatics software. These discrepancies may stem from many factors, including the composition of the isomiR profile, their abundance, length, and investigated species. In conclusion, in this study, we shed light on the bioinformatics aspects of miRNAome profiling, elucidating its complexity and pinpointing potential features influencing validation. Thus, qPCR validation results should be open to interpretation when not fully concordant with NGS results until further, additional analyses are conducted.

Integrative miRNA-mRNA network analysis to identify crucial pathways of salinity adaptation in brain transcriptome of Labeo rohita

Introduction: Brain being the master regulator of the physiology of animal, the current study focuses on the gene expression pattern of the brain tissue with special emphasis on regulation of growth, developmental process of an organism and cellular adaptation of Labeo rohita against unfavourable environmental conditions. Methods: RNA-seq study was performed on collected brain samples at 8ppt salt concentration and analyzed for differential gene expression, functional annotation and miRNA-mRNA regulatory network. Results: We found that 2450 genes were having significant differential up and down regulation. The study identified 20 hub genes based on maximal clique centrality algorithm. These hub genes were mainly involved in various signaling pathways, energy metabolism and ion transportation. Further, 326 up and 1214 down regulated genes were found to be targeted by 7 differentially expressed miRNAs i.e., oni-miR-10712, oni-miR-10736, ssa-miR-221-3p, ssa-miR-130d-1-5p, ssa-miR-144-5p and oni-miR-10628. Gene ontology analysis of these differentially expressed genes led to the finding that these genes were involved in signal transduction i.e., calcium, FOXO, PI3K-AKT, TGF-β, Wnt and p53 signalling pathways. Differentially expressed genes were also involved in regulation of immune response, environmental adaptation i.e., neuroactive ligand-receptor interaction, ECM-receptor interaction, cell adhesion molecules and circadian entrainment, osmoregulation and energy metabolism, which are critical for salinity adaptation. Discussion: The findings of whole transcriptomic study on brain deciphered the miRNA-mRNA interaction patterns and pathways associated with salinity adaptation of L. rohita.

cbx2 is a functional target of the let-7 family in the gonad of Japanese flounder (Paralichthys olivaceus)

As a key member of the miRNA family, the role and target gene of the let-7 family in the gonad of Japanese flounder (Paralichthys olivaceus) is unclear. Chromobox homolog 2 (CBX2) is one of the core components of the polycomb group complex (PcG) and significantly influences gonadal development. The deletion of CBX2 can lead to sex reversal in mammals. Therefore, exploring the relationship between the let-7 family and cbx2 is crucial to clarify the role played by the let-7 family in the gonad of Japanese flounder. We predicted and verified the target interaction between the let-7 family and cbx2. The results showed that cbx2 was a direct target of let-7d, let-7e, let-7g, let-7j, and let-7b. Among them, let-7d, let-7e, let-7g, and let-7j exhibited an extremely significant targeting relationship with cbx2 (p < 0.001). Taking let-7g as an example, we further investigated the regulatory role between let-7g and cbx2 in the gonad by miRNA overexpression and inhibition experiments in primary testis cells. The results revealed that let-7g could negatively regulate cbx2 at the level of primary testis cells. And the expression of sf1 (steroidogenic factor 1) was also significantly decreased after the interference of cbx2 siRNA. This suggests that the let-7 family may be involved in the Japanese flounder gonadal development via targeting cbx2.

Identification and Characterization of Sex-Biased miRNAs in the Golden Pompano (Trachinotus blochii)

The golden pompano (Trachinotus blochii) is a marine fish of considerable commercial importance in China. It shows notable sexual size dimorphism; the growth rate of females is faster than that of males. Therefore, sex-biased research is of great importance in T. blochii breeding. However, there have been few studies on sex differentiation and mechanisms underlying sex determination in T. blochii. MicroRNAs (miRNAs) play crucial roles in sex differentiation and determination in animals. However, limited miRNA data are available on fish. In this study, two small RNA libraries prepared from the gonads of T. blochii were constructed and sequenced. The RNA-seq analysis yielded 1366 known and 69 novel miRNAs with 289 significantly differentially expressed miRNAs (p < 0.05). Gene ontology (GO) analysis confirmed that the TFIIA transcription factor complex (GO: 0005672) was the most significantly enriched GO term. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the differentially expressed miRNAs and target genes were mainly related to sex determination and gonadal developmental signaling pathways, specifically the Wnt signaling pathway, MAPK signaling pathway, and steroid biosynthetic pathway. MiRNA-mRNA co-expression network analysis strongly suggested a role for sex-biased miRNAs in sex determination/differentiation and gonadal development. For example, gata4, foxo3, wt1, and sf1 genes were found to be regulated by bta-miR-2898; esr2 and foxo3 by novel_176, and ar by oar-let-7b. Quantitative real-time polymerase chain reaction analysis of selected mRNAs and miRNAs validated the integrated analysis. This study established a set of sex-biased miRNAs that are potential regulatory factors in gonadal development in T. blochii. These results provide new insight into the function of miRNAs in sex differentiation and determination in T. blochii and highlight some key miRNAs for future studies.

Dynamic Expression and Gene Regulation of MicroRNAs During Bighead Carp (Hypophthalmichthys nobilis) Early Development

The early development of fish is regulated through dynamic and complex mechanisms involving the regulation of various genes. Many genes are subjected to post-transcriptional regulation by microRNAs (miRNAs). In the Chinese aquaculture industry, the native species bighead carp (Hypophthalmichthys nobilis) is important. However, the genetic regulation related to the early development of bighead carp is unknown. Here, we generated developmental profiles by miRNA sequencing to study the dynamic regulation of miRNAs during bighead carp early development. This study identified 1 046 miRNAs, comprising 312 known miRNAs and 734 uncharacterized miRNAs. Changes in miRNA expression were identified in the six early development stages. An obviously increased expression trend was detected during the development process, with the main burst of activity occurring after the earliest stage (early blastula, DS1). Investigations revealed that several miRNAs were dominantly expressed during the development process, especially in the later stages (e.g., miR-10b-5p, miR-21, miR-92a-3p, miR-206-3p, and miR-430a-3p), suggesting that these miRNAs exerted important functions during embryonic development. The differentially expressed miRNAs (DEMs) and time-serial analysis (profiles) of DEMs were analyzed. A total of 372 miRNAs were identified as DEMs (fold-change &gt;2, and false discovery rate &lt;0.05), and three expression profiles of the DEMs were detected to have co-expression patterns (r &gt; 0.7, and p &lt; 0.05). The broad negative regulation of target genes by miRNAs was speculated, and many development-related biological processes and pathways were enriched for the targets of the DEMs, which might be associated with maternal genome degradation and embryogenesis processes. In conclusion, we revealed the repertoire of miRNAs that are active during early development of bighead carp. These findings will increase our understanding of the regulatory mechanisms of early development of fish.

Epigenetic mechanisms in sex determination and in the evolutionary transitions between sexual systems

The hypothesis that epigenetic mechanisms of gene expression regulation have two main roles in vertebrate sex is presented. First, and within a given generation, by contributing to the acquisition and maintenance of (i) the male or female function once during the lifetime in individuals of gonochoristic species; and (ii) the male and female function in the same individual, either at the same time in simultaneous hermaphrodites, or first as one sex and then as the other in sequential hermaphrodites. Second, if environmental conditions change, epigenetic mechanisms may have also a role across generations, by providing the necessary phenotypic plasticity to facilitate the transition: (i) from one sexual system to another, or (ii) from one sex-determining mechanism to another. Furthermore, if the environmental change lasts enough time, epimutations could facilitate assimilation into genetic changes that stabilize the new sexual system or sex-determining mechanism. Examples supporting these assertions are presented, caveats or difficulties and knowledge gaps identified, and possible ways to test this hypothesis suggested.

This article is part of the theme issue ‘Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)’.

Identification of sex differentiation-related microRNA and long non-coding RNA in Takifugu rubripes gonads

Although sex determination and differentiation are key developmental processes in animals, the involvement of non-coding RNA in the regulation of this process is still not clarified. The tiger pufferfish (Takifugu rubripes) is one of the most economically important marine cultured species in Asia, but analyses of miRNA and long non-coding RNA (lncRNA) at early sex differentiation stages have not been conducted yet. In our study, high-throughput sequencing technology was used to sequence transcriptome libraries from undifferentiated gonads of T. rubripes. In total, 231 (107 conserved, and 124 novel) miRNAs were obtained, while 2774 (523 conserved, and 2251 novel) lncRNAs were identified. Of these, several miRNAs and lncRNAs were predicted to be the regulators of the expression of sex-related genes (including fru-miR-15b/foxl2, novel-167, novel-318, and novel-538/dmrt1, novel-548/amh, lnc_000338, lnc_000690, lnc_000370, XLOC_021951, and XR_965485.1/gsdf). Analysis of differentially expressed miRNAs and lncRNAs showed that three mature miRNAs up-regulated and five mature miRNAs were down-regulated in male gonads compared to female gonads, while 79 lncRNAs were up-regulated and 51 were down-regulated. These findings could highlight a group of interesting miRNAs and lncRNAs for future studies and may reveal new insights into the function of miRNAs and lncRNAs in sex determination and differentiation.

Distinct expression patterns of seven crucial microRNAs during early embryonic development in medaka (Oryzias latipes)

MicroRNAs (i.e. miRNAs) are small non-coding RNAs that play essential modulation roles in embryonic development in vertebrates. Paternal and maternal miRNAs contribute to the development of post-fertilization embryo and zygotic genome activation. The pattern of expression and their roles in embryonic development of medaka are not clearly understood. The present study, therefore, examined a temporal expression of seven miRNAs, ola-let-7a, ola-miR-202-3p, ola-miR-126-3p, ola-miR-122, ola-miR-92a, ola-miR-125a-3p and ola-miR-430a in sperm, oocytes, and embryos during early developmental stages. Three unique expression patterns of miRNAs were observed. ola-let7a, ola-miR-202-3p and ola-miR-126-3p showed both paternal and maternal expression, and ola-miR-122, ola-miR-92a, ola-miR-125a-3p showed maternal expression only. The expression of six out of seven miRNAs significantly decreased after maternal-zygotic transition (MZT), whereas ola-miR-430a expression initiated only after MZT. The temporal dynamic expression of these miRNAs suggests their potential roles in early embryogenesis and genome-zygotic activation in medaka.

miR-430a regulates the development of left-right asymmetry by targeting sqt in the teleost

microRNAs (miRNAs) are short, endogenous non-coding RNAs that contain approximately 18-22 nucleotides. miRNAs are involved in gene regulation by recognizing and binding the 3'UTR of target gene. In our previous data, miR-430 family showed significant differential expression modes through metamorphosis in Japanese flounder. It was speculated that miR-430a plays a key role in left-right patterning. We predicted the targets of miR-430a and gene ontology (GO) was performed. We speculated miR-430a is involved in the basal molecular function and organ development. In Japanese flounder, sqt as a target of miR-430a was enriched into heart development term. Sqt has been reported to participate in mesendoderm formation and organ development. Cardiac morphogenesis is the first asymmetric development process, which breaks left-right symmetry in bilateria. It was used as a marker to detect L-R asymmetric effects of miR-430a. Overexpression and suppression of miR-430a resulted in abnormal KV (Kupffer's vesicles) development and disordered in nodal-related expression with consequent cardiac laterality. Squint mRNA of Japanese flounder (Posqt) as a target of miR-430a was overexpressed and caused similar phenotype with miR-430a suppression group, such as longer cilia in KV and high range of clmc2 and spaw ectopic expression. Moreover, rescue experiments were performed and suggested that cardiac and KV defections, induced by overexpressing miR-430a, could be rescued by injecting Posqt mRNA. These results suggested that miR-430a regulates the development of left-right asymmetry by targeting sqt in the teleost.

The role of miR-92 in regulating early development and metamorphosis of Japanese flounder Paralichthys olivaceus

MicroRNAs are a class of short non-coding RNAs that contain approximately 22 nucleotides and play a regulatory role in RNA silencing and translational repression. miR-92 belongs to the miR-17-92 family and has a regulatory effect on cell proliferation, apoptosis, and expression of proto-oncogenes and tumor suppressor genes. However, its function in flatfish is unclear. In this study, we used farmed Japanese flounder, Paralichthys olivaceus, and showed that gata5 is a target gene of miR-92. Experiments on miR-92 overexpression indicated that gata5 and sox17 were downregulated, while the transcription level of ntl increased. By contrast, depletion of miR-92 resulted in increased gata5 and sox17 levels and reduced ntl level. Moreover, thiourea treatment indicated that miR-92 may inhibit the metamorphic development of Japanese flounder. Our study suggests that miR-92 regulates the fate of endoderm and mesoderm by controlling gata5.

Inventory of European Sea Bass (Dicentrarchus labrax) sncRNAs Vital During Early Teleost Development

During early animal ontogenesis, a plethora of small non-coding RNAs (sncRNAs) are greatly expressed and have been shown to be involved in several regulatory pathways vital to proper development. The rapid advancements in sequencing and computing methodologies in the last decade have paved the way for the production of sequencing data in a broad range of organisms, including teleost species. Consequently, this has led to the discovery of sncRNAs as well as the potentially novel roles of sncRNA in gene regulation. Among the several classes of sncRNAs, microRNAs (miRNAs) have, in particular, been shown to play a key role in development. The present work aims to identify the miRNAs that play important roles during early European sea bass (Dicentrarchus labrax) development. The European sea bass is a species of high commercial impact in European and especially Mediterranean aquaculture. This study reports, for the first time, the identification and characterization of small RNAs that play a part in the 10 developmental stages (from morula to all fins) of the European sea bass. From 10 developmental stages, more than 135 million reads, generated by next-generation sequencing, were retrieved from publicly available databases as well as newly generated. The analysis resulted in about 2,000 sample grouped reads, and their subsequently annotation revealed that the majority of transcripts belonged to the class of miRNAs followed by small nuclear RNAs and small nucleolar RNAs. The analysis of small RNA expression among the developmental stages under study revealed that miRNAs are active throughout development, with the main activity occurring after the earlier stages (morula and 50% epiboly) and at the later stages (first feeding, flexion, and all fins). Furthermore, investigating miRNAs exclusively expressed in one of the stages unraveled five miRNAs with a higher abundance only in the morula stage (miR-155, miR-430a, d1, d2, and miR-458), indicating possible important key roles of those miRNAs in further embryonic development. An additional target search showed putative miRNA-mRNA interactions with possible direct and indirect regulatory functions of the identified miRNAs.

Induced androgenetic development in rainbow trout and transcriptome analysis of irradiated eggs

Ionizing radiation is administered to damage nuclear genome in fish eggs during induced androgenesis. In this study, we examined whether 350 Gy of X-ray applied to damage chromosomes in the rainbow trout eggs affects maternal RNA. Shortly after irradiation, we did not find any symptoms of RNA degradation in the treated eggs. Significant (p < 0.01) differences between non-irradiated and irradiated eggs concerned only a few transcripts including increased expression of immediate early response 2 (IER2) and early growth response 1 (EGR1) genes observed in the irradiated eggs. Both genes belong to the group of “immediate early genes” that respond quickly to the diverse extracellular stimuli. Elevated expression of these genes was accompanied by decreased level of ssa-miR-10b-5p and ssa-miR-21b-5p (p < 0.05), for which IER2 and EGR1 are target genes. The level of RNA in the fertilized irradiated eggs was highly significantly lower than in the non-irradiated eggs (p < 0.001) and in the unfertilized irradiated eggs (p < 0.0001). However, transcriptome profiles of fertilized non-irradiated eggs and fertilized irradiated eggs did not differ significantly. Thus, we assume that reduced abundance of mRNA in the fertilized irradiated eggs was associated with post-translational degradation and clearance of the maternal transcripts rather than from the irradiation of eggs.

Global microRNA and isomiR expression associated with liver metabolism is induced by organophosphorus flame retardant exposure in male Chinese rare minnow (Gobiocypris rarus)

To reveal the adverse effects of organophosphorus flame retardants (OPFRs) on aquatic organisms at the epigenetic level, male Chinese rare minnows were exposed to 0.24 mg/L tris(2‑butoxyethyl) phosphate (TBOEP), 0.04 mg/L tris(1,3‑dichloro‑2‑propyl) phosphate (TDCIPP), or 0.012 mg/L triphenyl phosphate (TPHP) for 14 days. The effects of sub-acute OPFR exposure on liver miRNA and the 3' isomiR expression profiles of Chinese rare minnows were investigated. Through small RNA sequencing and bioinformatics analysis, a total of 32, 84, and 19 differentially expressed miRNAs were detected for TBOEP, TDCIPP, and TPHP exposure, respectively (p < 0.05). Target prediction of the differentially expressed miRNAs and pathway enrichment analysis indicated that predicted altered mRNAs for all three OPFRs were associated with metabolic pathways, whereas base excision repair was only predicted to be perturbed by the TPHP treatment. In addition, 3' isomiR-Us were unexpectedly abundant in all groups (e.g., miR-143), and TDCIPP strongly increased the ratio of 3' isomiR-U expression. Finally, histological examination and metabolic enzyme activity analyses werein agreement with the predicted metabolic pathways. As such, our study indicates that the investigation of epigenetics changes in miRNA gene transcription is a considerable method for the assessment of aquatic toxicity.

Marine Environmental Epigenetics

Marine organisms' persistence hinges on the capacity for acclimatization and adaptation to the myriad of interacting environmental stressors associated with global climate change. In this context, epigenetics-mechanisms that facilitate phenotypic variation through genotype-environment interactions-are of great interest ecologically and evolutionarily. Our comprehensive review of marine environmental epigenetics guides our recommendations of four key areas for future research: the dynamics of wash-in and wash-out of epigenetic effects, the mechanistic understanding of the interplay of different epigenetic marks and the interaction with the microbiome, the capacity for and mechanisms of transgenerational epigenetic inheritance, and the evolutionary implications of the interaction of genetic and epigenetic features. Emerging insights in marine environmental epigenetics can be applied to critical issues such as aquaculture, biomonitoring, and biological invasions, thereby improving our ability to explain and predict the responses of marine taxa to global climate change.

Investigating the aetiology of adverse events following HPV vaccination with systems vaccinology

In contrast to the insidious and poorly immunogenic human papillomavirus (HPV) infections, vaccination with the HPV virus-like particles (vlps) is non-infectious and stimulates a strong neutralizing-antibody response that protects HPV-naïve vaccinees from viral infection and associated cancers. However, controversy about alleged adverse events following immunization (AEFI) with the vlps have led to extensive reductions in vaccine acceptance, with countries like Japan dropping it altogether. The AEFIs are grouped into chronic fatigue syndrome/myalgic encephalomyelitis (CFS/ME). In this review, we present a hypothesis that the AEFIs might arise from malfunctions within the immune system when confronted with the unusual antigen. In addition, we outline how the pathophysiology of the AEFIs can be cost-effectively investigated with the holistic principles of systems vaccinology in a two-step process. First, comprehensive immunological profiles of HPV vaccinees exhibiting the AEFIs are generated by integrating the data derived from serological profiling for prominent HPV antibodies and serum cytokines, with data from serum metabolomics, peripheral white blood cells transcriptomics and gut microbiome profiling. Next, the immunological profiles are compared with corresponding profiles generated for matched (a) HPV vaccinees without AEFIs; (b) non-HPV-vaccinated individuals with CFS/ME-like symptoms; and (c) non-HPV-vaccinated individuals without CFS/ME. In these comparisons, any causal links between HPV vaccine and the AEFIs, as well as the underlying molecular basis for the links will be revealed. Such a study should provide an objective basis for evaluating HPV vaccine safety and for identifying biomarkers for individuals at risk of developing AEFI with HPV vaccination.

Circulating miRNA measurements are reflective of cholesterol-based changes in rainbow trout (Oncorhynchus mykiss)

MicroRNAs (miRNAs) are a class of small non-coding RNAs which are known to posttranscriptionally regulate the expression of most genes in both animals and plants. Meanwhile, studies have shown that numbers of miRNAs are present in body fluids including the plasma. Despite the mode of action of these circulating miRNAs still remains unknown, they have been found to be promising biomarkers for disease diagnosis, prognosis and response to treatment. In order to evaluate the potential of miRNAs as non-invasive biomarkers in aquaculture, a time-course experiment was implemented to investigate the postprandial regulation of miRNAs levels in liver and plasma as well as the hepatic expression of genes involved in cholesterol metabolism. We showed that miR-1, miR-33a, miR-122, miR-128 and miR-223 were expressed in the liver of rainbow trout and present at detectable level in the plasma. We also demonstrated that hepatic expression of miR-1, miR-122 and miR-128 were regulated by feed intake and reached their highest levels 12 hours after the meal. Interestingly, we observed that circulating levels of miR-128 and miR-223 are subjected to postprandial regulations similar to that observed in their hepatic counterparts. Statistical correlations were observed between liver and plasma for miR-128 and miR-223 and between hepatic and circulating miR-122, miR-128 and miR-223 and expression of genes related to cholesterol synthesis and efflux or glucose phosphorylation. These results demonstrated that circulating miR-122, miR-128 and miR-223 are potential biomarkers of cholesterol metabolism in rainbow trout.

MiR-205 Mediated Cu-Induced Lipid Accumulation in Yellow Catfish Pelteobagrus fulvidraco

The present working hypothesis is that the Cu-induced changes in lipid metabolism may be mediated by miRNAs. Here, we describe the miRNA profile of the liver tissues of yellow catfish exposed to waterborne Cu, based on larger-scale sequencing of small RNA libraries. We identified a total of 172 distinct miRNAs. Among these miRNAs, compared to the control, mRNA expression levels of 16 miRNAs (miR-203a, 205, 1788-3p, 375, 31, 196a, 203b-3p, 2187-5p, 196d, 459-3p, 153a and miR-725, and two novel-miRNAs: chr4-1432, chr-7684) were down-regulated, and mRNA levels of miR-212 and chr20-5274 were up-regulated in Cu-exposed group. The functions of their target genes mainly involved ether lipid metabolism, glycerophospholipid metabolism, linoleic acid metabolism and α-linolenic acid metabolism. Cu exposure inhibited the expression of miR-205, whose predicted target genes were enriched in the pathway of lipid metabolism, including fas, lxrα, ddit3, lamp2, casp3a and baxa. These potential target genes were further verified by Dual-luciferase reporter gene assay. Using primary hepatocytes of yellow catfish, Cu incubation down-regulated miR-205 expression, and increased TG contents and FAS activity. LXR antagonist effectively ameliorate the Cu-induced change of TG content and FAS activity. These data suggest that down-regulation of the miRNA-205 may be an important step in Cu-induced changes in lipid metabolism in yellow catfish.

MiR-202 controls female fecundity by regulating medaka oogenesis

Female gamete production relies on coordinated molecular and cellular processes that occur in the ovary throughout oogenesis. In fish, as in other vertebrates, these processes have been extensively studied both in terms of endocrine/paracrine regulation and protein expression and activity. The role of small non-coding RNAs in the regulation of animal reproduction remains however largely unknown and poorly investigated, despite a growing interest for the importance of miRNAs in a wide variety of biological processes. Here, we analyzed the role of miR-202, a miRNA predominantly expressed in male and female gonads in several vertebrate species. We studied its expression in the medaka ovary and generated a mutant line (using CRISPR/Cas9 genome editing) to determine its importance for reproductive success with special interest for egg production. Our results show that miR-202-5p is the most abundant mature form of the miRNA and that it is expressed in granulosa cells and in the unfertilized egg. The knock out (KO) of mir-202 gene resulted in a strong phenotype both in terms of number and quality of eggs produced. Mutant females exhibited either no egg production or produced a dramatically reduced number of eggs that could not be fertilized, ultimately leading to no reproductive success. We quantified the size distribution of the oocytes in the ovary of KO females and performed a large-scale transcriptomic analysis approach to identified dysregulated molecular pathways. Together, cellular and molecular analyses indicate that the lack of miR-202 impairs the early steps of oogenesis/folliculogenesis and decreases the number of large (i.e. vitellogenic) follicles, ultimately leading to dramatically reduced female fecundity. This study sheds new light on the regulatory mechanisms that control the early steps of follicular development, including possible targets of miR-202-5p, and provides the first in vivo functional evidence that a gonad-predominant microRNA may have a major role in female reproduction.

The role of small non-coding RNAs in the regulation of animal reproduction remains poorly investigated, despite a growing interest for the importance of miRNAs in a wide variety of biological processes. Here, we analyzed the role of miR-202, a miRNA predominantly expressed in gonads in vertebrate. We studied its expression in the medaka ovary and knocked out the mir-202 gene to study its importance for reproductive success. We showed that the lack of miR-202 results in the sterility of both females and males. In particular, it led to a drastic reduction of both the number and the quality of eggs produced by females. Mutant females exhibited either no egg production or produced a drastically reduced number of eggs that could not be fertilized, ultimately leading to no reproductive success. Quantitative histological and molecular analyses indicated that mir-202 KO impairs oocyte development and is also associated with the dysregulation of many genes that are critical for reproduction. This study sheds new light on the regulatory mechanisms that control oogenesis, including possible targets of miR-202-5p, and provides the first in vivo functional evidence that a gonad-predominant microRNA may have a major role in female reproduction.

Dumbbell-PCR for Discriminative Quantification of a Small RNA Variant

Cellular RNAs are often expressed as multiple isoforms of complex heterogeneity in both length and terminal sequences. IsomiRs, the isoforms of microRNAs, are such an example. Distinct quantification of each RNA variant is necessary to unravel the biogenesis mechanism and biological significance of heterogenetic RNA expression. Here we describe Dumbbell-PCR (Db-PCR), a TaqMan RT-PCR-based method that distinctively quantifies a specific small RNA variant with single-nucleotide resolution at terminal sequences. Db-PCR enables the quantitative analysis of RNA terminal heterogeneity without performing Next-Generation Sequencing.

Differential miRNA expression in the three-spined stickleback, response to environmental changes

miRNAs play important role in the various physiological and evolutionary processes, however, there is no data allowing comparison of evolutionary differences between various ecotypes adapted to different environmental conditions and specimen demonstrating immediate physiological response to the environmental changes. We compared miRNA expression profiles between marine and freshwater stickleback populations of the three-spined stickleback to identify the evolutionary differences. To study the immediate physiological response to foreign environment, we explored the changes induced by transfer of marine sticklebacks into freshwater environment and vice versa. Comparative analysis of changes in miRNA expression suggested that they are driven by three independent factors: (1) non-specific changes in miRNA expression under different environmental conditions; (2) specific response to freshwater conditions in the marine stickleback ecotype; (3) specific response to extreme osmotic conditions for both marine and freshwater ecotypes during the contact with non-native environment. Gene Ontology enrichment analysis of differential expressed miRNA targets supports our current hypothesis.

Identification and profiling of growth-related microRNAs in Chinese perch (Siniperca chuatsi)

Background

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play important roles in the regulation of diverse biological processes in eukaryotes. Chinese perch (Siniperca chuatsi) is one of the most economically important fish species widely cultured in China. Growth is an extremely important characteristic in fish. Individual differences in body size are common in Siniperca chuatsi, which significantly influence the aquaculture production of Siniperca chuatsi. However, the underline growth-related regulatory factors, such as miRNAs, are still unknown.

Results

To investigate the growth-related miRNAs in Siniperca chuatsi, two RNA libraries from four growth-related tissues (brain, pituitary, liver, and muscle) of Siniperca chuatsi at 6-month stage with relatively high or low growth rates (big-size group or small-size group) were obtained and sequenced using Solexa sequencing. A total of 252 known miRNAs and 12 novel miRNAs were identified. The expression patterns of these miRNAs in big-size group and small-size group were compared, and the results showed that 31 known and 5 novel miRNAs were differently expressed (DE). Furthermore, to verify the Solexa sequencing, five DE miRNAs were randomly selected and quantified by quantitative reverse transcription polymerase chain reaction (qRT-PCR). The results showed that their expression patterns were consistent with those of Solexa sequencing. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of target genes of DE miRNAs was performed. It showed that the target genes were involved in multiple biological processes including metabolic process, suggesting that metabolic process played an important role in growth of fish.

Conclusions

Siniperca chuatsi is a popular and economically important species in aquaculture. In this study, miRNAs in Siniperca chuatsi with different growth rates were identified, and their expression profiles were compared. The data provides the first large-scale miRNA profiles related to growth of Siniperca chuatsi, which is expected to contribute to a better understanding of the role of miRNAs in regulating the biological processes of growth and possibly useful for Siniperca chuatsi breeding.

Electronic supplementary material

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

Comparative microRNAome analysis of the testis and ovary of the Chinese giant salamander

MicroRNAs (miRNAs) are 18-24 nucleotides non-coding RNAs that regulate gene expression by post-transcriptional suppression of mRNA. The Chinese giant salamander (CGS, Andrias davidianus), which is an endangered species, has become one of the important models of animal evolution; however, no miRNA studies on this species have been conducted. In this study, two small RNA libraries of CGS ovary and testis were constructed using deep sequencing technology. A bioinformatics pipeline was developed to distinguish miRNA sequences from other classes of small RNAs represented in the sequencing data. We found that many miRNAs and other small RNAs such as piRNA and tsRNA were abundant in CGS tissue. A total of 757 and 756 unique miRNAs were annotated as miRNA candidates in the ovary and testis respectively. We identified 145 miRNAs in CGS ovary and 155 miRNAs in CGS testis that were homologous to those in Xenopus laevis ovary and testis respectively. Forty-five miRNAs were more highly expressed in ovary than in testis and 21 miRNAs were more highly expressed in testis than in ovary. The expression profiles of the selected miRNAs (miR-451, miR-10c, miR-101, miR-202, miR-7a and miR-499) had their own different roles in other eight tissues and different development stages of testis and ovary, suggesting that these miRNAs play vital regulatory roles in sexual differentiation, gametogenesis and development in CGS. To our knowledge, this is the first study to reveal miRNA profiles that are related to male and female CGS gonads and provide insights into sex differences in miRNA expression in CGS.

Evolution of Fish Let-7 MicroRNAs and Their Expression Correlated to Growth Development in Blunt Snout Bream

The lethal-7 (let-7) miRNA, known as one of the first founding miRNAs, is present in multiple copies in a genome and has diverse functions in animals. In this study, comparative genomic analysis of let-7 miRNAs members in fish species indicated that let-7 miRNA is a sequence conserved family in fish, while different species have the variable gene copy numbers. Among the ten members including let-7a/b/c/d/e/f/g/h/i/j, the let-7a precursor sequence was more similar to ancestral sequences, whereas other let-7 miRNA members were separate from the late differentiation of let-7a. The mostly predicted target genes of let-7 miRNAs are involved in biological process, especially developmental process and growth through Gene Ontology (GO) enrichment analysis. In order to identify the possible different functions of these ten miRNAs in fish growth development, their expression levels were quantified in adult males and females of Megalobrama amblycephala, as well as in 3-, 6-, and 12-months-old individuals with relatively slow- and fast-growth rates. These ten miRNAs had similar tissue expression patterns between males and females, with higher expression levels in the brain and pituitary than that in other tissues (p < 0.05). Among these miRNAs, the relative expression level of let-7a was the highest among almost all the tested tissues, followed by let-7b, let-7d and let-7c/e/f/g/h/i/j. As to the groups with different growth rates, the expression levels of let-7 miRNAs in pituitary and brain from the slow-growth group were always significantly higher than that in the fast-growth group (p < 0.05). These results suggest that let-7 miRNA members could play an important role in the regulation of growth development in M. amblycephala through negatively regulating expression of their target genes.

Dynamics of miRNA transcriptome during gonadal development of zebrafish

Studies in non-teleost vertebrates have found microRNAs (miRNAs) to be essential for proper gonadal development. However, comparatively little is known about their role during gonadal development in teleost fishes. So far in zebrafish, a model teleost, transcript profiling throughout gonadal development has not been established because of a tiny size of an organ in juvenile stages and its poor distinguishability from surrounding tissues. We performed small RNA sequencing on isolated gonads of See-Thru-Gonad line, from the undifferentiated state at 3 weeks post fertilization (wpf) to fully mature adults at 24 wpf. We identified 520 gonadal mature miRNAs; 111 of them had significant changes in abundance over time, while 50 miRNAs were either testis- or ovary-enriched significantly in at least one developmental stage. We characterized patterns of miRNA abundance over time including isomiR variants. We identified putative germline versus gonadal somatic miRNAs through differential small RNA sequencing of isolated gametes versus the whole gonads. This report is the most comprehensive analysis of the miRNA repertoire in zebrafish gonads during the sexual development to date and provides an important database from which functional studies can be performed.

Identification and characterization of the expression profile of the microRNAs in the Amazon species Colossoma macropomum by next generation sequencing

Colossoma macropomum is a resistant species native of Amazonas and Orinoco river basins. It is regarded as the second largest finfish of Solimões and Amazon rivers, representing a major fishery resource in Amazonas and an important species in tropical aquaculture. MicroRNAs are non-coding endogenous riboregulators of nearly 22 nucleotides that play a key role in post-transcriptional gene regulation of several organisms. We analyzed samples of liver and skin from specimens of C. macropomum using next generation sequencing. The dataset was evaluated using computational programs to check the quality of sequences, identification of miRNAs, as well as to evaluate the expression levels of these microRNAs and interaction of target genes. We identified 279 conserved miRNAs, being 257 from liver and 272 from skin, with several miRNAs shared between tissues, with divergence in the number of reads. The strands miR-5p and miR-3p were observed in 72 miRNAs, some of them presenting a higher number of 3p reads. The functional annotation of the most expressed miRNAs resulted in 27 pathways for the liver and skin mainly related to the "biological processes" domain. Based on the identified pathways, we visualized a large gene network, suggesting the regulation of selected miRNA over this interactive dataset. We were able to identify and characterize the expression levels of miRNAs in two tissues of great activity in C. macropomum, which stands out as the beginning of several studies that can be carried out to elucidate the influence of miRNAs in this species and their applicability as biotechnological tools.

Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production

Flatfish have a high market acceptance thus representing a profitable aquaculture production. The main farmed species is the turbot (Scophthalmus maximus) followed by Japanese flounder (Paralichthys olivaceous) and tongue sole (Cynoglossus semilaevis), but other species like Atlantic halibut (Hippoglossus hippoglossus), Senegalese sole (Solea senegalensis) and common sole (Solea solea) also register an important production and are very promising for farming. Important genomic resources are available for most of these species including whole genome sequencing projects, genetic maps and transcriptomes. In this work, we integrate all available genomic information of these species within a common framework, taking as reference the whole assembled genomes of turbot and tongue sole (>210× coverage). New insights related to the genetic basis of productive traits and new data useful to understand the evolutionary origin and diversification of this group were obtained. Despite a general 1:1 chromosome syntenic relationship between species, the comparison of turbot and tongue sole genomes showed huge intrachromosomic reorganizations. The integration of available mapping information supported specific chromosome fusions along flatfish evolution and facilitated the comparison between species of previously reported genetic associations for productive traits. When comparing transcriptomic resources of the six species, a common set of ~2500 othologues and ~150 common miRNAs were identified, and specific sets of putative missing genes were detected in flatfish transcriptomes, likely reflecting their evolutionary diversification.

Illumina Sequencing Reveals Aberrant Expression of MicroRNAs and Their Variants in Whitefish (Coregonus lavaretus) Liver after Exposure to Microcystin-LR

Molecular analyses show that challenging fish with microcystin-LR (MC-LR) causes perturbations of microRNA (miRNA) signaling. However, the significance and scope of these alterations is currently unknown. To address this issue, we studied miRNA gene expression in the liver of juvenile whitefish, C. lavaretus, during 28 days of exposure to a subacute dose of MC-LR (100 μg·kg^-1 body mass). Using genomic resources of Atlantic salmon (AGKD03), the mature miRNA library of Atlantic salmon (miRBase-21) and bioinformatics tools (sRNAbench), we discovered and annotated a total of 377 distinct mature miRNAs belonging to 93 families of evolutionary conserved miRNAs, as well as 24 novel mature miRNA candidates that were mapped to 14 distinct S. salar miRNA precursors. miRNA-Seq transcriptome profiling of liver tissues revealed differential miRNA expression in control and treated fish at 14 days (73 miRNAs were modulated) and at 28 days (83 miRNAs) of the treatment, subsequently validated by qPCR for nine selected differentially expressed miRNAs. Additional qPCR study confirmed the miRNA-Seq data and revealed consistent, aberrant miRNAs expression profile in the later phase of MC-LR hepatotoxicity (7–28 d). Functional annotation analysis revealed that the aberrantly expressed miRNAs have target genes involved in cytoskeletal remodeling, cell metabolism, cell cycle regulation and apoptosis; dysregulation of these processes in liver cells leads to cirrhosis and hepatocellular carcinoma in humans. To enable deeper insight into the molecular responses of liver cells in fish exposed to MC-LR, we expanded the miRNAome analysis by inclusion of miRNA variants (isomiRs) profiles, and we showed that the isomiR profiles of liver specific MiR122, and a few other miRNAs, correlated with MC-LR treatment. Given the importance of isomiRs for disease biology in mammals, we believe that further research focused on the miRNA isoforms will bring us closer to better understanding the molecular mechanisms of MC-LR hepatotoxicity.

Identification of novel microRNA genes in freshwater and marine ecotypes of the three-spined stickleback (Gasterosteus aculeatus)

The three-spined stickleback (Gasterosteus aculeatus L.) is an important model organism for studying the molecular mechanisms of speciation and adaptation to salinity. Despite increased interest to microRNA discovery and recent publication on microRNA prediction in the three-spined stickleback using bioinformatics approaches, there is still a lack of experimental support for these data. In this paper, high-throughput sequencing technology was applied to identify microRNA genes in gills of the three-spined stickleback. In total, 595 miRNA genes were discovered; half of them were predicted in previous computational studies and were confirmed here as microRNAs expressed in gill tissue. Moreover, 298 novel microRNA genes were identified. The presence of miRNA genes in selected 'divergence islands' was analysed and 10 miRNA genes were identified as not randomly located in 'divergence islands'. Regulatory regions of miRNA genes were found enriched with selective SNPs that may play a role in freshwater adaptation.

Dysregulation of miRNA isoform level at 5' end in Alzheimer's disease

Alzheimer's disease (AD) is the most common form of dementia, whose mechanism is still not yet fully understood. A miRNA-based signature method, commonly according to the changes of expression levels, is widely used for AD analysis in previous studies. Recently, miRNA isoforms called as isomiR variants, which is considered to play important biological roles, have been demonstrated as the applications of high throughput sequencing platforms. Here, we presented an entropy-based model to detect the miRNA isoform level at the 5' end, and found many miRNAs with significant changes of isoform levels between the early stage and the late stage of AD by the application of this model to the public data. The statistical significance of the overlap between isoform-level changed miRNAs and AD related miRNAs extracted from HMDD2 supports that these miRNA isoforms are not degradation products. Based on the most common isomiR seed analysis of isoform-level changed AD related miRNAs, the predicted targets are also found to be enriched for genes involved in transcriptional regulation and the nervous system. After comparing with the expression level based method, we detected that changes of 5' isoform levels are more stable than those of expression levels for AD related miRNA detecting.

Comprehensive analysis of small RNA-seq data reveals that combination of miRNA with its isomiRs increase the accuracy of target prediction in Arabidopsis thaliana

Along with the canonical miRNA, distinct miRNA-like sequences called sibling miRNAs (sib-miRs) are generated from the same pre-miRNA. Among them, isomeric sequences featuring slight variations at the terminals, relative to the canonical miRNA, constitute a pool of isomeric sibling miRNAs (isomiRs). Despite the high prevalence of isomiRs in eukaryotes, their features and relevance remain elusive. In this study, we performed a comprehensive analysis of mature precursor miRNA (pre-miRNA) sequences from Arabidopsis to understand their features and regulatory targets. The influence of isomiR terminal heterogeneity in target binding was examined comprehensively. Our comprehensive analyses suggested a novel computational strategy that utilizes miRNA and its isomiRs to enhance the accuracy of their regulatory target prediction in Arabidopsis. A few targets are shared by several members of isomiRs; however, this phenomenon was not typical. Gene Ontology (GO) enrichment analysis showed that commonly targeted mRNAs were enriched for certain GO terms. Moreover, comparison of these commonly targeted genes with validated targets from published data demonstrated that the validated targets are bound by most isomiRs and not only the canonical miRNA. Furthermore, the biological role of isomiRs in target cleavage was supported by degradome data. Incorporating this finding, we predicted potential target genes of several miRNAs and confirmed them by experimental assays. This study proposes a novel strategy to improve the accuracy of predicting miRNA targets through combined use of miRNA with its isomiRs.

Micromanaging metabolism-a role for miRNAs in teleost energy metabolism

MicroRNAs (miRNAs) are small, non-protein coding RNA sequences, which are found in most eukaryotes. Since their initial discovery, miRNAs have emerged as important regulators of many biological processes. One of the most important processes profoundly regulated by miRNAs is energy metabolism. Traditionally, metabolic functions of miRNAs have been studied in genome-sequenced mammalian organisms, especially the mouse model. However, partially driven by commercial interest in aquaculture, increasingly feasible large-scale molecular techniques have resulted in the characterization of miRNA repertoires, and importantly, several genome sequences of several (commercially important) teleost species, which also hold important roles as research models in the comparative physiology of energy metabolism. This review aims to introduce the recent advances in miRNA research in teleost fish and to describe the current knowledge of miRNA function in teleost energy metabolism. The most pressing research needs and questions to determine metabolic roles of miRNAs in teleost models are presented, as well as applicable technical approaches and current bottlenecks. Rainbow trout, which possess the advantages of newly available molecular tools and a long history as comparative research model in teleost energy metabolism, are discussed as a promising research model to address these questions.

De novo sequencing and analysis of Lophophora williamsii transcriptome, and searching for putative genes involved in mescaline biosynthesis

Background

Lophophora williamsii (commonly named peyote) is a small, spineless cactus with psychoactive alkaloids, particularly mescaline. Peyote utilizes crassulacean acid metabolism (CAM), an alternative form of photosynthesis that exists in succulents such as cacti and other desert plants. Therefore, its transcriptome can be considered an important resource for future research focused on understanding how these plants make more efficient use of water in marginal environments and also for research focused on better understanding of the overall mechanisms leading to production of plant natural products and secondary metabolites.

Results

In this study, two cDNA libraries were generated from L. williamsii. These libraries, representing buttons (tops of stems) and roots were sequenced using different sequencing platforms (GS-FLX, GS-Junior and PGM, respectively). A total of 5,541,550 raw reads were generated, which were assembled into 63,704 unigenes with an average length of 564.04 bp. A total of 25,149 unigenes (62.19 %) was annotated using public databases. 681 unigenes were found to be differentially expressed when comparing the two libraries, where 400 were preferentially expressed in buttons and 281 in roots. Some of the major alkaloids, including mescaline, were identified by GC-MS and relevant metabolic pathways were reconstructed using the Kyoto encyclopedia of genes and genomes database (KEGG). Subsequently, the expression patterns of preferentially expressed genes putatively involved in mescaline production were examined and validated by qRT-PCR.

Conclusions

High throughput transcriptome sequencing (RNA-seq) analysis allowed us to efficiently identify candidate genes involved in mescaline biosynthetic pathway in L. williamsii; these included tyrosine/DOPA decarboxylase, hydroxylases, and O-methyltransferases. This study sets the theoretical foundation for bioassay design directed at confirming the participation of these genes in mescaline production.

Electronic supplementary material

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

Variability of miRNA expression during the differentiation of human embryonic stem cells into retinal pigment epithelial cells

Embryonic stem cells (ESCs) and induced pluripotent stem cells can be induced to differentiate into retinal pigment epithelium (RPE). MiRNAs have been characterized and found playing important roles in the differentiation process of ESCs, but their length and sequence heterogeneity (isomiRs), and their non-canonical forms of miRNAs are underestimated or ignored. In this report, we found some non-canonical miRNAs (dominant isomiRs) in all differentiation stages, and 27 statistically significant editing sites were identified in 24 different miRNAs. Moreover, we found marked major-to-minor arm-switching events in 14 pre-miRNAs during the hESC to RPE cell differentiation phases. Our study for the first time reports exploring the variability of miRNA expression during the differentiation of hESCs into RPE cells and the results show that miRNA variability is a ubiquitous phenomenon in the ESC differentiation.

Dynamics of gene expression patterns during early development of the European seabass (Dicentrarchus labrax)

Larval and embryonic stages are the most critical period in the life cycle of marine fish. Key developmental events occur early in development and are influenced by external parameters like stress, temperature, salinity, and photoperiodism. Any failure may cause malformations, developmental delays, poor growth, and massive mortalities. Advanced understanding of molecular processes underlying marine larval development may lead to superior larval rearing conditions. Today, the new sequencing and bioinformatic methods allow transcriptome screens comprising messenger (mRNA) and microRNA (miRNA) with the scope of detecting differential expression for any species of interest. In the present study, we applied Illumina technology to investigate the transcriptome of early developmental stages of the European seabass ( Dicentrarchus labrax). The European seabass, in its natural environment, is a euryhaline species and has shown high adaptation processes in early life phases. During its embryonic and larval phases the European seabass lives in a marine environment and as a juvenile it migrates to coastal zones, estuaries, and lagoons. Investigating the dynamics of gene expression in its early development may shed light on factors promoting phenotypic plasticity and may also contribute to the improvement and advancement of rearing methods of the European seabass, a species of high economic importance in European and Mediterranean aquaculture. We present the identification, characterization, and expression of mRNA and miRNA, comprising paralogous genes and differentially spliced transcripts from early developmental stages of the European seabass. We further investigated the detection of possible interactions of miRNA with mRNA.

Temperature during early development has long-term effects on microRNA expression in Atlantic cod

Background

Environmental temperature has serious implications in life cycle of aquatic ectotherms. Understanding the molecular mechanisms of temperature acclimation and adaptation of marine organisms is of the uttermost importance for ecology, fisheries, and aquaculture, as it allows modeling the effects of global warming on population dynamics. Regulatory molecules are major modulators of acclimation and adaptation; among them, microRNAs (miRNAs) are versatile and substantial contributors to regulatory networks of development and adaptive plasticity. However, their role in thermal plasticity is poorly known. We have asked whether the temperature and its shift during the early ontogeny (embryonic and larval development) affect the miRNA repertoire of Atlantic cod (Gadus morhua), and if thermal experience has long-term consequences in the miRNA profile.

Results

We characterized miRNA during different developmental stages and in juvenile tissues using next generation sequencing. We identified 389 putative miRNA precursor loci, 120 novel precursor miRNAs, and 281 mature miRNAs. Some miRNAs showed stage- or tissue-enriched expression and miRNAs, such as the miR-17 ~ 92 cluster, myomiRs (miR-206), neuromiRs (miR-9, miR-124), miR-130b, and miR-430 showed differential expression in different temperature regimes. Long-term effect of embryonic incubation temperature was revealed on expression of some miRNAs in juvenile pituitary (miR-449), gonad (miR-27c, miR-30c, and miR-200a), and liver (let-7 h, miR-7a, miR-22, miR-34c, miR-132a, miR-192, miR-221, miR-451, miR-2188, and miR-7550), but not in brain. Some of differentially expressed miRNAs in the liver were confirmed using LNA-based rt-qPCR. The effect of temperature on methylation status of selected miRNA promoter regions was mostly inconclusive.

Conclusions

Temperature elevation by several degrees during embryonic and larval developmental stages significantly alters the miRNA profile, both short-term and long-term. Our results suggest that a further rise in seas temperature might affect life history of Atlantic cod.

Electronic supplementary material

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

Entropy-based model for miRNA isoform analysis

MiRNAs have been widely studied due to their important post-transcriptional regulatory roles in gene expression. Many reports have demonstrated the evidence of miRNA isoform products (isomiRs) in high-throughput small RNA sequencing data. However, the biological function involved in these molecules is still not well investigated. Here, we developed a Shannon entropy-based model to estimate isomiR expression profiles of high-throughput small RNA sequencing data extracted from miRBase webserver. By using the Kolmogorov-Smirnov statistical test (KS test), we demonstrated that the 5p and 3p miRNAs present more variants than the single arm miRNAs. We also found that the isomiR variant, except the 3’ isomiR variant, is strongly correlated with Minimum Free Energy (MFE) of pre-miRNA, suggesting the intrinsic feature of pre-miRNA should be one of the important factors for the miRNA regulation. The functional enrichment analysis showed that the miRNAs with high variation, particularly the 5’ end variation, are enriched in a set of critical functions, supporting these molecules should not be randomly produced. Our results provide a probabilistic framework for miRNA isoforms analysis, and give functional insights into pre-miRNA processing.

Dumbbell-PCR: a method to quantify specific small RNA variants with a single nucleotide resolution at terminal sequences

Recent advances in next-generation sequencing technologies have revealed that cellular functional RNAs are not always expressed as single entities with fixed terminal sequences but as multiple isoforms bearing complex heterogeneity in both length and terminal sequences, such as isomiRs, the isoforms of microRNAs. Unraveling the biogenesis and biological significance of heterogenetic RNA expression requires distinctive analysis of each RNA variant. Here, we report the development of dumbbell PCR (Db-PCR), an efficient and convenient method to distinctively quantify a specific individual small RNA variant. In Db-PCR, 5′- and 3′-stem–loop adapters are specifically hybridized and ligated to the 5′- and 3′-ends of target RNAs, respectively, by T4 RNA ligase 2 (Rnl2). The resultant ligation products with ‘dumbbell-like’ structures are subsequently quantified by TaqMan RT-PCR. We confirmed that high specificity of Rnl2 ligation and TaqMan RT-PCR toward target RNAs assured both 5′- and 3′-terminal sequences of target RNAs with single nucleotide resolution so that Db-PCR specifically detected target RNAs but not their corresponding terminal variants. Db-PCR had broad applicability for the quantification of various small RNAs in different cell types, and the results were consistent with those from other quantification method. Therefore, Db-PCR provides a much-needed simple method for analyzing RNA terminal heterogeneity.

MicroRNAs in melanocyte and melanoma biology

The importance of microRNAs as key molecular components of cellular processes is now being recognized. Recent reports have shown that microRNAs regulate processes as diverse as protein expression and nuclear functions inside cells and are able to signal extracellularly, delivered via exosomes, to influence cell fate at a distance. The versatility of microRNAs as molecular tools inspires the design of novel strategies to control gene expression, protein stability, DNA repair and chromatin accessibility that may prove very useful for therapeutic approaches due to the extensive manageability of these small molecules. However, we still lack a comprehensive understanding of the microRNA network and its interactions with the other layers of regulatory elements in cellular and extracellular functions. This knowledge may be necessary before we exploit microRNA versatility in therapeutic settings. To identify rules of interactions between microRNAs and other regulatory systems, we begin by reviewing microRNA activities in a single cell type: the melanocyte, from development to disease.

Characterization and expression analysis of microRNAs in the tube foot of sea cucumber Apostichopus japonicus

MicroRNAs (miRNAs) are a class of endogenous non-coding small RNA with average length of 22 nucleotides, participating in the post-transcriptional regulation of gene expression. In this study, we report the identification and characterization of miRNAs in the tube foot of sea cucumber (Apostichopus japonicus) by next generation sequencing with Illumina HiSeq 2000 platform. Through the bioinformatic analysis, we identified 260 conserved miRNAs and six novel miRNAs from the tube foot small RNA transcriptome. Quantitative realtime PCR (qRT-PCR) was performed to characterize the specific expression in the tube foot. The results indicated that four miRNAs, including miR-29a, miR-29b, miR-2005 and miR-278-3p, were significantly up-regulated in the tube foot. The target genes of the four specifically expressed miRNAs were predicted in silico and validated by performing qRT-PCR. Gene ontology (GO) and KEGG pathway analyses with the target genes of these four miRNAs were conducted to further understand the regulatory function in the tube foot. This is the first study to profile the miRNA transcriptome of the tube foot in sea cucumber. This work will provide valuable genomic resources to understand the mechanisms of gene regulation in the tube foot, and will be useful to assist the molecular breeding in sea cucumber.

MicroRNA in teleost fish

MicroRNAs (miRNAs) are transcriptional and posttranscriptional regulators involved in nearly all known biological processes in distant eukaryotic clades. Their discovery and functional characterization have broadened our understanding of biological regulatory mechanisms in animals and plants. They show both evolutionary conserved and unique features across Metazoa. Here, we present the current status of the knowledge about the role of miRNA in development, growth, and physiology of teleost fishes, in comparison to other vertebrates. Infraclass Teleostei is the most abundant group among vertebrate lineage. Fish are an important component of aquatic ecosystems and human life, being the prolific source of animal proteins worldwide and a vertebrate model for biomedical research. We review miRNA biogenesis, regulation, modifications, and mechanisms of action. Specific sections are devoted to the role of miRNA in teleost development, organogenesis, tissue differentiation, growth, regeneration, reproduction, endocrine system, and responses to environmental stimuli. Each section discusses gaps in the current knowledge and pinpoints the future directions of research on miRNA in teleosts.

Identification and differential expression of hepatopancreas microRNAs in red swamp crayfish fed with emodin diet

Using high-throughput Illumina Solexa system, the differential miRNA expressions from hepatopancreas in red swamp crayfish (Procambarus clarkii) fed with diets containing 0 (control) and 75 mg emodin kg(-1) (trial) were identified, respectively. As a result, 13,335,928 raw reads from the control sample and 14,938,951 raw reads from the trial sample were obtained while 13,053,344 (98.77%) and 14,517,522 (98.34%) small RNA were identified, respectively. 106 mature miRNAs (belonging to 68 miRNA gene families) were identified. 35 miRNAs displayed significantly differential expressions between two libraries. Of these, comparing to the control library, 6 miRNAs were significantly up-regulated and 29 miRNAs were significantly down-regulated. Moreover, 5 novel miRNAs (2 from control sample, 3 from trial sample) and target genes were predicted. GO analysis suggested that these miRNAs might be involved in innate immune response, growth, metabolism, cellular process, biological regulation and stimulus response. Our knowledge from this study could contribute to a better understanding of the miRNAs roles in regulating innate immune response and the study of miRNA function in crayfish.

Thermal plasticity of the miRNA transcriptome during Senegalese sole development

Background

Several miRNAs are known to control myogenesis in vertebrates. Some of them are specifically expressed in muscle while others have a broader tissue expression but are still involved in establishing the muscle phenotype. In teleosts, water temperature markedly affects embryonic development and larval growth. It has been previously shown that higher embryonic temperatures promoted faster development and increased size of Senegalese sole (Solea senegalensis) larvae relatively to a lower temperature. The role of miRNAs in thermal-plasticity of growth is hitherto unknown. Hence, we have used high-throughput SOLiD sequencing to determine potential changes in the miRNA transcriptome in Senegalese sole embryos that were incubated at 15°C or 21°C until hatching and then reared at a common temperature of 21°C.

Results

We have identified 320 conserved miRNAs in Senegalese sole, of which 48 had not been previously described in teleosts. mir-17a-5p, mir-26a, mir-130c, mir-206-3p, mir-181a-5p, mir-181a-3p and mir-199a-5p expression levels were further validated by RT- qPCR. The majority of miRNAs were dynamically expressed during early development, with peaks of expression at pre-metamorphosis or metamorphosis. Also, a higher incubation temperature (21°C) was associated with expression of some miRNAs positively related with growth (e.g., miR-17a, miR-181-5p and miR-206) during segmentation and at hatching. Target prediction revealed that these miRNAs may regulate myogenesis through MAPK and mTOR pathways. Expression of miRNAs involved in lipid metabolism and energy production (e.g., miR-122) also differed between temperatures. A miRNA that can potentially target calpain (miR-181-3p), and therefore negatively regulate myogenesis, was preferentially expressed during segmentation at 15°C compared to 21°C.

Conclusions

Temperature has a strong influence on expression of miRNAs during embryonic and larval development in fish. Higher expression levels of miR-17a, miR-181-5p and miR-206-3p and down-regulation of miR-181a-3p at 21°C may promote myogenesis and are in agreement with previous studies in Senegalese sole, which reported enhanced growth at higher embryonic temperatures compared to 15°C. Moreover, miRNAs involved in lipid metabolism and energy production may also contribute to increased larval growth at 21°C compared to 15°C. Taken together, our data indicate that miRNAs may play a role in temperature-induced phenotypic plasticity of growth in teleosts.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-525) contains supplementary material, which is available to authorized users.

Differential expression profiling of spleen microRNAs in response to two distinct type II interferons in Tetraodon nigroviridis

MicroRNAs are endogenous, small non-coding RNAs approximately 18-26 nucleotides in length that regulate target gene expression at the post-transcription level. Interferon-γ (IFN-γ) is a Th1 cytokine that is involved in both the innate and adaptive immune responses. We previously identified two IFN-γ genes in green-spotted puffer fish (Tetraodon nigroviridis). To determine whether miRNAs participate in IFN-γ-related immune responses, T. nigroviridis spleen cells were treated with recombinant IFN-γ isoforms, and a Solexa high-throughput sequencing method was used to identify miRNAs. In total, 1,556, 1,538 and 1,573 miRNAs were found in the three samples, and differentially expressed miRNAs were determined. In total, 398 miRNAs were differentially expressed after rIFN-γ1 treatment, and 438 miRNAs were differentially expressed after rIFN-γ2 treatment; additionally, 403 miRNAs were differentially expressed between the treatment groups. Ten differentially expressed miRNAs were chosen for validation using qRT-PCR. Target genes for the differentially expressed miRNAs were predicted, and GO and KEGG analyses were performed. This study provides basic knowledge regarding fish IFN-γ-induced miRNAs and offers clues for further studies into the mechanisms underlying fish IFN-γ-mediated immune responses.

Performance comparison of digital microRNA profiling technologies applied on human breast cancer cell lines

MicroRNA profiling represents an important first-step in deducting individual RNA-based regulatory function in a cell, tissue, or at a specific developmental stage. Currently there are several different platforms to choose from in order to make the initial miRNA profiles. In this study we investigate recently developed digital microRNA high-throughput technologies. Four different platforms were compared including next generation SOLiD ligation sequencing and Illumina HiSeq sequencing, hybridization-based NanoString nCounter, and miRCURY locked nucleic acid RT-qPCR. For all four technologies, full microRNA profiles were generated from human cell lines that represent noninvasive and invasive tumorigenic breast cancer. This study reports the correlation between platforms, as well as a more extensive analysis of the accuracy and sensitivity of data generated when using different platforms and important consideration when verifying results by the use of additional technologies. We found all the platforms to be highly capable for microRNA analysis. Furthermore, the two NGS platforms and RT-qPCR all have equally high sensitivity, and the fold change accuracy is independent of individual miRNA concentration for NGS and RT-qPCR. Based on these findings we propose new guidelines and considerations when performing microRNA profiling.

Characterization of novel precursor miRNAs using next generation sequencing and prediction of miRNA targets in Atlantic halibut

BACKGROUND

microRNAs (miRNAs) are implicated in regulation of many cellular processes. miRNAs are processed to their mature functional form in a step-wise manner by multiple proteins and cofactors in the nucleus and cytoplasm. Many miRNAs are conserved across vertebrates. Mature miRNAs have recently been characterized in Atlantic halibut (Hippoglossus hippoglossus L.). The aim of this study was to identify and characterize precursor miRNA (pre-miRNAs) and miRNA targets in this non-model flatfish. Discovery of miRNA precursor forms and targets in non-model organisms is difficult because of limited source information available. Therefore, we have developed a methodology to overcome this limitation.

METHODS

Genomic DNA and small transcriptome of Atlantic halibut were sequenced using Roche 454 pyrosequencing and SOLiD next generation sequencing (NGS), respectively. Identified pre- miRNAs were further validated with reverse-transcription PCR. miRNA targets were identified using miRanda and RNAhybrid target prediction tools using sequences from public databases. Some of miRNA targets were also identified using RACE-PCR. miRNA binding sites were validated with luciferase assay using the RTS34st cell line.

RESULTS

We obtained more than 1.3 M and 92 M sequence reads from 454 genomic DNA sequencing and SOLiD small RNA sequencing, respectively. We identified 34 known and 9 novel pre-miRNAs. We predicted a number of miRNA target genes involved in various biological pathways. miR-24 binding to kisspeptin 1 receptor-2 (kiss1-r2) was confirmed using luciferase assay.

CONCLUSION

This study demonstrates that identification of conserved and novel pre-miRNAs in a non-model vertebrate lacking substantial genomic resources can be performed by combining different next generation sequencing technologies. Our results indicate a wide conservation of miRNA precursors and involvement of miRNA in multiple regulatory pathways, and provide resources for further research on miRNA in non-model animals.

A combined strategy involving Sanger and 454 pyrosequencing increases genomic resources to aid in the management of reproduction, disease control and genetic selection in the turbot (Scophthalmus maximus)

Background

Genomic resources for plant and animal species that are under exploitation primarily for human consumption are increasingly important, among other things, for understanding physiological processes and for establishing adequate genetic selection programs. Current available techniques for high-throughput sequencing have been implemented in a number of species, including fish, to obtain a proper description of the transcriptome. The objective of this study was to generate a comprehensive transcriptomic database in turbot, a highly priced farmed fish species in Europe, with potential expansion to other areas of the world, for which there are unsolved production bottlenecks, to understand better reproductive- and immune-related functions. This information is essential to implement marker assisted selection programs useful for the turbot industry.

Results

Expressed sequence tags were generated by Sanger sequencing of cDNA libraries from different immune-related tissues after several parasitic challenges. The resulting database (“Turbot 2 database”) was enlarged with sequences generated from a 454 sequencing run of brain-hypophysis-gonadal axis-derived RNA obtained from turbot at different development stages. The assembly of Sanger and 454 sequences generated 52,427 consensus sequences (“Turbot 3 database”), of which 23,661 were successfully annotated. A total of 1,410 sequences were confirmed to be related to reproduction and key genes involved in sex differentiation and maturation were identified for the first time in turbot (AR, AMH, SRY-related genes, CYP19A, ZPGs, STAR FSHR, etc.). Similarly, 2,241 sequences were related to the immune system and several novel key immune genes were identified (BCL, TRAF, NCK, CD28 and TOLLIP, among others). The number of genes of many relevant reproduction- and immune-related pathways present in the database was 50–90% of the total gene count of each pathway. In addition, 1,237 microsatellites and 7,362 single nucleotide polymorphisms (SNPs) were also compiled. Further, 2,976 putative natural antisense transcripts (NATs) including microRNAs were also identified.

Conclusions

The combined sequencing strategies employed here significantly increased the turbot genomic resources available, including 34,400 novel sequences. The generated database contains a larger number of genes relevant for reproduction- and immune-associated studies, with an excellent coverage of most genes present in many relevant physiological pathways. This database also allowed the identification of many microsatellites and SNP markers that will be very useful for population and genome screening and a valuable aid in marker assisted selection programs.