The miRNA biogenesis in marine bivalves

The Evolution and Characterization of the RNA Interference Pathways in Lophotrochozoa

In animals, three main RNA interference mechanisms have been described so far, which respectively maturate three types of small noncoding RNAs (sncRNAs): miRNAs, piRNAs, and endo-siRNAs. The diversification of these mechanisms is deeply linked with the evolution of the Argonaute gene superfamily since each type of sncRNA is typically loaded by a specific Argonaute homolog. Moreover, other protein families play pivotal roles in the maturation of sncRNAs, like the DICER ribonuclease family, whose DICER1 and DICER2 paralogs maturate respectively miRNAs and endo-siRNAs. Within Metazoa, the distribution of these families has been only studied in major groups, and there are very few data for clades like Lophotrochozoa. Thus, we here inferred the evolutionary history of the animal Argonaute and DICER families including 43 lophotrochozoan species. Phylogenetic analyses along with newly sequenced sncRNA libraries suggested that in all Trochozoa, the proteins related to the endo-siRNA pathway have been lost, a part of them in some phyla (i.e. Nemertea, Bryozoa, Entoprocta), while all of them in all the others. On the contrary, early diverging phyla, Platyhelminthes and Syndermata, showed a complete endo-siRNA pathway. On the other hand, miRNAs were revealed the most conserved and ubiquitous mechanism of the metazoan RNA interference machinery, confirming their pivotal role in animal cell regulation.

Comparative Analysis of Bivalve and Sea Urchin Genetics and Development: Investigating the Dichotomy in Bilateria

This comprehensive review presents a comparative analysis of early embryogenesis in Protostomia and Deuterostomia, the first of which exhibit a mosaic pattern of development, where cells are fated deterministically, while Deuterostomia display a regulatory pattern of development, where the fate of cells is indeterminate. Despite these fundamental differences, there are common transcriptional mechanisms that underline their evolutionary linkages, particularly in the field of functional genomics. By elucidating both conserved and unique regulatory strategies, this review provides essential insights into the comparative embryology and developmental dynamics of these groups. The objective of this review is to clarify the shared and distinctive characteristics of transcriptional regulatory mechanisms. This will contribute to the extensive areas of functional genomics, evolutionary biology and developmental biology, and possibly lay the foundation for future research and discussion on this seminal topic.

Synergistic modulation of neuroendocrine-inflammation pathway by microRNAs facilitates intertidal adaptation of molluscs

Neuroendocrine-immune system is an evolution-conserved regulatory network in maintaining the homeostasis of animals. While knowledge on the roles of neuroendocrine-immune system in the disease and stress responses of organisms is growing, the ecological roles of neuroendocrine-immune system, especially how it shapes the unique lifestyle of organisms remain insufficiently investigated. As an endemic and dominant mollusc in intertidal region, oysters have evolved with a primitive neuroendocrine-immune system and with a sessile lifestyle. Recently, a novel neuroendocrine-immune pathway, Ca2+/calmodulin (CaM)-nitrite oxide synthase (NOS)/nitrite oxide (NO)-tumor necrosis factor (TNF) pathway, is identified in oysters and found altered dynamically during aerial exposure, one common but challenging stresses for intertidal organisms and a decisive factor shaping their habitat. Since the pathway proves fatal in prolonged aerial exposure, we hypothesized that the activation/deactivation of pathway could be strictly modulated in adaptation to the sessile lifestyle of oysters. Here, a synergistic modulation on the Ca2+/CaM-NOS/NO-TNF pathway by four members of miR-92 family and two oyster-specific miRNAs was identified, which further hallmarks the resilience and survival strategy of oysters to aerial exposure. Briefly, these six miRNAs were down-regulating CgCaM24243 post-transcriptionally and deactivating the pathway during the early-stage of stress. However, a robust recession of these miRNAs occurred at the late-stage of stress, resulting in the reactivation of pathway and overwhelming accumulation of cytokines. These results demonstrated a complicated interaction between miRNAs and ancient neuroendocrine-immune system, which facilitates the environmental adaptation of intertidal oysters and provides novel insight on the function and evolution of neuroendocrine-immune system in ecological context.

Mussel antiviral transcriptome response and elimination of viral haemorrhagic septicaemia virus (VHSV)

