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

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.

Integrated mRNA and miRNA analysis reveals the regulatory network of oxidative stress and inflammation in Coilia nasus brains during air exposure and salinity mitigation

BACKGROUND

Air exposure is an inevitable source of stress that leads to significant mortality in Coilia nasus. Our previous research demonstrated that adding 10‰ NaCl to aquatic water could enhance survival rates, albeit the molecular mechanisms involved in air exposure and salinity mitigation remained unclear. Conversely, salinity mitigation resulted in decreased plasma glucose levels and improved antioxidative activity. To shed light on this phenomenon, we characterized the transcriptomic changes in the C. nasus brain upon air exposure and salinity mitigation by integrated miRNA-mRNA analysis.

RESULTS

The plasma glucose level was elevated during air exposure, whereas it decreased during salinity mitigation. Antioxidant activity was suppressed during air exposure, but was enhanced during salinity mitigation. A total of 629 differentially expressed miRNAs (DEMs) and 791 differentially expressed genes (DEGs) were detected during air exposure, while 429 DEMs and 1016 DEGs were identified during salinity mitigation. GO analysis revealed that the target genes of DEMs and DEGs were enriched in biological process and cellular component during air exposure and salinity mitigation. KEGG analysis revealed that the target genes of DEMs and DEGs were enriched in metabolism. Integrated analysis showed that 24 and 36 predicted miRNA-mRNA regulatory pairs participating in regulating glucose metabolism, Ca2+ transport, inflammation, and oxidative stress. Interestingly, most of these miRNAs were novel miRNAs.

CONCLUSION

In this study, substantial miRNA-mRNA regulation pairs were predicted via integrated analysis of small RNA sequencing and RNA-Seq. Based on predicted miRNA-mRNA regulation and potential function of DEGs, miRNA-mRNA regulatory network involved in glucose metabolism and Ca2+ transport, inflammation, and oxidative stress in C. nasus brain during air exposure and salinity mitigation. They regulated the increased/decreased plasma glucose and inhibited/promoted antioxidant activity during air exposure and salinity mitigation. Our findings would propose novel insights to the mechanisms underlying fish responses to air exposure and salinity mitigation.

Analysis of adaptive molecular mechanisms in response to low salinity in antennal gland of mud crab, Scylla paramamosain

As an important marine aquaculture species, the mud crab (Scylla paramamosain) is a good candidate for studying the osmoregulatory mechanism of crustaceans. While previous studies have focused on the osmoregulatory function of the gills, this study aims to explore the osmoregulatory function of the antennal glands. By the comparative transcriptomic analysis, we found the pathways of ion regulation including "proximal tubule bicarbonate reclamation" and "mineral absorption" were activated in the antennal glands of the crabs long-term dwelling in low salinity. The enhanced ionic reabsorption was associated with up-regulated ion transport genes such as NKA, CA-c, VPA, and NHE, and with energy metabolism genes such as MDH, SLC25, and PEPCK. The upregulation of NKA and CA-c was also verified by the increased enzyme activity. The lowered osmolality and ion concentration of the hemolymph and the enlarged labyrinth lumen and hemolymph capillary inside the antennal glands indicated the infiltration of external water and the responsively increase of urine excretion, which explained the requirement of enhanced ionic reabsorption. To further confirm these findings, we examined the change of gene expression, enzyme activity, internal ion concentration, and external ion concentration during a 96 h low salinity challenge with seven intervals. The results were basically consistent with the results as shown in the long-term low salinity adaptation. The present study provides valuable information on the osmoregulatory function of the antennal glands of S. paramamosain. The implication of this study in marine aquaculture is that it provides valuable information on the osmoregulatory mechanism of mud crabs, which can be used to improve their culture conditions and enhance their tolerance to salinity stress. The identified genes and pathways involved in osmoregulation can also be potential targets for genetic selection and breeding programs to develop more resilient mud crab strains for aquaculture.

Mechanisms of ammonotelism, epithelium damage, cellular apoptosis, and proliferation in gill of Litopenaeus vannamei under NH4Cl exposure

