Next Article in Journal
A Possible Novel Protective Effect of Piceatannol against Isoproterenol (ISO)-Induced Histopathological, Histochemical, and Immunohistochemical Changes in Male Wistar Rats
Next Article in Special Issue
TNF-α Suppresses Apelin Receptor Expression in Mouse Quadriceps Femoris-Derived Cells
Previous Article in Journal
Expression of GOT2 Is Epigenetically Regulated by DNA Methylation and Correlates with Immune Infiltrates in Clear-Cell Renal Cell Carcinoma
Previous Article in Special Issue
A Glance at the Molecules That Regulate Oligodendrocyte Myelination
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Insulin-like Peptide (ILP) and Its Potential Role in Ovarian Development of the Cuttlefish Sepiella japonica

1
National Engineering Laboratory of Marine Germplasm Resources Exploration and Utilization, College of Marine Science and Technology, Zhejiang Ocean University, Zhoushan 316022, China
2
National Engineering Research Center for Facilitated Marine Aquaculture, Zhejiang Ocean University, Zhoushan 316022, China
*
Author to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2022, 44(6), 2490-2504; https://doi.org/10.3390/cimb44060170
Submission received: 6 April 2022 / Revised: 19 May 2022 / Accepted: 24 May 2022 / Published: 27 May 2022

Abstract

:
The insulin-like peptide (ILP) family is well known for regulating reproduction in invertebrates, while its role in mollusks remains largely unknown. In this study, we first isolated and characterized the ILP gene in the cuttlefish Sepiella japonica. The full-length SjILP cDNA obtained was 926 bp and encoded a precursor protein of 161 amino acids. The precursor protein consisted of a signal peptide, a B chain, a C-peptide, and an A chain. It possessed the typical features of ILP proteins, including two cleavage sites (KR) and eight conserved cysteines. To define the function of SjILP, the expression of SjILP in different tissues and ovarian development stages were analyzed using qRT-PCR. SjILP was mainly expressed in the ovary, and its gene expression correlated with ovarian development. Furthermore, silencing SjILP using RNA interference (RNAi) dramatically decreased the expression levels of four ovarian-development-related genes (vitellogenin1, vitellogenin2, cathepsin L1-like, and follistatin). These data suggest the critical role of SjILP in the regulation of ovarian development in S. japonica.

1. Introduction

Insulin-like peptides (ILPs), the most widely distributed peptide among invertebrate and invertebrate species, are essential regulatory factors in animal growth, development, metabolism, and reproduction [1]. In mollusks, since the first identification of insulin-like peptides (MIPs) in the pond snail Lymnaea stagnalis [2], a variable number of ILPs have been identified in different species, including gastropods [3], bivalves [4,5], and cephalopods. Unlike their vertebrate counterparts, the primary source of ILPs in mollusks is the nervous system [6]. Beyond that, ILPs are also produced in other tissues, such as the gill [7], intestinal epithelial [8], digestive gland [9], adductor muscle [10], siphon, and foot [11]. These peptides are stored in the secretory granules in response to cellular or molecular stimuli. Once released, the peptides can diffuse to target cells, where they combine with their receptors to regulate multiple processes [12].
The biological functions of insulin observed in vertebrate species are generally conserved in invertebrates. In invertebrates such as mollusks, multiple ILPs are suggested to have functions other than just regulating animal growth [10], including the regulation of reproduction [13]. In Planorbarius corneus [14] and Lymnaea stagnalis [6], ILPs have been detected in the brain and are involved in the hormonal control of reproduction [15]. In isolated gonadal cells of Helix aspersa, Monnier and Bride [16] showed that vertebrate insulin stimulated protein synthesis, suggesting that ILPs may be involved in the regulation of germinal cell proliferation and maturation. In the Pacific oyster Crassostrea gigas, Cg-ILP is mainly expressed in the gonadal area, suggesting involvement in the control of sexual reproduction [5]. Combined with the finding that one receptor for C. gigas IRPs (CIR) is expressed in the gonadal area, these results suggest that ILPs are involved in reproduction via binding to the insulin receptor (IR) [13]. Genetic studies have indicated that multiple ILPs may activate this receptor and stimulate signal transduction through the PI3K/AKT pathway (Figure 1) [17].
The cuttlefish Sepiella japonica is one of the most important commercial species widely distributed in the Northwest Pacific Ocean and North Indian Ocean. However, due to the over-fishing of this species, the wild population has declined. Therefore, many artificial breeding programs have been established to increase the production of S. japonica [18]. However, it appeared to have sexual precocity in an artificial environment. Previous studies showed that the sexual precocity of female individuals was associated with ovarian development [19]. Furthermore, ovarian development is influenced by both genetic and environmental factors [20,21,22]. Among the genetic factors, some ovarian-development-related genes, such as ILP, have been indicated to be involved in this process in some molluscan species [13,14,17,23], while gene identification and functional analyses of it have not been conducted in cuttlefish. Only one ILP (QHI00148.1) in Sepia latimanus and two ILPs (ILP-7, XP_014782891.1, and ILP-3, XP_014772597.1) in Octopus bimaculoides have been identified, but their functions remain unclear. To the best of our knowledge, there is a paucity of relevant information regarding the presence of ILPs in cuttlefish. Moreover, these reports mainly concern the cloning and characterization of genes using RACE, expression profile using qRT-PCR, and cellular localization using situ hybridization. However, few studies have examined the functional characterization of ILPs in ovarian development via RNA interference.
Thus, considering the vital function of ILPs and the limited knowledge of cephalopods, this study aimed to investigate their potential involvement in ovarian development in S. japonica. In this study, we first identified the ILP gene from S. japonica. Then, we analyzed SjILP expression in different tissues and ovarian development stages using qRT-PCR. In addition, the expressions of several ovarian-development-related genes were determined after silencing SjILP. These data may provide valuable information for a better understanding of the molecular mechanisms of ovarian development in cephalopods.

2. Materials and Methods

2.1. Animals

