<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mar. Sci.</journal-id>
<journal-title>Frontiers in Marine Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mar. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-7745</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmars.2021.666763</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Marine Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>RNAi Reveals Role of Insulin-Like Androgenic Gland Hormone 2 (IAG2) in Sexual Differentiation and Growth in Hermaphrodite Shrimp</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Liu</surname> <given-names>Fang</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1253535/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Shi</surname> <given-names>Wenyuan</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Ye</surname> <given-names>Haihui</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/489038/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Liu</surname> <given-names>An</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhu</surname> <given-names>Zhihuang</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>College of Ocean and Earth Sciences, Xiamen University</institution>, <addr-line>Xiamen</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>College of Fisheries, Jimei University</institution>, <addr-line>Xiamen</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>Fisheries Research Institute of Fujian</institution>, <addr-line>Xiamen</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Tomomi Watanabe-Asaka, Tohoku Medical and Pharmaceutical University, Japan</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Amir Sagi, Ben-Gurion University of the Negev, Israel; Hidekazu Katayama, Tokai University, Japan</p></fn>
<corresp id="c001">&#x002A;Correspondence: Haihui Ye, <email>hhye@jmu.edu.cn</email></corresp>
<fn fn-type="other" id="fn004"><p>This article was submitted to Aquatic Physiology, a section of the journal Frontiers in Marine Science</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>06</day>
<month>04</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>8</volume>
<elocation-id>666763</elocation-id>
<history>
<date date-type="received">
<day>11</day>
<month>02</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>16</day>
<month>03</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2021 Liu, Shi, Ye, Liu and Zhu.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Liu, Shi, Ye, Liu and Zhu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Insulin-like androgenic gland hormone (IAG) is the most widely known hormone that regulates sexual differentiation in crustaceans. Previously, a transcriptome study described two transcripts of IAGs (<italic>Lvit-IAG1</italic> and <italic>Lvit-IAG2</italic>) in the peppermint shrimp <italic>Lysmata vittata</italic>, a species characterized by a rare reproductive system of protandric simultaneous hermaphroditism (PSH). Herein, we explored the function of <italic>Lvit</italic>-IAG2 via RNA interference (RNAi) experiments, and then compared the functional differences between the two IAGs. We demonstrated that <italic>Lvit</italic>-IAG2 positively regulated the expression of crustacean hyperglycemic hormone (<italic>Lvit-CHH</italic>) in eyestalk ganglion but exhibited no significant effect on the expression of gonad-inhibiting hormone (<italic>Lvit-GIHs</italic>) and crustacean female sex hormone (<italic>Lvit-CFSHs</italic>). Besides, <italic>Lvit-IAG2</italic> gene knockdown caused degeneration in appendices masculinae (AM) and suppressed germ cells at the secondary spermatocyte stage. Moreover, silencing the <italic>Lvit-IAG2</italic> gene impeded ovarian development, including smaller previtellogenic oocytes, and lower expression of <italic>vitellogenin</italic> (<italic>Lvit-Vg</italic>) gene in hepatopancreas and <italic>vitellogenin receptor</italic> (<italic>Lvit-VgR</italic>) gene in the ovarian region. Notably, silencing the <italic>Lvit-IAG2</italic> gene impeded individual growth of the species. Collectively, findings from this study demonstrate that <italic>Lvit</italic>-IAG2 and <italic>Lvit</italic>-IAG1 coordinatively function to modulate sexual differentiation in <italic>L. vittata</italic>; meanwhile, <italic>Lvit</italic>-IAG2 stimulates the growth of the PSH species.</p>
</abstract>
<kwd-group>
<kwd>sexual differentiation</kwd>
<kwd>IAG</kwd>
<kwd>PSH</kwd>
<kwd>reproductive endocrine</kwd>
<kwd>crustacean</kwd>
</kwd-group>
<counts>
<fig-count count="8"/>
<table-count count="3"/>
<equation-count count="0"/>
<ref-count count="49"/>
<page-count count="15"/>
<word-count count="0"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1">
<title>Introduction</title>
<p>Gonochorism is the most common reproductive strategy in decapod crustacean species (<xref ref-type="bibr" rid="B18">Juchault, 1999</xref>). Males and females have significantly different reproductive systems, as well as external features (<xref ref-type="bibr" rid="B35">Nagaraju, 2011</xref>). For instance, the reproductive system of male American lobster <italic>Homarus americanus</italic> consists of the paired testes, vas deferens, spermiducts, and gonopores located on the fifth pair of walking legs. Besides, the specialized first pair of pleopods serves as a gonopod during mating (<xref ref-type="bibr" rid="B8">Comeau and Benhalima, 2018b</xref>). However, the female reproductive system of <italic>H. americanus</italic> consists of the paired ovaries, oviducts, and gonopores located on the third pair of walking legs. The specialized organ, seminal receptacle, arises on the thoracic ventral side between the fourth and the fifth pairs of walking legs (<xref ref-type="bibr" rid="B7">Comeau and Benhalima, 2018a</xref>). Interestingly, a review of the literature shows that Caridean shrimps exhibit several other protandry sexual systems apart from gonochorism (<xref ref-type="bibr" rid="B3">Bauer, 2000</xref>). Among them, strictly sequential protandric hermaphroditism (SPH) and protandric simultaneous hermaphroditism (PSH) are well-studied. For example, Northern spot shrimp <italic>Pandalus platyceros</italic>, a SPH species, first mature as functional males with hermaphrodite gonads, and then go to a transitional phase followed by a functional female phase (<xref ref-type="bibr" rid="B25">Levy et al., 2020b</xref>). Meanwhile, all the <italic>Lysmata</italic> species with the PSH sexual system, eventually acquire reproductive functions for both males and females following the male phase and a short transitional phase (<xref ref-type="bibr" rid="B3">Bauer, 2000</xref>; <xref ref-type="bibr" rid="B5">Chen et al., 2019</xref>).</p>
<p>As with insect, sexual differentiation is controlled by key reproductive hormones in decapods. It is generally believed that insulin-like androgenic gland hormone (IAG) regulates male sexual differentiation in dioecious decapods (<xref ref-type="bibr" rid="B43">Ventura et al., 2009</xref>). In decapod species, IAG was first identified in the androgenic gland (AG) of the red-claw crayfish <italic>Cherax quadricarinatus</italic> in 2007 (<xref ref-type="bibr" rid="B30">Manor et al., 2007</xref>). Since then, numerous studies have found IAG in dioecious (<xref ref-type="bibr" rid="B43">Ventura et al., 2009</xref>; <xref ref-type="bibr" rid="B39">Rosen et al., 2010</xref>; <xref ref-type="bibr" rid="B44">Ventura et al., 2012</xref>; <xref ref-type="bibr" rid="B6">Chung, 2014</xref>; <xref ref-type="bibr" rid="B16">Huang et al., 2014</xref>), parthenogenetic (<xref ref-type="bibr" rid="B23">Levy et al., 2017</xref>), and hermaphroditic (<xref ref-type="bibr" rid="B47">Zhang et al., 2017</xref>; <xref ref-type="bibr" rid="B24">Levy et al., 2020a</xref>; <xref ref-type="bibr" rid="B29">Liu et al., 2020</xref>) species. Furthermore, its critical role in sexual differentiation has been properly studied via loss-of-function experiments (<xref ref-type="bibr" rid="B44">Ventura et al., 2012</xref>; <xref ref-type="bibr" rid="B49">Zmora and Chung, 2014</xref>; <xref ref-type="bibr" rid="B37">Priyadarshi et al., 2017</xref>; <xref ref-type="bibr" rid="B11">Fu et al., 2020</xref>). For instance, in the giant freshwater prawn, <italic>Macrobrachium rosenbergii</italic>, IAG gene silencing induced degeneration of male primary and secondary characteristics (<xref ref-type="bibr" rid="B43">Ventura et al., 2009</xref>), and feminization in young males (<xref ref-type="bibr" rid="B44">Ventura et al., 2012</xref>). Silencing IAG gene could also feminize male-related phenotypes in the intersex red-claw crayfish <italic>C. quadricarinatus</italic> (<xref ref-type="bibr" rid="B39">Rosen et al., 2010</xref>) and the male Chinese mitten crab <italic>Eriocheir sinensis</italic> (<xref ref-type="bibr" rid="B11">Fu et al., 2020</xref>). Moreover, augmentation of IAG hormone transcripts in adult <italic>M. rosenbergii</italic> males resulted in significantly higher (orange claw) OC to (blue claw) BC transformations and vice versa (<xref ref-type="bibr" rid="B37">Priyadarshi et al., 2017</xref>). Considering its universal and pivotal role as a master regulator of crustacean sexual development, IAG is regard as the sexual &#x201C;IAG-switch&#x201D; in decapod crustacean (<xref ref-type="bibr" rid="B22">Levy and Sagi, 2020</xref>).</p>
<p>Crustacean female sex hormone (CFSH) is a neurohormone that is closely related to female sexual differentiation in decapods. It was first isolated from the eyestalk ganglion in the Atlantic blue crab <italic>Callinectes sapidus</italic>, and was found to be pivotal in the development of the female&#x2019;s mating and egg brooding systems, including the gonopores and ovigerous setae (<xref ref-type="bibr" rid="B49">Zmora and Chung, 2014</xref>). Similar functions have also been established in the mud crab <italic>Scylla paramamosain</italic> (<xref ref-type="bibr" rid="B17">Jiang et al., 2020</xref>). Besides, CFSH was also reported as an inhibitor of IAG, as recombinant <italic>Sp</italic>-CFSH protein significantly suppressed <italic>Sp-IAG</italic> expression in the AG explants of the mud crab <italic>S. paramamosain in vitro</italic> (<xref ref-type="bibr" rid="B28">Liu et al., 2017</xref>).</p>
<p>In crustacean males, previous studies have revealed that the X-organ-sinus gland complex (XO-SG) in the eyestalk ganglion, a major source of neuropeptides, modulates reproductive functions via the eyestalk-AG-testicular axis (<xref ref-type="bibr" rid="B20">Khalaila et al., 2002</xref>; <xref ref-type="bibr" rid="B35">Nagaraju, 2011</xref>). Besides, its inhibitory effects on AG development and IAG expression have been intensively explored through eyestalk ablation. The active components suppressing AG development have been identified in XO-SG. Among them, researchers have comprehensively characterized the crustacean hyperglycemic hormone (CHH) superfamily neuropeptides as well as CFSH mentioned above. Most members of the CHH superfamily exerted inhibitory effects on either AG development or <italic>IAG</italic> expression (<xref ref-type="bibr" rid="B10">Das et al., 2015</xref>; <xref ref-type="bibr" rid="B26">Li et al., 2015</xref>; <xref ref-type="bibr" rid="B13">Guo et al., 2019</xref>). Following the silencing of <italic>gonad-inhibiting hormone</italic> (<italic>GIH</italic>), <italic>molt-inhibiting hormone</italic> (<italic>MIH</italic>), and <italic>CHH</italic> gene, <italic>IAG</italic> transcription level increased significantly in the black tiger prawn <italic>Penaeus monodon</italic> (<xref ref-type="bibr" rid="B10">Das et al., 2015</xref>), the oriental river prawn <italic>Macrobrachium nipponense</italic> (<xref ref-type="bibr" rid="B26">Li et al., 2015</xref>), and the Pacific white shrimp <italic>Penaeus</italic> (<italic>Litopenaeus</italic>) <italic>vannamei</italic> (<xref ref-type="bibr" rid="B13">Guo et al., 2019</xref>).</p>
<p>To date, however, there are rare reports elaborating sexual differentiation mechanism mediated by IAG in protandric crustaceans. In the SPH shrimp <italic>P. platyceros</italic>, <italic>Pnp-IAG</italic> knockdown elevated <italic>vitellogenin</italic> gene expression in the hepatopancreas and transformation of the gonad from ovotestis to ovary (<xref ref-type="bibr" rid="B24">Levy et al., 2020a</xref>). <xref ref-type="bibr" rid="B47">Zhang et al. (2017)</xref> identified an IAG transcript (<italic>Lw-IAG</italic>) in a PSH shrimp <italic>Lysmata wurdemanni</italic>, and examined its expression profile during gonadal development. The very low but detectable expression of <italic>Lw-IAG</italic> during the euhermaphrodite phases suggested that <italic>Lw-IAG</italic> was possibly responsible for the maintenance of the male reproductive activity in euhermaphrodite phase. Recently, a transcriptomic study identified two insulin-like peptide (ILP) transcripts (<italic>Lvit</italic>-IAG1 and <italic>Lvit</italic>-IAG2) in the peppermint shrimp <italic>Lysmata vittata</italic> (<xref ref-type="bibr" rid="B2">Bao et al., 2020</xref>), and biological functions of <italic>Lvit</italic>-IAG1 have been explored in detail (<xref ref-type="bibr" rid="B29">Liu et al., 2020</xref>). In <italic>L.vittata</italic>, silencing <italic>Lvit-IAG1</italic> impeded development of male-related phenotypes (appendices masculinae (AM) and male gonopores) while suppressing the germ cells at the primary spermatocyte stage, demonstrating that <italic>Lvit</italic>-IAG1 was indeed a functional IAG. Meanwhile, <italic>Lvit-</italic>IAG1 was also suggested to regulate the ovarian development by inhibiting <italic>Lvit-</italic>GIHs and <italic>Lvit-</italic>CFSHs expression in the eyestalk ganglion (<xref ref-type="bibr" rid="B29">Liu et al., 2020</xref>). In decapods, four types of insulin and related peptides were identified: insulin-like androgenic gland hormone (IAG), insulin, relaxins, and gonadulins (<xref ref-type="bibr" rid="B42">Veenstra, 2020</xref>). Besides, there is usually only one IAG gene identified in a species (<xref ref-type="bibr" rid="B27">Li et al., 2012</xref>; <xref ref-type="bibr" rid="B15">Huang et al., 2017b</xref>). Therefore, further in-depth studies are warranted to explore whether <italic>Lvit-</italic>IAG2 is another functional IAG or other insulin and related peptides.</p>
<p>The peppermint shrimp <italic>L. vittata</italic>, a relatively small caridean species (maximum carapace length &#x003C; 1 cm), like other species, displays a unique sexual system, PSH, whereby individuals first mature as males; however, with increasing age and size, they acquire reproductive functions for both males and females (<xref ref-type="bibr" rid="B3">Bauer, 2000</xref>; <xref ref-type="bibr" rid="B1">Alves et al., 2019</xref>). Given its impressive reproductive fecundity, short generation time, and relatively clear ontogenetic gonad development staging system, <italic>L. vittata</italic> is a good model organism for exploring the molecular mechanisms for endocrine regulation of the unique PSH sexual system in crustaceans (<xref ref-type="bibr" rid="B5">Chen et al., 2019</xref>).</p>
<p>In this study, we continuously explored the putative function of <italic>Lvit</italic>-IAG2 in <italic>L. vittata</italic>. Following its similar spatial and temporal expression profiles with <italic>Lvit-IAG1</italic>, we hypothesized that <italic>Lvit</italic>-IAG2 and <italic>Lvit</italic>-IAG1 exert similar biological functions and they coordinatively regulate sexual differentiation of the species. To validate this hypothesis, we performed both short-term and long-term gene knockdown via RNA interference (RNAi). In addition to morphological characteristics, expression levels of gonadal development-associated genes [<italic>GIH</italic>, <italic>vitellogenin</italic> (<italic>Vg</italic>), and <italic>vitellogenin Receptor</italic> (<italic>VgR</italic>)], carbohydrate metabolism (CHH), and sexual differentiation (<italic>CFSH</italic>) were assessed via qRT-PCR.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="S2.SS1">
<title>Animals</title>
<p>The experimental animals (<italic>L. vittata</italic>) were captive-bred at the Fisheries Research Institute of Fujian Province in Xiamen city, China. They were then acclimated in seawater aquaria for 2 days under these conditions: temperature of 25.5 &#x00B1; 0.5&#x00B0;C and salinity of 32 &#x00B1; 1 PSU. During this period, the animals were fed on a commercially formulated shrimp diet daily. A recent study from this laboratory had described two developmental phases covering four gonadal development stages for <italic>L. vittata</italic>, which we followed to define the gonadal stages in the present experiment. Three gonadal development stages (Stage I: ovarian region is smaller than testicular region and both regions are transparent; Stage II: both regions are cloudy white and they are of similar size; and Stage III: ovarian region becomes earth brown and bigger than the testicular region while testicular region is still cloudy white) are included in the male phase, during which the testicular part of the gonad becomes sequentially mature while the ovarian part is still immature. Both testicular region and ovarian region are fully developed and filled with mature germ cells in the euhermaphrodite phase (Stage IV), during which individuals acquire reproductive functions for both males and females (<xref ref-type="bibr" rid="B5">Chen et al., 2019</xref>). Animal handling and experimental procedures were performed with strict adherence to the guidelines approved by the Xiamen University Animal Care and Use Committee.</p>
</sec>
<sec id="S2.SS2">
<title>Fragment Cloning of <italic>Lvit-IAG2</italic></title>
<p>The total RNA was extracted from the androgenic gland at gonadal development stage I using HiPure Universal RNA Kit (Magen) following protocol stipulated by the manufacturer. The first-strand cDNA was generated from 1 &#x03BC;g total RNA using RevertAid First Strand cDNA Synthesis Kit (Fermentas). To verify the accuracy of the predicted open reading frame (ORF) polymerase chain reaction (PCR) was performed by specific primers to ascertain its suitability for the subsequent studies. The PCR reaction was prepared with Ex-Taq polymerase (TaKaRa) and run under the following conditions: 95&#x00B0;C for 3 min; 35 cycles of 95&#x00B0;C for 30 s, 60&#x00B0;C for 30 s, and 72&#x00B0;C for 30 s, followed by 72&#x00B0;C for 5 min final extension. The PCR products were visually examined in 1.0% agarose gel. Subsequently, we purified <italic>Lvit</italic>-IAG2 fragments and inserted them into the pMD19-T vector (TaKaRa) for sequencing. Specific primers used are listed in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap position="float" id="T1">
<label>TABLE 1</label>
<caption><p>Summary of primers used in this study.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left"><bold>Primer</bold></td>
<td valign="top" align="left"><bold>Sequence (5&#x2032;&#x2013;3&#x2032;)</bold></td>
<td valign="top" align="left"><bold>Application</bold></td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">IAG2F</td>
<td valign="top" align="left">TAACCAAGAAATTCACCGTGAAAATGG</td>
<td valign="top" align="left">Fragment validation</td>
</tr>
<tr>
<td valign="top" align="left">IAG2R</td>
<td valign="top" align="left">ACGCTGTAGTCAAGCCATTGGACC</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">CFSH1QF</td>
<td valign="top" align="left">ATCCACACCTCAGAACTCATC</td>
<td valign="top" align="left">RT-PCR/qRT-PCR</td>
</tr>
<tr>
<td valign="top" align="left">CFSH1QR</td>
<td valign="top" align="left">GCACAGGCTACGGTTATCT</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">CFSH2QF</td>
<td valign="top" align="left">CAAGGACGGCGATGATGA</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">CFSH2QR</td>
<td valign="top" align="left">GCGAAGGATCTGAGATGTGTA</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">IAG1QF</td>
<td valign="top" align="left">CTAATCTTGCTGCCTCATTCTAC</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">IAG1QR</td>
<td valign="top" align="left">GCGTCGTTCTCTGTAATAATCG</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">IAG2QF</td>
<td valign="top" align="left">TCAGTCTCAGCCATCTCCT</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">IAG2QR</td>
<td valign="top" align="left">TGAACCGACCACCTCTAATG</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">CHHQF</td>
<td valign="top" align="left">CATCTATGACCGTGAACTCTT</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">CHHQF</td>
<td valign="top" align="left">TACTTGCCGACCATCTGA</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">GIH1QF</td>
<td valign="top" align="left">GACTTCCTGTGGTGCGTGTA</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">GIH1QR</td>
<td valign="top" align="left">GCTCGCAGTATGCTCATGGA</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">GIH2QF</td>
<td valign="top" align="left">ATATGGCGTGTGGTTCTG</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">GIH2QR</td>
<td valign="top" align="left">GAAGTGAGCGGACTACATT</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">VgQF</td>
<td valign="top" align="left">GCAAAAGTGGGAGCCGAAAG</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">VgQR</td>
<td valign="top" align="left">ATCACCCGTAGAGGGTAGGG</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">VgRQF</td>
<td valign="top" align="left">CTGCGTCTCGGAACTCAA</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">VgRQR</td>
<td valign="top" align="left">GTGCTGGTGGTGAAGATGA</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">actinF</td>
<td valign="top" align="left">CGTGACCTGACTGATTACC</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">actinR</td>
<td valign="top" align="left">CGTTACCGATAGTGATTACCT</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">IAG2dsF</td>
<td valign="top" align="left">CTCTGTAAATCAGTCTCAGCCATCT</td>
<td valign="top" align="left">dsRNA synthesis</td>
</tr>
<tr>
<td valign="top" align="left">IAG2dsR</td>
<td valign="top" align="left">ATACCGTCTTGCAGAATTTCACA</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">GFPdsF</td>
<td valign="top" align="left">TGGGCGTGGATAGCGGTTTG</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">GFPdsR</td>
<td valign="top" align="left">GGTCGGGGTAGCGGCTGAAG</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">T7primer</td>
<td valign="top" align="left">TAATACGACTCACTATAGGG</td>
<td/>
</tr>
<tr>
<td valign="top" align="left">SP6primer</td>
<td valign="top" align="left">ATTTAGGTGACACTATAG</td>
<td/>
</tr>
</tbody>
</table></table-wrap>
</sec>
<sec id="S2.SS3">
<title>Bioinformatics Analyses</title>
<p>The primers used for fragment cloning and dsRNA preparation were designed via the Primer 5.0 software. Then we adopted the ORF Finder software<sup><xref ref-type="fn" rid="footnote1">1</xref></sup> to predict the open reading frame (ORF), whereas the SignalP-5.0 Server<sup><xref ref-type="fn" rid="footnote2">2</xref></sup> was used to predict the signal peptides. Further, cysteine residues and putative disulfide bonds were predicted via the DiANNA 1.1 web server<sup><xref ref-type="fn" rid="footnote3">3</xref></sup>. To align the deduced amino acid sequences with reported sequences, we used the Clustal Omega website<sup><xref ref-type="fn" rid="footnote4">4</xref></sup>. The N-glycosylation motif was predicted by the NetNGlyc 1.0 Server<sup><xref ref-type="fn" rid="footnote5">5</xref></sup>.</p>
<p>The Maximum Likelihood method with 1,000 bootstrap replicates based on the JTT matrix-based model in MEGA7 was applied to generate a phylogenetic tree entailing the deduced amino acid sequence alignments and excluding signal peptide of insulin-related peptides. IAG sequences were shown in <xref ref-type="table" rid="T2">Table 2</xref>. Insulin, relaxins, and gonadulins sequences were borrowed from previous works by <xref ref-type="bibr" rid="B42">Veenstra (2020)</xref>.</p>
<table-wrap position="float" id="T2">
<label>TABLE 2</label>
<caption><p>Summary of IAG sequences used in multiple sequence alignment and phylogenetic analysis.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left"><bold>Sequence</bold></td>
<td valign="top" align="left"><bold>Species</bold></td>
<td valign="top" align="left"><bold>GenBank accession number</bold></td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left"><italic>Pch-IAG1</italic></td>
<td valign="top" align="left"><italic>Penaeus chinensis</italic></td>
<td valign="top" align="left">AFU60548.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Pch-IAG2</italic></td>
<td valign="top" align="left"><italic>Penaeus chinensis</italic></td>
<td valign="top" align="left">AFU60549.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Pm-IAG</italic></td>
<td valign="top" align="left"><italic>Penaeus monodon</italic></td>
<td valign="top" align="left">ADA67878.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Pv-IAG</italic></td>
<td valign="top" align="left"><italic>Penaeus vannamei</italic></td>
<td valign="top" align="left">AIR09497.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Pj-IAG</italic></td>
<td valign="top" align="left"><italic>Penaeus japonicus</italic></td>
<td valign="top" align="left">BAK20460.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Je-IAG</italic></td>
<td valign="top" align="left"><italic>Jasus edwardsii</italic></td>
<td valign="top" align="left">AIM55892.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Sv-IAG</italic></td>
<td valign="top" align="left"><italic>Sagmariasus verreauxi</italic></td>
<td valign="top" align="left">AHY99679.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Cd-IAG</italic></td>
<td valign="top" align="left"><italic>Cherax destructor</italic></td>
<td valign="top" align="left">ACD91988.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Cqua-IAG</italic></td>
<td valign="top" align="left"><italic>Cherax quadricarinatus</italic></td>
<td valign="top" align="left">ABH07705.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Pc-IAG</italic></td>
<td valign="top" align="left"><italic>Procambarus clarkii</italic></td>
<td valign="top" align="left">ALX72789.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Pf-IAG</italic></td>
<td valign="top" align="left"><italic>Procambarus fallax</italic></td>
<td valign="top" align="left">ASM94213.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-IAG1</italic></td>
<td valign="top" align="left"><italic>Lysmata vittata</italic></td>
<td valign="top" align="left">MT114196</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-IAG2</italic></td>
<td valign="top" align="left"><italic>Lysmata vittata</italic></td>
<td valign="top" align="left">MT114197</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Es-IAG</italic></td>
<td valign="top" align="left"><italic>Eriocheir sinensis</italic></td>
<td valign="top" align="left">AVK43106.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Cqui-IAG</italic></td>
<td valign="top" align="left"><italic>Chaceon quinquedens</italic></td>
<td valign="top" align="left">ASA45642.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Sp-IAG</italic></td>
<td valign="top" align="left"><italic>Scylla paramamosain</italic></td>
<td valign="top" align="left">AFY09905.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Cs-IAG</italic></td>
<td valign="top" align="left"><italic>Callinectes sapidus</italic></td>
<td valign="top" align="left">AEI72263.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Pp-IAG</italic></td>
<td valign="top" align="left"><italic>Pandalus platyceros</italic></td>
<td valign="top" align="left">ASM94212.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Ppac-IAG</italic></td>
<td valign="top" align="left"><italic>Palaemon pacificus</italic></td>
<td valign="top" align="left">BAJ84109.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Ppau-IAG</italic></td>
<td valign="top" align="left"><italic>Palaemon paucidens</italic></td>
<td valign="top" align="left">BAJ84108.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Ml-IAG</italic></td>
<td valign="top" align="left"><italic>Macrobrachium lar</italic></td>
<td valign="top" align="left">BAJ78349.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Mv-IAG</italic></td>
<td valign="top" align="left"><italic>Macrobrachium vollenhovenii</italic></td>
<td valign="top" align="left">AHZ34725.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lw-IAG</italic></td>
<td valign="top" align="left"><italic>Lysmata wurdemanni</italic></td>
<td valign="top" align="left">10.1371/journal. pone.0172782 (<xref ref-type="bibr" rid="B47">Zhang et al., 2017</xref>)</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Mn-IAG</italic></td>
<td valign="top" align="left"><italic>Macrobrachium nipponense</italic></td>
<td valign="top" align="left">AHA33389.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Mr-IAG</italic></td>
<td valign="top" align="left"><italic>Macrobrachium rosenbergii</italic></td>
<td valign="top" align="left">AWU67706.1</td>
</tr>
</tbody>
</table></table-wrap>
</sec>
<sec id="S2.SS4">
<title>The qRT-PCR Assays</title>
<p>Primers used for quantitative real-time PCR (qRT-PCR) were designed by the Beacon Designer 8 software. RT-PCR products were sequenced as described in section &#x201C;Fragment cloning of <italic>Lvit-IAG2</italic>&#x201D; for accuracy. The melting curves were subjected to intensive analysis to ensure primer specificity. A standard curve was used to calculate the amplification efficiency of each primer pair. Notably, each selected primer pair exhibited a suitable PCR amplification efficiency (96.3&#x2013;104.5%, see <xref ref-type="table" rid="T3">Table 3</xref>). The first-strand cDNA was generated from 300 ng total RNA using TransScript<sup>&#x00AE;</sup> II One-Step gDNA Removal and cDNA short SuperMix Kit (TransGen). The cDNA was diluted to four folds using RNase-free water before it was utilized in qRT-PCR detection. Components including, 10 &#x03BC;l TB Green Premix Ex Taq II (2X) (TaKaRa), 2 &#x03BC;l diluted cDNA, 0.5 &#x03BC;l forward/reverse primer (1 mM), and 7 &#x03BC;l RNase-free water, were used for a 20 &#x03BC;l qRT-PCR reaction system. The reaction was performed using 7500 Real-Time PCR (Applied Biosystems) with the following steps: 95&#x00B0;C for 30 s, followed by 40 cycles of 95&#x00B0;C for 15 s, 58.5&#x00B0;C for 15 s, and 72&#x00B0;C for 30 s. The result was calculated using the 2<sup>&#x2013;&#x0394;&#x0394;<italic>Ct</italic></sup> method, <italic>Lvit-</italic>&#x03B2;<italic>-actin</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT114194">MT114194</ext-link>) as the reference gene.</p>
<table-wrap position="float" id="T3">
<label>TABLE 3</label>
<caption><p>Summary of qPCR efficiency of the primer pairs.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left"><bold>Gene</bold></td>
<td valign="top" align="left"><bold>Primer pairs</bold></td>
<td valign="top" align="center"><bold>Annealing temperature (&#x00B0;C)</bold></td>
<td valign="top" align="center"><bold>Amplification efficiency (%)</bold></td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left"><italic>Lvit-IAG1</italic></td>
<td valign="top" align="left">IAG1QF/IAG1QR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">98.9</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-IAG2</italic></td>
<td valign="top" align="left">IAG2QF/IAG2QR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">104.5</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-CFSH1</italic></td>
<td valign="top" align="left">CFSH1QF/CFSH1QR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">99.1</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-CFSH2</italic></td>
<td valign="top" align="left">CFSH2QF/CFSH2QR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">103.2</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-GIH1</italic></td>
<td valign="top" align="left">GIH1QF/GIH1QR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">98.8</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-GIH2</italic></td>
<td valign="top" align="left">GIH2QF/GIH2QR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">98.5</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-CHH</italic></td>
<td valign="top" align="left">CHHQF/CHHQR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">101.7</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-Vg</italic></td>
<td valign="top" align="left">VgQF/VgQR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">96.3</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-VgR</italic></td>
<td valign="top" align="left">VgRQF/VgRQR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">101.8</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Lvit-</italic>&#x03B2;<italic>-actin</italic></td>
<td valign="top" align="left">actinQF/actinQR</td>
<td valign="top" align="center">58.5</td>
<td valign="top" align="center">101.2</td>
</tr>
</tbody>
</table></table-wrap>
</sec>
<sec id="S2.SS5">
<title>Tissue Expression Profile of <italic>Lvit-IAG2</italic> in <italic>L. vittata</italic></title>
<p>As described in section &#x201C;Fragment cloning of <italic>Lvit-IAG2</italic>,&#x201D; total RNA was extracted from various tissues (eyestalk ganglion, brain, thoracic ganglion, abdominal ganglion, ovary, testis, AG, hepatopancreas, stomach, intestine, heart, gill, and muscle). The first-strand cDNA synthesis was performed as described in section &#x201C;The qRT-PCR assays.&#x201D; RT-PCR tissue expression profile detection was conducted under the following conditions: 95&#x00B0;C for 3 min; 35 cycles of 95&#x00B0;C for 30 s, 58.5&#x00B0;C for 30 s and 72&#x00B0;C for 30 s, followed by 72&#x00B0;C for 5 min final extension. Meanwhile, <italic>Lvit-</italic>&#x03B2;<italic>-actin</italic> was amplified as a positive control, with similar PCR conditions as described above. RT-PCR products were examined using 1.5% agarose gel, then images taken by a UV detector (Geldoc, Thermo Fisher Scientific).</p>
</sec>
<sec id="S2.SS6">
<title>Expression Profiles of <italic>Lvit-IAG2</italic> and <italic>Lvit-IAG1</italic> During Gonadal Development</title>
<p>The AGs of <italic>L. vittata</italic> at different gonadal development stages (I&#x2013;IV) (<italic>n</italic> = 5) were obtained. This was followed by RNA extraction, the first-strand cDNA synthesis, and qRT-PCR analysis as described in sections &#x201C;Fragment cloning of <italic>Lvit-IAG2&#x201D;</italic> and &#x201C;The qRT-PCR assays.&#x201D;</p>
</sec>
<sec id="S2.SS7">
<title>dsRNA Preparation</title>
<p><italic>Lvit-IAG2</italic> fragment encoding C peptide and A chain and green fluorescent protein gene (<italic>GFP</italic>) (exogenous gene control) fragment were cloned into pGEM-T Easy Vector (Promega). dsRNA synthesis was performed using T7 RNA Polymerase (Takara) and SP6 RNA Polymerase (TaKaRa) following the standard protocols. Finally, dsRNA was diluted with 10 mM phosphate-buffered saline (PBS, pH 7.4).</p>
</sec>
<sec id="S2.SS8">
<title>Short-Term Silencing Experiment <italic>in vivo</italic></title>
<p>The efficacy of gene knockdown via RNAi was evaluated by a short-term silencing experiment carried out with <italic>L. vittata</italic> at gonadal development stage I. A total of 15 shrimp (carapace length: 3.21 &#x00B1; 0.35 mm, body weight: 50.40 &#x00B1; 12.14 mg) were randomly and equally assigned to evaluate the efficacy of gene knockdown via RNAi, we prepared a short-term silencing experiment using <italic>L. vittata</italic> at gonadal development stage I. A total of 15 shrimp (carapace length: 3.21 &#x00B1; 0.35 mm, bodyweight: 50.40 &#x00B1; 12.14 mg) were randomly and equally assigned to the following 3 treatment groups (<italic>n</italic> = 5): dsRNA <italic>Lvit-IAG2</italic>-injected, dsRNA <italic>GFP</italic>-injected, and PBS. The delivery of dsRNA (2 &#x03BC;g/g) (<xref ref-type="bibr" rid="B29">Liu et al., 2020</xref>) was by intramuscular injection in the abdominal segment of shrimp, and the PBS-injected treatment received an equivalent volume of PBS. Sampling was performed 24 h after injection. AG and eyestalk ganglion were collected after the shrimp were anesthetized on ice for 5 min. Expression levels of <italic>Lvit-IAG1</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT114196">MT114196</ext-link>) and <italic>Lvit-IAG2</italic> in the AG were detected by qRT-PCR to test the efficacy and specificity of dsRNA-mediated silencing on <italic>Lvit-IAG2</italic>. Meanwhile, <italic>Lvit-GIH1</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT113121">MT113121</ext-link>), <italic>Lvit-GIH2</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT313290">MT313290</ext-link>), <italic>Lvit-CHH</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT701562">MT701562</ext-link>), <italic>Lvit-CFSH1</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT114198">MT114198</ext-link>) and <italic>Lvit-CFSH2</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT114199">MT114199</ext-link>) expression levels in the eyestalk ganglion were also detected. RNA extraction, qRT-PCR, and the first-strand cDNA were performed as described in section &#x201C;Expression Profile of <italic>Lvit-IAG2</italic> and <italic>Lvit-IAG1</italic> During Gonadal Development.&#x201D;</p>
</sec>
<sec id="S2.SS9">
<title>Long-Term Silencing Experiment <italic>in vivo</italic></title>
<p>A long-term silencing experiment was prepared to explore the potential role of <italic>Lvit-IAG2</italic> in sexual differentiation and gonadal development in <italic>L. vittata</italic>. Shrimp (carapace length 3.07 &#x00B1; 0.20 mm, bodyweight 45.98 &#x00B1; 7.43 mg) at stage I were randomly categorized into 3 treatment groups (<italic>n</italic> = 11) as described in section &#x201C;Short-Term Silencing Experiment <italic>in vivo</italic>.&#x201D; An equivalent dose of dsRNA (2 &#x03BC;g/g) or equivoluminal PBS was injected into the abdominal segment of shrimp once every 4 days (for a total of 8 injections in a 29 day duration); during which shrimp were kept in seawater aquaria under the following conditions: Temperature, 25.5 &#x00B1; 0.5&#x00B0;C; salinity, 32 &#x00B1; 1 PSU; a 12L:12D photoperiod. The shrimp were fed with a commercially formulated shrimp diet twice a day. On day 30 (24 h after the 8th injection), all shrimps were sampled after anesthetization. Measurements of carapace length and body weight were recorded. Male and female external sexual characteristics and gonad shape were photographed using a stereomicroscope (model M165FC; Leica Application Suite X). Hematoxylin&#x2013;eosin (H&#x0026;E) staining was applied to visualize morphological and histological changes in ovotestis. The long and short axis lengths of each vitellogenic oocyte were measured and averaged, yielding a mean diameter for each cell. Only cells with visible nucleus were measured. The vitellogenic oocyte diameter of each individual was then calculated from 5 cells/field and 3 fields/section. Samples of AG, eyestalk ganglion, the ovarian region of the gonad, and hepatopancreas were collected to examine the relative mRNA expression levels of <italic>Lvit-IAG2</italic>, <italic>Lvit-GIH1</italic>, <italic>Lvit-GIH2</italic>, <italic>Lvit-CHH</italic>, <italic>Lvit-CFSH1</italic>, <italic>Lvit-CFSH2</italic>, <italic>Lvit-Vg</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT113122">MT113122</ext-link>), and <italic>Lvit-VgR</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT114195">MT114195</ext-link>) by qRT-PCR. RNA extraction, the first-strand cDNA synthesis, and qRT-PCR analyses were performed as described in section &#x201C;Expression Profile of <italic>Lvit-IAG2</italic> and <italic>Lvit-IAG1</italic> During Gonadal Development.&#x201D;</p>
</sec>
<sec id="S2.SS10">
<title>Statistical Analyses</title>
<p>Normality of data was established by the Kolmogorov-Smirnov test. All the data were presented in a normal distribution and tested for variances homogeneity by the Levene&#x2019;s test. All statistical analyses were performed using the SPSS 18.0 software; statistical significance (<italic>p</italic> &#x003C; 0.05) of the data was determined using one-way ANOVA followed by Tukey&#x2019;s multiple range tests. All data were presented as mean &#x00B1; SEM (<italic>n</italic> = 4&#x2013;6).</p>
</sec>
</sec>
<sec id="S3">
<title>Results</title>
<sec id="S3.SS1">
<title>Sequence Analysis of <italic>Lvit-IAG2</italic></title>
<p>A schematic diagram of preproprotein of <italic>Lvit-IAG2</italic> was depicted in <xref ref-type="fig" rid="F1">Figure 1A</xref>. The <italic>Lvit-IAG2</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT114197">MT114197</ext-link>) coding region was 441-bp in length and encoded a 146-aa polypeptide, including a 27-aa signal peptide, 36-aa B chain, 45-aa C peptide, and 38-aa A chain. We predicted 6 conserved cysteine residues in the B chain (C<sub>B</sub><sub>12</sub> and C<sub>B</sub><sub>23</sub>) and A chain (C<sub>A</sub><sub>13</sub>, C<sub>A</sub><sub>14</sub>, C<sub>A</sub><sub>19</sub>, and C<sub>A</sub><sub>27</sub>). The predicted mature <italic>Lvit</italic>-IAG2 peptide comprised the B and A chains with two interchain disulfide bonds (between C<sub>B</sub><sub>12</sub> and C<sub>A</sub><sub>14</sub>, C<sub>B</sub><sub>23</sub>, and C<sub>A</sub><sub>27</sub>, respectively) and an intrachain disulfide bond (between C<sub>A</sub><sub>13</sub> and C<sub>A</sub><sub>19</sub>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Sequence analysis of <italic>Lvit</italic>-IAG2. <bold>(A)</bold> Schematic diagram of preproprotein of <italic>Lvit</italic>-IAG2 containing a signal peptide, B chain, C peptide, and A chain. The six cysteine residues of B chain and A chain were predicted to form three putative disulfide bridges connected with lines. <bold>(B)</bold> Multiple sequence alignment of putative B chain and A chain of decapod IAGs. The sequence information was given in <xref ref-type="table" rid="T2">Table 2</xref>. Six conserved cysteine residues were boxed in red, and three putative disulfide bridges were connected with lines. Predicted N-glycosylation sites were marked by gray shadow. <italic>Lvit</italic>-IAG1 was marked with solid arrows.</p></caption>
<graphic xlink:href="fmars-08-666763-g001.tif"/>
</fig>
</sec>
<sec id="S3.SS2">
<title>Homology and Phylogenetic Analysis</title>
<p>Multiple sequence alignment of the putative B chain and A chain of IAGs from decapod species was highlighted in <xref ref-type="fig" rid="F1">Figure 1B</xref>. Notably, 6 cysteine residues forming two interchain disulfide bonds (between C<sub>B</sub><sub>15</sub> and C<sub>A</sub><sub>18</sub>, C<sub>B</sub><sub>26</sub> and C<sub>A</sub><sub>37</sub>, respectively) and an intrachain disulfide bond (between C<sub>A</sub><sub>18</sub> and C<sub>A</sub><sub>28</sub>) were fully conserved among decapod species. Besides, <italic>Lvit</italic>-IAG2 shared the highest identity with <italic>Lvit</italic>-IAG1 (58.82% for B chain and 44.12% for A chain, respectively). Based on phylogenetic analysis, insulin and related peptides formed four major clades: IAG, insulin, gonadulin, and relaxin. Meanwhile, the IAGs in decapods formed three subclades: (i) The family Caridea; (ii) sequences from the family Brachyura; (iii) the families Astacidea, Achelata, and Penaeoidea (<xref ref-type="fig" rid="F2">Figure 2</xref>). We classified <italic>Lvit</italic>-IAG2 into the subclade containing the family Caridea.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p>Phylogenetic tree of amino acid sequences of the insulin and related peptide in decapods. Except for <italic>Lvit</italic>-IAG2 (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="MT114197">MT114197</ext-link>, <italic>Lysmata vittata</italic>) and <italic>Cs-</italic>IAG<italic>-ov</italic> (GenBank accession number: <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="APC42865.1">APC42865.1</ext-link>, <italic>Callinectes sapidus</italic>), IAG sequences shown in <xref ref-type="table" rid="T2">Table 2</xref>; other sequences were obtained from <xref ref-type="bibr" rid="B42">Veenstra (2020)</xref>. Deduced amino acid sequences of insulin and related peptide excluding signal peptide were used for phylogenic analysis. Phylogenetic analysis was conducted by Maximum Likelihood method based on the JTT matrix-based model in MEGA7. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1,000 replicates) was shown next to the branches. Four major clades were represented by different colors. Each subclade of IAGs was noted by a different color of dashed box. <italic>Lvit-</italic>IAG2 was marked by a solid arrow.</p></caption>
<graphic xlink:href="fmars-08-666763-g002.tif"/>
</fig>
</sec>
<sec id="S3.SS3">
<title>Spatial and Temporal Expression Profiles of <italic>Lvit-IAG2</italic></title>
<p>To explore the spatial distribution profiles of <italic>Lvit-IAG2</italic>, RT-PCR was performed on <italic>L. vittata</italic> at gonadal development stage I. Results demonstrated exclusive expression of <italic>Lvit-IAG2</italic> in the AG (<xref ref-type="fig" rid="F3">Figure 3A</xref>). The relative expression of <italic>Lvit-IAG2</italic> in the AG during gonadal development was also assessed through qRT-PCR (<xref ref-type="fig" rid="F3">Figure 3B</xref>). Notably, the expression levels of <italic>Lvit-IAG2</italic> reached a peak at stage I, decreased sharply at stage II, and were continuously maintained at low levels at stages III and IV [<italic>F</italic><sub>(</sub><sub>3</sub>, <sub>1</sub><sub>6</sub><sub>)</sub> = 19.797, <italic>p</italic> &#x003C; 0.05]. Similar trend was observed in <italic>Lvit-IAG1</italic> expression profiles except that the sharp decrease presented at stage III [<italic>F</italic><sub>(</sub><sub>3</sub>, <sub>1</sub><sub>6</sub><sub>)</sub> = 161.408, <italic>p</italic> &#x003C; 0.05], which was consistent with results of the former study (<xref ref-type="bibr" rid="B29">Liu et al., 2020</xref>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p>Spatial expression profile of <italic>Lvit-IAG2</italic> and temporal expression profiles of <italic>Lvit-IAG2</italic> and <italic>Lvit-IAG1.</italic> <bold>(A)</bold> Distribution of <italic>Lvit-IAG2</italic> in different tissues of <italic>L. vittata</italic>. The analysis was generated by PCR assays with cDNAs from various tissues of individuals at the gonadal development stage I. The <italic>Lvit-</italic>&#x03B2;<italic>-actin</italic> gene was used as a reference control gene. <bold>(B)</bold> Expression profiles of <italic>Lvit-IAG2</italic> and <italic>Lvit-IAG1</italic> in the androgenic gland during gonad development by qRT-PCR. The <italic>Lvit-IAGs</italic> expression levels standardized by <italic>Lvit-</italic>&#x03B2;<italic>-actin</italic> expression levels were represented as mean &#x00B1; SEM (&#x201C;a, b, and c,&#x201D; <italic>p</italic> &#x003C; 0.05; &#x201C;&#x03B1;, &#x03B2;, and &#x03B3;,&#x201D; <italic>p</italic> &#x003C; 0.05; one-way ANOVA followed by Tukey&#x2019;s multiple range tests; <italic>n</italic> = 5).</p></caption>
<graphic xlink:href="fmars-08-666763-g003.tif"/>
</fig>
</sec>
<sec id="S3.SS4">
<title>Short-Term Silencing Experiment <italic>in vivo</italic></title>
<p>Results demonstrated that the transcript level of <italic>Lvit-IAG2</italic> was explicitly inhibited up to 77% [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>12</sub><sub>)</sub> = 22.851, <italic>p</italic> &#x003C; 0.05]. Meanwhile, <italic>Lvit-IAG2</italic> knockdown significantly suppressed the expression of <italic>Lvit-CHH</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>12</sub><sub>)</sub> = 113.857, <italic>p</italic> &#x003C; 0.05]. On the contrary, no significant difference in <italic>Lvit-IAG1</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>12</sub><sub>)</sub> = 0.173, <italic>p</italic> &#x003E; 0.05], <italic>Lvit-CFSH1</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>12</sub><sub>)</sub> = 0.876, <italic>p</italic> &#x003E; 0.05], <italic>Lvit-CFSH2</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>12</sub><sub>)</sub> = 2.270, <italic>p</italic> &#x003E; 0.05], <italic>Lvit-GIH1</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>12</sub><sub>)</sub> = 0.687, <italic>p</italic> &#x003E; 0.05], and <italic>Lvit-GIH2</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>12</sub><sub>)</sub> = 1.378, <italic>p</italic> &#x003E; 0.05] expression was found (<xref ref-type="fig" rid="F4">Figure 4</xref>).</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p>Effects of short-term <italic>Lvit-IAG2</italic> silencing on gene expression of <italic>L. vittata</italic>. The effectiveness of gene knockdown in the short-term <italic>Lvit-IAG2</italic> silencing experiment was evaluated by qRT-PCR. The expression levels of <italic>Lvit-IAG2</italic>, <italic>Lvit-IAG1</italic>, <italic>Lvit-CFSH1</italic>, <italic>Lvit-CFSH2</italic>, <italic>Lvit-GIH1, Lvit-GIH2</italic>, and <italic>Lvit-CHH</italic> were detected following <italic>in vivo</italic> injection with PBS, dsRNA <italic>GFP</italic> or dsRNA <italic>Lvit-IAG2</italic>. The gene expression levels were standardized by <italic>Lvit-</italic>&#x03B2;<italic>-actin</italic> expression levels and represented as mean &#x00B1; SEM (&#x201C;a and b,&#x201D; <italic>p</italic> &#x003C; 0.05; one-way ANOVA followed by Tukey&#x2019;s multiple range tests; <italic>n</italic> = 5).</p></caption>
<graphic xlink:href="fmars-08-666763-g004.tif"/>
</fig>
</sec>
<sec id="S3.SS5">
<title>Long-Term Silencing Experiment <italic>in vivo</italic></title>
<sec id="S3.SS5.SSS1">
<title>Effects of <italic>Lvit-IAG2</italic> Silencing on Gene Expression</title>
<p>At the end of the 30 day long-term trial, the efficacy of gene knockdown was evaluated. Compared to the PBS treatment, <italic>Lvit-IAG2</italic> transcript was 95% inhibited [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 172.195, <italic>p</italic> &#x003C; 0.05]. Besides, <italic>Lvit-IAG2</italic> knockdown not only significantly inhibited the expression of <italic>Lvit-CHH</italic> in the eyestalk ganglion [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 38.197, <italic>p</italic> &#x003C; 0.05] but also repressed <italic>Lvit-Vg</italic> expression in the hepatopancreas [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 6.851, <italic>p</italic> &#x003C; 0.05] and <italic>Lvit-VgR</italic> expression in the ovarian region [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 9.605, <italic>p</italic> &#x003C; 0.05]. However, we found no significant difference in the expression levels of <italic>Lvit-CFSH1</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 0.953, <italic>p</italic> &#x003E; 0.05], <italic>Lvit-CFSH2</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 1.129, <italic>p</italic> &#x003E; 0.05], <italic>Lvit-GIH1</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 1.333, <italic>p</italic> &#x003E; 0.05], and <italic>Lvit-GIH2</italic> [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 1.249, <italic>p</italic> &#x003E; 0.05]; similar results were reported with short-term silencing experiment (<xref ref-type="fig" rid="F4">Figures 4</xref>,<xref ref-type="fig" rid="F5"> 5</xref>).</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>Effects of long-term <italic>Lvit-IAG2</italic> silencing on gene expression of <italic>L. vittata</italic>. The effectiveness of gene knockdown in the long-term <italic>Lvit-IAG2</italic> silencing experiment was evaluated by qRT-PCR. The expression levels of <italic>Lvit-CFSH1</italic>, <italic>Lvit-CFSH2</italic>, <italic>Lvit-GIH1</italic>, <italic>Lvit-GIH2</italic>, <italic>Lvit-CHH</italic>, <italic>Lvit-Vg</italic>, and <italic>Lvit-VgR</italic> were detected following <italic>in vivo</italic> injection with PBS, dsRNA <italic>GFP</italic> or dsRNA <italic>Lvit-IAG2</italic>. The gene expression levels were standardized by <italic>Lvit-</italic>&#x03B2;<italic>-actin</italic> expression levels and represented as mean &#x00B1; SEM (&#x201C;a, b, and c,&#x201D; <italic>p</italic> &#x003C; 0.05; one-way ANOVA followed by Tukey&#x2019;s multiple range tests; <italic>n</italic> = 4&#x2013;6).</p></caption>
<graphic xlink:href="fmars-08-666763-g005.tif"/>
</fig>
</sec>
<sec id="S3.SS5.SSS2">
<title>Effects of <italic>Lvit-IAG2</italic> Silencing on Growth and Development of Sexual Characteristics</title>
<p>After a 30 day long-term experiment, we recorded the average carapace length and bodyweight of the shrimps. Shrimps from the <italic>Lvit-IAG2</italic> silencing treatment group (4.27 &#x00B1; 0.05 mm, 90.05 &#x00B1; 3.08 mg) were much smaller compared to those of the PBS (4.83 &#x00B1; 0.05 mm, 144.78 &#x00B1; 3.04 mg) or dsRNA <italic>GFP</italic> (4.96 &#x00B1; 0.05 mm, 147.87 &#x00B1; 4.33 mg) treatment groups [carapace length: <italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 7.619, <italic>p</italic> &#x003C; 0.05; body weight: <italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 11.082, <italic>p</italic> &#x003C; 0.05] (<xref ref-type="fig" rid="F6">Figures 6A,B</xref>). Moreover, changes in male and female sexual characteristics were documented at the end of the experiment through photography (<xref ref-type="fig" rid="F7">Figure 7</xref>). Notably, <italic>Lvit-IAG2</italic> gene knockdown retarded appendices masculinae (AM) development [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 24.824, <italic>p</italic> &#x003C; 0.05] (<xref ref-type="fig" rid="F6">Figures 6C</xref>, <xref ref-type="fig" rid="F7">7</xref>). No significant difference in female characteristics (female gonophores) and other male sexual characteristics (male gonophores and cincinnuli) was observed (<xref ref-type="fig" rid="F7">Figure 7</xref>).</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption><p>Effects of <italic>Lvit-IAG2</italic> silencing on carapace length, body weight, AM length and oocyte diameter of <italic>L. vittata</italic>. <bold>(A)</bold> Carapace length; <bold>(B)</bold> body weight; <bold>(C)</bold> the normalized length of AM (AM/AI); <bold>(D)</bold> the oocyte diameter. Data were represented as mean &#x00B1; SEM (&#x201C;a and b,&#x201D; <italic>p</italic> &#x003C; 0.05; one-way ANOVA followed by Tukey&#x2019;s multiple range tests; <italic>n</italic> = 4&#x2013;6).</p></caption>
<graphic xlink:href="fmars-08-666763-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption><p>Effects of <italic>Lvit-IAG2</italic> silencing on development of male and female sexual characteristics. Both male (cincinnuli, AM, male gonopore) and female characteristics (female gonopore) were photographed at the end of the 30 day experiment. The AI was marked by white dashed lines, while the AM was marked by yellow dashed lines. Male and female gonopores were noted by blue and red dotted circle, respectively. AM, appendices masculinae; AI, appendix interna.</p></caption>
<graphic xlink:href="fmars-08-666763-g007.tif"/>
</fig>
</sec>
<sec id="S3.SS5.SSS3">
<title>Effects of <italic>Lvit-IAG2</italic> Silencing on Gonadal Development</title>
<p>We examined the morphological and histological characteristics of gonads following <italic>Lvit-IAG2</italic> knockdown; notably, gene silencing resulted in abnormal development of gonad (<xref ref-type="fig" rid="F8">Figure 8</xref>). In the testicular regions, different compositions of germ cell types were reported across the three treatments. For both <italic>GFP</italic>-injected and PBS-injected control, active spermatogenesis occurred (evident from the abundant spermatogonia, spermatid, and spermatozoa found in the testicular region) (<xref ref-type="fig" rid="F8">Figures 8D,E</xref>). However, with <italic>Lvit-IAG2</italic> silencing, a large number of secondary spermatocytes and little spermatid were found, demonstrating arrested spermatogenesis (<xref ref-type="fig" rid="F8">Figure 8F</xref>). Besides, the ovarian region was much smaller (<xref ref-type="fig" rid="F8">Figure 8C</xref>), and in terms of the arrangement, germ cells were tighter and more disordered compared to those of the controls (<xref ref-type="fig" rid="F8">Figure 8I</xref>). Additionally, the oocyte size significantly decreased following <italic>Lvit-IAG2</italic> knockdown [<italic>F</italic><sub>(</sub><sub>2</sub>, <sub>13</sub><sub>)</sub> = 67.576, <italic>p</italic> &#x003C; 0.05] (<xref ref-type="fig" rid="F6">Figures 6D</xref>,<xref ref-type="fig" rid="F8"> 8I</xref>).</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption><p>Effects of <italic>Lvit-IAG2</italic> silencing on gonadal development of <italic>L. vittata.</italic> Gonad morphology and characteristics were photographed at the end of the 30 day experiment following <italic>in vivo</italic> injection of PBS <bold>(A)</bold>, dsRNA <italic>GFP</italic> <bold>(B)</bold>, or dsRNA <italic>Lvit-IAG2</italic> <bold>(C)</bold>. Hematoxylin and eosin (H &#x0026; E)-stained sections was used for structure description. Testicular region mainly consisted of Sg and Sd/Sz in PBS <bold>(D)</bold> and dsRNA <italic>GFP</italic> <bold>(E)</bold> injection treatments while the majority of cells in the testicular region of dsRNA <italic>Lvit-IAG2</italic> injected treatment were Sc II <bold>(F)</bold>. Ovarian regions of the shrimp from the <italic>Lvit-IAG2</italic> silencing treatment were more transparent and smaller in relative size <bold>(C)</bold>. Compared with the PBS-injected <bold>(G)</bold> and dsRNA <italic>GFP</italic>-injected <bold>(H)</bold> group, oocytes were less developed and arranged more tightly and disorderly in the dsRNA <italic>Lvit-IAG2</italic>-injected group <bold>(I)</bold>. Ova, ovary; Tes, testis; Ovd, oviduct; Spd, sperm duct; Oog, oogonia; Ooc, oocytes; Fc, follicular cell; Sg, spermatogonia; Sc I, primary spermatocyte; Sc II, secondary spermatocyte; Sd, spermatid; Sz, spermatozoa.</p></caption>
<graphic xlink:href="fmars-08-666763-g008.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec id="S4">
<title>Discussion</title>
<p>Following a previous transcriptome study, two ILP transcripts, <italic>Lvit-IAG1</italic> and <italic>Lvit-IAG2</italic>, exist in the PSH shrimp, <italic>L. vittata</italic> (<xref ref-type="bibr" rid="B2">Bao et al., 2020</xref>). <italic>Lvit-</italic>IAG1 has also been identified as an IAG via functional study (<xref ref-type="bibr" rid="B29">Liu et al., 2020</xref>). In this study, we expounded on the biological functions of another ILP, <italic>Lvit-IAG2</italic>. <italic>Lvit-</italic>IAG2 contained features that were structurally conserved in IAGs: a preproprotein comprised of a signal peptide, B chain, C peptide, and A chain; six strictly conserved cysteine residues forming three disulfide bonds in the B chain and A chain among predicted mature peptides. Besides, based on phylogenetic tree analysis, <italic>Lvit</italic>-IAG2 was classified into the IAG clade and subclade of the family Caridea. These results confirmed that <italic>Lvit</italic>-IAG2 was identified structurally as an IAG. Meanwhile, we noticed that some studies reported more than one transcript of IAG in a decapod species (<xref ref-type="bibr" rid="B27">Li et al., 2012</xref>; <xref ref-type="bibr" rid="B15">Huang et al., 2017b</xref>). However, these transcripts were generated by alternative splicing of one IAG gene, and deduced mature amino acid sequences of these encoded IAG proteins were nearly identical (<xref ref-type="bibr" rid="B27">Li et al., 2012</xref>; <xref ref-type="bibr" rid="B15">Huang et al., 2017b</xref>). Herein, mature peptides of <italic>Lvit</italic>-IAG2 and <italic>Lvit</italic>-IAG1 shared only 58.82% similarity for B chain and 44.12% similarity for A chain, respectively. These results might suggest that two different IAG genes exist in the species, which has never been reported.</p>
<p>Earlier experiments demonstrated that <italic>IAG</italic> genes were mainly expressed in the AG in decapod crustaceans (<xref ref-type="bibr" rid="B30">Manor et al., 2007</xref>). More recently, studies detected <italic>IAG</italic> mRNA expression in the hepatopancreas and ovary (<xref ref-type="bibr" rid="B6">Chung, 2014</xref>; <xref ref-type="bibr" rid="B16">Huang et al., 2014</xref>). Herein, the <italic>Lvit-IAG2</italic> gene was exclusively expressed in the AG; its expression level was significantly higher in the male phase (Stage I-III) than that in the euhermaphrodite phase (Stage IV) during the life cycle of the PSH species. These reports implied that <italic>Lvit-IAG2</italic> potentially can initiate male sexual differentiation in the male phase and maintain male reproductive activity in the euhermaphrodite phase. Notably, <italic>Lvit-IAG2</italic> and <italic>Lvit-IAG1</italic> displayed slightly different expression trends during gonadal development. The expression levels of <italic>Lvit-IAG2</italic> decreased sharply at stage II, while significant downregulation of <italic>Lvit-IAG1</italic> expression appeared at stage III. The distinction suggested different functions of the two IAGs.</p>
<p>The present results demonstrated that gene silencing via RNAi successfully induced a specific knockdown of <italic>Lvit-IAG2</italic> transcripts levels by 77 and 95% in short-term and long-term silencing experiments, respectively. Furthermore, we assessed the influence of <italic>Lvit</italic>-IAG2 on the development of male features by comparing male external features and germ cell type composition in the testicular region. Notably, previous studies revealed that <italic>IAG</italic> silencing delayed the appearance of male sexual characteristics and induced testicular spermatogenesis arrest in the freshwater prawn, <italic>M. rosenbergii</italic> (<xref ref-type="bibr" rid="B43">Ventura et al., 2009</xref>) and the red swamp crayfish <italic>Procambarus clarkii</italic> (<xref ref-type="bibr" rid="B40">Shi et al., 2019</xref>). Similar results were observed for <italic>Lvit-IAG1</italic> in <italic>L. vittata</italic>. <italic>Lvit-IAG1</italic> knockdown induced abnormal development of both AM and male gonopores, and suppressed germ cells at the primary spermatocytes (<xref ref-type="bibr" rid="B29">Liu et al., 2020</xref>). However, we found that <italic>Lvit-</italic>IAG2 and <italic>Lvit-</italic>IAG1 exerted distinct functions in the male sexual differentiation process. <italic>Lvit-IAG2</italic> knockdown solely suppressed AM development and germ cells arrested at the secondary spermatocytes, implying its crucial contribution to secondary spermatocyte-to-spermatid transition. These results demonstrated that <italic>Lvit</italic>-IAG2 was indeed another functional IAG of <italic>L.vittata</italic>. On the other hand, either primary-to-secondary spermatocyte transition or secondary spermatocyte-to-spermatid transition was essential for the normal development of testicular region. Meanwhile, both of the two <italic>Lvit</italic>-IAGs might be involved in developmental processes specific to AM. Thus, <italic>Lvit</italic>-IAG2 and <italic>Lvit</italic>-IAG1 acted in concert, and both of them were indispensable in the male sexual differentiation process of the PSH species.</p>
<p>Simultaneously, we evaluated the effects of <italic>Lvit</italic>-IAG2 on female development; notably, no significant difference in <italic>Lvit-CFSHs</italic> gene expression among the three treatments in both short-term and long-term silencing experiments was reported. Further morphological characterization revealed similar results as <italic>Lvit-IAG2</italic> knockdown did not impact the development of female gonopores. Nonetheless, a significant inhibitory effect on the development of ovarian region was observed. Contrary to the controls, shrimp in the <italic>Lvit-IAG2</italic> silencing treatments exhibited decreased volume of the ovarian region, and germ cells were much smaller and arranged in a disorganized mass. It is interesting to note that silencing IAG gene in another hermaphrodite shrimp <italic>P. platyceros</italic> induced feminization of individuals, including <italic>Vg</italic> expression in the hepatopancreas and ovarian development (<xref ref-type="bibr" rid="B24">Levy et al., 2020a</xref>), which is consistent with studies in dioecious species. Different outcome might be induced by the difference between SPH and PSH. For SPH shirimp <italic>P. platyceros</italic>, individuals first mature as male and eventually acquire reproductive function of female, indicating that functional males and females have different reproductive systems (<xref ref-type="bibr" rid="B25">Levy et al., 2020b</xref>). This phase-specific gender pattern is similar, to some extent, with dioecious species, and so it is with reproductive endocrine system. However, For SPH shrimp <italic>L. vitata</italic>, individuals finally acquire the same euhermaphrodite reproductive system with both male and female functions. It is therefore reasonable to speculate that different reproductive endocrine system presents in PSH species, and IAG is likely to possess different functions in ovarian development. Following reports from previous experiments, <italic>Lvit</italic>-IAG1 could regulate ovarian development by inhibiting <italic>Lvit</italic>-GIHs expression (<xref ref-type="bibr" rid="B29">Liu et al., 2020</xref>). However, in the present experiment, we reported no significant change in <italic>Lvit-GIHs</italic> transcript levels. Remarkably, knockdown of <italic>Lvit-IAG1</italic> and <italic>Lvit-IAG2</italic> suppressed ovarian development. Thus, we could not eliminate the possibility that insulin and related peptides could directly regulate <italic>Vg</italic> expression in the hepatopancreas, and promote oocyte growth and development. Following a review by <xref ref-type="bibr" rid="B9">Das and Arur (2017)</xref>, insulin signaling played a pivotal role in the regulation of oocyte growth, development, and maturation in both vertebrates and invertebrates. For instance, in the rainbow trout <italic>Oncorhynchus mykiss</italic>, insulin was reported to induce the accumulation of vitellogenin in cultured oocytes (<xref ref-type="bibr" rid="B41">Shibata et al., 1993</xref>). Besides, ILPs could also control germline cyst growth and vitellogenesis in <italic>Drosophila</italic> (<xref ref-type="bibr" rid="B21">LaFever and Drummond-Barbosa, 2005</xref>). In the silkmoth <italic>Bombyx mori</italic>, knock-out of an insulin-like growth factor (IGF)-like peptide gene induced drastically lower ovary weights and number of eggs in female individuals (<xref ref-type="bibr" rid="B12">Fujinaga et al., 2019</xref>). In the kuruma shrimp <italic>Marsupenaeus japonicas</italic>, <italic>Mj</italic>-IAG was proved to have insulin-type disulfide, rather than AGH-type (<xref ref-type="bibr" rid="B19">Katayama et al., 2014</xref>), indicating that IAG potentially possessed similar functions as IGF or insulin. Moreover, a recent <italic>in vitro</italic> experiment revealed that bovine insulin induced the <italic>Vg</italic> expression in the hepatopancreas of the mud crab <italic>S. paramamosain</italic>, implying that insulin-like peptide could directly modulate vitellogenesis in the hepatopancreas (<xref ref-type="bibr" rid="B14">Huang et al., 2017a</xref>). Thus, <italic>Lvit</italic>-IAG2 potentially possessed similar functions as IGF or insulin in the PSH <italic>L. vittata</italic>, which might be explained by their structural similarities (<xref ref-type="bibr" rid="B31">Martin et al., 1998</xref>, <xref ref-type="bibr" rid="B32">1999</xref>; <xref ref-type="bibr" rid="B30">Manor et al., 2007</xref>). On the other hand, growth and reproduction are closely related (<xref ref-type="bibr" rid="B33">Michalakis et al., 2013</xref>). Herein, we noted that <italic>Lvit-IAG2</italic> silencing induced significantly slower growth in the treated shrimp, implying that <italic>Lvit</italic>-IAG2 is potentially a growth regulator. Similar results have been reported in the freshwater prawn <italic>M. rosenbergii</italic> (<xref ref-type="bibr" rid="B43">Ventura et al., 2009</xref>) and the mud crab <italic>S. paramamosain</italic> (<xref ref-type="bibr" rid="B16">Huang et al., 2014</xref>), and with the function of insulin and related peptide in mammals (<xref ref-type="bibr" rid="B4">Bruning et al., 2000</xref>; <xref ref-type="bibr" rid="B36">Neirijnck et al., 2019</xref>), fishes (<xref ref-type="bibr" rid="B38">Reinecke, 2010</xref>), and insects (<xref ref-type="bibr" rid="B46">Wu and Brown, 2006</xref>). Moreover, transcript levels of <italic>Lvit-CHH</italic> significantly decreased following knock down of <italic>Lvit-IAG2</italic> in both short-term and long-term experiment. CHH is defined as a metabolism associated hormone as its participation in multiple metabolic processes, including glycometabolism (<xref ref-type="bibr" rid="B45">Webster et al., 2012</xref>), lipid metabolism (<xref ref-type="bibr" rid="B34">Montiel-Arzate et al., 2020</xref>) and ammonia metabolism (<xref ref-type="bibr" rid="B48">Zhang et al., 2020</xref>). Of note, the low expression of <italic>Lvit-CHH</italic> implied the growth and developmental abnormalities of individuals. Hence, it could also be reasonably hypothesized that retardation in the development of the ovarian region in <italic>Lvit-IAG2</italic> silencing shrimp potentially arose from the abnormity in the growth and development of the shrimp. Further <italic>in vitro</italic> experiments are required to test the two hypotheses.</p>
<p>Taken together, the present work characterized an <italic>IAG</italic> gene, <italic>Lvit-IAG2</italic>, from the AG of the peppermint shrimp <italic>L. vittata</italic> and explored its crucial functions in regulating primary and secondary sexual characteristics of the PSH species. Moreover, this study implicated <italic>Lvit</italic>-IAG2 and <italic>Lvit</italic>-IAG1 to act in concert, aimed at regulating male sexual differentiation in the PSH shrimp. Also, <italic>Lvit</italic>-IAG2 may directly contribute to ovarian development, or via the promotion of growth in <italic>L. vittata</italic>. These findings provide new insights into molecular mechanisms of sexual differentiation in PSH crustaceans, expounding our current knowledge on crustacean reproductive endocrinology.</p>
</sec>
<sec id="S5">
<title>Data Availability Statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.</p>
</sec>
<sec id="S6">
<title>Author Contributions</title>
<p>FL contributed to conceptualization, methodology, software, validation, formal analysis, investigation, data curation, visualization and writing&#x2014;original draft preparation of the study. HY contributed to conceptualization, methodology, validation, data curation, writing&#x2014;review and editing, supervision, project administration and funding acquisition. WS contributed to investigation. AL contributed to validation and writing&#x2014;review and editing of the study. ZZ contributed to funding acquisition and provided resources. All authors contributed to manuscript revision, read, and approved the submitted version.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</body>
<back>
<fn-group>
<fn fn-type="financial-disclosure">
<p><bold>Funding.</bold> The work was supported by the special fund of marine and fishery structure adjustment in Fujian (2020HYJG01 and 2020HYJG08).</p>
</fn>
</fn-group>
<ack>
<p>We also thank all laboratory members for their constructive suggestions and discussions.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Alves</surname> <given-names>D. F. R.</given-names></name> <name><surname>L&#x00F3;pez Greco</surname> <given-names>L. S.</given-names></name> <name><surname>Barros-Alves</surname> <given-names>S.</given-names></name> <name><surname>de</surname> <given-names>P.</given-names></name> <name><surname>Hirose</surname> <given-names>G. L.</given-names></name></person-group> (<year>2019</year>). <article-title>Sexual system, reproductive cycle and embryonic development of the red-striped shrimp <italic>Lysmata vittata</italic>, an invader in the western Atlantic Ocean.</article-title> <source><italic>PLoS One</italic></source> <volume>14</volume>:<issue>e0210723</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0210723</pub-id> <pub-id pub-id-type="pmid">30645636</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bao</surname> <given-names>C.</given-names></name> <name><surname>Liu</surname> <given-names>F.</given-names></name> <name><surname>Yang</surname> <given-names>Y.</given-names></name> <name><surname>Lin</surname> <given-names>Q.</given-names></name> <name><surname>Ye</surname> <given-names>H.</given-names></name></person-group> (<year>2020</year>). <article-title>Identification of peptides and pheir GPCRs in the peppermint shrimp <italic>Lysmata vittata</italic>, a protandric simultaneous hermaphrodite species.</article-title> <source><italic>Front. Endocrinol.</italic></source> <volume>11</volume>:<issue>229</issue>. <pub-id pub-id-type="doi">10.3389/fendo.2020.00226</pub-id> <pub-id pub-id-type="pmid">32425883</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bauer</surname> <given-names>R. T.</given-names></name></person-group> (<year>2000</year>). <article-title>Simultaneous hermaphroditism in Caridean shrimps: a unique and puzzling sexual system in the Decapoda.</article-title> <source><italic>J. Crustac. Biol.</italic></source> <volume>20</volume> <fpage>116</fpage>&#x2013;<lpage>128</lpage>. <pub-id pub-id-type="doi">10.1163/1937240X-90000014</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bruning</surname> <given-names>J. C.</given-names></name> <name><surname>Gautam</surname> <given-names>D.</given-names></name> <name><surname>Burks</surname> <given-names>D. J.</given-names></name> <name><surname>Gillette</surname> <given-names>J.</given-names></name> <name><surname>Schubert</surname> <given-names>M.</given-names></name> <name><surname>Orban</surname> <given-names>P. C.</given-names></name><etal/></person-group> (<year>2000</year>). <article-title>Role of brain insulin receptor in control of body weight and reproduction.</article-title> <source><italic>Science</italic></source> <volume>289</volume> <fpage>2122</fpage>&#x2013;<lpage>2125</lpage>. <pub-id pub-id-type="doi">10.1126/science.289.5487.2122</pub-id> <pub-id pub-id-type="pmid">11000114</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>D.</given-names></name> <name><surname>Liu</surname> <given-names>F.</given-names></name> <name><surname>Zhu</surname> <given-names>Z.</given-names></name> <name><surname>Lin</surname> <given-names>Q.</given-names></name> <name><surname>Zeng</surname> <given-names>C.</given-names></name> <name><surname>Ye</surname> <given-names>H.</given-names></name></person-group> (<year>2019</year>). <article-title>Ontogenetic development of gonads and external sexual characters of the protandric simultaneous hermaphrodite peppermint shrimp, <italic>Lysmata vittata</italic> (Caridea: Hippolytidae).</article-title> <source><italic>PLoS One</italic></source> <volume>14</volume>:<issue>e0215406</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0215406</pub-id> <pub-id pub-id-type="pmid">31002693</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chung</surname> <given-names>J. S.</given-names></name></person-group> (<year>2014</year>). <article-title>An insulin-like growth factor found in hepatopancreas implicates carbohydrate metabolism of the blue crab <italic>Callinectes sapidus</italic>.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>199</volume> <fpage>56</fpage>&#x2013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2014.01.012</pub-id> <pub-id pub-id-type="pmid">24503150</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Comeau</surname> <given-names>M.</given-names></name> <name><surname>Benhalima</surname> <given-names>K.</given-names></name></person-group> (<year>2018a</year>). <article-title>Functional anatomy of the female reproductive system of the American lobster (<italic>Homarus americanus</italic>).</article-title> <source><italic>J. Morphol.</italic></source> <volume>279</volume> <fpage>1603</fpage>&#x2013;<lpage>1614</lpage>. <pub-id pub-id-type="doi">10.1002/jmor.20889</pub-id> <pub-id pub-id-type="pmid">30397936</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Comeau</surname> <given-names>M.</given-names></name> <name><surname>Benhalima</surname> <given-names>K.</given-names></name></person-group> (<year>2018b</year>). <article-title>Functional anatomy of the male reproductive system of the American lobster (<italic>Homarus americanus</italic>).</article-title> <source><italic>J. Morphol.</italic></source> <volume>279</volume> <fpage>1431</fpage>&#x2013;<lpage>1443</lpage>. <pub-id pub-id-type="doi">10.1002/jmor.20878</pub-id> <pub-id pub-id-type="pmid">30192993</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Das</surname> <given-names>D.</given-names></name> <name><surname>Arur</surname> <given-names>S.</given-names></name></person-group> (<year>2017</year>). <article-title>Conserved insulin signaling in the regulation of oocyte growth, development, and maturation.</article-title> <source><italic>Mol. Reprod. Dev.</italic></source> <volume>84</volume> <fpage>444</fpage>&#x2013;<lpage>459</lpage>. <pub-id pub-id-type="doi">10.1002/mrd.22806</pub-id> <pub-id pub-id-type="pmid">28379636</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Das</surname> <given-names>R.</given-names></name> <name><surname>Krishna</surname> <given-names>G.</given-names></name> <name><surname>Priyadarshi</surname> <given-names>H. P.</given-names></name> <name><surname>Babu</surname> <given-names>G.</given-names></name> <name><surname>Pavan-Kumar</surname> <given-names>A.</given-names></name> <name><surname>Rajendran</surname> <given-names>K. V.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Captive maturation studies in <italic>Penaeus monodon</italic> by GIH silencing using constitutively expressed long hairpin RNA.</article-title> <source><italic>Aquaculture</italic></source> <volume>448</volume> <fpage>512</fpage>&#x2013;<lpage>520</lpage>. <pub-id pub-id-type="doi">10.1016/j.aquaculture.2015.06.036</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fu</surname> <given-names>C.</given-names></name> <name><surname>Li</surname> <given-names>F.</given-names></name> <name><surname>Wang</surname> <given-names>L.</given-names></name> <name><surname>Wu</surname> <given-names>F.</given-names></name> <name><surname>Wang</surname> <given-names>J.</given-names></name> <name><surname>Fan</surname> <given-names>X.</given-names></name><etal/></person-group> (<year>2020</year>). <article-title>Molecular characteristics and abundance of insulin-like androgenic gland hormone and effects of RNA interference in <italic>Eriocheir sinensis</italic>.</article-title> <source><italic>Anim. Reprod. Sci.</italic></source> <volume>215</volume>:<issue>106332</issue>. <pub-id pub-id-type="doi">10.1016/j.anireprosci.2020.106332</pub-id> <pub-id pub-id-type="pmid">32216938</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fujinaga</surname> <given-names>D.</given-names></name> <name><surname>Shiomi</surname> <given-names>K.</given-names></name> <name><surname>Yagi</surname> <given-names>Y.</given-names></name> <name><surname>Kataoka</surname> <given-names>H.</given-names></name> <name><surname>Mizoguchi</surname> <given-names>A.</given-names></name></person-group> (<year>2019</year>). <article-title>An insulin-like growth factor-like peptide promotes ovarian development in the silkmoth <italic>Bombyx mori</italic>.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>9</volume>:<issue>18446</issue>. <pub-id pub-id-type="doi">10.1038/s41598-019-54962-w</pub-id> <pub-id pub-id-type="pmid">31804598</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guo</surname> <given-names>Q.</given-names></name> <name><surname>Li</surname> <given-names>S.</given-names></name> <name><surname>Lv</surname> <given-names>X.</given-names></name> <name><surname>Xiang</surname> <given-names>J.</given-names></name> <name><surname>Manor</surname> <given-names>R.</given-names></name> <name><surname>Sagi</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2019</year>). <article-title>Sex-biased CHHs and their putative receptor regulate the expression of IAG gene in the Shrimp <italic>Litopenaeus vannamei</italic>.</article-title> <source><italic>Front. Physiol.</italic></source> <volume>10</volume>:<issue>1525</issue>. <pub-id pub-id-type="doi">10.3389/fphys.2019.01525</pub-id> <pub-id pub-id-type="pmid">31920723</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>X.</given-names></name> <name><surname>Feng</surname> <given-names>B.</given-names></name> <name><surname>Huang</surname> <given-names>H.</given-names></name> <name><surname>Ye</surname> <given-names>H.</given-names></name></person-group> (<year>2017a</year>). <article-title>In vitro stimulation of vitellogenin expression by insulin in the mud crab, <italic>Scylla paramamosain</italic>, mediated through PI3K/Akt/TOR pathway.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>250</volume> <fpage>175</fpage>&#x2013;<lpage>180</lpage>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2017.06.013</pub-id> <pub-id pub-id-type="pmid">28647319</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>X.</given-names></name> <name><surname>Ye</surname> <given-names>H.</given-names></name> <name><surname>Chung</surname> <given-names>J. S.</given-names></name></person-group> (<year>2017b</year>). <article-title>The presence of an insulin-like androgenic gland factor (IAG) and insulin-like peptide binding protein (ILPBP) in the ovary of the blue crab, <italic>Callinectes sapidus</italic> and their roles in ovarian development.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>249</volume> <fpage>64</fpage>&#x2013;<lpage>70</lpage>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2017.05.001</pub-id> <pub-id pub-id-type="pmid">28479084</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>X.</given-names></name> <name><surname>Ye</surname> <given-names>H.</given-names></name> <name><surname>Huang</surname> <given-names>H.</given-names></name> <name><surname>Yang</surname> <given-names>Y.</given-names></name> <name><surname>Gong</surname> <given-names>J.</given-names></name></person-group> (<year>2014</year>). <article-title>An insulin-like androgenic gland hormone gene in the mud crab, <italic>Scylla paramamosain</italic>, extensively expressed and involved in the processes of growth and female reproduction.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>204</volume> <fpage>229</fpage>&#x2013;<lpage>238</lpage>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2014.06.002</pub-id> <pub-id pub-id-type="pmid">24929228</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jiang</surname> <given-names>Q.</given-names></name> <name><surname>Lu</surname> <given-names>B.</given-names></name> <name><surname>Lin</surname> <given-names>D.</given-names></name> <name><surname>Huang</surname> <given-names>H.</given-names></name> <name><surname>Chen</surname> <given-names>X.</given-names></name> <name><surname>Ye</surname> <given-names>H.</given-names></name></person-group> (<year>2020</year>). <article-title>Role of crustacean female sex hormone (CFSH) in sex differentiation in early juvenile mud crabs, <italic>Scylla paramamosain</italic>.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>289</volume>:<issue>113383</issue>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2019.113383</pub-id> <pub-id pub-id-type="pmid">31904358</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Juchault</surname> <given-names>P.</given-names></name></person-group> (<year>1999</year>). <article-title>Hermaphroditism and gonochorism. a new hypothesis on the evolution of sexuality in Crustacea.</article-title> <source><italic>C. R. Acad. Sci. III</italic></source> <volume>322</volume> <fpage>423</fpage>&#x2013;<lpage>427</lpage>. <pub-id pub-id-type="doi">10.1016/S0764-4469(99)80078-X</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Katayama</surname> <given-names>H.</given-names></name> <name><surname>Kubota</surname> <given-names>N.</given-names></name> <name><surname>Hojo</surname> <given-names>H.</given-names></name> <name><surname>Okada</surname> <given-names>A.</given-names></name> <name><surname>Kotaka</surname> <given-names>S.</given-names></name> <name><surname>Tsutsui</surname> <given-names>N.</given-names></name><etal/></person-group> (<year>2014</year>). <article-title>Direct evidence for the function of crustacean insulin-like androgenic gland factor (IAG): total chemical synthesis of IAG.</article-title> <source><italic>Bioorg. Med. Chem.</italic></source> <volume>22</volume> <fpage>5783</fpage>&#x2013;<lpage>5789</lpage>. <pub-id pub-id-type="doi">10.1016/j.bmc.2014.09.031</pub-id> <pub-id pub-id-type="pmid">25270404</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Khalaila</surname> <given-names>I.</given-names></name> <name><surname>Manor</surname> <given-names>R.</given-names></name> <name><surname>Weil</surname> <given-names>S.</given-names></name> <name><surname>Granot</surname> <given-names>Y.</given-names></name> <name><surname>Keller</surname> <given-names>R.</given-names></name> <name><surname>Sagi</surname> <given-names>A.</given-names></name></person-group> (<year>2002</year>). <article-title>The eyestalk-androgenic gland-testis endocrine axis in the crayfish <italic>Cherax quadricarinatus</italic>.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>127</volume> <fpage>147</fpage>&#x2013;<lpage>156</lpage>. <pub-id pub-id-type="doi">10.1016/S0016-6480(02)00031-X</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>LaFever</surname> <given-names>L.</given-names></name> <name><surname>Drummond-Barbosa</surname> <given-names>D.</given-names></name></person-group> (<year>2005</year>). <article-title>Direct control of germline stem cell division and cyst growth by neural insulin in <italic>Drosophila</italic>.</article-title> <source><italic>Science</italic></source> <volume>309</volume> <fpage>1071</fpage>&#x2013;<lpage>1073</lpage>. <pub-id pub-id-type="doi">10.1126/science.1111410</pub-id> <pub-id pub-id-type="pmid">16099985</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Levy</surname> <given-names>T.</given-names></name> <name><surname>Sagi</surname> <given-names>A.</given-names></name></person-group> (<year>2020</year>). <article-title>The &#x201C;IAG-Switch&#x201D;-a key controlling element in decapod crustacean sex differentiation.</article-title> <source><italic>Front. Endocrinol.</italic></source> <volume>11</volume>:<issue>651</issue>. <pub-id pub-id-type="doi">10.3389/fendo.2020.00651</pub-id> <pub-id pub-id-type="pmid">33013714</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Levy</surname> <given-names>T.</given-names></name> <name><surname>Rosen</surname> <given-names>O.</given-names></name> <name><surname>Simons</surname> <given-names>O.</given-names></name> <name><surname>Savaya Alkalay</surname> <given-names>A.</given-names></name> <name><surname>Sagi</surname> <given-names>A.</given-names></name></person-group> (<year>2017</year>). <article-title>The gene encoding the insulin-like androgenic gland hormone in an all-female parthenogenetic crayfish.</article-title> <source><italic>PLoS One</italic></source> <volume>12</volume>:<issue>e0189982</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0189982</pub-id> <pub-id pub-id-type="pmid">29261765</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Levy</surname> <given-names>T.</given-names></name> <name><surname>Tamone</surname> <given-names>S. L.</given-names></name> <name><surname>Manor</surname> <given-names>R.</given-names></name> <name><surname>Aflalo</surname> <given-names>E. D.</given-names></name> <name><surname>Sklarz</surname> <given-names>M. Y.</given-names></name> <name><surname>Chalifa-Caspi</surname> <given-names>V.</given-names></name><etal/></person-group> (<year>2020a</year>). <article-title>The IAG-switch and further transcriptomic insights into sexual differentiation of a protandric shrimp.</article-title> <source><italic>Front. Mar. Sci.</italic></source> <volume>7</volume>:<issue>587454</issue>. <pub-id pub-id-type="doi">10.3389/fmars.2020.587454</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Levy</surname> <given-names>T.</given-names></name> <name><surname>Tamone</surname> <given-names>S. L.</given-names></name> <name><surname>Manor</surname> <given-names>R.</given-names></name> <name><surname>Bower</surname> <given-names>E. D.</given-names></name> <name><surname>Sagi</surname> <given-names>A.</given-names></name></person-group> (<year>2020b</year>). <article-title>The protandric life history of the Northern spot shrimp <italic>Pandalus platyceros</italic>: molecular insights and implications for fishery management.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>10</volume>:<issue>1287</issue>. <pub-id pub-id-type="doi">10.1038/s41598-020-58262-6</pub-id> <pub-id pub-id-type="pmid">31992795</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>F.</given-names></name> <name><surname>Bai</surname> <given-names>H.</given-names></name> <name><surname>Zhang</surname> <given-names>W.</given-names></name> <name><surname>Fu</surname> <given-names>H.</given-names></name> <name><surname>Jiang</surname> <given-names>F.</given-names></name> <name><surname>Liang</surname> <given-names>G.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Cloning of genomic sequences of three crustacean hyperglycemic hormone superfamily genes and elucidation of their roles of regulating insulin-like androgenic gland hormone gene.</article-title> <source><italic>Gene</italic></source> <volume>561</volume> <fpage>68</fpage>&#x2013;<lpage>75</lpage>. <pub-id pub-id-type="doi">10.1016/j.gene.2015.02.012</pub-id> <pub-id pub-id-type="pmid">25680292</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>S.</given-names></name> <name><surname>Li</surname> <given-names>F.</given-names></name> <name><surname>Sun</surname> <given-names>Z.</given-names></name> <name><surname>Xiang</surname> <given-names>J.</given-names></name></person-group> (<year>2012</year>). <article-title>Two spliced variants of insulin-like androgenic gland hormone gene in the Chinese shrimp, <italic>Fenneropenaeus chinensis</italic>.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>117</volume> <fpage>246</fpage>&#x2013;<lpage>255</lpage>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2012.04.010</pub-id> <pub-id pub-id-type="pmid">22561290</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>A.</given-names></name> <name><surname>Liu</surname> <given-names>J.</given-names></name> <name><surname>Liu</surname> <given-names>F.</given-names></name> <name><surname>Huang</surname> <given-names>Y.</given-names></name> <name><surname>Wang</surname> <given-names>G.</given-names></name> <name><surname>Ye</surname> <given-names>H.</given-names></name></person-group> (<year>2017</year>). <article-title>Crustacean female sex hormone from the mud crab <italic>Scylla paramamosain</italic> is highly expressed in prepubertal males and inhibits the development of androgenic gland.</article-title> <source><italic>Front. Physiol.</italic></source> <volume>9</volume>:<issue>924</issue>. <pub-id pub-id-type="doi">10.3389/fphys.2018.00924</pub-id> <pub-id pub-id-type="pmid">30065661</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>F.</given-names></name> <name><surname>Shi</surname> <given-names>W.</given-names></name> <name><surname>Ye</surname> <given-names>H.</given-names></name> <name><surname>Zeng</surname> <given-names>C.</given-names></name> <name><surname>Zhu</surname> <given-names>Z.</given-names></name></person-group> (<year>2020</year>). <article-title>Insulin-like androgenic gland hormone 1 (IAG1) regulates sexual differentiation in a hermaphrodite shrimp through feedback to neuroendocrine factors.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>303</volume>:<issue>113706</issue>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2020.113706</pub-id> <pub-id pub-id-type="pmid">33359802</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Manor</surname> <given-names>R.</given-names></name> <name><surname>Weil</surname> <given-names>S.</given-names></name> <name><surname>Oren</surname> <given-names>S.</given-names></name> <name><surname>Glazer</surname> <given-names>L.</given-names></name> <name><surname>Aflalo</surname> <given-names>E. D.</given-names></name> <name><surname>Ventura</surname> <given-names>T.</given-names></name><etal/></person-group> (<year>2007</year>). <article-title>Insulin and gender: an insulin-like gene expressed exclusively in the androgenic gland of the male crayfish.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>150</volume> <fpage>326</fpage>&#x2013;<lpage>336</lpage>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2006.09.006</pub-id> <pub-id pub-id-type="pmid">17094989</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Martin</surname> <given-names>G.</given-names></name> <name><surname>Sorokine</surname> <given-names>O.</given-names></name> <name><surname>Moniatte</surname> <given-names>M.</given-names></name> <name><surname>Vandorsselaer</surname> <given-names>A.</given-names></name></person-group> (<year>1998</year>). <article-title>The androgenic hormone of the crustacean isopod <italic>Armadillidium vulgare</italic>.</article-title> <source><italic>Ann. NY Acad. Sci.</italic></source> <volume>839</volume> <fpage>111</fpage>&#x2013;<lpage>117</lpage>. <pub-id pub-id-type="doi">10.1111/j.1749-6632.1998.tb10741.x</pub-id></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Martin</surname> <given-names>G.</given-names></name> <name><surname>Sorokine</surname> <given-names>O.</given-names></name> <name><surname>Moniatte</surname> <given-names>M.</given-names></name> <name><surname>Bulet</surname> <given-names>P.</given-names></name> <name><surname>Hetru</surname> <given-names>C.</given-names></name> <name><surname>Van Dorsselaer</surname> <given-names>A.</given-names></name></person-group> (<year>1999</year>). <article-title>The structure of a glycosylated protein hormone responsible for sex determination in the isopod, <italic>Armadillidium vulgare</italic>.</article-title> <source><italic>Eur. J. Biochem.</italic></source> <volume>262</volume> <fpage>727</fpage>&#x2013;<lpage>736</lpage>. <pub-id pub-id-type="doi">10.1046/j.1432-1327.1999.00442.x</pub-id> <pub-id pub-id-type="pmid">10411634</pub-id></citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Michalakis</surname> <given-names>K.</given-names></name> <name><surname>Mintziori</surname> <given-names>G.</given-names></name> <name><surname>Kaprara</surname> <given-names>A.</given-names></name> <name><surname>Tarlatzis</surname> <given-names>B. C.</given-names></name> <name><surname>Goulis</surname> <given-names>D. G.</given-names></name></person-group> (<year>2013</year>). <article-title>The complex interaction between obesity, metabolic syndrome and reproductive axis: a narrative review.</article-title> <source><italic>Metabolism</italic></source> <volume>62</volume> <fpage>457</fpage>&#x2013;<lpage>478</lpage>. <pub-id pub-id-type="doi">10.1016/j.metabol.2012.08.012</pub-id> <pub-id pub-id-type="pmid">22999785</pub-id></citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Montiel-Arzate</surname> <given-names>A.</given-names></name> <name><surname>S&#x00E1;nchez-Castrej&#x00F3;n</surname> <given-names>E.</given-names></name> <name><surname>Camacho-Jim&#x00E9;nez</surname> <given-names>L.</given-names></name> <name><surname>D&#x00ED;az</surname> <given-names>F.</given-names></name> <name><surname>Ponce-Rivas</surname> <given-names>E.</given-names></name></person-group> (<year>2020</year>). <article-title>Effect of recombinant crustacean hyperglycemic hormones rCHH-B1 and rCHH-B2 on lipid metabolism in the Pacific white shrimp <italic>Litopenaeus vannamei</italic>.</article-title> <source><italic>Aquac. Res.</italic></source> <volume>51</volume> <fpage>4267</fpage>&#x2013;<lpage>4278</lpage>. <pub-id pub-id-type="doi">10.1111/are.14769</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nagaraju</surname> <given-names>G. P. C.</given-names></name></person-group> (<year>2011</year>). <article-title>Reproductive regulators in decapod crustaceans: an overview.</article-title> <source><italic>J. Exp. Biol.</italic></source> <volume>214</volume> <fpage>3</fpage>&#x2013;<lpage>16</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.047183</pub-id> <pub-id pub-id-type="pmid">21147963</pub-id></citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Neirijnck</surname> <given-names>Y.</given-names></name> <name><surname>Papaioannou</surname> <given-names>M. D.</given-names></name> <name><surname>Nef</surname> <given-names>S.</given-names></name></person-group> (<year>2019</year>). <article-title>The Insulin/IGF system in mammalian sexual development and reproduction.</article-title> <source><italic>Int. J. Mol. Sci.</italic></source> <volume>20</volume>:<issue>4440</issue>. <pub-id pub-id-type="doi">10.3390/ijms20184440</pub-id> <pub-id pub-id-type="pmid">31505893</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Priyadarshi</surname> <given-names>H.</given-names></name> <name><surname>Das</surname> <given-names>R.</given-names></name> <name><surname>Pavan-Kumar</surname> <given-names>A.</given-names></name> <name><surname>Gireesh-Babu</surname> <given-names>P.</given-names></name> <name><surname>Javed</surname> <given-names>H.</given-names></name> <name><surname>Kumar</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2017</year>). <article-title>Silencing and augmentation of IAG hormone transcripts in adult <italic>Macrobrachium rosenbergii</italic> males affects morphotype transformation.</article-title> <source><italic>J. Exp. Biol.</italic></source> <volume>220</volume> <fpage>4101</fpage>&#x2013;<lpage>4108</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.163410</pub-id> <pub-id pub-id-type="pmid">28851817</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reinecke</surname> <given-names>M.</given-names></name></person-group> (<year>2010</year>). <article-title>Insulin-like growth factors and fish reproduction.</article-title> <source><italic>Biol. Reprod.</italic></source> <volume>82</volume> <fpage>656</fpage>&#x2013;<lpage>661</lpage>. <pub-id pub-id-type="doi">10.1095/biolreprod.109.080093</pub-id> <pub-id pub-id-type="pmid">19864315</pub-id></citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rosen</surname> <given-names>O.</given-names></name> <name><surname>Manor</surname> <given-names>R.</given-names></name> <name><surname>Weil</surname> <given-names>S.</given-names></name> <name><surname>Gafni</surname> <given-names>O.</given-names></name> <name><surname>Linial</surname> <given-names>A.</given-names></name> <name><surname>Aflalo</surname> <given-names>E. D.</given-names></name><etal/></person-group> (<year>2010</year>). <article-title>A sexual shift induced by silencing of a single insulin-like gene in crayfish: ovarian upregulation and testicular degeneration.</article-title> <source><italic>PLoS One</italic></source> <volume>5</volume>:<issue>e15281</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0015281</pub-id> <pub-id pub-id-type="pmid">21151555</pub-id></citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shi</surname> <given-names>L.</given-names></name> <name><surname>Han</surname> <given-names>S.</given-names></name> <name><surname>Fei</surname> <given-names>J.</given-names></name> <name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Ray</surname> <given-names>J. W.</given-names></name> <name><surname>Wang</surname> <given-names>W.</given-names></name><etal/></person-group> (<year>2019</year>). <article-title>Molecular characterization and functional study of insulin-like androgenic gland hormone gene in the red swamp crayfish, <italic>Procambarus clarkii</italic>.</article-title> <source><italic>Genes</italic></source> <volume>10</volume>:<issue>645</issue>. <pub-id pub-id-type="doi">10.3390/genes10090645</pub-id> <pub-id pub-id-type="pmid">31455039</pub-id></citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shibata</surname> <given-names>N.</given-names></name> <name><surname>Yoshikuni</surname> <given-names>M.</given-names></name> <name><surname>Nagahama</surname> <given-names>Y.</given-names></name></person-group> (<year>1993</year>). <article-title>Vitellogenin incorporation into oocytes of rainbow trout, <italic>Oncorhynchus mykiss</italic>, in vitro: effect of hormones on denuded oocytes.</article-title> <source><italic>Dev. Growth Diff.</italic></source> <volume>35</volume> <fpage>115</fpage>&#x2013;<lpage>121</lpage>. <pub-id pub-id-type="doi">10.1111/j.1440-169X.1993.00115.x</pub-id></citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Veenstra</surname> <given-names>J. A.</given-names></name></person-group> (<year>2020</year>). <article-title>Gonadulins, the fourth type of insulin-related peptides in decapods.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>296</volume>:<issue>113528</issue>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2020.113528</pub-id> <pub-id pub-id-type="pmid">32526328</pub-id></citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ventura</surname> <given-names>T.</given-names></name> <name><surname>Manor</surname> <given-names>R.</given-names></name> <name><surname>Aflalo</surname> <given-names>E. D.</given-names></name> <name><surname>Weil</surname> <given-names>S.</given-names></name> <name><surname>Raviv</surname> <given-names>S.</given-names></name> <name><surname>Glazer</surname> <given-names>L.</given-names></name><etal/></person-group> (<year>2009</year>). <article-title>Temporal silencing of an androgenic gland-specific insulin-like gene affecting phenotypical gender differences and spermatogenesis.</article-title> <source><italic>Endocrinology</italic></source> <volume>150</volume> <fpage>1278</fpage>&#x2013;<lpage>1286</lpage>. <pub-id pub-id-type="doi">10.1210/en.2008-0906</pub-id> <pub-id pub-id-type="pmid">18988670</pub-id></citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ventura</surname> <given-names>T.</given-names></name> <name><surname>Manor</surname> <given-names>R.</given-names></name> <name><surname>Aflalo</surname> <given-names>E. D.</given-names></name> <name><surname>Weil</surname> <given-names>S.</given-names></name> <name><surname>Rosen</surname> <given-names>O.</given-names></name> <name><surname>Sagi</surname> <given-names>A.</given-names></name></person-group> (<year>2012</year>). <article-title>Timing sexual differentiation: full functional sex reversal achieved through silencing of a single insulin-like gene in the prawn, <italic>Macrobrachium rosenbergii</italic>1.</article-title> <source><italic>Biol. Reprod.</italic></source> <volume>86</volume>:<issue>90</issue>. <pub-id pub-id-type="doi">10.1095/biolreprod.111.097261</pub-id> <pub-id pub-id-type="pmid">22133694</pub-id></citation></ref>
<ref id="B45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Webster</surname> <given-names>S. G.</given-names></name> <name><surname>Keller</surname> <given-names>R.</given-names></name> <name><surname>Dircksen</surname> <given-names>H.</given-names></name></person-group> (<year>2012</year>). <article-title>The CHH-superfamily of multifunctional peptide hormones controlling crustacean metabolism, osmoregulation, moulting, and reproduction.</article-title> <source><italic>Gen. Comp. Endocrinol.</italic></source> <volume>175</volume> <fpage>217</fpage>&#x2013;<lpage>233</lpage>. <pub-id pub-id-type="doi">10.1016/j.ygcen.2011.11.035</pub-id> <pub-id pub-id-type="pmid">22146796</pub-id></citation></ref>
<ref id="B46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wu</surname> <given-names>Q.</given-names></name> <name><surname>Brown</surname> <given-names>M. R.</given-names></name></person-group> (<year>2006</year>). <article-title>Signaling and function of insulin-like peptides in insects.</article-title> <source><italic>Annu. Rev. Entomol.</italic></source> <volume>51</volume> <fpage>1</fpage>&#x2013;<lpage>24</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.ento.51.110104.151011</pub-id> <pub-id pub-id-type="pmid">16332201</pub-id></citation></ref>
<ref id="B47"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>D.</given-names></name> <name><surname>Sun</surname> <given-names>M.</given-names></name> <name><surname>Liu</surname> <given-names>X.</given-names></name></person-group> (<year>2017</year>). <article-title>Phase-specific expression of an insulin-like androgenic gland factor in a marine shrimp <italic>Lysmata wurdemanni</italic>: implication for maintaining protandric simultaneous hermaphroditism.</article-title> <source><italic>PLoS One</italic></source> <volume>12</volume>:<issue>e0172782</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0172782</pub-id> <pub-id pub-id-type="pmid">28235100</pub-id></citation></ref>
<ref id="B48"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>X.</given-names></name> <name><surname>Pan</surname> <given-names>L.</given-names></name> <name><surname>Wei</surname> <given-names>C.</given-names></name> <name><surname>Tong</surname> <given-names>R.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Ding</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2020</year>). <article-title>Crustacean hyperglycemic hormone (CHH) regulates the ammonia excretion and metabolism in white shrimp, <italic>Litopenaeus vannamei</italic> under ammonia-N stress.</article-title> <source><italic>Sci. Total Environ.</italic></source> <volume>723</volume>:<issue>138128</issue>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2020.138128</pub-id> <pub-id pub-id-type="pmid">32222513</pub-id></citation></ref>
<ref id="B49"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zmora</surname> <given-names>N.</given-names></name> <name><surname>Chung</surname> <given-names>J. S.</given-names></name></person-group> (<year>2014</year>). <article-title>A novel hormone is required for the development of reproductive phenotypes in adult female crabs.</article-title> <source><italic>Endocrinology</italic></source> <volume>155</volume> <fpage>230</fpage>&#x2013;<lpage>239</lpage>. <pub-id pub-id-type="doi">10.1210/en.2013-1603</pub-id> <pub-id pub-id-type="pmid">24280057</pub-id></citation></ref>
</ref-list><fn-group>
<fn id="footnote1">
<label>1</label>
<p><ext-link ext-link-type="uri" xlink:href="https://www.ncbi.nlm.nih.gov/orffinder/">https://www.ncbi.nlm.nih.gov/orffinder/</ext-link></p></fn>
<fn id="footnote2">
<label>2</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.cbs.dtu.dk/services/SignalP/">http://www.cbs.dtu.dk/services/SignalP/</ext-link></p></fn>
<fn id="footnote3">
<label>3</label>
<p><ext-link ext-link-type="uri" xlink:href="http://clavius.bc.edu/&#x007e;clotelab/DiANNA/">http://clavius.bc.edu/&#x007e;clotelab/DiANNA/</ext-link></p></fn>
<fn id="footnote4">
<label>4</label>
<p><ext-link ext-link-type="uri" xlink:href="https://www.ebi.ac.uk/Tools/msa/clustalo/">https://www.ebi.ac.uk/Tools/msa/clustalo/</ext-link></p></fn>
<fn id="footnote5">
<label>5</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.cbs.dtu.dk/services/NetNGlyc/">http://www.cbs.dtu.dk/services/NetNGlyc/</ext-link></p></fn>
</fn-group>
</back>
</article>