De novo assembly and annotation of the whole transcriptome of Oratosquilla oratoria

Full-length transcriptome reveals the circularly polarized light response-related molecular genetic characteristics of Oratosquilla oratoria

The mantis shrimp is the only animal that can recognize circularly polarized light (CPL), but its molecular genetic characteristics are unclear. Multi-tissue level full-length (FL) transcriptome sequencing of Oratosquilla oratoria, a representative widely distributed mantis shrimp, was performed in the present study. We used comparative transcriptomics to explore the critical genes of O. oratoria selected by CPL and the GNβ gene associated with CPL signal transduction was hypothesized to be positively selected. Furthermore, the FL transcriptomes of O. oratoria compound eyes under five light conditions were sequenced and used to detect alternative splicing (AS). The ASs associated with CPL recognition mainly occurred in the LWS, ARR and TRPC regions. The number of FL transcripts with AS events and annotation information also provided evidence that O. oratoria could recognize LCPL. Additionally, 51 sequences belonging to the LWS, UV and Peropsin gene families were identified based on conserved 7tm domains. The LWS, UV and Peropsin opsins have similar 3D structures with seven domains across the cell membrane and conserved KSLRTPSN, DRY, and QAKK motifs. In conclusion, these results are undoubtedly valuable for perfecting the vision theory of O. oratoria and other mantis shrimp.

How Oratosquilla oratoria compound eye response to the polarization of light: In the perspective of vision genes and related proteins

The special rhabdom structure of the mid-band ommatidium in compound eye contributes to the mantis shrimp being the only animal species known to science that can recognize circularly polarized light (CPL). Although the number of mid-band ommatidium of Oratosquilla oratoria is reduced, the mid-band ommatidium still has orthogonal geometric interleaved rhabdom and short oval distal rhabdom, which may mean that the O. oratoria has weakened circular polarized light vision (CPLV). Here we explored the molecular mechanisms of how O. oratoria response to the polarization of light. Based on the specific expression patterns of vision-related functional genes and proteins, we suggest that the order of light response by O. oratoria compound eye was first natural light, then left-circularly polarized light (LCPL), linearly polarized light, right-circularly polarized light (RCPL) and dark. Meanwhile, we found that the expression levels of vision-related functional genes and proteins in O. oratoria compound eye under RCPL were not significantly different from those in DL, which may imply that O. oratoria cannot respond to RCPL. Furthermore, the response of LCPL is likely facilitated by the differential expression of opsin and microvilli - related functional genes and proteins (arrestin and sodium-coupled neutral amino acid transporter). In conclusion, this study systematically illustrated for the first time how O. oratoria compound eye response to the polarization of light at the genetic level, and it can improve the visual ecological theory behind polarized light vision evolution.

Exploitation pattern assessment of the Japanese mantis shrimp (Oratosquilla oratoria) resource in the coastal waters of the Shandong Peninsula

The exploitation pattern of the invertebrate resource (Oratosquilla oratoria) was investigated in the coastal waters of the Shandong Peninsula, along with the seasonal variation in body length structure, spatial distribution in abundance, and interannual stock status. Results showed that the model with only catch data suggested both stocks in the north and south were suffering from extreme fishing pressure without explicit recovery (North: B2019/BMSY = 0.468 and F2019/FMSY = 1.88 in CMSY. South: B2019/BMSY = 0.349 and F2019/FMSY = 2.59 in CMSY). However, the other two assessment models indicated that the northern stock began to gradually recover as the fishing pressure dropped to an appropriate level after the original overfished status (North: B2019/BMSY = 0.738 and F2019/FMSY = 0.882 in AMSY, B2019/BMSY = 0.831 and F2019/FMSY = 0.774 in BSM. South: B2019/BMSY = 0.164 and F2019/FMSY = 1.44 in AMSY, B2019/BMSY = 0.384 and F2019/FMSY = 1.76 in BSM). Overall, the stock status in the north was better than that in the south. This study suggested that spatial exploitation pattern and quarterly differences should be considered in fishery management process.

Metabolomics integrated with transcriptomics reveals the changes during developmental stages in Shiraia bambusicola

