Bioinformatics analysis of hemocyte miRNAs of scallop Chlamys farreri against acute viral necrobiotic virus (AVNV)

miRNA-seq provides novel insight into the response to hyper- and hypo- salinity acclimation in Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is a significant bivalve species with economic importance. It primarily inhabits the estuarine intertidal zones in southern China, making it susceptible to salinity fluctuations. Consequently, investigating the molecular mechanisms governing salinity regulation in C. hongkongensis is essential. In this study, we conducted miRNA-seq on C. hongkongensis to compare miRNA expression differences under varying salinities (5‰, 25‰, and 35‰). The miRNA sequencing revealed 51 known miRNAs and 95 novel miRNAs across nine small RNA libraries (S5, S25, and S35). Among these miRNAs, we identified 6 down-regulated differentially expressed (DE) miRNAs in response to hypo-salinity stress (5‰), while 1 up-regulated DE miRNA and 5 down-regulated DE miRNAs were associated with hyper-salinity stress (35‰). Additionally, we predicted 931 and 768 potential target genes for hypo- and hyper-salinity stress, respectively. Functional gene annotation indicated that the target genes under hypo-salinity stress were linked to vesicle-mediated transport and metal ion binding. Conversely, those under hyper-salinity stress were primarily involved in signal transduction and metabolic processes. These findings have provided insights into the regulatory role of miRNAs, their potential target genes and associated pathways in oyster hypo- and hyper-salinity stress, which establish a foundation for future studies on the roles of miRNAs in salinity acclimation mechanisms in C. hongkongensis.

Insights into ADAR gene complement, expression patterns, and RNA editing landscape in Chlamys farreri

Adenosine Deaminases Acting on RNA (ADARs) are evolutionarily conserved enzymes known to convert adenosine to inosine in double-stranded RNAs and participate in host-virus interactions. Conducting a meta-analysis of available transcriptome data, we identified and characterised eight ADAR transcripts in Chlamys farreri, a farmed marine scallop susceptible to Acute viral necrosis virus (AVNV) infections and mortality outbreaks. Accordingly, we identified six ADAR genes in the Zhikong scallop genome, revised previous gene annotations, and traced alternative splicing variants. In detail, each ADAR gene encodes a unique combination of functional domains, always including the Adenosine deaminase domain, RNA binding domains and, in one case, two copies of a Z-DNA binding domain. After phylogenetic analysis, five C. farreri ADARs clustered in the ADAR1 clade along with sequences from diverse animal phyla. Gene expression analysis indicated CF051320 as the most expressed ADAR, especially in the eye and male gonad. The other four ADAR1 genes and one ADAR2 gene exhibited variable expression levels, with CF105370 and CF051320 significantly increasing during early scallop development. ADAR-mediated single-base editing, evaluated across adult C. farreri tissues and developmental stages, was mainly detectable in intergenic regions (83 % and 85 %, respectively). Overall, the expression patterns of the six ADAR genes together with the editing and hyper-editing values computed on scallops RNA-seq samples support the adaptive value of ADAR1-mediated editing, particularly in the pre-settling larval stages.

MicroRNA transcriptome analysis reveals the immune regulatory mechanism of Crassostrea hongkongesis against Vibrio harveyi infection

MicroRNAs (miRNAs) are small non-coding RNA molecules that modulate target-genes expression and play crucial roles in post-transcriptional regulation and immune system regulation. The Hong Kong oyster (Crassostrea hongkongesis), as the main marine aquaculture shellfish in the South China Sea, not only has high economic and ecological value, but also is an ideal model for conducting research on pathogen host interaction. Vibrio harveyi, a Gram negative luminescent marine bacterium, is widely distributed in coastal water environments and can cause large-scale death of C. hongkongesis. However, little in formation is available on the immune regulatory mechanisms of C. hongkongesis infected with V. harveyi. Therefore, we performed microRNA transcriptome analysis for elucidating the immunoregulation mechanism of C. hongkongesis infected with V. harveyi. The results show that a total of 308468208 clean reads and 288371159 clean tags were obtained. 222 differentially expressed miRNAs were identified. A total of 388 target genes that were differentially expressed and negatively correlated with miRNA expression were predicted by 222 DEmiRs. GO enrichment analysis of 388 DETGs showed that they were mainly enriched in the immune-related term of membrane-bounded vesicle, endocytic vesicle lumen, antigen processing and presentation of exogenous peptide antigen via MHC class I, antigen processing and presentation of peptide antigen via MHC class I, and other immune-related term. KEGG enrichment analysis showed that DETGs were mainly enriched in the Complement and coagulation cascades, Herpes simplex virus 1 infection, Bacterial invasion of epithelial cells, Antigen processing and presentation and NOD-like receptor signaling pathway. The 16 key DEmiRs and their target genes form a regulatory network for seven immune-related pathways. These results suggest that V. harveyi infection induces a complex miRNA response with wide-ranging effects on immune gene expression in the C. hongkongesis. This study explored the immune response of C. hongkongesis to V. harveyi infection at the level of miRNAs, which provides new ideas for the healthy culture and selective breeding of C. hongkongesis.

