Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2011; 6(9): e24427.
Published online 2011 September 7. doi:  10.1371/journal.pone.0024427
PMCID: PMC3168472

Insights into Eyestalk Ablation Mechanism to Induce Ovarian Maturation in the Black Tiger Shrimp

Irina Agoulnik, Editor


Eyestalk ablation is commonly practiced in crustacean to induce ovarian maturation in captivity. The molecular mechanism of the ablation has not been well understood, preventing a search for alternative measures to induce ovarian maturation in aquaculture. This is the first study to employ cDNA microarray to examine effects of eyestalk ablation at the transcriptomic level and pathway mapping analysis to identify potentially affected biological pathways in the black tiger shrimp (Penaeus monodon). Microarray analysis comparing between gene expression levels of ovaries from eyestalk-intact and eyestalk-ablated brooders revealed 682 differentially expressed transcripts. Based on Hierarchical clustering of gene expression patterns, Gene Ontology annotation, and relevant functions of these differentially expressed genes, several gene groups were further examined by pathway mapping analysis. Reverse-transcriptase quantitative PCR analysis for some representative transcripts confirmed microarray data. Known reproductive genes involved in vitellogenesis were dramatically increased during the ablation. Besides these transcripts expected to be induced by the ablation, transcripts whose functions involved in electron transfer mechanism, immune responses and calcium signal transduction were significantly altered following the ablation. Pathway mapping analysis revealed that the activation of gonadotropin-releasing hormone signaling, calcium signaling, and progesterone-mediated oocyte maturation pathways were putatively crucial to ovarian maturation induced by the ablation. These findings shed light on several possible molecular mechanisms of the eyestalk ablation effect and allow more focused investigation for an ultimate goal of finding alternative methods to replace the undesirable practice of the eyestalk ablation in the future.


Eyestalk ablation has been employed to induce reproductive maturation in crustacean [1]. While the ablation can induce ovarian maturation, it also jeopardizes growth, shortens molting cycle, increases energetic demands [2], and resulting in an eventual loss in egg quality and high mortality [3]. Accordingly, predictable maturation and spawning in captive shrimp without the eyestalk ablation is a long–term goal for the shrimp industry [4].

Albeit undesirable effects, the eyestalk ablation is currently unavoidable to solve a serious problem of poor reproductive maturation in captivity of an economically important black tiger shrimp (Penaeus monodon). The reproductive biology of both male and female of this animal still remains poorly understood. In female, according to the gonadosomatic index (GSI; ratio of gonad weight to body weight) and physiological changes [5], the ovarian development can be separated into four stages: pre-vitellogenic, vitellogenic, early cortical rod, and spent (late cortical rod). Unlike ovarian development in natural habitat, the ovarian maturation of P. kerathurus [6] and P. monodon [7] does not naturally reach the later stages of ovarian maturation in captivity.

Crustacean eyestalk is a location for the X-organ sinus gland which is important for neuroendocrine system [8], [9]. In the spiny lobster Panulirus argus, the sinus gland seemed to contain a heat-stable factor inhibiting gonadal maturation [10]. In P. monodon, one of the important hormones secreted from this gland is gonad-inhibiting hormone (GIH) which is known to inhibit synthesis of a yolk protein precursor, vitellogenin [11]. Besides the GIH evidence, another hypothesis suggests that removing eyestalk also reduces light intensity and thereby inducing ovarian maturation. In the banana prawn P. merguiensis, dim light was reported to favor ovarian maturation and spawning [12]. To date, the exact mechanism of eyestalk ablation on the ovarian maturation is not conclusive.

To avoid the use of eyestalk ablation, understanding molecular effects of the ablation is essential for finding alternative measures to trigger ovarian maturation without the undesirable effects from the ablation. DNA microarray is a transcriptomic approach to study responses from an organism to stimuli. This high-throughput technique provides a global picture of gene expression patterns. Previously, we have successfully constructed and employed cDNA microarrays to study ovarian and testicular development in P. monodon [13], [14], [15]. In this study, we employed the cDNA microarrays to shed light on molecular mechanism of the eyestalk ablation effects by comparing gene expression levels of ovaries from non-ablated and ablated female broodstock over the course of seven days after the ablation. In addition, biological pathway information was used to identify putatively important pathways affected by the eyestalk ablation. The reverse-transcriptase quantitative PCR analysis was performed to validate the gene expression results and emphasize the potential involvement of these genes in the mechanism of ovarian maturation induced by the eyestalk ablation.


Induction of ovarian maturation by eyestalk ablation

To examine effects of eyestalk ablation on reproductive maturation of domesticated female black tiger shrimp broodstock at a molecular level, an eyestalk ablation experiment was conducted (Fig. 1A). Growth rates (length and body weight), ovary weight and gonadosomatic index (GSI) of the female shrimp significantly increased after the eyestalk ablation (Table 1). Length and body weight increased steadily over the eyestalk ablation time course and significantly increased from 17.3±0.5 cm and 62.6±5.5 g on Day 0 to 18.6±0.8 cm and 74.2±8.9 g on Day 7, respectively. Ovary weight increased from 0.7±0.2 to 3.6±1.9 g after seven days of the ablation. The significant increase in ovary weight during Day 4 to Day 7 indicated rapid maturation during that period after the ablation. In addition, GSI values, a known indicator of ovarian maturation in the female black tiger shrimp, significantly increased from 1.1±0.2% to 4.7±2.3% by Day 7. In addition to GSI values, the physiological changes and colors of ovaries indicated rapid maturation after the eyestalk ablation (Fig. 1B). The size and color of ovaries from the non-ablated shrimp were small and white, whereas those from the ablated shrimp increased in size and turned from white to yellow at Day 1, to light green at Day 4, and to dark green at Day 7 (Fig. 1B, lower panel).

Figure 1
Diagram of the eyestalk ablation experiment and physiology of shrimp and ovary.
Table 1
Summary of shrimp samples used in this study.

Global gene expression changes due to eyestalk ablation