Shiraia bambusicola is a fungus with high economic value widely used in medicine, agriculture, and food. We wished to understand the genes and metabolites changes involved in the different developmental stages of S. bambusicola. So, to reveal key genes and metabolites in the main active metabolite, the  were analyzed in different developmental stages of S. bambusicola fruiting body. A total of 29,137 Unigenes were annotated. In the whole growth process, differentially expressed genes were involved in the pathways of cytochrome P450, transcription factors, transporters, and so on, while in the early stage of growth, genes enriching to synthesis pathways of basic substances. In the middle stage of growth, genes with more prominent changes were involved in the pathways of the cell cycle, cancer mechanisms, and aminobenzoate degradation; in the later stage of growth, differentially expressed genes that enriched synthesis pathways of secondary metabolites. A total of 612 metabolites were detected from different growth stages of S. bambusicola. Among them, coumarins, alkaloids, rutin, liquiritigenin, quercetin, and other medically relevant metabolites were detected for the first time. We have identified 31 secondary metabolites,  relevantly only accumulated in the early and middle stage, but not detected in the later stage, such as flavonols, coumarins, nucleotides and its derivates and hydroxycinnamoyl derivatives. The differential genes and metabolites of the same group were enriched in 127 pathways, and more significantly in ubiquinone and other terpenoid quinone biosynthesis, phenylalanine metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, and phenylpropanoid biosynthesis. The correlation networks of several significantly enriched pathways were analyzed, and the relationships within and between these pathways, genes, and metabolites, were analyzed. The synthetic pathway of hypocrellin has been speculated upon. We believe that hypocrellin is synthesized in S. bambusicola via the shikimic acid pathway followed by phenylalanine, tyrosine, and tryptophan biosynthesis pathway, then the ubiquinone and other terpenoid quinone biosynthesis pathway, and finally a series of polymerization and modification reactions. Several genes and metabolites involved in the biosynthesis of hypocrellin have been identified. This study provides a reference for further research on S. bambusicola, by providing a basis for its use and development.

Metabolomic Alterations in the Digestive System of the Mantis Shrimp Oratosquilla oratoria Following Short-Term Exposure to Cadmium

Mantis shrimp Oratosquilla oratoria is an economically critical aquatic species along the coast of China but strongly accumulates marine pollutant cadmium (Cd) in its digestive system. It is necessary to characterize the toxicity of Cd in the digestive system of mantis shrimp. The metabolic process is an essential target of Cd toxicity response. In this work, we used ultra-performance liquid chromatography coupled with time-of-flight mass spectrometry (UPLC-TOF-MS) for untargeted metabolomics to characterize the metabolic changes in the digestive system of O. oratoria, exposed to 0.05 mg/L for 96 h. The aim of this study was to further investigate the effect of O. oratoria on Cd response to toxicity and develop biomarkers. Metabolomics analysis showed the alteration of metabolism in the digestive system of mantis shrimp under Cd stress. A total of 91 metabolites were differentially expressed and their main functions were classified into amino acids, phospholipids, and fatty acid esters. The enrichment results of differential metabolite functional pathways showed that biological processes such as amino acid metabolism, transmembrane transport, energy metabolism, and signal transduction are significantly affected. Based on the above results, the Cd-induced oxidative stress and energy metabolism disorders were characterized by the differential expression of amino acids and ADP in mantis shrimp, while the interference of transmembrane transport and signal transduction was due to the differential expression of phospholipids. Overall, this work initially discussed the toxicological response of Cd stress to O. oratoria from the metabolic level and provided new insights into the mechanism.

Identification of putative key genes for thermal adaptation in the Japanese mantis shrimp (Oratosquilla oratoria) through population genomic analysis

Little is known about the mechanisms underlying the relationship between genetic variation and the adaptation of Oratosquilla oratoria populations to different habitat temperature. Here, the genome-wide genetic information of three O. oratoria populations were obtained by IIB restriction site-associated DNA (2b-RAD) sequencing and 2403 single-nucleotide polymorphisms (SNPs) were identified. Based on the 2403 SNPs, we found a remarkable genetic differentiation between the Yellow Sea and the East China Sea groups of O. oratoria. Furthermore, 63 SNPs are thought to be associated with different sea temperatures. Based on the 63 SNPs, it is hypothesised that the long-term temperature differences may contribute to the variation of genes associated with multiple biological functions, such as material metabolism, cytoskeleton, cellular processes, inflammatory response and hormonal regulation. This study provides new information for elucidating the molecular mechanisms underlying the relationship between genetic variation and the adaptation of Oratosquilla oratoria populations to different temperature.

Local adaptation of a dominant coastal tree to freshwater availability and solar radiation suggested by genomic and ecophysiological approaches

Local adaptation is often a product of environmental variations in geographical space and has implications for biodiversity conservation. We investigated the role of latitudinal heterogeneity in climate on the organization of genetic and phenotypic variation in the dominant coastal tree Avicennia schaueriana. In a common garden experiment, samples from an equatorial region, with pronounced seasonality in precipitation, accumulated less biomass, and showed lower stomatal conductance and transpiration, narrower xylem vessels, smaller leaves and higher reflectance of long wavelengths by the stem epidermis than samples from a subtropical region, with seasonality in temperature and no dry season. Transcriptomic differences identified between trees sampled under field conditions at equatorial and subtropical sites, were enriched in functional categories such as responses to temperature, solar radiation, water deficit, photosynthesis and cell wall biosynthesis. Remarkably, the diversity based on genome-wide SNPs revealed a north-south genetic structure and signatures of selection were identified for loci associated with photosynthesis, anthocyanin accumulation and the responses to osmotic and hypoxia stresses. Our results suggest the existence of divergence in key resource-use characteristics, likely driven by seasonality in water deficit and solar radiation. These findings provide a basis for conservation plans and for predicting coastal plants responses to climate change.

The Single-molecule long-read sequencing of Scylla paramamosain