Excessive ammonia-N in coastal environment and aquaculture threatens the health of marine organisms. To explore the mechanism of gill damage induced by ammonia-N, transcriptome of Litopenaeus vannamei 's gill was carried out under 20 mg/L NH4Cl for 0, 6, and 48 h. K-means clustering analysis suggested that ammonia excretion and metabolism-related genes were elevated. GO and KEGG enrichment analysis suggested that glycosyltransferase activity and amino acid metabolism were affected by ammonia. Moreover, histological observation via three staining methods gave clues on the changes of gill after ammonia-N exposure. Increased mucus, hemocyte infiltration, and lifting of the lamellar epithelium suggested that gill epithelium was suffering damage under ammonia-N stress. Meanwhile, the composition of extracellular matrix (ECM) in connective tissue changed. Based on the findings of transcriptomic and histological analysis, we further investigated the molecular mechanism of gill damage under multiple concentrations of NH4Cl (0, 2, 10, 20 mg/L) for multiple timepoints (0, 3, 6, 12, 24, 48, 72 h). First, ammonia excretion was elevated via ion channel, transporter, and exocytosis pathways, but hemolymph ammonia still kept at a high level under 20 mg/L NH4Cl exposure. Second, we focused on glycosaminoglycan metabolism which was related to the dynamics of ECM. It turned out that the degradation and biosynthesis of chondroitin sulfate (CS) were elevated, suggesting that the structure of CS might be destructed under ammonia-N stress and CS played an important role in maintaining gill structure. It was enlightening that the destructions occurred in extracellular regions were vital to gill damage. Third, ammonia-N stress induced a series of cellular responses including enhanced apoptosis, active inflammation, and inhibited proliferation which were closely linked and jointly led to the impairment of gill. Our results provided some insights into the physiological changes induced by ammonia-N and enriched the understandings of gill damage under environmental stress.

Whole Transcriptome Analysis Reveals the Global Molecular Responses of mRNAs, lncRNAs, miRNAs, circRNAs, and Their ceRNA Networks to Salinity Stress in Hong Kong Oysters, Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is an estuarine bivalve with remarkable commercial value in South China, and the increase of salinity in estuaries during the dry season has posed a major threat to the oyster farming. To explore the global transcriptional response to salinity stress, a whole-transcriptome analysis was performed with the gills of oysters in 6‰, 18‰, and 30‰ filtered seawater. Overall, 2243, 194, 371, and 167 differentially expressed mRNAs (DEmRNAs), differentially expressed long non-coding RNAs (DElncRNAs), differentially expressed circular RNAs (DEcircRNAs), and differentially expressed microRNAs (DEmiRNAs) were identified, respectively. Based on GO enrichment and KEGG pathway analysis, these important DEmRNAs, DElncRNAs, DEcircRNAs, and DEmiRNAs were predicted to be mainly involved in amino acids metabolism, microtubule movement, and immune defense. This demonstrated the complexity of dynamic transcriptomic profiles of C. hongkongensis in response to salinity fluctuation. The regulatory relationships of DEmiRNAs-DEmRNAs, DElncRNAs-DEmiRNAs, and DEcircRNAs-DEmiRNAs were also predicted, and finally, a circRNA-associated competing endogenous RNA (ceRNA) network was constructed, consisting of six DEcircRNAs, eight DEmiRNAs, and five DEmRNAs. The key roles of taurine and hypotaurine metabolism and phenylalanine metabolism were highlighted in this ceRNA network, which was consistent with the major contribution of free amino acids to intracellular osmolality and cell volume regulation. Collectively, this study provides comprehensive data, contributing to the exploration of coding and non-coding RNAs in C. hongkongensis salinity response. The results would benefit the understanding of the response mechanism of bivalves against salinity fluctuation, and provide clues for genetic improvement of C. hongkongensis with hyper-salinity tolerance.

Topologically associating domains in the POLLED region are the same for Angus- and Brahman-specific Hi-C reads from F1 hybrid fetal tissue

Horns, a form of headgear carried by Bovidae, have ethical and economic implications for ruminant production species such as cattle and goats. Hornless (polled) individuals are preferred. In cattle, four genetic variants (Celtic, Friesian, Mongolian and Guarani) are associated with the polled phenotype, which are clustered in a 300-kb region on chromosome 1. As the variants are intergenic, the functional effect is unknown. The aim of this study was to determine if the POLLED variants affect chromatin structure or disrupt enhancers using publicly available data. Topologically associating domains (TADs) were analyzed using Angus- and Brahman-specific Hi-C reads from lung tissue of an Angus (Celtic allele) cross Brahman (horned) fetus. Predicted bovine enhancers and chromatin immunoprecipitation sequencing peaks for histone modifications associated with enhancers (H3K27ac and H3K4me1) were mapped to the POLLED region. TADs analyzed from Angus- and Brahman-specific Hi-C reads were the same, therefore, the Celtic variant does not appear to affect this level of chromatin structure. The Celtic variant is located in a different TAD from the Friesian, Mongolian, and Guarani variants. Predicted enhancers and histone modifications overlapped with the Guarani and Friesian variants but not the Celtic or Mongolian variants. This study provides insight into the mechanisms of the POLLED variants for disrupting horn development. These results should be validated using data produced from the horn bud region of horned and polled bovine fetuses.