To investigate transcriptomic changes as a result of the eyestalk ablation, a cDNA microarray, UniShrimpChip, was used to compare expression levels before and after the ablation in domesticated female black tiger shrimp broodstock (Fig. 2). From the microarray experiment, 2,743 transcripts present in at least 7 from 8 microarrays with expression levels ≥ median value±1SD (~1.9-fold change) in at least 1 microarray were selected for further analysis (Fig. 2A). From these, 682 features with >2-fold change in at least 4 from 8 microarrays were considered differentially expressed and categorized based on Gene Ontology (GO; Fig. 3).

Figure 2
Gene expression analysis by cDNA microarray.
Figure 3
Gene Ontology (GO) annotation.

GO distribution of the 682 differentially expressed transcripts was compared to that of all 5,568 features on the UniShrimpChip (Fig. 3). According to the three GO categories, the highest numbers of the features belonged to the “biological process” (45.8% and 46.6% for the all array features and for the differentially expressed transcripts, respectively) followed by the “cellular component” (31.6% and 27.1%, respectively) and the “molecular function” (22.6% and 26.3%, respectively).

In the “biological process” category, various subgroups, such as cellular process, metabolic process, localization, biological regulation and biological adhesion, contained noticeably different distributions of the features between the array features and the differentially expressed genes (Fig. 3C).

Although the overall GO distribution of the differentially expressed transcripts (27.1%) was lower than of the array features (31.6%) for the “cellular components” category, the distributions among the subgroups in this category for the array features and the differentially expressed transcripts were mostly similar (Fig. 3D). In this category, transcripts involved in the cell subgroup were predominant, followed by those involved in organelle.

In the “molecular function” category, transcripts involved in binding activity were predominant, followed by those involved in catalytic activity (Fig. 3E). Among the subgroups in this category, distributions of the array features and the differentially expressed transcripts belonged to the transporter activity and electron carrier activity subgroups were the most distinct.

In addition to the GO analysis, the differentially expressed transcripts were grouped according to the similarity in gene expression profiles using Hierarchical cluster method (Fig. 2B). Three expression patterns of late induced (140 features), repressed (272 features) and early induced (49 features) transcripts were further examined as they include some transcripts with relevant functions (Fig. 2C). Cluster I contained genes whose expression levels were higher after the eyestalk ablation at the later days and was so called “Late induced group.” This transcript group includes a known marker for vitellogenesis such as vitellogenin and genes involved in energy metabolism mechanism such as NADH dehydrogenase subunit4l, ATP synthase F0 subunit6, cytochrome C oxidase subunitI, II, III, cytochromeB, and acyl-CoA oxidase. In contrast, genes in Cluster II were expressed lower after the ablation and categorized as “Repressed group.” This group includes genes involved in immune responses such as HLA-B, crustin-like antimicrobial peptide, C-type lectin, antilipopolysaccharide factor, and hemolymph proteinase5. Cluster III or “Early induced group” contained transcripts that were induced within a day after the eyestalk ablation. This group includes many known reproductive and hormonal genes (ovarian lipoprotein receptor, thrombospondin, female steile m3 and ovarian peritrophin), genes involved in calcium signal transduction (calmodulin and calnexin isoform2) and genes involved in developmental processes (innexin2 and thyroid hormone receptor associated protein).

Candidate pathways relevant to ovarian maturation

The differentially expressed transcripts were used to identify pathways in which the genes could be active after the eyestalk ablation. Of the 682 differentially expressed transcripts identified by the microarray analysis, 190 transcripts were automatically mapped on 117 different predicted pathways from the KEGG database, which features large numbers of biological pathways. The numbers of participated genes per pathway ranged from 1 to 16 (Fig. 4). A majority of the predicted pathways (63 pathways) were mapped with only a single participated gene, which cannot provide sufficient information in gene regulatory mechanism. Thus, the pathways with at least two participated genes were considered. According to the number of participated gene in each pathway, the pathways were ranked in the following order: ribosome, oxidative phosphorylation, RNA transport, protein processing in endoplasmic reticulum, glycolysis, fatty acid metabolism, ribosome biogenesis in eukaryotes, focal adhesion, insulin signaling pathway, tight junction pathway, GnRH signaling pathway, calcium signaling pathway, and progesterone-mediated oocyte maturation (Table 2).

Figure 4
Distribution of number of genes assigned to predicted pathways from the KEGG database.
Table 2
The pathways with the highest numbers of participated genes.

Considering the pathway with the highest number of participated genes, the ribosome pathway with 16 differentially expressed genes participated exhibited down-regulated expression patterns during the eyestalk ablation. The example genes in this pathway included ribosomal proteins and elongation factor 1 alpha (EF-1α). The oxidative phosphorylation pathway ranked second with 15 differentially expressed genes exhibiting later induction pattern at Days 4–7 after the ablation, and the participated genes included those encoding for protein subunits of the five enzyme complexes from electron transport chain such as NADH dehydrogenase subunit4l, ATP synthase F0 subunit6, cytochrome C oxidase subunitI, II, III, VII, and cytochromeB (Table 3).

Table 3
List of genes participated in the putative pathways analyzed by KASS.

In addition to the top ranked pathways, three pathways (GnRH signaling pathway, calcium signaling pathway, and progesterone-mediated oocyte maturation pathway) were selected for mapping analysis using GenMAPP software because the early induced expression patterns of their participated genes after the ablation suggesting that these pathways might be potentially important as molecular responses to the eyestalk ablation. Moreover, these pathways have previously been known as a potential signaling pathway active during ovarian maturation in rats and crustaceans [16], [17].

Therefore, the biological relationships among the three pathways were investigated together by visualizing the expression patterns of the participated gene expression (Fig. 5). A total of 11 genes involved in these signaling pathways and with expression data from microarray were guanine nucleotide-binding protein G(q) subunit alpha (G q/11), calreticulin (CALR), calmodulin (CaM), calcineurinB (CaN), Ca2+ transporting ATPase (PMCA), guanine nucleotide-binding protein G(i) subunit alpha (G), adenylate cyclase (AC), cell division cycle25 (Cdc25), cell division cycle2 (Cdc2), cyclinB (CycB), and mitogen-activated protein kinase (MAPK). All of these transcripts in the putative pathways were up-regulated one day after the ablation.