Healthy adult S. japonica was collected from the aquaculture station of Marine Fisheries Research Institute of Zhejiang on Xishan island (Zhoushan, Zhejiang, China). Following a previous method [20], ovarian development was classified into the oogonium production period (stage I), protoplasmic growth period (stage II), interstitial growth period (stage III), and trophoplasmic growth period (stage IV) according to the color and size of the oocyte. To evaluate the expression of SjILP during the reproductive cycle of female S. japonica, animals in four ovarian development stages were collected. Before tissue sample collection, cuttlefish were anesthetized in a solution of magnesium chloride solution (27.0 g/L). Tissues were dissected, frozen immediately in liquid nitrogen, and then stored at −80 °C until RNA extraction. In Figure 2, we lay out a roadmap for this paper.

2.2. Cloning of SjILP

Total RNA was extracted from the ovary of adult S. japonica using Trizol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s protocol. After RNA quality check using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), 10µg of total RNA was reverse-transcribed into first-strand cDNA using M-MLV reverse transcriptase (RNase H) (TaKaRa Bio Inc., Shiga, Japan).
A partial fragment of ILP was obtained from the transcriptome data of S. japonica [19]. Based on the partial fragment, specific primers were designed using the Primer 6.0 program (Table S4). The full-length cDNA of SjILP was determined using the rapid amplification of cDNA ends (RACE). PCR reaction was performed as described in Pang et al. [21]. The primers used in this study are shown in Table S1. PCR products were separated by electrophoresis on 1.5% agarose gel. Bands with the expected size were excised and cloned into the pGEM-Teasy vector (Promega) and then sequenced.

2.3. Bioinformatics Analysis

The open reading frame (ORF) of ILP cDNA was predicted using the ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder) (accessed on 10 February 2022) and checked using the BLAST program (https://blast.ncbi.nlm.nih.gov/Blast.cgi) (accessed on 10 February 2022) to confirm the identity of the sequence. The deduced amino acid sequence was subsequently used to predict the domain architecture using SMART (http://smart.embl.de/) (accessed on 10 February 2022). The signal peptide of the amino acid was examined using the SignalP 4.0 Server (http://www.cbs.dtu.dk/services/SignalP) (accessed on 10 February 2022), and the cleavage sites were predicted using Neuro pred (http://stagbeetle.animal.uiuc.edu/cgi-bin/neuropred.py) (accessed on 10 February 2022). The three-dimensional (3D) structure of SjILP was generated using I-TASSER software [24,25] and edited using Swiss-Pdb viewer (v4.1.0) [26]. Multiple sequence alignments of ILP amino acid sequences were performed using Clustal W2 software (http://www.ebi.ac.uk/Tools/msa/clustalw2/) (accessed on 10 February 2022). The phylogenetic tree was constructed using the neighbor-joining (NJ) method (1000 bootstrap replicates) in MEGA5.0 [27]. The detailed information of these sequences is provided in Supplementary Table S1.

2.4. Subcellular Localization of SjILP

To investigate the subcellular localization of SjILP, HEK293 T cells were used as the in vitro model according to our previous reports [21,22]. HEK293 cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM), supplemented with 10% FBS and 2% P/S, and maintained in a humidified incubator (37 °C with 5% CO2). The recombinant plasmids were constructed as described in a previous study. Cells were transfected with SjILP-pEGFP or pEGFP-N1 as a negative control using lipofectamine 2000 (Invitrogen) according to the manufacturer’s protocol. After transfection for 4 h, cells were fixed with 4% formaldehyde for 20 min and washed twice with PBS. Cytomembrane was stained with Dil for 20 min in the dark. Then, nuclei were stained with DAPI for 10 min at 25 °C. Subcellular localization of SjILP was observed using a digital confocal microscope (Leica TCSSP5, Germany).

2.5. Expression of SjILP in Different Tissues

Total RNA was isolated from eight tissues, namely, the brain, liver, pancreas, muscle, heart, ovary, gill, and intestine. The RNA extraction and the cDNA synthesis process were performed as described in the present study. The expression of SjILP in different tissues was analyzed using quantitative Real-Time PCR (qRT-PCR). β-actin was applied as the internal reference gene [28]. Gene-specific primers and reference genes are listed in Table S4. The system and procedure of qRT-PCR were as described in our previous study [21]. Briefly, PCR reactions were performed in a medium with a volume of 10 μL, containing 5 μL TB Green™ Premix Ex Taq™ II, 0.4 μL of each primer (10 μM), 0.4 μL cDNA template, 0.2 μL ROX Reference Dye II, and 3.6 μL ultrapure H2O. The cycling protocol was as follows: initial incubation at 94 °C for 10 min, followed by 40 cycles of 95 °C for 15 s, 60 °C for 45 s, and 72 °C for 2 min. The relative gene expression was calculated using the 2–ΔΔCT method [29]. Each sample was tested in three replicates.

2.6. Expression of SjILP in Different Ovarian Development Stages

To evaluate the expression of SjILP during the reproductive cycle of female S. japonica, tissues of the pancreas, ovary, and liver from four ovarian development stages were dissected. The RNA extraction, cDNA synthesis, and qRT-PCR process were performed as described above in the present study. Each sample was tested in three replicates.

2.7. RNA Interference In Vivo

2.7.1. Synthesis and Injection of dsRNA

Based on the complete sequence of SjILP, two dsRNAs and one negative control dsRNA were designed and synthesized using RT-PCR with primers containing the T7 polymerase promoter sequence (Table S4). The dsRNAs were dissolved in DEPC water to obtain a final concentration of 2 μmol/mL.
About 32 cuttlefish were randomly divided into four groups (n = 8 each group): two treatment groups (dsRNA1 and dsRNA2 injected) and two control groups (scramble dsRNA and normal saline injected). Then, 100 μL dsRNAs (dsRNA1 or dsRNA2) of SjILP, scramble dsRNA, or normal saline (NS) was injected according to the procedure described by Gore et al. [30]. At 48 h after injection, the ovary of each group was dissected.

2.7.2. Functional Analysis of SjILP with Gene Silencing

To evaluate the possible functions of SjILP in S. japonica reproduction, a highly effective dsRNA was selected to knockdown ILPs in females to examine their impact on the expression of ovarian-development-related genes. Four genes were selected, namely, vitellogenin (Vg1 and Vg2), cathepsin L1-like (CtsL1-like), and follistatin (FS). The expressions of these genes were detected using qRT-PCR as described above. Primers for these genes (Table S1) were designed based on sequencing information from the transcriptome data of S. japonica [19].