Transcriptome Analysis of Marbled Rockfish Sebastiscus marmoratus under Salinity Stress

The marbled rockfish, Sebastiscus marmoratus, belongs to the euryhaline fish and is an oviparous scleractinian fish. There are few studies on the adaptation mechanism, functional genes, and related pathways of S. marmoratus and salinity. The results showed that a total of 72.1 GB of clean reads were obtained and all clean reads annotated a total of 25,278 Unigenes, of which 2,160 were novel genes. Compared to 20‱, 479 and 520 differential genes were obtained for 35‱ and 10‱, respectively. Gene ontology (GO) enrichment analysis revealed significant enrichment in protein binding, ion binding, ATP binding, and catalytic activity. Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that differentially expressed genes significantly expressed under salinity stress were mainly involved in the pathways of the cytochrome P450 metabolism of xenobiotics, tryptophan metabolism, cellular senescence, and calcium signaling pathways. Among them, pik3r6b, cPLA2γ-like, and WSB1 were differentially expressed in all three groups, and they were associated with apoptosis, inflammation, DNA damage, immune regulation, and other physiological processes. Six differentially expressed genes were randomly selected for qRT-PCR validation, and the results showed that the transcriptomic data were of high confidence.

Comprehensive Analysis of Whole-Transcriptome Profiles in Response to Acute Hypersaline Challenge in Chinese Razor Clam Sinonovacula constricta

The Chinese razor clam (Sinonovacula constricta) is an important for Chinese aquaculture marine bivalve that naturally occurs across intertidal and estuarine areas subjected to significant changes in salinity level. However, the information on the molecular mechanisms related to high salinity stress in the species remain limited. In this study, nine gill samples of S. constricta treated with 20, 30, and 40 ppt salinity for 24 h were used for whole-transcriptome RNA sequencing, and a regulatory network of competing endogenous RNAs (ceRNAs) was constructed to better understand the mechanisms responsible for adaptation of the species to high salinity. A total of 83,262 lncRNAs, 52,422 mRNAs, 2890 circRNAs, and 498 miRNAs were identified, and 4175 of them displayed differential expression pattern among the three groups examined. The KEGG analyses of differentially expressed RNAs evidenced that amino acid synthesis and membrane transport were the dominant factors involved in the adaptation of the Chinese razor clam to acute salinity increase, while lipid metabolism and signaling played only a supporting role. In addition, lncRNA/circRNA-miRNA-mRNA regulatory networks (ceRNA network) showed clearly regulatory relationships among different RNAs. Moreover, the expression of four candidate genes, including tyrosine aminotransferase (TAT), hyaluronidase 4 (HYAL4), cysteine sulfinic acid decarboxylase (CSAD), and ∆¹-pyrroline-5-carboxylate synthase (P5CS) at different challenge time were detected by qRT-PCR. The expression trend of TAT and HYAL4 was consistent with that of the ceRNA network, supporting the reliability of established network. The expression of TAT, CSAD, and P5CS were upregulated in response to increased salinity. This might be associated with increased amino acid synthesis rate, which seems to play an essential role in adaptation of the species to high salinity stress. In contrast, the expression level of HYAL4 gene decreased in response to elevated salinity level, which is associated with reduction Hyaluronan hydrolysis to help maintain water in the cell. Our findings provide a very rich reference for understanding the important role of ncRNAs in the salinity adaptation of shellfish. Moreover, the acquired information may be useful for optimization of the artificial breeding of the Chinese razor clam under aquaculture conditions.

Global Metabolomics of Fireflies (Coleoptera: Lampyridae) Explore Metabolic Adaptation to Fresh Water in Insects