Figure 5
Putative pathways affected by the eyestalk ablation.

By integrating gene expression data and retrieving biological pathway information affected by eyestalk ablation, major cellular pathway interactions in ovarian maturation of P. monodon were proposed (Fig. 6). The previtellogenic oocytes were induced early by activation of the gonadotropin-releasing hormone (GnRH) signaling pathway and the calcium signaling pathway, but inhibition of gonad- or vitellogenesis- inhibiting hormones (GIH or VIH) upon eyestalk ablation. Changes in signaling molecules of these biological pathways stimulated the release of progesterone leading oocyte maturation. The increased levels of transcripts including vitellogenin and energy-production genes seemed to drive the development of the previtellogenic to matured oocytes with germinal vesicle breakdown (GVBD) before spawning.

Figure 6
The proposed scheme to induce ovarian maturation by eyestalk ablation in P. monodon.

Validation of gene expression by reverse-transcriptase quantitative PCR (RT-qPCR)

Based on the different gene expression patterns of microarray analysis, pathway mapping and gene functions potentially involved in reproductive maturation, gene expression patterns of 12 transcripts were analyzed using reverse-transcriptase quantitative PCR (RT-qPCR) (Table 4). Of these, nine genes were selected to represent the three gene expression patterns identified by microarray (Fig. 2C): Late induced (NADH dehydrogenase, cytochrome C oxidase, ATP synthase F0, and vitellogenin), Repressed (crustin-like antimicrobial peptide, C-type lectin, and HLA-B), and Early induced (innexin2, and calmodulin). Since calmodulin was one of the differentially expressed transcripts with the early induced pattern and was mapped into calcium signaling pathway (Fig. 5), two additional genes (calcineurinB and calmodulin kinase) found in the same pathway were further examined to confirm the involvement of this pathway. Additionally, due to potentially relevant functions in reproductive maturation and the early induced expression pattern from the microarray, thyroid hormone receptor associated protein was also included in the RT-qPCR analysis.

Table 4
Primers used in Reverse-transcriptase Quantitative PCR (RT-qPCR) analysis.

For the late induced genes, the NADH dehydogenase, cytochrome C oxidase, ATP synthase F0, and vitellogenin transcripts exhibited significant increases in gene expression levels at Day 7 when compared to earlier time points and before eyestalk ablation (Fig. 7A). Among the four examined genes, at Day 7, the vitellogenin transcript level was dramatically increased by 240-fold from its initial level before the ablation (D0). Similarly, the cytochrome C oxidase and ATP synthase F0 transcript levels were significantly expressed highest at Day 7 after the ablation with ~58- and 67-fold increases, respectively. Although the induction of the NADH dehydrogenase transcript was ~2 fold which was not as high as the other genes examined from this group, it was significantly induced at Day 7 (P<0.05).

Figure 7
Reverse-transcriptase quantitative PCR (RT-qPCR) analysis.

For the repressed group, the crustin-like antimicrobial peptide, C-type lectin, and HLA-B transcripts were selected for RT-qPCR analysis. All three genes were significantly repressed at either Day 4 or Day 7 (P<0.05). The decreasing pattern of expression levels was in agreement with the microarray data. At Day 7, the expression levels of these transcripts were ~2–5 folds lower than those before the eyestalk ablation (Fig. 7B).

The transcript levels of innexin2, thyroid hormone receptor associated protein, calmodulin, calmodulin kinase, and calcineurinB from the RT-qPCR analysis confirmed their early induction patterns determined by the microarray data as well (Fig. 7C). The expression levels were significantly induced to ~3-, 8-, 4-, 17- and 12-fold, respectively, by one day after the eyestalk ablation (P<0.05). After the significant induction of these genes at Day 1, their expression levels were gradually reduced on the later days. At the end of the ablation experiment, some of the expression levels of these transcripts (innexin2, calmodulin, and calmodulin kinase) were even lower than their initial levels (D0).


Eyestalk ablation is commonly practiced to induce ovarian maturation in the black tiger shrimp farming [1]. This study aimed to understand effects of the eyestalk ablation in female domesticated broodstock at a molecular level by using cDNA microarray to rapidly screen for transcripts differentially expressed and using biological pathway information to reveal pathways potentially involved in the molecular response to the ablation.

In this study, based on the increased gonadosomatic index (GSI; ratio of gonad weight to body weight) and physiological changes from the bigger size and greener color of the ovaries, the unilateral eyestalk ablation evidently induced ovarian maturation within seven days. The GSI value is a good indicator of female reproductive maturation [18] in penaeid shrimp such as in Penaeus setiferus [19], P. duorarum [20], and P. monodon [5]. Based on the GSI values, physiological changes such as ovary color, and morphology of ovarian cells, ovarian maturation in penaeid shrimp can be categorized into four stages: pre-vitellogenic (undeveloped), vitellogenesis (developing), cortical rod (nearly ripe), and late cortical rod (ripe) [5]. From the result, by seven days after the eyestalk ablation, the ovaries were induced to mature to the late cortical rod stage.

From the microarray result, the eyestalk ablation in the female brooders altered expression levels of several important genes group. First and foremost, after the eyestalk ablation, the significant increase in the vitellogenin transcript level at the time when the color of the ovaries turned green and the GSI value reached the mature stage of ovary confirmed that the ablation led to vitellogenesis, which is the accumulation of a major yolk protein marking the start of ovarian development [21]. Vitellogenin has also been known as one of the most reliable markers correlated to ovarian maturation as its gene expression positively correlated to the ovarian maturation levels in Marsupenaeus japonicus [22] and in P. monodon [13].