Transcriptional landscape of small non-coding RNAs reveals diversity of categories and functions in molluscs

Small non-coding RNAs (sncRNAs) are non-coding RNA molecules that play various roles in metazoans. Among the sncRNAs, microRNAs (miRNAs) guide post-translational gene regulation during cellular development, proliferation, apoptosis, and differentiation, while PIWI-interacting RNAs (piRNAs) suppress transposon activity to safeguard the genome from detrimental insertion mutagenesis. While an increasing number of piRNAs are being identified in the soma and germlines of various organisms, they are scarcely reported in molluscs. To unravel the small RNA (sRNA) expression patterns and genomic function in molluscs, we generated a comprehensive sRNA dataset by sRNA sequencing (sRNA-seq) of eight mollusc species. Abundant miRNAs were identified and characterized in all investigated molluscs, and ubiquitous piRNAs were discovered in both somatic and gonadal tissues in six of the investigated molluscs, which are more closely associated with transposon silencing. Tens of piRNA clusters were also identified based on the genomic mapping results, which varied among different tissues and species. Our dataset serves as important reference data for future genomic and genetic studies on sRNAs in these molluscs and related species, especially in elucidating the ancestral state of piRNAs in bilaterians.

Immunomodulatory effects of one novel microRNA miR-63 in pearl oyster Pinctada fucata martensii

Novel microRNA miR-63 (novel-miR-63) from pearl oyster Pinctada fucata martensii (Pm-novel-miR-63) is a species-specific miRNA. Our previous research has shown that the expression of Pm-novel-miR-63 was significantly downregulated at 24 h after nucleus transplantation. In this study, we analyzed the function and regulatory role of Pm-novel-miR-63 in the immune response of pearl oysters. The results showed that Pm-novel-miR-63 expression increased after the stimulation of pathogen associated molecular patterns at 6-12 h, and the activity of immune and antioxidant enzymes in the serum decreased after Pm-novel-miR-63 overexpression. Transcriptome analysis revealed that Pm-novel-miR-63 participated in regulating transplantation immunity through the Notch and mRNA surveillance signaling pathways. Target prediction and dual luciferase analysis revealed that Pm-GDP-FucTP, Pm-CysLTR2, and Pm-RLR were the target genes of Pm-novel-miR-63. These results suggested that Pm-novel-miR-63 participated in regulating the immune response in pearl oysters and can serve as a new interference target to reasonably control excessive immune rejection in pearl culture.

Identification and comparative analysis of miRNA transcriptomes after allograft and xenograft transplantation in Pinctada fucata martensii

Effective immune regulation after transplantation during pearl production is crucial for the cultivation of high-quality pearls. MicroRNAs (miRNAs) play an important role in a variety of physiological processes. To understand the regulatory rules of miRNAs after transplantation in Pinctada funcata martensii, we constructed 13 miRNA transcriptomes, including the control group (Con), allograft (Al), and xenograft (Xe) transplantation at six time points (6, 12, and 24 h and 3, 6, and 12 days), in which the xenografted mantle tissue was from Pinctada maxima. We identified 159 differentially expressed miRNAs (DEMs) and found that these DEMs showed high expression at 12 h, 24 h, and 3 days after transplantation. A total of 130 DEMs, such as Let-7, were present in the Al and Xe groups; miR-34 and 16 other DEMs were specifically present in the Al group; miR-216b and 13 other DEMs were specifically present in the Xe group. Compared with the Con group, the target genes of DEMs in the Al group were significantly enriched in protein complex, cytoskeleton, and macromolecular complex, and the Xe group was significantly enriched in ribonucleoside metabolic process, nucleoside binding, and cell division. Compared with the Al group, the target genes in the Xe group were significantly enriched in response to DNA damage stimulation. Overall, multiple pathways associated with cellular activity were enriched in higher numbers of genes in the Xe group than in the Al group. These findings enriched the information on immune regulatory mechanisms at the expression level of miRNAs in P. f. martensii after transplantation.

Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture

The European flat oyster (Ostrea edulis L.) is a native bivalve of the European coasts. Harvest of this species has declined during the last decades because of the appearance of two parasites that have led to the collapse of the stocks and the loss of the natural oyster beds. O. edulis has been the subject of numerous studies in population genetics and on the detection of the parasites Bonamia ostreae and Marteilia refringens. These studies investigated immune responses to these parasites at the molecular and cellular levels. Several genetic improvement programs have been initiated especially for parasite resistance. Within the framework of a European project (PERLE 2) that aims to produce genetic lines of O. edulis with hardiness traits (growth, survival, resistance) for the purpose of repopulating natural oyster beds in Brittany and reviving the culture of this species in the foreshore, obtaining a reference genome becomes essential as done recently in many bivalve species of aquaculture interest. Here, we present a chromosome-level genome assembly and annotation for the European flat oyster, generated by combining PacBio, Illumina, 10X linked, and Hi-C sequencing. The finished assembly is 887.2 Mb with a scaffold-N50 of 97.1 Mb scaffolded on the expected 10 pseudochromosomes. Annotation of the genome revealed the presence of 35,962 protein-coding genes. We analyzed in detail the transposable element (TE) diversity in the flat oyster genome, highlighted some specificities in tRNA and miRNA composition, and provided the first insight into the molecular response of O. edulis to M. refringens. This genome provides a reference for genomic studies on O. edulis to better understand its basic physiology and as a useful resource for genetic breeding in support of aquaculture and natural reef restoration.

MicroRNA expression signature in the striated and smooth adductor muscles of Yesso scallop Patinopecten yessoensis

Increasing evidences point to the potential role of microRNAs (miRNAs) in muscle growth and development in animals. However, knowledge on the identity of miRNAs and their targets in molluscs remains largely unknown. Scallops have one large adductor muscle, composed of fast (striated) and slow (smooth) muscle types, which display great differences in muscle fibers, meat quality, cell types and molecular components. In the present study, we performed a comprehensive investigation of miRNA transcriptomes in fast and slow adductor muscles of Yesso scallop Patinopecten yessoensis. As a result, 47 differentially expressed miRNAs representing ten miRNA families were identified between the striated and smooth adductor muscles. The KEGG enrichment analysis of their target genes were mainly associated with amino acid metabolism, energy metabolism and glycan biosynthesis. The target genes of miR-133 and miR-71 were validated by the dual-luciferase reporter assays and miRNA antagomir treatment in vivo. The identification and functional validation of these different miRNAs in scallops will greatly help our understanding of miRNA regulatory mechanism that achieves the unique muscle phenotypes in scallops. The present findings provide the direct evidences for muscle-specific miRNAs involved in muscle growth and differentiation in molluscs.

Molecular allocation of PC4s provides implications for deciphering thermal response in Zhikong scallop (Chlamys farreri)

The increasing sea temperature caused by global warming has led to serious death of Zhikong scallop (Chlamys farreri) and improving its thermal tolerance has become an active research area in scallop aquaculture industry. Gene transcriptional coactivator p15 (PC4) plays pivotally multi-faced roles in most vertebrates and some invertebrates, but the systematic identification and characterization of PC4 genes have less been reported in scallops. In this study, 15 PC4 genes (CfPC4s) were identified in Zhikong scallop through whole-genome scanning, including two pairs of tandem duplicate genes located in the same scaffold (CF-19495.9 and CF-19495.10, CF-6819.1 and CF-6819.2). Protein structural and phylogenetic analyses were performed to verify identities and evolutionary relationships of these genes. Spatiotemporal expression patterns were determined at different development stages and in healthy adult tissues, as well as expression regulations in selected tissues (mantles, gills, hemocytes and hearts) after high temperatures challenge (27 °C) with different durations (3 h, 6 h, 12 h, 24 h, 3 d, 6 d, 15 d and 30 d). Spatiotemporal expressions of CfPC4s were ubiquitous but exhibited different patterns, suggesting the functional roles of CfPC4s in all stages of growth and development of the scallop. Expression regulations of CfPC4s and their functional related factors (TFIIA, TFIID, TFIIH and RNAPII) in pre-initiation complex (PIC) in various tissues displayed up- and/or down-regulated responses at different time points, showing time- and/or tissue-dependent expression patterns with function allocation upon different thermal durations. Collectively, this study demonstrated that gene allocation of CfPC4s provided implications for deciphering thermal response in Zhikong scallop and potentially helped in developing strategies for long-term healthy sustainable Zhikong scallop culture.

Paired miRNA and RNA sequencing provides a first insight into molecular defense mechanisms of Scapharca broughtonii during ostreid herpesvirus-1 infection