Aquatic insects are well-adapted to freshwater environments, but metabolic mechanisms of such adaptations, particularly to primary environmental factors (e.g., hypoxia, water pressure, dark light, and abundant microbes), are poorly known. Most firefly species (Coleoptera: Lampyridae) are terrestrial, but the larvae of a few species are aquatic. We generated 24 global metabolomic profiles of larvae and adults of Aquatica leii (freshwater) and Lychnuris praetexta (terrestrial) to identify freshwater adaptation-related metabolites (AARMs). We identified 110 differentially abundant metabolites (DAMs) in A. leii (adults vs. aquatic larvae) and 183 DAMs in L. praetexta (adults vs. terrestrial larvae). Furthermore, 100 DAMs specific to aquatic A. leii larvae were screened as AARMs via interspecific comparisons (A. leii vs. L. praetexta), which were primarily involved in antioxidant activity, immune response, energy production and metabolism, and chitin biosynthesis. They were assigned to six categories/superclasses (e.g., lipids and lipid-like molecules, organic acids and derivatives, and organoheterocyclic compound). Finally, ten metabolic pathways shared between KEGG terms specific to aquatic fireflies and enriched by AARMs were screened as aquatic adaptation-related pathways (AARPs). These AARPs were primarily involved in energy metabolism, xenobiotic biodegradation, protection of oxidative/immune damage, oxidative stress response, and sense function (e.g., glycine, serine and threonine metabolism, drug metabolism-cytochrome P450, and taste transduction), and certain aspects of morphology (e.g., steroid hormone biosynthesis). These results provide evidence suggesting that abundance changes in metabolomes contribute to freshwater adaptation of fireflies. The metabolites identified here may be vital targets for future work to determine the mechanism of freshwater adaptation in insects.

Sexual reproduction contributes to the evolution of resistance-breaking isolates of the spinach pathogen Peronospora effusa

Peronospora effusa causes downy mildew, the economically most important disease of cultivated spinach worldwide. To date, 19 P. effusa races have been denominated based on their capacity to break spinach resistances, but their genetic diversity and the evolutionary processes that contribute to race emergence are unknown. Here, we performed the first systematic analysis of P. effusa races showing that those emerge by both asexual and sexual reproduction. Specifically, we studied the diversity of 26 P. effusa isolates from 16 denominated races based on mitochondrial and nuclear comparative genomics. Mitochondrial genomes based on long-read sequencing coupled with diversity assessment based on short-read sequencing uncovered two mitochondrial haplogroups, each with distinct genome organization. Nuclear genome-wide comparisons of the 26 isolates revealed that 10 isolates from six races could clearly be divided into three asexually evolving groups, in concordance with their mitochondrial phylogeny. The remaining isolates showed signals of reticulated evolution and discordance between nuclear and mitochondrial phylogenies, suggesting that these evolved through sexual reproduction. Increased understanding of this pathogen's reproductive modes will provide the framework for future studies into the molecular mechanisms underlying race emergence and into the P. effusa-spinach interaction, thus assisting in sustainable production of spinach through knowledge-driven resistance breeding.

Application of Transcriptome Analysis to Understand the Adverse Effects of Hypotonic Stress on Different Development Stages in the Giant Freshwater Prawn Macrobrachium rosenbergii Post-Larvae

Salinity is one of the important environmental factors affecting survival and growth of aquatic animals. However, the impact of low-salinity stress on M. rosenbergii post-larvae at different development stages remains elusive. Therefore, the aim of this study was to explore the underlying mechanisms of hypotonic stress at different development stages of M. rosenbergii post-larvae through transcriptome analysis and antioxidant parameters detection. The salinity of the control group was 15 psu (S15) and the hypotonic stress group was 6 psu (S6). Samples were collected at 7 days-post-hatch (dph), 14 dph and 21 dph larvae. The results showed that hypotonic stress caused oxidative damage in post-larvae evidenced by decreased glutathione peroxidase (GSH-Px); superoxide dismutase (SOD); anti-superoxide anion free radical (ASAFR); and increased malondialdehyde (MDA); nitric oxide (NO); and inducible nitric oxide synthase (iNOS) levels. Transcriptome analysis showed that there were 1428, 1187, 132 DEGs including 301, 366, 4 up-regulated genes and 1127, 821, 128 down-regulated genes at 7 dph, 14 dph and 21 dph larvae under hypotonic stress, respectively. Furthermore, GO and KEGG enrichment indicated that hypotonic stress led to dysregulation of immune signals including lysosome and autophagy in the 7 dph larvae. The autophagy-related genes including beclin 1-associated autophagy-related key regulator (Barkor); ubiquitin-like modifier-activating enzyme ATG7 (ATG7); Beclin; autophagy-related protein 13 (ATG13); nuclear receptor-binding factor 2 (Nrbf2); ubiquitin-like-conjugating enzyme ATG3 (ATG3); vacuole membrane protein 1 (VMP1); and autophagy-related protein 2 (ATG2) decreased at 7 dph, and 14 dph larvae, and then increased at 21 dph larvae under hypotonic stress. In the 14 dph and 21 dph larvae, the renin-angiotensin system was activated. In conclusion, our data indicated that hypotonic stress reduced the antioxidant capacity and impaired the immune system in post-larvae, but as development progresses, the adaptability of post-larvae to hypotonic stress gradually increased, and might reach a new homeostasis through the RAS signaling pathway.