As filter-feeding bivalves, mussels have been traditionally studied as possible vectors of different bacterial or viral pathogens. The absence of a known viral pathogen in these bivalves makes it particularly interesting to study the interaction of the mussel innate immune system with a virus of interest. In the present work, mussels were challenged with viral haemorrhagic septicaemia virus (VHSV), which is a pathogen in several fish species. The viral load was eliminated after 24 h and mussels evidenced antiviral activity towards VHSV, demonstrating that the virus was recognized and eliminated by the immune system of the host and confirming that mussels are not VHSV vectors in the marine environment. The transcriptome activating the antiviral response was studied, revealing the involvement of cytoplasmic viral sensors with the subsequent activation of the JAK-STAT pathway and several downstream antiviral effectors. The inflammatory response was inhibited with the profound downregulation of MyD88, shifting the immune balance towards antiviral functions. High modulation of retrotransposon activity was observed, revealing a mechanism that facilitates the antiviral response and that had not been previously observed in these species. The expression of several inhibitors of apoptosis and apoptosis-promoting genes was modulated, although clear inhibition of apoptosis in bivalves after severe viral infection and subsequent disease was not observed in this study. Finally, the modulated expression of several long noncoding RNAs that were correlated with genes involved in the immune response was detected.

The Expression of miRNAs Involved in Long-Term Memory Formation in the CNS of the Mollusk Helix lucorum

Mollusks are unique animals with a relatively simple central nervous system (CNS) containing giant neurons with identified functions. With such simple CNS, mollusks yet display sufficiently complex behavior, thus ideal for various studies of behavioral processes, including long-term memory (LTM) formation. For our research, we use the formation of the fear avoidance reflex in the terrestrial mollusk Helix lucorum as a learning model. We have shown previously that LTM formation in Helix requires epigenetic modifications of histones leading to both activation and inactivation of the specific genes. It is known that microRNAs (miRNAs) negatively regulate the expression of genes; however, the role of miRNAs in behavioral regulation has been poorly investigated. Currently, there is no miRNAs sequencing data being published on Helix lucorum, which makes it impossible to investigate the role of miRNAs in the memory formation of this mollusk. In this study, we have performed sequencing and comparative bioinformatics analysis of the miRNAs from the CNS of Helix lucorum. We have identified 95 different microRNAs, including microRNAs belonging to the MIR-9, MIR-10, MIR-22, MIR-124, MIR-137, and MIR-153 families, known to be involved in various CNS processes of vertebrates and other species, particularly, in the fear behavior and LTM. We have shown that in the CNS of Helix lucorum MIR-10 family (26 miRNAs) is the most representative one, including Hlu-Mir-10-S5-5p and Hlu-Mir-10-S9-5p as top hits. Moreover, we have shown the involvement of the MIR-10 family in LTM formation in Helix. The expression of 17 representatives of MIR-10 differentially changes during different periods of LTM consolidation in the CNS of Helix. In addition, using comparative analysis of microRNA expression upon learning in normal snails and snails with deficient learning abilities with dysfunction of the serotonergic system, we identified a number of microRNAs from several families, including MIR-10, which expression changes only in normal animals. The obtained data can be used for further fundamental and applied behavioral research.

Small RNA analysis of Perna viridis after exposure to Prorocentrum lima, a DSP toxins-producing dinoflagellate

Diarrheic shellfish poisoning toxins (DSP toxins) are a set of the most important phycotoxins produced by some dinoflagellates. Studies have shown that DSP toxins have various toxicities such as genotoxicity, cytotoxicity, and immunotoxicity to bivalve mollusks. However, these toxicities appear decreasing with exposure time and concentration of DSP toxins. The underlying mechanism involved remains unclear. In this study, small RNA sequencing was performed in the digestive gland of the mussel Perna viridis after exposure to DSP toxins-producing dinoflagellate Prorocentrum lima for different time periods. The potential roles of miRNAs in response and detoxification to DSP toxins in the mussel were analyzed. Small RNA sequencing of 12 samples from 72 individuals was conducted by BGISEQ-500. A total of 123 mature miRNAs were identified, including 90 conserved miRNAs and 33 potential novel miRNAs. After exposure to P. lima, multiple important miRNAs displayed some alterations. Further miRNA target prediction revealed some important genes involved in cytoskeleton, apoptosis, complement system and immune stress. qPCR demonstrated that miR-71_5, miR-750_1 and novel_mir4 were significantly up-regulated at 6 h after exposure to P. lima, while miR-100_2 was significantly down-regulated after 96 h of exposure. Accordingly, putative target genes of these differentially expressed miRNAs experienced some changes. After 6 h of DSP toxins exposure, NHLRC2 and C1q-like were significantly down-regulated. After 96 h of DSP toxins exposure, NHLRC2 was significantly up-regulated. It is reasonable to speculate that the mussel P. viridis might respond to DSP toxins through miR-750_1, novel_mir4 and miR-71_5 regulating the expression of relevant target genes involved in apoptosis, cytoskeleton, and immune response, etc. This study might provide new clues to uncover the toxic response of bivalve to DSP toxins and lay a foundation for revealing the roles of miRNAs in the environmental adaptation in shellfish.