Ostreid herpesvirus 1 (OsHV-1) infection caused mortalities with relevant economic losses in bivalve aquaculture industry worldwide. Initially described as an oyster pathogen, OsHV-1 can infect other bivalve species, like the blood clam Scapharca broughtonii. However, at present, little is known about the molecular interactions during OsHV-1 infection in the blood clam. We produced paired miRNA and total RNA-seq data to investigate the blood clam transcriptional changes from 0 to 72 h after experimental infection with OsHV-1. High-throughput miRNA sequencing of 24 libraries revealed 580 conserved and 270 new blood clam miRNAs, whereas no genuine miRNA was identified for OsHV-1. Total 88-203 differently expressed miRNAs were identified per time point, mostly up-regulated and mainly targeting metabolic pathways. Most of the blood clam mRNAs, in contrast, were down-regulated up to 60 h post-injection, with the trend analysis revealing the activation of immune genes only when comparing the early and latest stage of infection. Taken together, paired short and long RNA data suggested a miRNA-mediated down-regulation of host metabolic and energetic processes as a possible antiviral strategy during early infection stages, whereas antiviral pathways appeared upregulated only at late infection.

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.

Differential Expression of Long Non-Coding RNA (lncRNA) in Mediterranean Mussel (Mytilus galloprovincialis) Hemocytes under Immune Stimuli

The Mediterranean mussel is one of the most economically relevant bivalve mollusk species in Europe and China. The absence of massive mortalities and their resistance to pathogens affecting other cultured bivalves has been under study in recent years. The transcriptome response of this species to different immune stimuli has been extensively studied, and even the complexity of its genome, which has recently been sequenced, has been suggested as one of the factors contributing to this resistance. However, studies concerning the non-coding RNA profiles remain practically unexplored-especially those corresponding to the lncRNAs. To the best of our knowledge, this is the second characterization and study of lncRNAs in this bivalve species. In this work, we identified the potential repertoire of lncRNAs expressed in mussel hemocytes, and using RNA-Seq we analyzed the lncRNA profile of mussel hemocytes stimulated in vitro with three different immune stimuli: LPS, poly I:C, and β-glucans. Compared to unstimulated hemocytes, LPS induced the highest modulation of lncRNAs, whereas poly I:C and β-glucans induced a similar discrete response. Based on the potential cis-regulatory activity of the lncRNAs, we identified the neighboring protein-coding genes of the regulated lncRNAs to estimate-at least partially-the processes in which they are implicated. After applying correlation analyses, it seems that-especially for LPS-the lncRNAs could participate in the regulation of gene expression, and substantially contribute to the immune response.

Identification and study of differentially expressed miRNAs in aged NAFLD rats based on high-throughput sequencing

INTRODUCTION AND OBJECTIVES

Hepatic microRNA (miR) expression profiles were explored in aged rats with NAFLD, in order to clarify the molecular mechanisms underlying the pathophysiological processes of aging-related NAFLD.

PATIENTS OR MATERIALS AND METHODS

24 aged rats (18-month-old) and 24 young rats (2-month-old) were randomly divided into two subgroups according to diet, control group and NAFLD group. After 8 weeks of administering 45% high-fat diet or normal diet, total hepatic RNA was extracted from liver tissues of the aged rats. Differentially expressed microRNAs (DE-miRs) in aged NAFLD group were detected and screened out using high-throughput sequencing technology. The data were subjected to Gene Ontology functional enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses using a bioinformatics approach. The sequencing results were further verified by RT-qPCR.

RESULTS

Compared with the aged control liver tissues, 6 significantly upregulated miRs (miR-881-3p, miR-871-3p, miR-335, miR-223-3p, miR-155-5p, miR-146b-5p) and 4 significantly downregulated miRs (miR-182, miR-193-3p, miR-31a-5p and miR-96-5p) were identified in the aged NAFLD liver tissues. These DE-miRs were found to be involved in the regulation of cell signaling transduction and metabolism processes, probably affecting signaling pathways relevant to insulin secretion and some senile diseases. RT-qPCR results corroborated the sequencing results and demonstrated that 6 significantly upregulated miRs were not identified in the young group.

CONCLUSIONS

A total of 10 DE-miRs identified in the aged NAFLD rats were involved in some certain insulin secretion and age-related functional pathways, which may serve as novel candidate targets for the diagnosis and treatment of aging-associated NAFLD.

Digging into bivalve miRNAomes: between conservation and innovation