Distribution of Membrane Trafficking System Components Across Ciliate Diversity Highlights Heterogenous Organelle-Associated Machinery

The ciliate phylum is a group of protists noted for their unusual membrane trafficking system and apparent environmental ubiquity; as highly successful microbial predators, they are found in all manner of environments and the ability for specific species to adapt to extremely challenging conditions makes them valued as bioindicators. Ciliates have also been used for many years as cell biological models due to their large cell size and ease of culturing, and for many fundamental cell structures, particularly membrane-bound organelles, ciliates were some of the earliest organisms in which these were observed via microscopy. In this study, we carried out a comparative genomic survey of selected membrane trafficking proteins in a pan-ciliate transcriptome and genome dataset. We observed considerable loss of membrane trafficking system (MTS) proteins that would indicate a loss of machinery that is generally conserved across eukaryotic diversity, even after controlling for potentially incomplete genome representation. In particular, the DSL1 complex was missing in all surveyed ciliates. This protein complex has been shown as involved in peroxisome biogenesis in some model systems, and a paucity of DSL1 components has been indicative of degenerate peroxisome. However, Tetrahymena thermophila (formerly Tetrahymena pyroformis) was one of the original models for visualising peroxisomes. Conversely, the AP3 complex essential for mucocyst maturation in T. thermophila, is poorly conserved despite the presence of secretory lysosome-related organelles across ciliate diversity. We discuss potential resolutions for these apparent paradoxes in the context of the heterogenous distribution of MTS machinery across the diversity of ciliates.

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.

Novel lncRNA lncRNA001074 participates in the low salinity-induced response in the sea cucumber Apostichopus japonicus by targeting the let-7/NKAα axis

Salinity fluctuations have severe impacts on sea cucumbers and therefore important consequences in sea cucumber farming. The responses of sea cucumbers to salinity changes are reflected in the expression profiles of multiple genes and non-coding RNAs (ncRNAs). The microRNA (let-7) which is a developmental regulator, the ion transporter gene sodium potassium ATPase gene (NKAα), and the long ncRNA lncRNA001074 were previously shown to be involved in responses to salinity changes in various marine species. To better understand the relationship between ncRNAs and target genes, the let-7/NKAα/lncRNA001074 predicted interaction was investigated in this study using luciferase reporter assays and gene knockdowns in the sea cucumber Apostichopus japonicus. The results showed that NKAα was the target gene of let-7 and NKAα expression levels were inversely correlated with let-7 expression based on the luciferase reporter assays and western blots. The let-7 abundance was negatively regulated by lncRNA001074 and NKAα both in vitro and in vivo. Knockdown of lncRNA001074 led to let-7 overexpression. These results demonstrated that lncRNA001074 binds to the 3'-UTR binding site of let-7 in a regulatory manner. Furthermore, the expression profiles of let-7, NKAα, and lncRNA001074 were analyzed in sea cucumbers after the knockdown of each of these genes. The results found that lncRNA001074 competitively bound let-7 to suppress NKAα expression under low salinity conditions. The downregulation of let-7, in conjunction with the upregulation of lncRNA001074 and NKAα, may be essential for the response to low salinity change in sea cucumbers. Therefore, the dynamic balance of the lncRNA001074, NKAα, and let-7 network might be a potential response mechanism to salinity change in sea cucumbers.

Comparative transcriptome analysis of the gills of Cardisoma armatum provides novel insights into the terrestrial adaptive related mechanism of air exposure stress

Cardisoma armatum is a typical member of the Gecarcinidae which show significant behavioral, morphological, physiological, and/or biochemical adaptations permitting extended activities on the land. The special gills (branchiostegal lung) of C. armatum play an important role in maintaining osmotic pressure balance and obtaining oxygen to adapt to the terrestrial environment. However, adaptive molecular mechanisms responding to air exposure in C. armatum are still poorly understood. In this study, transcriptomic analysis and histological analysis were conducted on the gills to test adaptive capabilities over 8 h between the aerial exposure (AE) and the water immersion (WI) group. Differentially expressed genes (DEGs) related to terrestrial adaptation were categorized into four broad categories: ion transport, acid-base balance, energy metabolism and immune response. This is the first research to reveal the molecular mechanism of terrestrial adaptation in C. armatum, and will provide new insight into the molecular genetic basis of terrestrial adaptation in crabs.

A chromosome-level genome of the mud crab (Scylla paramamosain estampador) provides insights into the evolution of chemical and light perception in this crustacean