Integrated analysis of miRNAome and transcriptome reveals miRNA-mRNA network regulation in Vibrio alginolyticus infected thick shell mussel Mytilus coruscus

The thick shell mussel Mytilus coruscus has developed into a model species for studying the interaction between molluscs and environmental stimuli. Herein, integrated analysis of miRNAome and transcriptome was performed to reveal miRNA-mRNA network regulation in Vibrio alginolyticus infected M. coruscus. There have detected some histological abnormalities in digestive gland and gills of V. alginolyticus challenged mussels, ascertaining the effective irritation by the present bacterial strain. A total of 265 novel miRNAs were finally predicted, of which 26 were differentially expressed miRNAs (DEMs). Additionally, 667 differentially expressed genes (DEGs) were detected, which may be potentially associated with innate immune response to V. alginolyticus infection. A regulatory network linked to 22 important pathways and 16 DEMs and 34 OGs was constructed. Some traditional immune-related signaling pathways such as toll-like receptor signaling pathway (TLR) signaling pathway, transforming growth factor-beta (TGF-beta) signaling pathway, peroxisome, phagosome, lysosome, mammalian target of rapamyoin (mTOR) signaling pathway were linked to specific miRNAs and genes in this network. Further, interactional relationship between certain miRNAs and TLR pathway was dissected, which the results predicted that a number of TLRs and TLR-associated signaling genes including TLR1, TLR2, TLR4, TLR6, IRAK1, TRAF6, MAPK, and IL-17 were negatively regulated by novel_miR_11, novel_miR_145, novel_miR_196, novel_miR_5, novel_miR_163 and novel_miR_217 in the TLR pathway. Additionally, interactional relationship between novel_miR_145 and TLR2 was validated by laboratory experiment. The integrated analysis of mRNA and microRNA deep sequencing data exhibited a sophisticated miRNA-mRNA regulation network in M. coruscus in response to V. alginolyticus challenge, which shed a new light on the underlying mechanism of molluscan confronting bacterial infection.

Discovery and functional understanding of MiRNAs in molluscs: a genome-wide profiling approach

Small non-coding RNAs play a pivotal role in gene regulation, repression of transposable element and viral activity in various organisms. Among the various categories of these small non-coding RNAs, microRNAs (miRNAs) guide post-translational gene regulation in cellular development, proliferation, apoptosis, oncogenesis, and differentiation. Here, we performed a genome-wide computational prediction of miRNAs to improve the understanding of miRNA observation and function in molluscs. As an initial step, hundreds of conserved miRNAs were predicted in 35 species of molluscs through genome scanning. Afterwards, the miRNAs' population, isoforms, organization, and function were characterized in detail. Furthermore, the key miRNA biogenesis factors, including AGO2, DGCR8, DICER, DROSHA, TRABP2, RAN, and XPO5, were elucidated based on homologue sequence searching. We also summarized the miRNAs' function in biomineralization, immune and stress response, as well as growth and development in molluscs. Because miRNAs play a vital role in various lifeforms, this study will provide insight into miRNA biogenesis and function in molluscs, as well as other invertebrates.

Stimulation of Mytilus galloprovincialis Hemocytes With Different Immune Challenges Induces Differential Transcriptomic, miRNomic, and Functional Responses