To find potential biomarkers whose expression levels correlated to ovarian maturation levels, genes with similar expression patterns to that of the vitellogenin transcript were considered. It is interesting to find that genes whose functions relevant to energy metabolism, such as NADH dehydrogenase, cytochrome C oxidase, and ATP synthase F0, share similar gene expression patterns to that of the vitellogenin transcript. The induced expression of energy-related genes may imply that there is an increased flux of energy during the ablation. A previous study revealed that eyestalk ablation resulted in lower glucose levels in hemolymph of female Litopenaeus vannamei [23] and increasing energy requirement to maintain glucose level in hemolymph during gonadal maturity in P. notialis [24]. The increasing energy required after the ablation supports the up-regulation of genes important to the energy production process in our study. The elevated expression levels of these energy metabolism transcripts in this study may help improving farming practice by emphasizing the important of high energy source feed for shrimp broodstock during vitellogenesis. Vitellogenesis undoubtedly requires high energy to synthesize a female-specific lipoprotein called vitellogenin, which carries lipids and proteins for egg formation [25]. As most of lipid classes have a crucial role for energy metabolism [26], it is not surprising to see an increase in ovarian lipids level in maturing ovaries in P. japonicus [27]. Better ovarian development was also observed when feed was supplemented with lipid rich components [28], [29], such as polychaete [30], [31], [32], [33], [34]. Additionally, some previous reports showed that lipid composition of P. monodon is influenced by the lipid composition of its food [35], and phospholipid and PUFA are required for reproduction in P. vannemei [28], [36].

To see the immediate effect of eyestalk ablation, transcripts exhibiting early induction pattern after the ablation were examined through pathway mapping analysis. These transcripts participated in the GnRH signaling, calcium signaling, and progesterone-mediated oocyte maturation pathways. These activated pathways from our results agree with previous studies of induced oocyte maturation mechanism in penaeid shrimp [37], [38], [39], in human [40], and in African clawed frog [41]. For instance, it has been suggested that GnRH-like peptides were found in P. monodon [42] and other crustaceans such as Marsupenaeus japonicus [43] and Macrobrachium rosenbergii [44]. Moreover, eyestalk ablation has previously been found to induce ovarian maturation by inhibiting the GIH in P. notialis [24], and double-stranded RNA of GIH has been successfully used to reduce the GIH transcript level resulting in induction of ovarian maturation in P. monodon [11]. Therefore, the GnRH signaling pathway might exist and involve in the P. monodon ovarian maturation.

In other organisms and crustaceans, the GnRH pathway controlled through calcium levels via calmodulin and calnexin can regulate hormone release, and in turn, trigger lipoprotein transporters for vitellogenesis [16], [17]. Calnexin, similar to calrecticulin, is one of the calcium binding chaperones [45] involved in oogenesis in mammalian cells [46] by regulating GnRH receptor in human [16]. Calmodulin is a calcium binding protein involved in the regulation of steroidogenesis and oocyte maturation in rat [47], in amphibian [48], [49], in fish [50], and in crustacean [51] by activating calmodulin kinase [52] leading to the release of gonadotropins which in turn activate folliculogenesis, estradiol and progesterone production. Then, high level of progesterone synthesized and released by follicular cell is postulated to induce oocyte maturation through progesterone-mediated oocyte maturation pathway resulting in activation of maturation-promoting factor (MPF) consisting of cell division cycle2 and cyclinB. MPF has an important role in promoting oocyte maturation into previtellogenic stage. This molecular change presented in early induced transcripts could be supported by physiological growth of oocyte in P. monodon [53]. Although the actual molecular mechanism of the involvement of the GnRH signaling pathway in the ovarian maturation is not known, the results from this study may suggest that the eyestalk ablation perhaps induce ovarian maturation through the peptide hormone, GnRH, which activates various molecules such as G q/11, calmodulin, calmodulin kinase, and calcineurinB to increase intracellular calcium ion (Ca2+), and then progesterone is released to induce oocyte maturation in P. monodon. However, this hypothesis must be further examined to reveal the exact molecular mechanism to induce ovarian maturation.

Besides the genes participated in these pathways, ovarian lipoprotein receptor, thrombospondin (TSP), ovarian peritrophin, innexin2 and thyroid hormone receptor associated protein (TRAP) were induced early after the ablation. Their highly induced expression levels and relevant function potentially serve as biomarkers for ovarian maturation. The highly induced ovarian lipoprotein receptor transcript level agrees with its high expression level during vitellogenesis in P. monodon [13], [54] and in P. japonicus [55]. In reproductive maturation, lipoproteins transport lipid source to gonad [25], [56]. Thrombospondin (TSP) and ovarian peritrophin involved in ovarian cell development as their expression levels were induced during oocyte maturation in M. japonicus [57], [58], in P. semisulcatus [59] and in P. monodon [13]. Innexin2 encodes for a transmembrane protein found in epithelial cells [60] and functions as an intercellular communicator between soma cells and germ-line cells in oogenesis of Drosophila malanogaster [61]. TRAP encodes for a transcriptional co-activator binding to the nuclear receptor family that functionally involved in several mechanisms including androgen receptor signaling pathway, thyroid hormone receptor binding and vitamin D receptor binding [62], [63], which is essential for mammary gland development [64] and the progression and maintenance of embryonic development [65]. Therefore, it would not be surprising if TRAP would have a role in developmental processes during ovarian development in P. monodon.

Not only did the ablation significantly induced subsets of transcripts, it also repressed an important group of genes of immune-relevant genes. Indeed, eyestalk ablation has been previously reported to affect immune responses in P. vannamei when by altering the concentrations of several hemolymph metabolites important in modulating the immune response [23]. Similarly, in the Indian spiny lobster (Panulirus homarus), total hemocyte count (THC) and prophenoloxidase (PO) activity were increased within a short period of 2 hours and then decreased at 1 week after the ablation [66]. The increasing hormone gene expression levels in pre-mating processes resulted in the decreasing immune gene expression levels in D. melanogaster [67]. The changes of the immune gene expression levels in various organisms from the previous and our studies suggested a linkage between the immune response with the eyestalk ablation.

In summary, besides the well-known effects of the eyestalk ablation on physiological changes leading to ovarian maturation in the female black tiger shrimp, the ablation leads to several crucial transcriptomic changes. Several transcripts were identified to be potential biomarkers for ovarian maturation. With the power of high-throughput gene expression analysis using cDNA microarray and biological pathway information, several groups of transcripts and pathways were identified as putative mechanism resulted from the eyestalk ablation. The removal of eyestalk results in increasing of cell localization, transporter activity responded to hormonal control and calcium signaling pathways in an early stage. Induction processes stimulate vitellogenesis to occur using transporter and electron carrier molecules to increase energy production. These findings pinpoint several potentially important molecular mechanisms of the eyestalk ablation to be further examined. The results also help guiding farming practice to include high energy source in maturation feed formulation.