Mud crabs, found throughout the Indo-Pacific region, are coastal species that are important fisheries resources in many tropical and subtropical Asian countries. Here, we present a chromosome-level genome assembly of a mud crab (Scylla paramamosain). The genome is 1.55 Gb (contig N50 191 kb) in length and encodes 17,821 proteins. The heterozygosity of the assembled genome was estimated to be 0.47%. Effective population size analysis suggested that an initial large population size of this species was maintained until 200 thousand years ago. The contraction of cuticle protein and opsin genes compared with Litopenaeus vannamei is assumed to be correlated with shell hardness and light perception ability, respectively. Furthermore, the analysis of three chemoreceptor gene families, the odorant receptor (OR), gustatory receptor (GR) and ionotropic receptor (IR) families, suggested that the mud crab has no OR genes and shows a contraction of GR genes and expansion of IR genes. The numbers of the three gene families were similar to those in three other decapods but different from those in two nondecapods and insects. In addition, IRs were more diversified in decapods than in nondecapod crustaceans, and most of the expanded IRs in the mud crab genome were clustered with the antennal IR clades. These findings suggested that IRs might exhibit more diverse functions in decapods than in nondecapods, which may compensate for the smaller number of GR genes. Decoding the S. paramamosain genome not only provides insight into the genetic changes underpinning ecological traits but also provides valuable information for improving the breeding and aquaculture of this species.

De novo transcriptome assembly for Pachygrapsus marmoratus, an intertidal brachyuran crab

The marble crab Pachygrapsus marmoratus inhabits the rocky shores of the Mediterranean Sea, Black Sea and East Atlantic Ocean. As other intertidal species, it is considered a model species to study the effects of environmental stressors on natural populations. In this study, we performed Illumina next-generation sequencing on eleven P. marmoratus specimens with the aims to (i) reconstruct their whole transcriptome, (ii) perform a functional annotation of the assembled transcriptome and (iii) develop gene-based markers for future genetic and genomic studies on this as well as other brachyuran species. We obtained a transcriptome assembly constituted by 56,308 unigenes and covering about 60.3 Mbp. We detected 43,915 Simple Sequence Repeats (SSRs) and 192,631 high-quality Single Nucleotide Polymorphisms (SNPs). Due to the scarcity of genomic resources in decapods, and crabs in particular, our results constitute a valuable resource for future studies on brachyuran crabs. The present data also represent a sound resource to investigate biological responses to pollution in intertidal and marine populations.

Gill Transcriptome Sequencing and De Novo Annotation of Acanthogobius ommaturus in Response to Salinity Stress

Acanthogobius ommaturus is a euryhaline fish widely distributed in coastal, bay and estuarine areas, showing a strong tolerance to salinity. In order to understand the mechanism of adaptation to salinity stress, RNA-seq was used to compare the transcriptome responses of Acanthogobius ommaturus to the changes of salinity. Four salinity gradients, 0 psu, 15 psu (control), 30 psu and 45 psu were set to conduct the experiment. In total, 131,225 unigenes were obtained from the gill tissue of A. ommaturus using the Illumina HiSeq 2000 platform (San Diego, USA). Compared with the gene expression profile of the control group, 572 differentially expressed genes (DEGs) were screened, with 150 at 0 psu, 170 at 30 psu, and 252 at 45 psu. Additionally, among these DEGs, Gene Ontology (GO) analysis indicated that binding, metabolic processes and cellular processes were significantly enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis detected 3, 5 and 8 pathways related to signal transduction, metabolism, digestive and endocrine systems at 0 psu, 30 psu and 45 psu, respectively. Based on GO enrichment analysis and manual literature searches, the results of the present study indicated that A. ommaturus mainly responded to energy metabolism, ion transport and signal transduction to resist the damage caused by salinity stress. Eight DEGs were randomly selected for further validation by quantitative real-time PCR (qRT-PCR) and the results were consistent with the RNA-seq data.

The miRNAs profiling revealed by high-throughput sequencing upon WSSV infection in mud crab Scylla paramamosain