Bivalves are a diverse mollusc group of economic and ecological importance. An evident resilience to pollution, parasites and extreme environments makes some bivalve species important models for studying adaptation and immunity. Despite substantial progress in sequencing projects of bivalves, information on non-coding genes and gene-regulatory aspects is still lacking. Here, we review the current repertoire of bivalve microRNAs (miRNAs), important regulators of gene expression in Metazoa. We exploited available short non-coding RNA (sncRNA) data for Pinctada martensii, Crassostrea gigas, Corbicula fluminea, Tegillarca granosa and Ruditapes philippinarum, and we produced new sncRNA data for two additional bivalves, the Mediterranean mussel Mytilus galloprovincialis and the blood clam Scapharca broughtonii. We found substantial heterogeneity and incorrect annotations of miRNAs; hence, we reannotated conserved miRNA families using recently established criteria for bona fide microRNA annotation. We found 106 miRNA families missing in the previously published bivalve datasets and 89 and 87 miRNA complements were identified in the two additional species. The overall results provide a homogeneous and evolutionarily consistent picture of miRNAs in bivalves and enable future comparative studies. The identification of two bivalve-specific miRNA families sheds further light on the complexity of transcription and its regulation in bivalve molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Stress responses in expressions of microRNAs in mussel Mytilus galloprovincialis exposed to cadmium

MicroRNAs (miRNAs) are known to have complicated functions in aquatic species, but little is known about the role of miRNAs in mollusk species under environmental stress. In this study, we performed small RNA sequencing to characterize the differentially expressed miRNAs in different tissues (whole tissues, digestive glands, gills, and gonads) of blue mussels (Mytilus galloprovincialis) exposed to cadmium (Cd). In summary, 107 known miRNAs and 32 novel miRNAs were significantly (p < 0.01) differentially expressed after Cd exposure. The peak size of miRNAs was 22 nucleotides. Target genes of these differentially expressions of miRNAs related to immune defense, apoptosis, lipid and xenobiotic metabolism showed significant changes under Cd stress. These findings provide the first characterization of miRNAs in mussel M. galloprovincialis and expressions of many target genes in response to Cd stress.

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.

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.

Parallel analysis of miRNAs and mRNAs suggests distinct regulatory networks in Crassostrea gigas infected by Ostreid herpesvirus 1

Background

Since 2008, the aquaculture production of Crassostrea gigas was heavily affected by mass mortalities associated to Ostreid herpesvirus 1 (OsHV-1) microvariants worldwide. Transcriptomic studies revealed the major antiviral pathways of the oyster immune response while other findings suggested that also small non-coding RNAs (sncRNA) such as microRNAs might act as key regulators of the oyster response against OsHV-1. To explore the explicit connection between small non-coding and protein-coding transcripts, we performed paired whole transcriptome analysis of sncRNA and messenger RNA (mRNA) in six oysters selected for different intensities of OsHV-1 infection.

Results

The mRNA profiles of the naturally infected oysters were mostly governed by the transcriptional activity of OsHV-1, with several differentially expressed genes mapping to the interferon, toll, apoptosis, and pro-PO pathways. In contrast, miRNA profiles suggested more complex regulatory mechanisms, with 15 differentially expressed miRNAs (DE-miRNA) pointing to a possible modulation of the host response during OsHV-1 infection. We predicted 68 interactions between DE-miRNAs and oyster 3′-UTRs, but only few of them involved antiviral genes. The sncRNA reads assigned to OsHV-1 rather resembled mRNA degradation products, suggesting the absence of genuine viral miRNAs.

Conclusions

We provided data describing the miRNAome during OsHV-1 infection in C. gigas. This information can be used to understand the role of miRNAs in healthy and diseased oysters, to identify new targets for functional studies and, eventually to disentangle cause and effect relationships during viral infections in marine mollusks.

Identification of candidate microRNAs from Ostreid herpesvirus-1 (OsHV-1) and their potential role in the infection of Pacific oysters (Crassostrea gigas)

Oyster production is an economic activity of great interest worldwide. Recently, oysters have been suffering significant mortalities from OsHV-1infection, which has resulted in substantial economic loses in several countries around the world. Understanding viral pathogenicity mechanisms is of central importance for the establishment of disease control measures. Thus, the present work aimed to identify and characterize miRNAs from OsHV-1 as well as to predict their target transcripts in the virus and the host. OsHV-1 genome was used for the in silico discovery of pre-miRNAs. Subsequently, viral and host target transcripts of the OsHV-1 miRNAs were predicted according to the base pairing interaction between mature miRNAs and mRNA 3' untranslated regions (UTRs). Six unique pre-miRNAs were found in different regions of the viral genome, ranging in length from 85 to 172 nucleotides. A complex network of self-regulation of viral gene expression mediated by the miRNAs was identified. These sequences also seem to have a broad ability to regulate the expression of host immune-related genes, especially those associated with pathogen recognition. Our results suggest that OsHV-1 encodes miRNAs with important functions in the infection process, inducing self-regulation of viral transcripts, as well as affecting the regulation of Pacific oyster transcripts related to immunity. Understanding the molecular basis of host-pathogen interactions can help mitigate the recurrent events of oyster mass mortalities by OsHV-1 observed worldwide.

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.