Materials and Methods

Animal collection and eyestalk ablation

Female domesticated P. monodon broodstock (14-month-old) were collected from the Shrimp Genetic Improvement Center (SGIC, Surat Thani, Thailand). Brooders were maintained at a biosecure station containing seawater pumped from the Gulf of Thailand at a salinity of 30 ppt with water temperature at around 27°C, and acclimated in tanks with 2 m in diameter for 3 weeks. Ovaries of the brooders were collected before (D0; n = 7) and after unilateral eyestalk ablation for 1, 4, and 7 days (n = 8 for each time point; Table 1 and Fig. 1A). All ovary samples were quickly frozen in liquid nitrogen and stored at −80°C until use.

RNA extraction and cDNA labeling for hybridization

RNA samples were extracted from frozen ovaries using TRI-REAGENT according to manufacturer's instruction (Molecular Research Center, USA). Contaminating genomic DNA was removed by treatment with DNase I at 0.15 U/µg total RNA at 37°C for 30 min. The quality and quantity of the RNA were assessed on an agarose gel electrophoresis and NanoDrop (ND-8000) before subsequent experiments. DNase-treated RNA samples were converted to the first-strand cDNAs and labeled with aminoallyl-dUTP (aa-dUTP; Sigma) using a LabelStar Array kit (Qiagen). The aa-cDNA samples were cleaned up using a Microcon YM-30 filter (Millipore) and resuspended in 6 µl of 0.1 M sodium borate buffer (pH 8.7).

cDNA microarray analysis

Construction of the P. monodon cDNA microarray (UniShrimpChip) used in these experiments was described previously [14], [15]. Briefly, the spotted cDNA microarray contains 5,568 cDNA features from 13 EST libraries of different P. monodon organs. UniShrimpChip was used to investigate effects of eyestalk ablation by comparing expression levels of ovaries before and after the ablation (Table 5). A reference sample (Cy3) was from ovary RNA samples pooled from eyestalk-intact domesticated broodstock (pooled from 7 individuals, D0, GSI = 1.1±0.2%), whereas Cy5 sample was from ovary RNA of individual broodstock after eyestalk ablation at Day 1 (n = 2, GSI = 1.2, 1.2%), Day 4 (n = 4, GSI = 1.5, 1.7, 3.8, 4.2%), and Day 7 (n = 2, GSI = 6.8, 7.4%).

Table 5
Summary of microarray experiment.

The slides were post-processed to block unused surfaced using a BSA-based method and microarray procedures were performed as previously described [13]. Briefly, the aa-cDNA was fluorescently labeled with Cy3- or Cy5- dyes (GE Healthcare) at room temperature for 1 h. The unincorporated dye was removed using a Microcon YM-30 filter and the purified probe was resuspended in 6 µl of TE. The Cy3- and Cy5- samples were mixed together in a hybridization buffer (50% formamide, 0.1% SDS, 5× SSC, 0.2 µg/µl of yeast tRNA and 0.2 µg/µl of salmon sperm DNA) in 30 µL of total reaction volume and denatured at 95°C in the dark for 3 min before applying onto the post-processed slide with a glass cover slip (22×32 mm, Deckglasser) for each experiment. Hybridization was performed overnight (~16 hr) in a sealed hybridization chamber (CORNING) at 42°C. The hybridized arrays were immersed into a pre-warmed washing buffer I (2× SSC and 0.1% SDS) to remove a cover slip before transferred to wash in the same buffer for 5 min with gently shaking. The arrays were washed in a washing buffer II (0.2× SSC) and a washing buffer III (0.05× SSC) for 5 min each at a room temperature with gently shaking. The arrays were then dried by centrifugation at 1,000 rpm for 5 min at room temperature (Allegra X-22 R, Beckman coulter) and scanned on an AXON GenePix 4000B microarray scanner (Axon Instruments, Molecular Devices, Sunnyvale, CA, USA).

Analysis of microarray data

Microarray data were extracted using GenePix Pro version 6.1 software. Relative gene expression levels were obtained by comparing the amounts of transcript present in the experimental samples (Cy5-labeled sample) to a reference sample (Cy3-labeled sample). Only spots with intensities greater than one standard deviation (SD) above the background intensity were further analyzed.