microRNAs (miRNAs) are known to regulate various immune functions by silencing the target genes in both vertebrates and invertebrates. However, in mud crab Scylla paramamosain, the role of miRNAs during the response to virus invasion remains unclear. To investigate the roles of miRNAs in S. paramamosain during virus infection, the mud crab was challenged with white spot syndrome virus (WSSV) and then subjected to the transcriptional analysis at different conditions. The results of high-throughput sequencing revealed that 940,379 and 1,306,023 high-quality mappable reads were detected in the hemocyte of normal and WSSV-infected mud crabs, respectively. Besides, the total number of 261 unique miRNAs were identified. Among them, 131 miRNAs were specifically expressed in the hemocytes of normal mud crabs, 46 miRNAs were specifically transcribed in those of WSSV-infected individuals, the other 84 miRNAs were expressed in both normal and WSSV-infected individuals. Furthermore, a number of 152 (89 down-regulated and 63 up-regulated) miRNAs were found to be differentially expressed in the WSSV-infected hemocytes, normalized to the controls. The identified miRNAs were subjected to GO analysis and target gene prediction and the results suggested that the differentially regulated miRNAs were mainly correlated with the changes of the immune responses of the hemocytes, including phagocytosis, melanism, and apoptosis as well. Taken together, the results demonstrated that the expressed miRNAs during the virus infection were mainly involved in the regulation of immunological pathways in mud crabs. Our findings not only enrich the understanding of the functions of miRNAs in the innate immune system but also provide some novel potential targets for the prevention of WSSV infection in crustaceans.

miR-10 involved in salinity-induced stress responses and targets TBC1D5 in the sea cucumber, Apostichopus japonicas

The sea cucumber is an economically important aquaculture species in China, where it encounter hypo-saline conditions caused by freshwater outflow from rivers and rainfall. MicroRNAs (miRNA) are small noncoding RNAs of about 22 nucleotides, which are crucial regulators of gene expression at the post-transcriptional level and are involved in a variety of physiological and pathophysiological processes. miR-10 is differentially expressed in salinity acclimation, and has a seed-region match with TBC1D5. The expression profiles of miR-10 and TBC1D5 indicate that miR-10 negatively regulates the expression of TBC1D5 in coelomocytes and sea cucumbers with a miR-10 agomir or antagomir. During salinity acclimation, up-regulation of miR-10 was induced after transfection in coelomocytes with a miR-10 inhibitor, while down-regulation of TBC1D5 was induced. The miR-10 expression maximum in coelomocytes appeared at 48 h post-transfection with a miR-10 inhibitor, was later than that of in sea cucumbers, which appeared 24 h after miR-10 antagomir injection. There was no longer a negative relationship between miR-10 and TBC1D5 expression in coelomocytes and sea cucumbers with miR-10 mimics or agomir during salinity acclimation. The miR-10 antagomir or agomir only affected sodium and NKA enzyme activities, and has little effect on other chloride and potassium ions. Our results demonstrate miR-10 directly regulates TBC1D5 by targeting its 3'-UTR, and that miR-10 suppression substantially increases TBC1D5 mRNA levels in vivo and in vitro. Furthermore, miR-10 and TBC1D5 fluctuating expression patterns after treatment with a miR-10 inhibitor or mimics during salinity acclimation may indicate that they are required for adaptation to salinity stress caused by environmental change. Especially, the miR-10 up-regulation in coelomocytes with miR-10 inhibitor during salinity acclimation indicated that they are required for adaptation to salinity stress caused by environmental change. We propose that miR-10 participates in a regulatory circuit that allows for rapid gene program transitions in response to osmotic stress.

Integrative Proteomic and MicroRNA Analysis: Insights Into Mechanisms of Eyestalk Ablation-Induced Ovarian Maturation in the Swimming Crab Portunus trituberculatus

Eyestalk ablation is the most common method to induce ovarian maturation in decapod crustacean aquaculture, but it jeopardizes broodstock survival and larvae production. It is important to understand the molecular basis underlying the maturation triggered by ablation and thereby develop an alternative measure for maturation manipulation. In this study, we investigate alterations of ovarian proteome and miRNA profile after ablation in a commercially important marine crab Portunus trituberculatus. Quantitative proteomic analysis using iTRAQ reveals that 163 proteins are differentially expressed following ablation, and modulation of methyl farnesoate metabolism and activation of calcium signaling may play important roles in the ovarian maturation induced by ablation. miRNA expression profiling identifies 31 miRNAs that show statistically significant changes. Integration of miRNA and proteome expression data with miRNA target prediction algorithms generates a potential regulatory network consisting of 26 miRNAs and 30 proteins linked by 71 possible functional associations. The miRNA-protein network analysis suggests that miRNAs are involved in promoting ovarian maturation by controlling expression of proteins related to methyl farnesoate synthesis, calcium signals, and energy metabolism. Experimental validation and temporal expression analysis indicate multiple miRNAs can act synergistically to regulate expression of Farnesoic acid O-methyltransferase and Calmodulin. Our findings provide new insights for elucidating the mechanisms underlying eyestalk ablation-induced ovarian maturation and could be useful for devising an alternative technique for manipulating reproduction in P. trituberculatus and other decapods.