Profiling of microRNAs and mRNAs in marine mussel Mytilus galloprovincialis

MicroRNAs (miRNAs) are a class of noncoding RNA molecules containing 18-24 nucleotides, and those with conserved structures are able to regulate the expression of eukaryotic genes by inhibition or enhancement of mRNA translation. However, miRNAs of the blue mussel, Mytilus galloprovincialis have not been reported. M. galloprovincialis is a primary species distributed along coastal zones worldwide. To reveal the repertoire of miRNAs in M. galloprovincialis, we constructed small RNA libraries prepared from three different mussels, which were then sequenced by Solexa deep sequencing technology. A total of 32,836,817, 33,359,113 and 33,093,562 clean reads from the tissues of the three M. galloprovincialis were obtained. Based on sequence similarities and hairpin structure predictions, 137 M. galloprovincialis miRNAs (mg-miRNA) were identified. Among the mg-miRNAs, 104 were conserved across species, whereas 33 might be novel and specific for M. galloprovincialis. Some of the mg-miRNAs, such as let-7 and the miR-100 family are playing key roles in many metabolic pathways and are worthy of further study. By performing a whole genome-scale characterization of mg-miRNAs and proposing their potential functions, these results provide a foundation for understanding the biological processes of the blue mussel, M. galloprovincialis.

Investigating the natural resistance of blackfoot p-a%%KERN_ERR%%ua Haliotis iris to abalone viral ganglioneuritis using whole transcriptome analysis

The natural resistance of New Zealand blackfoot p-a%%%%%%%%%%%%%%KERN_ERR%%KERN_ERR%%KERN_ERR%%KERN_ERR%%KERN_ERR%%KERN_ERR%%KERN_ERR%%ua Haliotis iris to infection by haliotid herpesvirus-1 (HaHV-1) and to the disease abalone viral ganglioneuritis was investigated in experimentally challenged p-aua using high throughput RNA-sequencing. HaHV-1-challenged p-aua up-regulated broad classes of genes that contained chitin-binding peritrophin-A domains, which seem to play diverse roles in the p-aua immune response. The p-aua also up-regulated vascular adhesion protein-1 (VAP-1), an important adhesion molecule for lymphocytes, and chitotriosidase-1 (CHIT-1), an immunologically important gene in mammalian immune systems. Moreover, several blood coagulation pathways were dysregulated in the p-aua, possibly indicating viral modulation. We also saw several indications that neurological tissues were specifically affected by HaHV-1, including the dysregulation of beta-1,4-N-acetylgalactosaminyltransferase (B4GALNT), GM2 ganglioside, neuroligin-4 and the Notch signalling pathway. This research may support the development of molecular therapeutics useful to control and/or manage viral outbreaks in abalone culture.

Dual Analysis of Virus-Host Interactions: The Case of Ostreid herpesvirus 1 and the Cupped Oyster Crassostrea gigas

Dual analyses of the interactions between Ostreid herpesvirus 1 (OsHV-1) and the bivalve Crassostrea gigas during infection can unveil events critical to the onset and progression of this viral disease and can provide novel strategies for mitigating and preventing oyster mortality. Among the currently used “omics” technologies, dual transcriptomics (dual RNA-seq) coupled with the analysis of viral DNA in the host tissues has greatly advanced the knowledge of genes and pathways mostly contributing to host defense responses, expression profiles of annotated and unknown OsHV-1 open reading frames (ORFs), and viral genome variability. In addition to dual RNA-seq, proteomics and metabolomics analyses have the potential to add complementary information, needed to understand how a malacoherpesvirus can redirect and exploit the vital processes of its host. This review explores our current knowledge of “omics” technologies in the study of host-pathogen interactions and highlights relevant applications of these fields of expertise to the complex case of C gigas infections by OsHV-1, which currently threaten the mollusk production sector worldwide.

Identification and characterization of miRNAs in the gills of the mud crab (Scylla paramamosain) in response to a sudden drop in salinity

BACKGROUND

The mud crab (Scylla paramamosain) is a euryhaline and commercially important species. MiRNAs participate in the regulation of many physiological activities.

RESULTS