Mediterranean mussels (Mytilus galloprovincialis) are marine bivalve molluscs with high resilience to biotic and abiotic stress. This resilience is one of the reasons why this species is such an interesting model for studying processes such as the immune response. In this work, we stimulated mussel hemocytes with poly I:C, β-glucans, and LPS and then sequenced hemocyte mRNAs (transcriptome) and microRNAs (miRNome) to investigate the molecular basis of the innate immune responses against these pathogen-associated molecular patterns (PAMPs). An immune transcriptome comprising 219,765 transcripts and an overview of the mussel miRNome based on 5,175,567 non-redundant miRNA reads were obtained. The expression analyses showed opposite results in the transcriptome and miRNome; LPS was the stimulus that triggered the highest transcriptomic response, with 648 differentially expressed genes (DEGs), while poly I:C was the stimulus that triggered the highest miRNA response, with 240 DE miRNAs. Our results reveal a powerful immune response to LPS as well as activation of certain immunometabolism- and ageing/senescence-related processes in response to all the immune challenges. Poly I:C exhibited powerful stimulating properties in mussels, since it triggered the highest miRNomic response and modulated important genes related to energy demand; these effects could be related to the stronger activation of these hemocytes (increased phagocytosis, increased NO synthesis, and increased velocity and accumulated distance). The transcriptome results suggest that after LPS stimulation, pathogen recognition, homeostasis and cell survival processes were activated, and phagocytosis was induced by LPS. β-glucans elicited a response related to cholesterol metabolism, which is important during the immune response, and it was the only stimulus that induced the synthesis of ROS. These results suggest a specific and distinct response of hemocytes to each stimulus from a transcriptomic, miRNomic, and functional point of view.

The functional roles of the non-coding RNAs in molluscs

This review focus on the current knowledge of non-coding RNAs (ncRNAs) in molluscs. In this review, we provide an overview of long ncRNAs (lncRNA), microRNAs (miRNA) and piwi-interacting RNAs (piRNA), followed by evidence for the regulation of ncRNAs in variety of biological process in molluscs, including development, biomineralization and innate immune response. This review advances our understanding on the roles of ncRNAs in molluscs and suggest the future direction to fully understand the epigenetic regulatory network of molluscs.

A species-specific miRNA participates in biomineralization by targeting CDS regions of Prisilkin-39 and ACCBP in Pinctada fucata

Biomineralization is a sophisticated biological process precisely regulated by multiple molecules and pathways. Accumulating miRNAs have been identified in invertebrates but their functions in biomineralization are poorly studied. Here, an oyster species-specific miRNA, novel_miR_1 was found to regulate biomineralization in Pinctada fucata. Target prediction showed that novel_miR_1 could target Prisilkin-39 and ACCBP by binding to their coding sequences (CDS). Tissue distribution analysis revealed that the expression level of novel_miR_1 was highest in the mantle, which was a key tissue participating in biomineralization. Gain-of-function assay in vivo showed that biomineralization-related genes including Prisilkin-39 and ACCBP were down-regulated and shell inner surfaces of both prismatic and nacreous layer were disrupted after the over-expression of novel_miR_1, indicating its dual roles in biomineralization. Furthermore, the shell notching results indicated that novel_miR_1 was involved in shell regeneration. Dual-luciferase reporter assay in vitro demonstrated that novel_miR_1 directly suppressed Prisilkin-39 and ACCBP genes by binding to the CDS regions. Taken together, these results suggest that novel_miR_1 is a direct negative regulator to Prisilkin-39 and ACCBP and plays an indispensable and important role in biomineralization in both prismatic and nacreous layer of P. fucata.

Transgenerational acclimation to changes in ocean acidification in marine invertebrates

The rapid pace of increasing oceanic acidity poses a major threat to the fitness of the marine ecosystem, as well as the buffering capacity of the oceans. Disruption in chemical equilibrium in the ocean leads to decreased carbonate ion precipitation, resulting in calcium carbonate saturation. If these trends continue, calcifying invertebrates will experience difficultly maintaining their calcium carbonate exoskeleton and shells. Because malfunction of exoskeleton formation by calcifiers in response to ocean acidification (OA) will have non-canonical biological cascading results in the marine ecosystem, many studies have investigated the direct and indirect consequences of OA on ecosystem- and physiology-related traits of marine invertebrates. Considering that evolutionary adaptation to OA depends on the duration of OA effects, long-term exposure to OA stress over multi-generations may result in adaptive mechanisms that increase the potential fitness of marine invertebrates in response to OA. Transgenerational studies have the potential to elucidate the roles of acclimation, carryover effects, and evolutionary adaptation within and over generations in response to OA. In particular, understanding mechanisms of transgenerational responses (e.g., antioxidant responses, metabolic changes, epigenetic reprogramming) to changes in OA will enhance our understanding of marine invertebrate in response to rapid climate change.