2.8. Statistic Analysis

Data are shown as mean ± standard deviations. All values were analyzed using one-way analysis of variance (ANOVA), and then followed by a multiple comparison in SPSS 14.0. The statistical significance was set at p < 0.05.

3. Results

3.1. Characterization of SjILP Sequence

The full-length ILP cDNA obtained from S. japonica was 926 bp. It contained a 310-bp 5′-untranslated region (UTR), a 486 bp open reading frame (ORF), and a 130 bp 3′-UTR (Figure 3). Bioinformatics analysis revealed that the ORF encoding a 161 amino acid (aa) polypeptide contained an insulin-like growth factor (IIGF) domain (aa 60–148, Figure 4A). Further analyses of this sequence identified two putative cleavage sites (arginine and lysine, KR) and eight conserved cysteine residues. The signal peptide and C-peptide of the ILP prepropeptide cleaved off to form a mature ILP comprising B and A chains (Figure 4C). These eight cysteine residues formed three inter-chain disulfide bonds (Cys48-Cys132, Cys63-Cys135, and Cys75-Cys148) and one intra-chain disulfide bond (Cys134-Cys139) (Figure 4C). In addition, its A chain presented a signature motif of MIPs: C(1X)CC(3X)C(8X)C. Based on these key features, the sequence obtained from S. japonica was termed SjILP (GenBank: MK611805).
Based on Homo sapiens insulin (PDB 7KD6) (sequence identity of 35.0%), the putative 3D structure of the mature protein was obtained, and it showed a similar structure with a TM score of 0.485 (Figure 4D). This similarity also demonstrates the conserved structure of SjILP at the tertiary level.

3.2. Multiple Alignment and Phylogenetic Analysis

The predicted mature peptide was further analyzed for sequence similarity with published sequences of insulins and ILPs. The result revealed that their primary sequences were poorly conserved, except for cysteine residues (Figure S1). BLAST analysis revealed the amino acid sequence shared a high identity with the ILP for Sepia latimanus (98.14%). However, the identity between S. japonica and other mollusk orders was relatively low (Table S2). The sequence similarity results were further reflected in a phylogenetic tree constructed with the neighbor-joining method. As shown in Figure 5, insulin and ILPs were clustered together according to phylum and formed two clades: vertebrate insulins and invertebrate ILPs. All molluskan ILPs formed a subgroup in which SjILP was most similar to the Sepia latimanus ILP.

3.3. Subcellular Localization of SjILP

HEK293T cells were transfected with pEGFP-N1/SjILP to analyze the possible subcellular localization of SjILP. In contrast to the universal distribution of green fluorescence in the EGFP control, SjILP was localized in the cytoplasm and membrane, indicating the typical features of the secretory protein of SjILP (Figure 6).

3.4. Expression Profile of SjILP Transcript in Different Tissues

As shown in Figure 7A, SjILP was expressed in all investigated tissues. However, the highest expression level was detected in the ovary, followed by the pancreas, liver, and brain. Moreover, its expression in the intestine, gill, muscle, and heart was relatively low (p < 0.05). Based on the tissue-specific expression profiles, the ovary, pancreas, and liver were used in expression detection for SjILP mRNA at different ovarian developmental stages.

3.5. Expression Profile of SjILP Transcript at Different Ovarian Developmental Stages

As shown in Figure 7B, low expression levels were detected at stage I, the expression level started to increase from stage II to stage III, the expression level reached the maximum at stage III, and then it obviously declined at stage IV.

3.6. Effects of SjILP Silencing on Ovarian Development

As shown in Figure 8, SjILP expression was effectively reduced by SjILP knockdown, with siRNA2 having the greater inhibitory effect (p < 0.05). Therefore, the dsRNA2-injected group was used as the treatment group to detect the expression of ovarian-development-related genes.
As expected, the expression levels of Vg1, Vg2, CtsL1-like, and FS in the ovary were also significantly reduced by RNAi treatment (Figure 9), further indicating the role of SjILP in the ovarian development of S. japonica.

4. Discussion