The miRNA transcriptome of the gills of S. paramamosain was used to investigate the expression profiles of miRNAs in response to a sudden drop in salinity. In total, seven known miRNAs and 43 novel miRNAs were identified, with 18 differentially expressed small RNAs. Fourteen thousand nine hundred fifty-one differentially expressed miRNAs target genes were screened by prediction. GO analysis of differentially expressed miRNAs target genes indicated that 578 genes associated with cellular processes, 523 associated with metabolic processes, and 422 associated with single-organism processes were the most strongly affected by a sudden drop in salinity from 23‰ to 3‰. KEGG pathway analysis showed 14 pathways were related to amino acid metabolism, which plays an important role in osmoregulation. Besides, several pathways were associated with starch and sucrose metabolism (ko00500), glycosaminoglycan degradation (ko00531), and galactose metabolism (ko00052).

CONCLUSIONS

S. paramamosain regulated osmotic pressure and energy balance by regulating target genes to adapt to a sudden changes in salinity. These results provided a basis for further investigations of miRNA-modulating networks underlying the osmoregulation of S. paramamosain.

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.

Advances in the application of high-throughput sequencing in invertebrate virology

Over the last decade, advances in high-throughput sequencing technologies have revolutionised biological research, making it possible for DNA/RNA sequencing of any organism of interest to be undertaken. Sequencing approaches are now routinely used in the detection and characterisation of (novel) viruses, investigation of host-pathogen interactions, and effective development of disease treatment strategies. For the sequencing and identification of viruses of interest, metagenomics approaches using infected host tissue are frequently used, as it is not always possible to culture and isolate these pathogens. High-throughput sequencing can also be used to investigate host-pathogen interactions by investigating (temporal) transcriptomic responses of both the host and virus, potentially leading to the discovery of novel opportunities for treatment and drug targets. In addition, viruses in environmental samples (e.g. water or soil samples) can be identified using eDNA/metagenomics approaches. The promise that recent developments in sequencing brings to the field of invertebrate virology are not devoid of technical challenges, including the need for better laboratory and bioinformatics strategies to sequence and assemble virus genomes within complex tissue or environmental samples, and the difficulties associated with the annotation of the large number of novel viruses being discovered.

Viruses infecting marine molluscs

Although a wide range of viruses have been reported in marine molluscs, most of these reports rely on ultrastructural examination and few of these viruses have been fully characterized. The lack of marine mollusc cell lines restricts virus isolation capacities and subsequent characterization works. Our current knowledge is mostly restricted to viruses affecting farmed species such as oysters Crassostrea gigas, abalone Haliotis diversicolor supertexta or the scallop Chlamys farreri. Molecular approaches which are needed to identify virus affiliation have been carried out for a small number of viruses, most of them belonging to the Herpesviridae and birnaviridae families. These last years, the use of New Generation Sequencing approach has allowed increasing the number of sequenced viral genomes and has improved our capacity to investigate the diversity of viruses infecting marine molluscs. This new information has in turn allowed designing more efficient diagnostic tools. Moreover, the development of experimental infection protocols has answered some questions regarding the pathogenesis of these viruses and their interactions with their hosts. Control and management of viral diseases in molluscs mostly involve active surveillance, implementation of effective bio security measures and development of breeding programs. However factors triggering pathogen development and the life cycle and status of the viruses outside their mollusc hosts still need further investigations.

Integrated analysis neurimmiRs of tilapia (Oreochromis niloticus) involved in immune response to Streptococcus agalactiae, a pathogen causing meningoencephalitis in teleosts

MicroRNAs (miRNAs) are a class of noncoding RNA molecules and play important roles in a wide spectrum of biological processes, including in immune response. Recent years have witnessed considerable amount of research interest in studies on miRNA-mediated modulation gene function during neuroinflammation. Here, we evaluated Streptococcus agalactiae infected tilapia (Oreochromis niloticus) brain for the expression profile of miRNAs, potential functions and their correlation with genes involved in inflammatory pathways. A total of 1981 miRNAs were identified, including in 486 miRNAs which have homologues in the currently available databases and 1945 novel miRNAs. The expression levels of 547 miRNAs were significantly altered at 6 h-48 h post-bacterial infection, and these miRNAs were therefore classified as differentially expressed tilapia miRNAs. Real-time PCR were implemented for 14 miRNAs co-expressed in five samples, and agreement was confirmed between the high-throughput sequencing and real-time PCR data. For the 486 differentially expressed miRNAs target 41,820 genes. GO and KEGG enrichment analysis revealed that some target genes of miRNAs were grouped mainly into the categories of apoptotic, signal pathwayand immune response. This is the first report of comprehensive identification of teleost miRNAs being differentially regulated in brain in normal conditions relating to bacterial infection.