Scylla paramamosain is an important aquaculture crab, which has great economical and nutritional value. To the best of our knowledge, few full-length crab transcriptomes are available. In this study, a library composed of 12 different tissues including gill, hepatopancreas, muscle, cerebral ganglion, eyestalk, thoracic ganglia, intestine, heart, testis, ovary, sperm reservoir, and hemocyte was constructed and sequenced using Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology. A total of 284803 full-length non-chimeric reads were obtained, from which 79005 high-quality unique transcripts were obtained after error correction and sequence clustering and redundant. Additionally, a total of 52544 transcripts were annotated against protein database (NCBI nonredundant, Swiss-Prot, KOG, and KEGG database). A total of 23644 long non-coding RNAs (lncRNAs) and 131561 simple sequence repeats (SSRs) were identified. Meanwhile, the isoforms of many genes were also identified in this study. Our study provides a rich set of full-length cDNA sequences for S. paramamosain, which will greatly facilitate S. paramamosain research.

Transcriptome analyses reveal alterations in muscle metabolism, immune responses and reproductive behavior of Japanese mantis shrimp (Oratosquilla oratoria) at different cold temperature

Low temperature reduction is thought to cause widespread effects on the physical and behavioral traits of marine organisms, which include metabolic processes, immune responses, and reproductive behavior. Crustaceans are generally considered sensitive to temperature reduction due to the lack of efficient regulators. To better understand the molecular regulatory mechanisms of crustacean exposure to cold stress, Japanese mantis shrimp (Oratosquilla oratoria) was chosen as a representative crustacean. Transcriptomic responses in O. oratoria from five temperatures (25 °C, 22 °C, 19 °C, 16 °C, and 13 °C) were studied using RNA-seq. A total of 64.91 Gb of clean transcriptomic data were generated in 10 libraries and then spliced into 52,107 unigenes with an average length of 1089 bp and an N50 length of 1872 bp. A total of 14,841 unigenes was annotated in at least one database using Blastx alignment. Compared with the control temperature (25 °C), 7, 21, 58, and 236 unigenes were significantly differentially expressed at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. GO analysis showed that 6, 20, 27, and 35 terms were significantly enriched at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. In addition, 2, 5, 2, and 10 significant pathways were presented at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. Combining NR, GO, and KEGG annotation information, many genes significantly differentially expressed at low temperatures may be associated with metabolic processes, immune response, and reproductive behavior. Additionally, we reconstructed the phylogenetic relationship based on 366 orthologous genes and the predicted differentiation time of O. oratoria and P. vannamei range from 212.82 to 365.30 Mya. Furthermore, 16 orthologous genes were identified as PSGs and 30 orthologous genes were identified as FEGs and these adaptive genes were associated with energy metabolism, stress response and immunity, and multiple cellular processing. These results provide fundamental information about molecular mechanisms regulating cold stress response of O. oratoria.

The Gastric Ganglion of Octopus vulgaris: Preliminary Characterization of Gene- and Putative Neurochemical-Complexity, and the Effect of Aggregata octopiana Digestive Tract Infection on Gene Expression

The gastric ganglion is the largest visceral ganglion in cephalopods. It is connected to the brain and is implicated in regulation of digestive tract functions. Here we have investigated the neurochemical complexity (through in silico gene expression analysis and immunohistochemistry) of the gastric ganglion in Octopus vulgaris and tested whether the expression of a selected number of genes was influenced by the magnitude of digestive tract parasitic infection by Aggregata octopiana. Novel evidence was obtained for putative peptide and non-peptide neurotransmitters in the gastric ganglion: cephalotocin, corticotrophin releasing factor, FMRFamide, gamma amino butyric acid, 5-hydroxytryptamine, molluscan insulin-related peptide 3, peptide PRQFV-amide, and tachykinin-related peptide. Receptors for cholecystokininA and cholecystokininB, and orexin2 were also identified in this context for the first time. We report evidence for acetylcholine, dopamine, noradrenaline, octopamine, small cardioactive peptide related peptide, and receptors for cephalotocin and octopressin, confirming previous publications. The effects of Aggregata observed here extend those previously described by showing effects on the gastric ganglion; in animals with a higher level of infection, genes implicated in inflammation (NFκB, fascin, serpinB10 and the toll-like 3 receptor) increased their relative expression, but TNF-α gene expression was lower as was expression of other genes implicated in oxidative stress (i.e., superoxide dismutase, peroxiredoxin 6, and glutathione peroxidase). Elevated Aggregata levels in the octopuses corresponded to an increase in the expression of the cholecystokininA receptor and the small cardioactive peptide-related peptide. In contrast, we observed decreased relative expression of cephalotocin, dopamine β-hydroxylase, peptide PRQFV-amide, and tachykinin-related peptide genes. A discussion is provided on (i) potential roles of the various molecules in food intake regulation and digestive tract motility control and (ii) the difference in relative gene expression in the gastric ganglion in octopus with relatively high and low parasitic loads and the similarities to changes in the enteric innervation of mammals with digestive tract parasites. Our results provide additional data to the described neurochemical complexity of O. vulgaris gastric ganglion.