In the present study, a putative ILP gene (SjILP) was identified from the cuttlefish S. japonica, and it shared a high similarity in gene organization with other molluskan ILPs. SjILP had two proteolytic cleavage sites (KR), which appeared to be highly conserved in lepidopteran ILP precursors (Figure 3) [31]. The C-peptide is removed from these two dibasic sites to form the mature ILP in the same way as vertebrate insulin [32]. Proteolytic cleavage of the C-peptide is a key post-translational modification for attaining hormone activity [33], which is a typical feature of the insulin family [34]. Combined with the finding that SjILP was localized in the cytoplasm and the membrane (Figure 6), these results suggest that SjILP belongs to the secretory protein. As is common to other molluskan ILPs, SjILP possessed eight cysteine residues constructing four disulfide bonds, which are essential for tertiary folding (Figure 4B,C) [5,6,35]. These cysteine residues also contribute to the hydrophobic core of insulins [4,35].
Although the structural organization of SjILP was similar to that of molluskan ILPs, it showed a low sequence identity (Table S1), possibly indicating the evolutionary divergence of ILPs in mollusks. This result is consistent with most previous reports about invertebrate ILPs [4,36]. Interestingly, in the various MIPs, not all the amino acid residues that are important for the maintenance of the basic conformation (i.e., the two-chain organization, proteolytic cleavage sites, and cysteine residues) are well-conserved; they are, in fact, only substitutions of surface residues that account for the differences between these peptides [37].
Gene expression is generally related to its function. In order to investigate the function of SjILP, we analyzed its expression in various tissues. In this study, it was found to be expressed in all tissues examined, which indicates a wide distribution and multiple functions of SjILP. The highest expression level of SjILP was found in the ovary, followed by the pancreas and liver (Figure 7A). This particular expression pattern of SjILP most likely reflects the specific demands of each tissue. The highest expression level in the ovary is in line with that reported by Iwami et al. [38], who found that the bombyxin transcript was localized in the B. mori ovary. Along with the report that bombyxin induced meiosis in Bombyx ovarian cells in vitro [39], our results suggest that SjILP plays an important role in ovarian development. Analogous findings have been reported for C. gigas [5]. In most vertebrates, insulin is primarily synthesized in endocrine pancreatic β cells, and it is mainly involved in growth control and carbohydrate metabolism [40]. The involvement of ILPs in the regulation of carbohydrate metabolism has also been reported in invertebrates, such as Lymnaea stagnalis [41] and Aplysia californica [3]. In this study, the relatively greater abundance of SjILP mRNA transcript in the pancreas suggests the paracrine function of insulin-like peptides during digestion (Figure 7A) [16]. In mollusks, neurosecretory cells in the cerebral ganglia are the main source of ILPs [42,43,44]. Moreover, in these animals, ILPs are not supposed to act as neurotransmitters on target neurons in the neural ganglia [37]. In this study, SjILP was also found to be expressed in the brain of S. japonica (Figure 7A), which may indicate a possible relationship between ILPs and neuroendocrine functional control in the brain to regulate reproduction. Other tissues, such as the intestine, heart, muscle, and gill, were also potential targets of insulin-like signaling, and they might be subjected to constant growth or cellular renewal [45].
Our analysis of SjILP expression during different ovarian development stages indicated that SjILP may regulate the reproduction of S. japonica (Figure 7B). The fluctuation of the SjILP expression level correlated with ovarian development, which suggests that ILP might be involved in ovarian development by stimulating oogenesis and/or vitellogenesis [46]. Combined with the finding that the insulin receptor (CIR) is expressed in the gonadal area of oysters and showed differential expression during embryogenesis and larval developmental stages [13,47], these results suggest that ILP might regulate reproduction through binding to its receptor.
To further investigate the function of SjILP in the reproduction of S. japonica, we detected the expressions of ovarian-development-related genes (Vg1, Vg2, CtsL1-like, and FS) after knockdown of the SjILP gene in vivo. In squid, vitellogenin is the major precursor of egg yolk protein expressed predominantly in the ovary during ovary vitellogenesis/maturation [48]. Studies in insects, nematodes, and fish have indicated that ILPs positively regulate vitellogenesis [49,50,51]. Our data show that silencing the SjILP gene resulted in a decrease in Vg1 and Vg2 expression levels (Figure 9), which is in line with previous observations in Schistocerca gregaria [52]. These findings indicate that SjILP may indirectly control reproduction by affecting vitellogenesis. CtsL1 is an important member of the cysteine protease family and plays a crucial role in oocyte maturation [53]. In shrimp, Mn-CTS L1 has been verified to be involved in hydrolyzing vitellogenin before ovarian maturity to provide nutrients for the maturation of oocytes [54]. Furthermore, silencing cathepsin L using RNAi led to the accumulation of vitellogenin [55]. Here, the expression of CtsL1 mRNA markedly reduced after SjILP silencing, which suggests that SjILP possibly involves vitellogenesis in the control of reproduction in S. japonica. FS is a monomeric glycosylated protein that plays a key role in regulating folliculogenesis and ovarian development [56]. Similar to CtsL1, SjILP gene silencing also caused the down-regulation of FS in S. japonica (Figure 9). Although a direct influence of ILP on FS expression has not been proven in other species, the down-regulation of FS expression in SjILP-dsRNA-injected females suggests that SjILP possibly regulates reproduction by affecting the expression of ovarian-development-related genes in S. japonica. These findings show that an exogenous injection of SjILP dsRNA can effectively inhibit ovarian development, thereby providing a novel approach to reduce the rapid sexual maturation of female cuttlefish.

5. Conclusions

In summary, this study characterized the SjILP gene from S. japonica, and it further investigated the ovarian development mechanism. The highest expression level of SjILP mRNA was found in the ovary, and its gene expression correlated with ovarian development. In addition, the silencing of SjILP markedly down-regulated the expression of the ovarian-development-related genes. Our data add valuable information in investigating the molecular mechanism of ovarian development in cephalopods. However, to understand the regulative mechanism of SjILP in cephalopod reproduction, more detailed studies on the network of SjILP signaling will be necessary.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb44060170/s1.

Author Contributions

Conceptualization, L.L. and Z.L.; data curation, S.Z. and Y.Z.; writing—original draft preparation, L.L. and S.Z.; Formal analysis, C.Y.; writing—review and editing, L.G. and B.L.; project administration, C.Y., S.Z. and Y.Z.; funding acquisition, L.L. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Chinese National Natural Science Foundation (41406138, and 41976121), Natural Science Foundation of Zhejiang Province (LY130190001), and the Open Project of Zhejiang Provincial Top Key Discipline of Biotechnology (KF2016008).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors are grateful to the aquaculture station of Marine Fisheries Research Institute of Zhejiang for providing the animals used in the experiments.

Conflicts of Interest

All authors declare no conflict of interest.