Genomics and immunity of the Mediterranean mussel Mytilus galloprovincialis in a changing environment

The Mediterranean mussel (Mytilus galloprovincialis) is a marine invasive species cultured all over the world. Mussels are an appreciated resource in local aquaculture enterprises because of their robust production and resilience that translates into a reliable economic value. So far, no massive mortalities have been reported in natural or cultured populations of this species. In the last years, the knowledge about its immune system has greatly improved but there are still many questions to be answered. One of them is why mussels, with their high filtering activity, are able to be exposed to a high number of potential pathogens without getting infected and without developing an elevated inflammatory response. The sequencing of the mussel genome has revealed a very complex organization with high heterozygosity, abundance of repetitive sequences and extreme intraspecific sequence diversity among individuals, mainly in immune related genes. Among those genes, antimicrobial peptides are the most expressed gene families in mussels, highly polymorphic and with antimicrobial effect against molluscs pathogens, but also against pathogens of lower vertebrates and humans. The combination of a complex genome with the adaptation of mussel immune system to a changing environment could explain this high variability, not only in immune-related genes, but also in the functional response among individuals sampled in the same location and date.

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.

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.

Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes

RNA interference (RNAi)-related pathways target viruses and transposable element (TE) transcripts in plants, fungi, and ecdysozoans (nematodes and arthropods), giving protection against infection and transmission. In each case, this produces abundant TE and virus-derived 20-30nt small RNAs, which provide a characteristic signature of RNAi-mediated defence. The broad phylogenetic distribution of the Argonaute and Dicer-family genes that mediate these pathways suggests that defensive RNAi is ancient, and probably shared by most animal (metazoan) phyla. Indeed, while vertebrates had been thought an exception, it has recently been argued that mammals also possess an antiviral RNAi pathway, although its immunological relevance is currently uncertain and the viral small RNAs (viRNAs) are not easily detectable. Here we use a metagenomic approach to test for the presence of viRNAs in five species from divergent animal phyla (Porifera, Cnidaria, Echinodermata, Mollusca, and Annelida), and in a brown alga-which represents an independent origin of multicellularity from plants, fungi, and animals. We use metagenomic RNA sequencing to identify around 80 virus-like contigs in these lineages, and small RNA sequencing to identify viRNAs derived from those viruses. We identified 21U small RNAs derived from an RNA virus in the brown alga, reminiscent of plant and fungal viRNAs, despite the deep divergence between these lineages. However, contrary to our expectations, we were unable to identify canonical (i.e. Drosophila- or nematode-like) viRNAs in any of the animals, despite the widespread presence of abundant micro-RNAs, and somatic transposon-derived piwi-interacting RNAs. We did identify a distinctive group of small RNAs derived from RNA viruses in the mollusc. However, unlike ecdysozoan viRNAs, these had a piRNA-like length distribution but lacked key signatures of piRNA biogenesis. We also identified primary piRNAs derived from putatively endogenous copies of DNA viruses in the cnidarian and the echinoderm, and an endogenous RNA virus in the mollusc. The absence of canonical virus-derived small RNAs from our samples may suggest that the majority of animal phyla lack an antiviral RNAi response. Alternatively, these phyla could possess an antiviral RNAi response resembling that reported for vertebrates, with cryptic viRNAs not detectable through simple metagenomic sequencing of wild-type individuals. In either case, our findings show that the antiviral RNAi responses of arthropods and nematodes, which are highly divergent from each other and from that of plants and fungi, are also highly diverged from the most likely ancestral metazoan state.

Genetic and molecular basis of the immune system in the brachiopod Lingula anatina

The extension of comparative immunology to non-model systems, such as mollusks and annelids, has revealed an unexpected diversity in the complement of immune receptors and effectors among evolutionary lineages. However, several lophotrochozoan phyla remain unexplored mainly due to the lack of genomic resources. The increasing accessibility of high-throughput sequencing technologies offers unique opportunities for extending genome-wide studies to non-model systems. As a result, the genome-based study of the immune system in brachiopods allows a better understanding of the alternative survival strategies developed by these immunologically neglected phyla. Here we present a detailed overview of the molecular components of the immune system identified in the genome of the brachiopod Lingula anatina. Our findings reveal conserved intracellular signaling pathways as well as unique strategies for pathogen detection and killing in brachiopods.

A potential microRNA regulation of immune-related genes in invertebrate haemocytes