The processed data were normalized within each array by global normalization method, and across arrays, using the Aroma package [68] (available from, run in R project environment ( The microarray data is MIAME compliant and the raw data has been deposited in a MIAME compliant database (NCBIs Gene Expression Omnibus,, with GEO accession number GSE 29025 at [69]). From all features on the microarray, only transcripts present in 7 from 8 microarrays and with expression changes ≥ the median expression values±1SD in one of the arrays were considered further. The differentially expressed transcripts were selected from those with ≥2-fold changes than the median expression±1SD in 4 of 8 microarrays. The differential expressed transcripts were grouped together according to their expression patterns by a Hierarchical clustering analysis using the Cluster 3.0 software and visualizing gene expression profile using the TreeView 1.6 software [70].

Mapping of candidate genes to pathways

The gene expression data from the microarray analysis were used to examine potential roles of these transcripts in physiological and biochemical processes by linking them with the Kyoto Encyclopaedia of Genes and Genomes (KEGG) database using two software programs called KEGG Automatic Annotation Server (KAAS) ( [71] and Gene Map Annotator and Pathway Profiler (GenMAPP) [72]. Briefly, the KASS was used to automatically link the microarray data to existing pathways from the KEGG database. The parameter on KASS was set for ESTs annotation with the bi-directional best hit (BBH) method. Later, the expression values of the differentially expressed genes and the retrieved information of pathways were integrated and visualized by using GenMAPP software comparing expression levels across the eyestalk ablation time course from Days 1–7.

Reverse-transcriptase quantitative PCR (RT-qPCR)

RT-qPCR experiment was performed for 12 transcripts (NADH dehydrogenase, cytochrome C oxidase, ATP synthase F0, vitellogenin, crustin-like antimicrobial peptide, C-type lectin, HLA-B, innexin2, calmodulin, thyroid hormone receptor-associated protein, calcineurinB, and calmodulin kinase) (Table 4). Elongation factor 1 alpha (EF-1α), a housekeeping gene, was used as an internal control. The expression study was performed using templates from individual eyestalk-intact broodstock (D0, n = 5), and eyestalk-ablated broodstock at Days 1, 4 and 7 (n = 5 for each group). Each reaction was performed in a 20-µl total volume containing 2× iQ™ SYBR® Green Supermix (Bio-Rad), 200 ng of first strand cDNA template, and 0.2 µM of primer pair. Cycling parameters were 95°C for 2.5 min; followed by 40 cycles of 95°C for 30 sec, 58°C for 20 sec, and 72°C for 30 sec. The specificity of PCR products was confirmed by melting curve analysis performed from 55°C-95°C with a continuous fluorescent reading with every 0.5°C increment. The relative expression levels of copy number of the target gene to that of a housekeeping gene EF-1α in different sample groups were statistically tested by ANOVA followed by Duncan's new multiple range test or Tukey test (P<0.05) [15].


We would like to thank Ms. Nuchjaree Pisamai and Ms. Sirawan Chatsapsin for their generous helps facilitating the project. We also would like to express our gratitude to Ms. Somjai Wongtripop, Mr. Panomkorn Chaitongsakul and Mr. Annop Kudech (SGIC, Thailand) for the shrimp sample and facility to conduct the eyestalk ablation experiment. We are grateful to Prof. Morakot Tanticharoen and Dr. Kanyawim Kirtikara for their mentorship.


Competing Interests: The authors have declared that no competing interests exist.

Funding: This research is supported by the National Center for Genetic Engineering and Biotechnology (BIOTEC) and National Science and Technology Development Agency (NSTDA). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Bray WA, Lawrence AL, editors. Reproduction of Penaeus species in captivity. Amsterdam, The Netherlands: Elsevier; 1992. pp. 93–170.
2. Browdy CL, Samocha TM. The effect of eyestalk ablation on spawning, molting and mating of Penaeus semisulcatus de Haan. Aquaculture. 1985;49:19–29.
3. Benzie JAH. Penaeid genetics and biotechnology. Aquaculture. 1998;164:23–47.
4. Quackenbush LS. Regulation of vitellogenesis in penaeid shrimp; In: DeLoach P, Dougherty WJ, Davidson MA, editors. Amsterdam: Elsevier; 1991. pp. 125–140.
5. Tan-Fermin JD, Pudadera RA. Ovarian maturation stages of the wild giant tiger prawn, Penaeus monodon Fabricius. Aquaculture. 1989;77:229–242.
6. Rodrguez A. A comparative study of the ovarian development in wild and pond-reared shrimp, Penaeus kerathurus (Forskål, 1775). Aquaculture. 1996;148:63–75.
7. Santiago J, Alfredo C. Successful spawning of cultured Penaeus monodon Fabricius after eyestalk ablation. Aquaculture. 1977;11:185–196.
8. Beltz BS. Crustacean neurohormones. Endocrinology of Selected Invertebrate Types. Invertebrate Endocrinology. 1988;2:235–258.
9. Chang ES, editor. Endocrinology. Amsterdam: Elsevier; 1992. pp. 53–91.
10. Quackenbush LS, Herrnkind WF. Regulation of molt and gonadal development in the spiny lobster, Panulirus argus (Crustacea: Palinuridae): Effect of eyestalk ablation. Comparative Biochemistry and Physiology Part A: Physiology. 1981;69:523–527.
11. Treerattrakool S, Panyim S, Chan SM, Withyachumnarnkul B, Udomkit A. Molecular characterization of gonad-inhibiting hormone of Penaeus monodon and elucidation of its inhibitory role in vitellogenin expression by RNA interference. FEBS J. 2008;275:970–980. [PubMed]
12. Hoang T, Lee SY, Keenan CP, Marsden GE. Ovarian maturation of the banana prawn, Penaeus merguiensis de Man under different light intensities. Aquaculture. 2002;208:159–168.
13. Karoonuthaisiri N, Sittikankeaw K, Preechaphol R, Kalachikov S, Wongsurawat T, et al. ReproArray(GTS): a cDNA microarray for identification of reproduction-related genes in the giant tiger shrimp Penaeus monodon and characterization of a novel nuclear autoantigenic sperm protein (NASP) gene. Comp Biochem Physiol Part D Genomics Proteomics. 2009;4:90–99. [PubMed]
14. Wongsurawat T, Leelatanawit R, Thamniemdee N, Uawisetwathana U, Karoonuthaisiri N, et al. Identification of testis-relevant genes using in silico analysis from testis ESTs and cDNA microarray in the black tiger shrimp (Penaeus monodon). BMC Mol Biol. 2010;11:55. [PMC free article] [PubMed]
15. Leelatanawit R, Uawisetwathana U, Klinbunga S, Karoonuthaisiri N. A cDNA microarray, UniShrimpChip, for identification of genes relevant to testicular development in the black tiger shrimp (Penaeus monodon). BMC Mol Biol. 2011;12:15. [PMC free article] [PubMed]
16. Brothers SP, Janovick JA, Conn PM. Calnexin regulated gonadotropin-releasing hormone receptor plasma membrane expression. J Mol Endocrinol. 2006;37:479–488. [PubMed]
17. Subramoniam T. Mechanisms and control of vitellogenesis in crustaceans. Fisheries Science. 2011;77:1–21.
18. Chu KH, Tam YK, Chung CK, Ng WL. Morphometric relationships and reproductive maturation of the shrimp, Metapenaeus ensis, from commercial catches in Hong Kong. Fisheries Research. 1993;18:187–197.
19. King JE. A study of the reproductive organs of the common marine shrimp, Penaeus setiferus (Linnaeus). Biol Bull. 1948;94:244–262. [PubMed]
20. Cummings WC. Maturation and Spawning of the Pink Shrimp, Penaeus duorarum Burkenroad. Transactions of the American Fisheries Society. 1961;90:462–468.
21. Meusy JJ, Payen GG. Female reproduction in malacostracan Crustacea. Zool Sci. 1988;5:217–265.
22. Okumura T, Kim YK, Kawazoe I, Yamano K, Tsutsui N, et al. Expression of vitellogenin and cortical rod proteins during induced ovarian development by eyestalk ablation in the kuruma prawn, Marsupenaeus japonicus. Comparative Biochemistry and Physiology - Part A: Molecular & Integrative Physiology. 2006;143:246–253. [PubMed]
23. Sainz-Hernández JC, Racotta IS, Dumas S, Hernández-López J. Effect of unilateral and bilateral eyestalk ablation in Litopenaeus vannamei male and female on several metabolic and immunologic variables. Aquaculture. 2008;283:188–193.
24. Rosas C, Fernandez I, Brito R, Diaz-Iglesia E. The effect of eyestalk ablation on the energy balance of the pink shrimp, Penaeus notialis. Comparative Biochemistry and Physiology Part A: Physiology. 1993;104:183–187.
25. Yepiz-Plascencia G, Vargas-Albores F, Higuera-Ciapara I. Penaeid shrimp hemolymph lipoproteins. Aquaculture. 2000;191:177–189.
26. Garrett R, Grisham CM. Biochemistry. Boston, MA, USA: Mary Finch; 2010.
27. Teshima S-I, Kanazawa A, Sasada H. Nutritional value of dietary cholesterol and other sterols to larval prawn, Penaeus japonicus Bate. Aquaculture. 1983;31:159–167.
28. Lytle JS, Lytle TF, Ogle JT. Polyunsaturated fatty acid profiles as a comparative tool in assessing maturation diets of Penaeus vannamei. Aquaculture. 1990;89:287–299.
29. Middleditch BS, Missler SR, Ward DG, McVey JB, Brown A, et al. MATURATION OF PENAEID SHRIMP: DIETARY FATTY ACIDS. Proceedings of the World Mariculture Society. 1979;10:472–476.
30. Benzie JAH. A review of the effect of genetics and environment on the maturation and larval quality of the giant tiger prawn Penaeus monodon. 1997. pp. 69–85. Proceedings of the fish and shellfish Larviculture Symposium Aquaculture.
31. Meunpol O, Iam-Pai S, Suthikrai W, Piyatiratitivorakul S. Identification of progesterone and 17[alpha]-hydroxyprogesterone in polychaetes (Perinereis sp.) and the effects of hormone extracts on penaeid oocyte development in vitro. Aquaculture. 2007;270:485–492.
32. Wouters R, Molina C, Lavens P, Calderón J. Lipid composition and vitamin content of wild female Litopenaeus vannamei in different stages of sexual maturation. Aquaculture. 2001;198:307–323.
33. Harrison KE. The role of nutrition in maturation, reproduction and embryonic development of decapod crustaceans: a review. J Shellfish Res. 1990;9:1–28.
34. Primavera JH. Broodstock of Surgo, Penaeus monodon Fabricius. 1983. Southeast Asian Fisheries Development Center, Extension Manual 7: Aquaculture Department SEAFDEC, Iloilo, Philippines.
35. O'Leary CD, Matthews AD. Lipid class distribution and fatty acid composition of wild and farmed prawn, Penaeus monodon (Fabricius). Aquaculture. 1990;89:65–81.
36. Cahu C, Guillaume JC, Stéphan G, Chim L. Influence of phospholipid and highly unsaturated fatty acids on spawning rate and egg and tissue composition in Penaeus vannamei fed semi-purified diets. Aquaculture. 1994;126:159–170.
37. Yano I. Final oocyte maturation, spawning and maturration in penaeid shrimp. J Exp Mar Biol Ecol. 1995;193:113–118.
38. Ngernsoungnern P, Ngernsoungnern A, Kavanaugh S, Sobhon P, Sower SA, et al. The presence and distribution of gonadotropin-releasing hormone-liked factor in the central nervous system of the black tiger shrimp, Penaeus monodon. General and Comparative Endocrinology. 2008;155:613–622. [PubMed]
39. Subramaniam T, Keller R. A new look at the endocrine regulation of egg maturation in the decapod crustacean. Curr Sci. 1993;65:619–623.
40. Seachord CL, VandeVoort CA, Duffy DM. Adipose differentiation-related protein: a gonadotropin- and prostaglandin-regulated protein in primate periovulatory follicles. Biol Reprod. 2005;72:1305–1314. [PubMed]
41. Schmitt A, Nebreda AR. Signalling pathways in oocyte meiotic maturation. J Cell Sci. 2002;115:2457–2459. [PubMed]
42. Ngernsoungnern A, Ngernsoungnern P, Weerachatyanukul W, Chavadej J, Sobhon P, et al. The existence of gonadotropin-releasing hormone (GnRH) immunoreactivity in the ovary and the effects of GnRHs on the ovarian maturation in the black tiger shrimp Penaeus monodon. Aquaculture. 2008;279:197–203.
43. Amano M, Okumura T, Okubo K, Amiya N, Takahashi A, et al. Biochemical Analysis and lmmunohistochemical Examination of a GnRH-like Immunoreactive Peptide in the Central Nervous System of a Decapod Crustacean, the Kuruma Prawn (Marsupenaeus japonicus). Zoolog Sci. 2009;26:840. [PubMed]
44. Ngernsoungnern P, Ngernsoungnern A, Sobhon P, Sretarugsa P. Gonadotropin-releasing hormone (GnRH) and a GnRH analog induce ovarian maturation in the giant freshwater prawn, Macrobrachium rosenbergii. Invertebrate Reproduction & Development. 2009;53:125–135.
45. Corbett EF, Michalak M. Calcium, a signaling molecule in the endoplasmic reticulum? Trends in Biochemical Sciences. 2000;25:307–311. [PubMed]
46. Okazaki Y, Ohno H, Takase K, Ochiai T, Saito T. Cell surface expression of calnexin, a molecular chaperone in the endoplasmic reticulum. J Biol Chem. 2000;275:35751–35758. [PubMed]
47. Carnegie JA, Tsang BK. The calcium-calmodulin system: participation in the regulation of steroidogenesis at different stages of granulosa cell differentiation. Biol Reprod. 1984;30:515–522. [PubMed]
48. Cicirelli MF, Smith LD. Do calcium and calmodulin trigger maturation in amphibian oocytes? Developmental Biology. 1987;121:48–57. [PubMed]
49. Francisco SK-S, Schuetz AW. Sources of calcium and the involvement of calmodulin during steroidogenesis and oocyte maturation in follicles of Rana pipiens. J Exp Zool. 1987;244:133–143. [PubMed]
50. Benninghoff AD, Thomas P. Involvement of calcium and calmodulin in the regulation of ovarian steroidogenesis in Atlantic croaker (Micropogonias undulatus) and modulation by Aroclor 1254. Gen Comp Endocrinol. 2005;144:211–223. [PubMed]
51. Nakatsuji T, Lee CY, Watson RD. Crustacean molt-inhibiting hormone: structure, function, and cellular mode of action. Comp Biochem Physiol A Mol Integr Physiol. 2009;152:139–148. [PubMed]
52. Kugu K, Dharmarajan AM, Preutthipan S, Wallach EE. Role of calcium/calmodulin-dependent protein kinase II in gonadotrophin-induced ovulation in in vitro perfused rabbit ovaries. J Reprod Fertil. 1995;103:273–278. [PubMed]
53. Kruevaisayawan H, Vanichviriyakit R, Weerachatyanukul W, Withyachumnarnkul B, Chavadej J, et al. Oogenesis and formation of cortical rods in the black tiger shrimp, Penaeus monodon. Aquaculture. 2010;301:91–98.
54. Preechaphol R, Klinbunga S, Yamano K, Menasveta P. Molecular cloning and expression of progestin membrane receptor component 1 (Pgmrc1) of the giant tiger shrimp Penaeus monodon. Gen Comp Endocrinol. 2010;168:440–449. [PubMed]
55. Tsutsui N, Kawazoe I, Ohira T, Jasmani S, Yang WJ, et al. Molecular Characterization of a cDNA Encoding Vitellogenin and Its Expression in the Hepatopancreas and Ovary during Vitellogenesis in the Kuruma Prawn, Penaeus japonicus. Zool Sci. 2000;17:651–660. [PubMed]
56. Ravid T, Tietz A, Khayat M, Boehm E, Michelis R, et al. Lipid Accumulation in the Ovaries of a Marine Shrimp Penaeus Semisulcatus (DE HAAN). J Exp Biol. 1999;202:1819–1829. [PubMed]
57. Qiu G-F, Yamano K, Unuma T. Cathepsin C transcripts are differentially expressed in the final stages of oocyte maturation in kuruma prawn Marsupenaeus japonicus. Comparative Biochemistry and Physiology Part B: Biochemistry and Molecular Biology. 2005;140:171–181. [PubMed]
58. Yamano K, Qiu GF, Unuma T. Molecular cloning and ovarian expression profiles of thrombospondin, a major component of cortical rods in mature oocytes of penaeid shrimp, Marsupenaeus japonicus. Biol Reprod. 2004;70:1670–1678. [PubMed]
59. Khayat M, Babin PJ, Funkenstein B, Sammar M, Nagasawa H, et al. Molecular characterization and high expression during oocyte development of a shrimp ovarian cortical rod protein homologous to insect intestinal peritrophins. Biol Reprod. 2001;64:1090–1099. [PubMed]
60. Bauer R, Lehmann C, Martini J, Eckardt F, Hoch M. Gap junction channel protein innexin 2 is essential for epithelial morphogenesis in the Drosophila embryo. Mol Biol Cell. 2004;15:2992–3004. [PMC free article] [PubMed]
61. Bohrmann J, Zimmermann J. Gap junctions in the ovary of Drosophila melanogaster: localization of innexins 1, 2, 3 and 4 and evidence for intercellular communication via innexin-2 containing channels. BMC Dev Biol. 2008;8:111. [PMC free article] [PubMed]
62. Ito M, Yuan CX, Malik S, Gu W, Fondell JD, et al. Identity between TRAP and SMCC complexes indicates novel pathways for the function of nuclear receptors and diverse mammalian activators. Mol Cell. 1999;3:361–370. [PubMed]
63. Yuan CX, Ito M, Fondell JD, Fu ZY, Roeder RG. The TRAP220 component of a thyroid hormone receptor- associated protein (TRAP) coactivator complex interacts directly with nuclear receptors in a ligand-dependent fashion. Proc Natl Acad Sci USA. 1998;95:7939–7944. [PubMed]
64. Zhu X, Cheng SY. New insights into regulation of lipid metabolism by thyroid hormone. Curr Opin Endocrinol Diabetes Obes. 2010;17:408–413. [PMC free article] [PubMed]
65. Ito M, Yuan CX, Okano HJ, Darnell RB, Roeder RG. Involvement of the TRAP220 component of the TRAP/SMCC coactivator complex in embryonic development and thyroid hormone action. Mol Cell. 2000;5:683–693. [PubMed]
66. Verghese B, Radhakrishnan EV, Padhi A. Effect of moulting, eyestalk ablation, starvation and transportation on the immune response of the Indian spiny lobster, Panulirus homarus. Aquaculture Research. 2008;39:1009–1013.
67. Lawniczak MKN, Barnes AI, Linklater JR, Boone JM, Wigby S, et al. Mating and immunity in invertebrates. Trends in Ecology & Evolution. 2007;22:48–55. [PubMed]
68. Bengtsson H, editor. 2004. {aroma} - {A}n {R} {O}bject-oriented {M}icroarray {A}nalysis environment}: Mathematical Statistics, Centre for Mathematical Sciences, Lund University, Sweden.
69. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–210. [PMC free article] [PubMed]
70. Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95:14863–14868. [PubMed]
71. Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35:W182–185. [PMC free article] [PubMed]
72. Salomonis N, Hanspers K, Zambon AC, Vranizan K, Lawlor SC, et al. GenMAPP 2: new features and resources for pathway analysis. BMC Bioinformatics. 2007;8:217. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science