References

  1. Okamoto, N. Molluscan insulin-related peptides. In Handbook of Hormones; Academic Press: Cambridge, MA, USA, 2021; pp. 681–683. [Google Scholar]
  2. Smit, A.B.; Vreugdenhil, E.; Ebberink, R.H.; Geraerts, W.P.; Klootwijk, J.; Joosse, J. Growth-controlling molluscan neurons produce the precursor of an insulin-related peptide. Nature 1988, 331, 535–538. [Google Scholar] [CrossRef] [PubMed]
  3. Floyd, P.D.; Li, L.J.; Rubakhin, S.S.; Sweedler, J.V.; Horn, C.C.; Kupfermann, I.; Alexeeva, V.; Ellis, T.A.; Dembrow, N.C. Insulin prohormone processing, distribution, and relation to metabolism in Aplysia Californica. J. Neurosci. 1999, 19, 7732–7741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hamano, K.; Awaji, M.; Usuki, H. cDNA structure of an insulin-related peptide in the Pacific oyster and seasonal changes in the gene expression. J. Endocrinol. 2005, 187, 55–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cherif-Feildel, M.; Berthelin, C.H.; Adeline, B.; Rivière, G.; Favrel, P.; Kellner, K. Molecular evolution and functional characterisation of insulin related peptides in molluscs: Contributions of Crassostrea gigas genomic and transcriptomic-wide screening. Gen. Comp. Endocrinol. 2019, 271, 15–29. [Google Scholar] [CrossRef]
  6. Smit, A.B.; Spijker, S.; Van Minnen, J.; Burke, J.F.; De Winter, F.; Van Elk, R.; Geraerts, W.P. Expression and characterization of molluscan insulin-related peptide VII from the mollusc Lymnaea stagnalis. Neuroscience 1996, 70, 589–596. [Google Scholar] [CrossRef]
  7. Zhang, H.; He, M.X. The role of a new insulin-like peptide in the pearl oyster pinctada fucata martensii. Sci. Rep. 2020, 10, 433. [Google Scholar] [CrossRef] [Green Version]
  8. Plisetskaya, E.V.K.; Kazakov, V.K.; Soltitskaya, L.; Leibson, L.G. Insulin producing cells in the gut of freshwater bivalve Molluscs Anodonfa Cygnaea and Unio pictorum and the role of insulin in the regulation of their carbohydrate metabolism. Gen. Camp. Endocrinol. 1978, 35, 133–145. [Google Scholar] [CrossRef]
  9. Abdraba, A.M.; Saleuddin, A.S.M. Localization and immunological characterization of insulin-like peptide(s) in the land snail Otala lactea (Mollusca:Pulmonata). Can. J. Zool. 2000, 78, 1515–1526. [Google Scholar] [CrossRef]
  10. Li, Y.J.; Fu, H.R.; Zhang, F.Q.; Ren, L.T.; Tian, J.; Li, Q.; Liu, S.K. Identification, characterization, and expression profiles of insulin-like peptides suggest their critical roles in growth regulation of the pacific oyster, crassostrea gigas. Gene 2020, 769, 145244. [Google Scholar] [CrossRef]
  11. Niu, D.; Wang, F.; Zhao, H.; Wang, Z.; Xie, S.; Li, J. Identification, expression, and innate immune responses of two insulin-like peptide genes in the razor clam sinonovacula constricta. Fish Shellfish. Immun. 2016, 51, 401–404. [Google Scholar] [CrossRef]
  12. Zhang, L.H.; Wang, G.D. Insulin-like peptides in invertebrates and their signaling pathways-take insects. Biotechnol. Bull. 2014, 20, 33–42. [Google Scholar]
  13. Gricourt, L.; Mathieu, M.; Kellner, K. An insulin-like system involved in the control of pacific oyster Crassostrea gigas reproduction: hrIGF-1 effect on germinal cell proliferation and maturation associated with expression of an homologous insulin receptor-related receptor. Aquaculture 2006, 251, 85–98. [Google Scholar] [CrossRef]
  14. Sonetti, D.; Van Heumen, W.R.A.; Roubos, E.W. Light-and electro-microscopic immunocytochemistry of a molluscan insulin-related peptide in the central nervous system of Planorbarius corneus. Cell Tissue Res. 1992, 267, 473–481. [Google Scholar] [CrossRef] [PubMed]
  15. Geraerts, W.P.M. The role of the lateral lobes in the control of growth and reproduction in the hermaphrodite freshwater snail Lymnaea stagnalis. Gen. Comp. Endocrinol 1976, 29, 97–108. [Google Scholar] [CrossRef]
  16. Monnier, Z.; Bride, M. In vitro effects of methionine-en-kephalin, somatostatin and insulin on cultured gonadal cells of the snail Helix aspersa. Experientia 1995, 51, 824–829. [Google Scholar] [CrossRef]
  17. Brown, M.R.; Clark, K.D.; Gulia, M.; Zhao, Z.W.; Garczynski, S.F.; Crim, J.W.; Suderman, R.J.; Strand, M.R. An insulin-like peptide regulates egg maturation and metabolism in the mosquito Aedes aegypti. Proc. Natl. Acad. Sci. USA 2008, 105, 5716–5721. [Google Scholar] [CrossRef] [Green Version]
  18. Dong, Z.Y.; Wu, C.W.; Ye, D.Y. A Fluorescent Method for Marking the Cuttlefish, Sepiella mainddroni de Rochebrune. J. Zhejiang Ocean. Univ. (Nat. Sci.) 2010, 29, 120–128. [Google Scholar]
  19. Lü, Z.M.; Liu, W.; Liu, L.Q.; Shi, H.L.; Ping, H.L.; Wang, T.M.; Chi, C.F.; Wu, C.W.; Chen, C.H.; Shen, K.N.; et al. De novo assembly and comparison of the ovarian transcriptomes of the common Chinese cuttlefish (Sepiella japonica) with different gonadal development. Genom. Data 2016, 7, 155–158. [Google Scholar] [CrossRef]
  20. Lü, Z.M.; Liu, W.; Liu, L.Q.; Wang, T.M.; Shi, H.L.; Ping, H.L.; Chi, C.F.; Yang, J.W.; Wu, C.W. Cloning, characterization, and expression profile of estrogen receptor in common Chinese cuttlefish, Sepiella japonica. J. Exp. Zool. Part A 2016, 325A, 181–193. [Google Scholar] [CrossRef]
  21. Pang, Z.; Lü, Z.M.; Wang, M.T.; Gong, L.; Liu, B.J.; Liu, L.Q. Characterization, relative abundances of mrna transcripts, and subcellular localization of two forms of membrane progestin receptors (MPRs) in the common chinese cuttlefish, sepiella japonica. Anim. Reprod. Sci. 2019, 208, 106107. [Google Scholar] [CrossRef]
  22. Lü, Z.M.; Zhu, K.H.; Pang, Z.; Liu, L.Q.; Jiang, H.; Liu, B.J.; Shi, H.L.; Ping, H.L.; Chi, C.F.; Gong, L. Identification, characterization and mrna transcript abundance profiles of estrogen related receptor (ERR) in Sepiella japonica imply its possible involvement in female reproduction. Anim. Reprod. Sci. 2019, 211, 106231. [Google Scholar] [CrossRef] [PubMed]
  23. Gomot, P. Studies on the control of spermatogenic DNA synthesis by the mesocerebrum in the snail Helix aspersa. Cell. Mol. Neurobiol. 1993, 13, 517–527. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, C.X.; Freddolino, P.L.; Zhang, Y. COFACTOR: Improved protein function prediction by combining structure, sequence and protein-protein interaction information. Nucleic Acids Res. 2017, 45, 291–299. [Google Scholar] [CrossRef]
  25. Zheng, W.; Zhang, C.X.; Li, Y.; Pearce, R.; Bell, E.W.; Zhang, Y. Folding non-homology proteins by coupling deep-learning contact maps with I-TASSER assembly simulations. Cell Rep. Methods 2021, 1, 100014. [Google Scholar] [CrossRef]
  26. Guex, N.; Peitsch, M.C. SWISS-MODEL and the Swiss-PdbViewer: An environment for comparative protein modeling. Electrophoresis 1997, 18, 2714–2723. [Google Scholar] [CrossRef]
  27. Kumar, S.; Dudley, J.; Nei, M.; Tamura, K. MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008, 9, 299–306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Olsvik, P.A.; Lie, K.K.; Jordal, A.E.O.; Nilsen, T.O.; Hordvik, I. Evaluation of potential reference genes in real-time RT-PCR studies of Atlantic salmon. BMC Mol. Biol. 2005, 6, 1–9. [Google Scholar] [CrossRef] [Green Version]
  29. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-DDCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  30. Gore, S.R.; Harms, C.A.; Kukanich, B.; Forsythe, J.; Lewbart, G.A.; Papich, M.G. Enrofloxacin pharmacokinetics in the European cuttlefish, Sepia officinalis, after a single iv injection and bath administration. J. Vet. Pharmacol. Ther. 2005, 28, 433–439. [Google Scholar] [CrossRef]
  31. Kimura-Kawakami, M.; Iwami, M.; Kawakami, A.; Nagasawa, H.; Suzuki, A.; Ishizaki, H. Structure and expression of bombyxin-related peptide genes of the moth Samia cynthia ricini. Gen. Comp. Endocr. 1992, 86, 257–268. [Google Scholar] [CrossRef]
  32. Krieger, M.J.B.; Jahan, N.; Riehle, M.A.; Cao, C.; Brown, M.R. Molecular characterization of insulin-like peptide genes and their expression in the african malaria mosquito, Anopheles gambiae. Insect Mol. Biol. 2010, 13, 305–315. [Google Scholar] [CrossRef] [PubMed]
  33. Okuno, A.; Hasegawa, Y.; Nishiyama, M.; Ohira, T.; Ko, R.; Kurihara, M.; Matsumoto, S.; Nagasawa, H. Preparation of an active recombinant peptide of crustacean androgenic gland hormone. Peptides 2002, 23, 567–572. [Google Scholar] [CrossRef]
  34. Mareddy, V.R.; Rosen, O.; Thaggard, H.B.; Manor, R.; Kuballa, A.V.; Aflalo, E.D.; Sagi, A.; Paterson, B.; Elizur, A. Isolation and characterization of the complete cdna sequence encoding a putative insulin-like peptide from the androgenic gland of penaeus monodon. Aquaculture 2011, 318, 364–370. [Google Scholar] [CrossRef]
  35. Velde, S.; Badisco, L.; Claeys, I.; Verleyen, P.; Chen, X.; Bosch, L.V.; Broeck, J.V.; Smagghe, G. Insulin-like peptides in Spodoptera littoralis (Lepidoptera): Detection, localization and identification. Gen. Comp. Endocr. 2007, 153, 72–79. [Google Scholar] [CrossRef] [PubMed]
  36. Ma, K.Y.; Lin, J.Y.; Guo, S.Z.; Chen, Y.; Li, J.L.; Qiu, G.F. Molecular characterization and expression analysis of an insulin-like gene from the androgenic gland of the oriental river prawn, macrobrachium nipponense. Gen. Comp. Endocr. 2013, 185, 90–96. [Google Scholar] [CrossRef]
  37. Smit, A.B.; Kesteren, R.E.V.; Li, K.W.; Minnen, J.V.; Spijker, S.; Heerikhuizen, H.V.; Geraerts, W.P.M. Towards understanding the role of insulin in the brain: Lessons from insulin-related signaling systems in the invertebrate brain. Prog. Neurobiol. 1998, 54, 35–54. [Google Scholar] [CrossRef]
  38. Iwami, M.; Furuya, I.; Kataoka, H. Bombyxin-related peptides: cDNA structure and expression in the brain of the hornworm Agrius convolvuli. Insect. Biochem. Mol. Biol. 1996, 26, 25–32. [Google Scholar] [CrossRef]
  39. Orikasa, C.; Yamauchi, H.; Nagasawa, H.; Suzuki, A.; Nagata, M. Induction of oocyte-nurse cell differentiation in the ovary by the brain during the initial stage of oogenesis in the silkworm, Bombyx mori (Lepidoptera: Bombycidae). Appl. Entomol. Zool. 1993, 28, 303–311. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, S.; Luo, X.N.; Zhang, S.H.; Yin, C.; Dou, Y.X.; Cai, X.P. Identification of putative insulin-like peptides and components of insulin signaling pathways in parasitic platyhelminths by the use of genome-wide screening. Febs J. 2014, 2818, 77–893. [Google Scholar] [CrossRef]
  41. Geraerts, W.P.M. Neurohormonal control of growth and carbohydrate metabolism by the light green cells in Lymnaea stagnalis. Gen. Comp. Endocrinol. 1992, 86, 433–444. [Google Scholar] [CrossRef]
  42. Meester, I.; Ramkema, M.D.; Van Minnen, J.; Boer, H.H. Differential expression of four genes encoding molluscan insu-lin-related peptides in the central nervous system of the pond snail Lymmea stagnalis. Cell Tissue Res. 1992, 269, 183–188. [Google Scholar] [CrossRef] [PubMed]
  43. Kellner-Cousin, K.M.E.; Mathieu, M. Identification of insulin-like peptides in cerebral ganglia neurosecretory cells in the mussel Mytilus Edulis. Tissue Cell. 1994, 26, 891–899. [Google Scholar] [CrossRef]
  44. Danton, E.; Kiyomoto, M.; Komaru, A.; Wada, K.T.; Awaji, M.; Mathieu, M. Comparative analysis of storage tissue and insulin-like neurosecretion in diploid and triploid mussels Mytilus galloprovincialis LMK in relation to their gametogenesis cycle. Invertebr. Reprod. Dev. 1996, 29, 37–46. [Google Scholar] [CrossRef]
  45. Gómez-Mendikute, A.; Elizondo, M.; Venier, P.; Cajaraville, M.P. Characterization of mussel gill cells in vivo and in vitro. Cell. Tissue Res. 2005, 321, 131–140. [Google Scholar] [CrossRef] [PubMed]
  46. Richard, D.S.; Rybczynski, R.; Wilson, T.G.; Wang, Y.; Wayne, M.L.; Zhou, Y.; Partridge, L.; Harshman, L.G. Insulin signaling is necessary for vitellogenesis in Drosophila melanogaster independent of the roles of juvenile hormone and ecdysteroids: Female sterility of the chico1 insulin signaling mutation is autonomous to the ovary. J. Insect Physiol. 2005, 51, 455–464. [Google Scholar] [CrossRef] [Green Version]
  47. Gricourt, L.; Bonnec, G.; Boujard, D.; Mathieu, M.; Kellner, K. Insulin-like system and growth regulation in the Pacific oyster Crassostrea gigas: hrIGF-1 effect on protein synthesis of mantle edge cells and expression of an homologous insulin receptor-related receptor. Gen. Comp. Endocrinol. 2003, 134, 44–56. [Google Scholar] [CrossRef]
  48. Chen, C.; Li, H.W.; Ku, W.L.; Lin, C.J.; Chang, C.F.; Wu, G.C. Two distinct vitellogenin genes are similar in function and expression in the bigfin reef squid sepioteuthis lessoniana. Biol. Reprod. 2018, 5, 5. [Google Scholar] [CrossRef]
  49. Sheng, Z.; Xu, J.; Bai, H.; Zhu, F.; Palli, S.R. Juvenile hormone regulates vitellogenin gene expression through insulin-like peptide signaling pathway in the red flour beetle, tribolium castaneum. J. Biol. Chem. 2011, 286, 41924–41936. [Google Scholar] [CrossRef] [Green Version]
  50. Murphy, C.T.; Mccarroll, S.A.; Bargmann, C.I.; Fraser, A.; Kamath, R.S.; Ahringer, J.; Hao, L.; Cynthia, K. Genes that act downstream of daf-16 to influence the lifespan of caenorhabditis elegans. Nature 2003, 424, 277. [Google Scholar] [CrossRef]
  51. Wuertz, S.; Nitsche, A.; Jastroch, M.; Gessner, J.; Klingenspor, M.; Kirschbaum, F.; Kloas, W. The role of the IGF-i system for vitellogenesis in maturing female sterlet, 1758. Acipenser ruthenus linnaeus. Gen. Comp. Endocrinol. 2007, 150, 140–150. [Google Scholar] [CrossRef]
  52. Badisco, L.; Marchal, E.; Wielendaele, P.V.; Verlinden, H.; Vleugels, R.; Broeck, J.V. RNA interference of insulin-related peptide and neuroparsins affects vitellogenesis in the desert locust schistocerca gregaria. Peptides 2011, 32, 573–580. [Google Scholar] [CrossRef] [PubMed]
  53. Tingaud-Sequeira, A.; Cerdà, J. Phylogenetic relationships and gene expression pattern of three different cathepsin L. (Ctsl) isoforms in zebrafish: Ctsla is the putative yolk processing enzyme. Gene 2007, 386, 98–106. [Google Scholar] [CrossRef] [PubMed]
  54. Zhu, J.; Fu, H.; Qiao, H.; Jin, S.; Xiong, Y. Expression and functional analysis of cathepsin L1 in ovarian development of the oriental river prawn, macrobrachium nipponense. Aquacult. Rep. 2021, 20, 100724. [Google Scholar] [CrossRef]
  55. Britton, C.; Murray, L. Cathepsin l protease (CPL-1) is essential for yolk processing during embryogenesis in caenorhabditis elegans. J. Cell. Sci. 2004, 117, 5133–5143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Kashimada, K.; Pelosi, E.; Chen, H.; Schlessinger, D.; Wilhelm, D.; Koopman, P. FOXL2 and BMP2 act cooperatively to regulate follistatin gene expression during ovarian development. Endocrinology 2011, 152, 272–280. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Proposed ILP signal transduction pathway in S. japonica. Abbreviations: IRS, insulin receptor substrate; PI3K, phosphatidylinositol-2,4,5-triphosphate; PIP2, phosphatidylinositol 4,5-bisphosphate; PIP3, phosphatidylinositol-3,4,5-triphosphate; PKB, protein kinase B; mTOR, mammalian target of rapamycin; S6K, ribosomal S6 protein kinase.