Identification and profiling of microRNAs responsive to cadmium toxicity in hepatopancreas of the freshwater crab Sinopotamon henanense

Background

Cadmium (Cd) is a ubiquitous environmental toxicant for aquatic animals. The freshwater crab, Sinopotamon henanense (S. henanense), is a useful model for monitoring Cd exposure since it is widely distributed in sediments whereby it tends to accumulate several toxicants, including Cd. In the recent years, the toxic effects of Cd in the hepatopancreas of S. henanense have been demonstrated by a series of biochemical analysis and ultrastructural observations as well as the deep sequencing approaches and gene expression profile analysis. However, the post-transcriptional regulatory network underlying Cd toxicity in S.henanense is still largely unknown.

Results

The miRNA transcriptional profile of the hepatopancreas of S. henanense was used to investigate the expression levels of miRNAs in response to Cd toxicity. In total, 464 known miRNAs and 191 novel miRNAs were identified. Among these 656 miRNAs, 126 known miRNAs could be matched with the miRNAs of Portunus trituberculatus, Eriocheir sinensis and Scylla paramamosain. Furthermore, a total of 24 conserved miRNAs were detected in these four crab species. Fifty-one differentially expressed miRNAs were identified in the Cd-exposed group, with 31 up-regulated and 20 down-regulated. Eight of the differentially expressed miRNAs were randomly selected and verified by the quantitative real-time PCR (qRT-PCR), and there was a general consistency (87.25%) between the qRT-PCR and miRNA transcriptome data. A total of 5258 target genes were screened by bioinformatics prediction. GO term analysis showed that, 17 GO terms were significantly enriched, which were mainly related to the regulation of oxidoreductase activity. KEGG pathway analysis showed that 18 pathways were significantly enriched, which were mainly associated with the biosynthesis, modification and degradation of proteins.

Conclusion

In response to Cd toxicity, in the hepatopancreas of S. henanense, the expressions of significant amount of miRNAs were altered, which may be an adaptation to resist the oxidative stress induced by Cd. These results provide a basis for further studies of miRNA-mediated functional adaptation of the animal to combat Cd toxicity.

High-throughput sequencing reveals the core gut microbiota of the mud crab (Scylla paramamosain) in different coastal regions of southern China

Background

Scylla paramamosain is a commercially important mud crab. The microbiota is a community that inhabits the crab intestine, and is important for physiological functional and host health.

Results

Proteobacteria, Firmicutes, Bacteroidetes, Tenericutes, Spirochaetae and Fusobacteria were the dominant phyla of the 36 representative phyla. Eleven genera of the 820 representative genera were considered as core gut microbiota and were distributed in the five dominant phyla. The core genus of the Proteobacteria included Arcobacter, Photobacterium, Vibrio, Shewanella and Desulfovibrio. The other four phyla contained one or two genera. Male and female crab samples had two different core genera, (male samples: Psychrilyobacter & Lactococcus; female samples: Clostridium_sensu_stricto_11 and Candidatus_Bacilloplasma).

Conclusions

This is the first time core intestinal microbiota have been identified in crab from nine coastal regions of southern China. This study provides sequencing data related to the gut microbiota of S. paramamosain, and may contribute to probiotic development for S. paramamosain aquaculture industries.

Gonadal microRNA Expression Profiles and Their Potential Role in Sex Differentiation and Gonadal Maturation of Mud Crab Scylla paramamosain

Although the sexual dimorphism in terms of gonadal development and gametogenesis of mud crab has been described, the internal regulating mechanism and sex differentiation process remain unclear. A comparative gonadal miRNA transcriptomic study was conducted to identify miRNAs that are differentially expressed between testes and ovaries, and potentially uncover miRNAs that might be involved in sex differentiation and gonadal maturation mechanisms of mud crabs (Scylla paramamosain). A total of 10 known miRNAs and 130 novel miRNAs were identified, among which 54 were differentially expressed. Target gene prediction revealed a significant enrichment in 30 KEGG pathways, including some reproduction-related pathways, e.g. phosphatidylinositol signalling system and inositol phosphate metabolism pathways. Further analysis on six differentially expressed known miRNAs, six differentially expressed novel miRNAs and their reproduction-related putative target genes shows that both miRNAs and putative target genes showed stage-specific expression during gonadal maturation, suggesting their potential regulatory roles in sex differentiation and reproductive development. This study reveals the sex-biased miRNA profile and establishes a solid foundation for understanding the sex differentiation and gonadal maturation mechanisms of S. paramamosain.