The miRNA biogenesis in marine bivalves

Small non-coding RNAs include powerful regulators of gene expression, transposon mobility and virus activity. Among the various categories, mature microRNAs (miRNAs) guide the translational repression and decay of several targeted mRNAs. The biogenesis of miRNAs depends on few gene products, essentially conserved from basal to higher metazoans, whose protein domains allow specific interactions with dsRNA. Here, we report the identification of key genes responsible of the miRNA biogenesis in 32 bivalves, with particular attention to the aquaculture species Mytilus galloprovincialis and Crassostrea gigas. In detail, we have identified and phylogenetically compared eight evolutionary conserved proteins: DROSHA, DGCR8, EXP5, RAN, DICER TARBP2, AGO and PIWI. In mussels, we recognized several other proteins participating in the miRNA biogenesis or in the subsequent RNA silencing. According to digital expression analysis, these genes display low and not inducible expression levels in adult mussels and oysters whereas they are considerably expressed during development. As miRNAs play an important role also in the antiviral responses, knowledge on their production and regulative effects can shed light on essential molecular processes and provide new hints for disease prevention in bivalves.

Molecular characterization and immune response expression of the QM gene from the scallop Chlamys farreri

The scallop Chlamys farreri is an important aquaculture species in northern China. However, the sustainable development of the scallop industry is currently threatened by several pathogens that cause mass mortality of this mollusk. Therefore, a complete understanding of the immune response mechanisms involved in host-virus interactions is necessary. This study reports a novel QM gene from C. farreri. This gene was first identified as a putative tumor suppressor gene from human and then confirmed to participate in several functions, including immune response. The QM gene from C. farreri (CfQM) was identified by suppression subtractive hybridization, and its full-length (763 bp) cDNA was obtained through rapid amplification of cDNA ends. The cDNA of CfQM contained a short 5'-UTR of 22 bp and a 3'-UTR of 84 bp. Its ORF comprised 657 nucleotides that encode 218 amino acids with a molecular weight of approximately 28.3 kDa and an isoelectric point of 10.06. The deduced amino acid sequence of CfQM contained a series of conserved functional motifs that belong to the QM family. Phylogenetic analysis revealed that CfQM was closely related to other mollusk QM proteins, and altogether they form a mollusk QM protein subfamily that displays evolutionary conservation from yeast to human. The tissue-specific expression and transcriptional regulation of CfQM were investigated by quantitative real-time PCR in response to bacterial (Vibrio anguillarum) and viral (acute viral necrobiotic virus) challenges. The transcript level of CfQM was high in all of the examined tissues in a constitutive manner. The highest and lowest expression levels of CfQM were measured in the hepatopancreas and hemocyte, respectively. Upon bacterial and viral challenges, the relative mRNA expression of CfQM sharply increased at 6 h post-infection (hpi) and then normalized at 48 hpi. These findings suggest that CfQM can respond to and protect against pathogen challenge. To the best of our knowledge, this study is the first report of the QM gene from scallop. The results presented herein provided new insights into the molecular basis of host-pathogen interactions in C. farreri.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

MicroRNA profile analysis of Epithelioma papulosum cyprini cell line before and after SVCV infection

MicroRNAs (miRNAs) play significant roles in regulating almost all of the biological processes in eukaryotes. An accumulating body of evidence shows that miRNAs are associated with cellular changes following viral infection. Spring viremia of carp virus (SVCV) is the pathogen of Spring viremia of carp (SVC), which results in heavy losses in the cultured common carp (Cyprinus carpio) industry in many countries. To study the involvement of miRNAs during SVCV infection, we adopted the Solexa sequencing technology to sequence small RNA libraries from the Epithelioma papulosum cyprini (EPC) cell line before and after infection with SVCV. In this study, a total of 161 conserved and 26 novel miRNAs were identified. Subsequently, the expression patterns of these miRNAs were compared between the uninfected (control library, M) and SVCV-infected (infection library, E) libraries. In addition, to verify the Solexa sequencing results, the expression patterns of 14 randomly selected miRNAs were validated by qRT-PCR. The targets of the significantly differentially expressed miRNAs were then predicted, and the miRNAs that could directly target the SVCV genome were also predicted. No miRNA encoded by SVCV itself was detected. To the best of our knowledge, this study presents the first miRNA profiling assessment in association with fish rhabdovirus infection, and the data presented lay a foundation for further investigations to determine the roles of miRNAs in regulating the molecular mechanism during SVCV infection.