Bivalve mollusks have been employed as sentinel organisms in environmental health programs due to their sedentary lifestyle, filter-feeding behavior and their ability to accumulate pathogens or toxin molecules inside tissues. Endocrine disrupting chemicals (EDCs) can be up taken and bioaccumulated, and due to sensibility of mollusks to these EDCs, being able to cause immune alterations. Recently, microRNAs (miRNAs) were shown to be involved in modulation and buffering developmental processes against the effects of environmental alterations and pathogenic microorganisms. Moreover, it is suggested that this miRNAs are incorporated into the estrogen-controlled immune network, regulating mechanism of immune gene expression at the posttranscriptional level, modulating immune responses as phagocytosis, redox reaction and apoptosis in bivalve haemocytes. Thus, miRNAs can be used as biomarkers that specifically elucidate immunotoxic effects caused by exogenous biotic or abiotic factors, and can act as useful tools in integrated monitoring environmental health programs. In this review, we aim to describe the investigations that have been carried out on miRNAs in bivalve mollusks, especially those associated with immune responses against infectious agents and xenobiotic exposure.

miRNAome expression profiles in the gonads of adult Melopsittacus undulatus

The budgerigar (Melopsittacus undulatus) is one of the most widely studied parrot species, serving as an excellent animal model for behavior and neuroscience research. Until recently, it was unknown how sexual differences in the behavior, physiology, and development of organisms are regulated by differential gene expression. MicroRNAs (miRNAs) are endogenous short non-coding RNA molecules that can post-transcriptionally regulate gene expression and play a critical role in gonadal differentiation as well as early development of animals. However, very little is known about the role gonadal miRNAs play in the early development of birds. Research on the sex-biased expression of miRNAs in avian gonads are limited, and little is known about M. undulatus. In the current study, we sequenced two small non-coding RNA libraries made from the gonads of adult male and female budgerigars using Illumina paired-end sequencing technology. We obtained 254 known and 141 novel miRNAs, and randomly validated five miRNAs. Of these, three miRNAs were differentially expressed miRNAs and 18 miRNAs involved in sexual differentiation as determined by functional analysis with GO annotation and KEGG pathway analysis. In conclusion, this work is the first report of sex-biased miRNAs expression in the budgerigar, and provides additional sequences to the avian miRNAome database which will foster further functional genomic research.

Epigenetic considerations in aquaculture

Epigenetics has attracted considerable attention with respect to its potential value in many areas of agricultural production, particularly under conditions where the environment can be manipulated or natural variation exists. Here we introduce key concepts and definitions of epigenetic mechanisms, including DNA methylation, histone modifications and non-coding RNA, review the current understanding of epigenetics in both fish and shellfish, and propose key areas of aquaculture where epigenetics could be applied. The first key area is environmental manipulation, where the intention is to induce an ‘epigenetic memory’ either within or between generations to produce a desired phenotype. The second key area is epigenetic selection, which, alone or combined with genetic selection, may increase the reliability of producing animals with desired phenotypes. Based on aspects of life history and husbandry practices in aquaculture species, the application of epigenetic knowledge could significantly affect the productivity and sustainability of aquaculture practices. Conversely, clarifying the role of epigenetic mechanisms in aquaculture species may upend traditional assumptions about selection practices. Ultimately, there are still many unanswered questions regarding how epigenetic mechanisms might be leveraged in aquaculture.

Bioinformatics Study of Structural Patterns in Plant MicroRNA Precursors

According to the RNA world theory, RNAs which stored genetic information and catalyzed chemical reactions had their contribution in the formation of current living organisms. In recent years, researchers studied this molecule diversity, i.a. focusing on small non-coding regulatory RNAs. Among them, of particular interest is evolutionarily ancient, 19–24 nt molecule of microRNA (miRNA). It has been already recognized as a regulator of gene expression in eukaryotes. In plants, miRNA plays a key role in the response to stress conditions and it participates in the process of growth and development. MicroRNAs originate from primary transcripts (pri-miRNA) encoded in the nuclear genome. They are processed from single-stranded stem-loop RNA precursors containing hairpin structures. While the mechanism of mature miRNA production in animals is better understood, its biogenesis in plants remains less clear. Herein, we present the results of bioinformatics analysis aimed at discovering how plant microRNAs are recognized within their precursors (pre-miRNAs). The study has been focused on sequential and structural motif identification in the neighbourhood of microRNA.