Figure 1. Proposed ILP signal transduction pathway in S. japonica. Abbreviations: IRS, insulin receptor substrate; PI3K, phosphatidylinositol-2,4,5-triphosphate; PIP2, phosphatidylinositol 4,5-bisphosphate; PIP3, phosphatidylinositol-3,4,5-triphosphate; PKB, protein kinase B; mTOR, mammalian target of rapamycin; S6K, ribosomal S6 protein kinase.
Cimb 44 00170 g001
Figure 2. Roadmap of the experimental program.
Figure 2. Roadmap of the experimental program.
Cimb 44 00170 g002
Figure 3. The nucleotide and deduced amino acid sequences of S. japonica ILP. The amino acid sequence of ORF consisted of a signal peptide (marked with yellow squares), B chain (marked with blue squares), C-peptide (marked with green squares), and A chain (marked with pink squares). Asterisk marks the stop codon. The putative cleavage sites and eight conserved cysteine residues are marked with red and black boxes, respectively. An insulin-like growth factor (IlGF) domain (aa 60–148) is indicated by a gray background.
Figure 3. The nucleotide and deduced amino acid sequences of S. japonica ILP. The amino acid sequence of ORF consisted of a signal peptide (marked with yellow squares), B chain (marked with blue squares), C-peptide (marked with green squares), and A chain (marked with pink squares). Asterisk marks the stop codon. The putative cleavage sites and eight conserved cysteine residues are marked with red and black boxes, respectively. An insulin-like growth factor (IlGF) domain (aa 60–148) is indicated by a gray background.
Cimb 44 00170 g003
Figure 4. SjILP putative structure. (A) Domain architecture of prepro-SjILP. SMART-defined architecture of SjILP: a putative signal peptide followed by an insulin-like growth factor (IIGF) domain (aa 60–148). (B) Linear model of the 161 aa SjILP pre-prohormone, composed of a signal peptide of 42 aa (yellow), B chain of 72 aa (blue), C-peptide of 9 aa (green), and A chain of 38 aa (pink); eight cysteine residue positions are marked. (C) Linear model of mature SjILP. Disulfide bridges between C48-C132, C63-C135, C75-C148, and C134-C139 allow functional tertiary structure. (D) Suggested 3D model of SjILP, structurally based on Homo sapiens insulin.
Figure 4. SjILP putative structure. (A) Domain architecture of prepro-SjILP. SMART-defined architecture of SjILP: a putative signal peptide followed by an insulin-like growth factor (IIGF) domain (aa 60–148). (B) Linear model of the 161 aa SjILP pre-prohormone, composed of a signal peptide of 42 aa (yellow), B chain of 72 aa (blue), C-peptide of 9 aa (green), and A chain of 38 aa (pink); eight cysteine residue positions are marked. (C) Linear model of mature SjILP. Disulfide bridges between C48-C132, C63-C135, C75-C148, and C134-C139 allow functional tertiary structure. (D) Suggested 3D model of SjILP, structurally based on Homo sapiens insulin.
Cimb 44 00170 g004
Figure 5. Phylogenetic tree based on amino acid sequences of mollusk ILPs and vertebrate INSs deposited in GeneBank. Neighbor-joining phylogram was constructed with 1000 bootstrap replicates in MEGA 6. Sepiella japonica ILP is indicated by a triangle.
Figure 5. Phylogenetic tree based on amino acid sequences of mollusk ILPs and vertebrate INSs deposited in GeneBank. Neighbor-joining phylogram was constructed with 1000 bootstrap replicates in MEGA 6. Sepiella japonica ILP is indicated by a triangle.
Cimb 44 00170 g005
Figure 6. SjILP-EGFP fusion protein in HEK293 cells. (A) Nucleus of cells stained blue with DAPI. (B) Plasma membranes of cells stained red with DiI. (C) Cytoplasm and membrane of cells stained green with ERR-EGFP. (D) ILP was localized in the cytoplasm and the membrane of HEK293 cells; nucleus and the plasma membrane are indicated with white arrows in (AD).
Figure 6. SjILP-EGFP fusion protein in HEK293 cells. (A) Nucleus of cells stained blue with DAPI. (B) Plasma membranes of cells stained red with DiI. (C) Cytoplasm and membrane of cells stained green with ERR-EGFP. (D) ILP was localized in the cytoplasm and the membrane of HEK293 cells; nucleus and the plasma membrane are indicated with white arrows in (AD).
Cimb 44 00170 g006
Figure 7. The expression of SjILP in different tissues (A) and ovarian development stages (B). Data are shown as means ± SD. Significant difference in SjILP mRNA levels is represented by different lowercase letters (p < 0.05).
Figure 7. The expression of SjILP in different tissues (A) and ovarian development stages (B). Data are shown as means ± SD. Significant difference in SjILP mRNA levels is represented by different lowercase letters (p < 0.05).
Cimb 44 00170 g007
Figure 8. Effect of SjILP knockdown on the expression of SjILP. Data are expressed as mean ± SD (n = 8). Significant difference among different groups is indicated by the different lowercase letters (p < 0.05).
Figure 8. Effect of SjILP knockdown on the expression of SjILP. Data are expressed as mean ± SD (n = 8). Significant difference among different groups is indicated by the different lowercase letters (p < 0.05).
Cimb 44 00170 g008
Figure 9. Effect of SjILP knockdown on the expression of reproduction-related genes. Vg1: vitellogenin1, Vg2: vitellogenin2, CtsL1-like: cathepsin L1-like, and FS: follistatin. Data are expressed as mean ± SD (n = 8). Significant differences (p < 0.05) are indicated by an asterisk (*).
Figure 9. Effect of SjILP knockdown on the expression of reproduction-related genes. Vg1: vitellogenin1, Vg2: vitellogenin2, CtsL1-like: cathepsin L1-like, and FS: follistatin. Data are expressed as mean ± SD (n = 8). Significant differences (p < 0.05) are indicated by an asterisk (*).
Cimb 44 00170 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lü, Z.; Yao, C.; Zhao, S.; Zhang, Y.; Gong, L.; Liu, B.; Liu, L. Characterization of Insulin-like Peptide (ILP) and Its Potential Role in Ovarian Development of the Cuttlefish Sepiella japonica. Curr. Issues Mol. Biol. 2022, 44, 2490-2504. https://doi.org/10.3390/cimb44060170

AMA Style

Lü Z, Yao C, Zhao S, Zhang Y, Gong L, Liu B, Liu L. Characterization of Insulin-like Peptide (ILP) and Its Potential Role in Ovarian Development of the Cuttlefish Sepiella japonica. Current Issues in Molecular Biology. 2022; 44(6):2490-2504. https://doi.org/10.3390/cimb44060170

Chicago/Turabian Style

Lü, Zhenming, Chenghao Yao, Shijie Zhao, Yao Zhang, Li Gong, Bingjian Liu, and Liqin Liu. 2022. "Characterization of Insulin-like Peptide (ILP) and Its Potential Role in Ovarian Development of the Cuttlefish Sepiella japonica" Current Issues in Molecular Biology 44, no. 6: 2490-2504. https://doi.org/10.3390/cimb44060170

APA Style

Lü, Z., Yao, C., Zhao, S., Zhang, Y., Gong, L., Liu, B., & Liu, L. (2022). Characterization of Insulin-like Peptide (ILP) and Its Potential Role in Ovarian Development of the Cuttlefish Sepiella japonica. Current Issues in Molecular Biology, 44(6), 2490-2504. https://doi.org/10.3390/cimb44060170

Article Metrics

Back to TopTop