Improved genome assembly of the whiteleg shrimp Penaeus (Litopenaeus) vannamei using long- and short-read sequences from public databases.

The Pacific whiteleg shrimp Penaeus (Litopenaeus) vannamei is a highly relevant species for the world's aquaculture development, for which an incomplete genome is available in public databases. In this work, PacBio long-reads from 14 publicly available genomic libraries (131.2 Gb) were mined to improve the reference genome assembly. The libraries were assembled, polished using Illumina short-reads, and scaffolded with P. vannamei, Feneropenaeus chinensis, and Penaeus monodon genomes. The reference-guided assembly, organized into 44 pseudo-chromosomes and 15,682 scaffolds, showed an improvement from previous reference genomes with a genome size of 2.055 Gb, N50 of 40.14 Mb, L50 of 21, and the longest scaffold of 65.79 Mb. Most orthologous genes (92.6%) of the Arthropoda_odb10 database were detected as "complete," and BRAKER predicted 21,816 gene models; from these, we detected 1,814 single-copy orthologues conserved across the genomic references for Marsupenaeus japonicus, F. chinensis, and P. monodon. Transcriptomic-assembly data aligned in more than 99% to the new reference-guided assembly. The collinearity analysis of the assembled pseudo-chromosomes against the P. vannamei and P. monodon reference genomes showed high conservation in different sets of pseudo-chromosomes. In addition, more than 21,000 publicly available genetic marker sequences were mapped to single-site positions. This new assembly represents a step forward to previously reported P. vannamei assemblies. It will be helpful as a reference genome for future studies on the evolutionary history of the species, the genetic architecture of physiological and sex-determination traits, and the analysis of the changes in genetic diversity and composition of cultivated stocks.

Identification and expression analysis of transcripts involved in taurine biosynthesis during early ontogeny of tropical gar Atractosteus tropicus.

In fishes, the availability of taurine is regulated during ontogenetic development, where its endogenous synthesis capacity is species dependent. Thus, different pathways and involved enzymes have been described: pathway I (cysteine sulfinate-dependent pathway), cysteine dioxygenase type 1 (cdo1) and cysteine sulfinic acid decarboxylase (csad); pathway II (cysteic acid pathway), cdo1 and glutamic acid decarboxylase (gad); and pathway III (cysteamine pathway), 2-aminoethanethiol dioxygenase (ado); whereas taurine transporter (taut) is responsible for taurine entry into cells on the cell membrane and the mitochondria. This study determined if the tropical gar (Atractosteus tropicus), an ancient holostean fish model, has the molecular mechanism to synthesize taurine through the identification and analysis expression of transcripts coding for proteins involved in its biosynthesis and transportation, at different embryo-larvae stages and in different organs of juveniles (31 dah). We observed a fluctuating expression of all transcripts involved in the three pathways at all analyzed stages. All transcripts are expressed during the beginning of larval development; however, ado and taut show a peak expression at 9 dah, and all transcripts but csad decreased at 23 dah, when the organism ended the larval period. Furthermore, at 31 dah, we observed taut expression in all examined organs. The transcripts involved in pathways I and III are expressed differently across all organs, whereas pathway II was only observed in the brain, eye, and skin. The results suggested that taurine biosynthesis in tropical gar is regulated during its early development before first feeding, and the pathway might also be organ-type dependent.

De novo transcriptome assembly and identification of G-Protein-Coupled-Receptors (GPCRs) in two species of monogenean parasites of fish.

Genomic resources for Platyhelminthes of the class Monogenea are scarce, despite the diversity of these parasites, some species of which are highly pathogenic to their fish hosts. This work aimed to generate de novo-assembled transcriptomes of two monogenean species, Scutogyrus longicornis (Dactylogyridae) and Rhabdosynochus viridisi (Diplectanidae), providing a protocol for cDNA library preparation with low input samples used in single cell transcriptomics. This allowed us to work with sub-microgram amounts of total RNA with success. These transcriptomes consist of 25,696 and 47,187 putative proteins, respectively, which were further annotated according to the Swiss-Prot, Pfam, GO, KEGG, and COG databases. The completeness values of these transcriptomes evaluated with BUSCO against Metazoa databases were 54.1% and 73%, respectively, which is in the range of other monogenean species. Among the annotations, a large number of terms related to G-protein-coupled receptors (GPCRs) were found. We identified 109 GPCR-like sequences in R. viridisi, and 102 in S. longicornis, including family members specific for Platyhelminthes. Rhodopsin was the largest family according to GRAFS classification. Two putative melatonin receptors found in S. longicornis represent the first record of this group of proteins in parasitic Platyhelminthes. Forty GPCRs of R. viridisi and 32 of S. longicornis that were absent in Vertebrata might be potential drug targets. The present study provides the first publicly available transcriptomes for monogeneans of the subclass Monopisthocotylea, which can serve as useful genomic datasets for functional genomic research of this important group of parasites.

Title: Assemblage de novo du transcriptome et identification des récepteurs couplés aux protéines G (RCPG) chez deux espèces de Monogènes parasites de poissons. Abstract: Les ressources génomiques pour les Plathelminthes de la classe Monogenea sont rares, malgré la diversité de ces parasites dont certaines espèces sont hautement pathogènes pour leurs hôtes poissons. Ce travail visait à générer des transcriptomes assemblés de novo pour deux espèces de monogènes, Scutogyrus longicornis (Dactylogyridae) et Rhabdosynochus viridisi (Diplectanidae), fournissant un protocole pour la préparation de la bibliothèque d'ADNc avec des échantillons à faible apport utilisés en transcriptomique unicellulaire, ce qui a permis de travailler avec des quantités inférieures au microgramme d'ARN total avec succès. Ces transcriptomes se composent de 25 696 et 47 187 protéines putatives, respectivement, qui ont ensuite été annotées selon les bases de données Swiss-Prot, Pfam, GO, KEGG et COG. L'exhaustivité de ces transcriptomes évaluée avec BUSCO par rapport aux bases de données des Métazoaires était respectivement de 54,1 % et 73 %, ce qui est dans la gamme des autres espèces de monogènes. Parmi les annotations, un grand nombre de termes liés aux récepteurs couplés aux protéines G (RCPG) ont été trouvés. Nous avons identifié 109 séquences de type RCPG chez R. viridisi et 102 chez S. longicornis, y compris des membres de la famille spécifiques de Platyhelminthes. La rhodopsine était la plus grande famille selon la classification GRAFS. Deux récepteurs putatifs de la mélatonine trouvés chez S. longicornis représentent le premier signalement de ce groupe de protéines chez les Plathelminthes parasites. Quarante RCPG de R. viridisi et 32 de S. longicornis, qui sont absents chez les Vertébrés, pourraient être des cibles médicamenteuses potentielles. La présente sont fournit les premiers transcriptomes accessibles au public pour les monogènes de la sous-classe Monopisthocotylea, qui peuvent servir d'ensembles de données génomiques utiles pour la recherche génomique fonctionnelle de cet important groupe de parasites.

Transcriptomic signaling in zebrafish (Danio rerio) embryos exposed to environmental concentrations of glyphosate.

Glyphosate [N-(phosphonomethyl)glycine] is one of the most popular herbicides worldwide. Globally, the use of glyphosate is increasing, and its residues have been found in drinking water and food products. The data regarding the possible toxic effects of this herbicide are controversial. Therefore, the aim of this study was to evaluate the effects of glyphosate at environmental concentrations in zebrafish (Danio rerio) embryos. Embryos were exposed to 0, 1, 100, and 1,000 µg/L glyphosate for 96 h, and mortality, heart rate, and hatching rate were evaluated. After the experiment, RNA was extracted from the embryos for transcriptional analysis. No mortality was recorded, and exposure to 100 µg/L and 1,000 µg/L of glyphosate resulted in lower heart rates at 48 h. In addition, RNA-seq analysis revealed that glyphosate exposure induced subtle changes in gene transcription profiles. We found 30 differentially expressed genes; however, the highest glyphosate concentration (1,000 µg/L) induced the greatest number of differentially expressed genes involved in oocyte maturation, metabolic processes, histone deacetylation, and nervous system development.

Epigenetic Transgenerational Modifications Induced by Xenobiotic Exposure in Zebrafish.

Zebrafish (Danio rerio) is a well-established vertebrate model in ecotoxicology research that responds to a wide range of xenobiotics such as pesticides, drugs, and endocrine-disrupting compounds. The epigenome can interact with the environment and transform internal and/or external signals into phenotypic responses through changes in gene transcription. Environmental exposures can also generate epigenetic variations in offspring even by indirect exposure. In this review, we address the advantages of using zebrafish as an experimental animal model to study transgenerational epigenetic processes upon exposure to xenobiotics. We focused mostly on DNA methylation, although studies on post-translational modifications of histones, and non-coding RNAs related to xenobiotic exposure in zebrafish are also discussed. A revision of the methods used to study epigenetic changes in zebrafish revealed the relevance and reproducibility for epigenetics-related research. PubMed and Google Scholar databases were consulted for original research articles published from 2013 to date, by using six keywords: zebrafish, epigenetics, exposure, parental, transgenerational, and F2. From 499 articles identified, 92 were considered, of which 14 were selected as included F2 and epigenetic mechanisms. Current knowledge regarding the effect of xenobiotics on DNA methylation, histone modifications, and changes in non-coding RNAs expressed in F2 is summarized, along with key experimental design considerations to characterize transgenerational effects.

Expression of long-chain polyunsaturated fatty acids biosynthesis genes during the early life-cycle stages of the tropical gar Atractosteus tropicus.

Long-chain (≥C20) polyunsaturated fatty acids (LC-PUFA), including eicosapentaenoic acid (EPA, 20:5n-3), arachidonic acid (ARA, 20:4n-6) and docosahexaenoic acid (DHA, 22:6n-3), are essential in multiple physiological processes, especially during early development of vertebrates. LC-PUFA biosynthesis is achieved by two key families of enzymes, fatty acyl desaturases (Fads) and elongation of very long-chain fatty acid (Elovl). The present study determined the expression patterns of genes encoding desaturases (fads1 and fads2) and elongases (elovl2 and elovl5) involved in the LC-PUFA biosynthesis during early life-stages of the tropical gar Atractosteus tropicus. We further analyzed the fatty acid profiles during early development of A. tropicus to evaluate the impact of Fads and Elovl enzymatic activities. Specific oligonucleotides were designed from A. tropicus transcriptome to perform qPCR (quantitative polymerase chain reaction) on embryonic and larval stages, along with several organs (intestine, white muscle, brain, liver, heart, mesenteric adipose, kidney, gill, swim bladder, stomach, and spleen) collected from juvenile specimens. Fatty acid content of feeds and embryonic and larval stages were analyzed. Results show that fads1, fads2, elovl2 and elovl5 expression was detected from embryonic stages with expression peaks from day 15 post hatching, which could be related to transcriptional and dietary factors. Moreover, fads1, fads2 and elovl2 showed a higher expression in intestine, while elovl5 showed a higher expression in liver, suggesting that the tropical gar activates its LC-PUFA biosynthetic machinery to produce ARA, EPA and DHA to satisfy physiological demands at crucial developmental milestones during early development.

De novo transcriptome assembly for two color types of the marine sponge Mycale (Carmia) cecilia.

Despite the ecological importance and the potential pharmacological application of the sponge Mycale (Carmia) cecilia, it is uncertain whether the body-color variation, even in individuals coexisting in the same area, is due to intraspecific phenotypic plasticity or corresponds to taxonomic divergence. This uncertainty is relatively common in several Porifera groups, which lack the resolution of morphological diagnostic characters and slow-evolving mitochondrial genomes as occurs in early splitting lineages. We sequenced the RNA of six individuals with two different body-color (green-morphotype and red-morphotype) collected at the same time side by side. High-throughput sequencing of cDNA libraries produced ~ 129 million reads with a length of 150 bp. Each morphotype was assembled separately owing to the low overlapping in the global assembly. Metatranscriptome de novo assembly of the trimmed and normalized reads produced 461 thousand transcripts for the green-morphotype and 342 thousand for the red-morphotype (respectively). Over 30% of the transcripts contained Open Reading Frames (ORFs) with functional significance. BUSCO analysis of the ORFs of putative poriferan origin (31.3% green or 30.4% red) indicated that our assemblies are 60% complete. This is the first attempt to evaluate the morphological diversity in the species M. (C.) cecilia and the phylum Porifera at the transcriptomic level. Due to the minimum overlap of the assembly and that, the red-morphotype diverged significantly from the green-morphotype (original color of M. (C.) cecilia). Therefore, we suggest that the red-morphotype should undergo a complete taxonomic investigation and its taxonomic status be reviewed. We expect that the transcriptome assembly metrics can be useful for comparing other transcriptome assemblies of non-model organisms.

The complete mitogenome of excavating sponge Thoosa mismalolli (Demospongiae, Tetractinellida, Thoosidae) from Northeastern Tropical Pacific.

The complete mitogenome of Thoosa mismalolli Carballo, Cruz-Barraza & Gómez, 2004 (Tetractinellida, Thoosidae) was sequenced. This is the first complete mitogenome of the suborden Thoosina and the third Tetractinellid so far. The mitochondrial genome of T. mismalolli was assembled based on reads obtained with the Illumina HiSeq platform. The length of complete mitogenome is 19,019 bp long and contained 14 protein-coding genes and 23 tRNA, with two tRNA genes. Phylogenetic reconstruction (maximum-likelihood) based on mitogenome of Tetractinellids, supports T. mismalolli as a sister group. This result is congruent with those obtained with molecular markers (CO1, 18S, and 28S), supporting the monophyletic status of Thoosa and providing additional molecular data in favor of the suborder Thoosina.

Effects of feed restriction on growth performance, lipid mobilization, and gene expression in rose spotted snapper (Lutjanus guttatus).

The effects of feed deprivation were evaluated for 1 week and 2 weeks in Lutjanus guttatus juveniles. A significant reduction in body mass was observed in both feed deprivation schemes, as well as in hepatosomatic, viscerosomatic and mesenteric fat indexes. The composition of fasted fish was characterized by a decrease in lipid content; the liver displayed an intense reduction of lipid reserves in both fasted groups, and increased expression of the lysosomal acid lipase. 1 week after re-feeding, both experimental groups showed an increase in specific growth rate, feed intake and feed conversion ratio. A recovery in hepatic lipid reserves was observed, and the expression of the lysosomal acid lipase decreased, although lipid content in both groups was still significantly lower than in control groups. Hepatic expression of the growth hormone receptor decreased after fasting, and remained low even after the fish were fed, whereas the expression of insulin-like growth factor 1 and 2 increased after fasting and decreased in both groups when fish were fed again. Altogether, these results showed a partial compensatory growth response in L. guttatus juveniles after fasting, with enhanced growth rates and improved feed efficiency. Fish used stored lipid reserves as the main energy source, and the expression of growth-related and lipid mobilization marker genes in the liver showed similar patterns in both fasting schemes. Based on the results, we suggest as much as 1-week fasting intervals during grow-out programmes to reduce visceral fat and increase growth rate in this species.

Post-transcriptional silencing of myostatin-1 in the spotted rose snapper (Lutjanus guttatus) promotes muscle hypertrophy.

Muscle growth is regulated by several factors including the growth differentiation factor 8, known as myostatin, an inhibitor of myocyte differentiation and proliferation. Research on myostatin regulation was already conducted to improve growth rates in farmed animals, including aquatic species. To explore the effects of myostatin inactivation in a commercial marine fish (spotted rose snapper, Lutjanus guttatus) in vivo, we induced post-transcriptional silencing (knockdown) of myostatin-1 (mstn-1) by injecting dsiRNA directly into the muscle of juvenile fish (87 days post-hatch) using a commercial polymer as vehicle. Results show a significant decrease in mstn-1 expression starting at 2 days after injection and for up to 5 days. Knockdown of mstn-1 caused muscle fiber hypertrophy (but not hyperplasia); however, there were no significant changes in weight or length. Although still experimental, this study provides evidence that temporary knockdown of mstn-1 in a commercial marine fish in vivo promotes fiber hypertrophy and therefore could potentially help grow-out programmes in fish aquaculture.

Genome-wide identification of ABC transporters in monogeneans.

ATP-Binding Cassette (ABC) transporters are proteins that actively mediate the transport of a wide variety of molecules, including drugs. Thus, in parasitology, ABC transporters have gained attention as potential targets for therapeutic drugs. Among the parasitic Platyhelminthes, ABC transporters have been identified and classified in a few species of Trematoda and Cestoda but not in Monogenea. Monogeneans are mainly ectoparasites of marine and freshwater fish, although they can also be found on other aquatic organisms. Severe epizootics caused by monogeneans have been reported around the world, mainly in confined and/or overcrowded fish. The purpose of this study was to identify the ABC transporters in four species of monogeneans (Gyrodactylus salaris, Protopolystoma xenopodis, Eudiplozoon nipponicum and Neobenedenia melleni) for which genomic resources are publicly available. For comparative purposes, ABC transporters were also identified in endoparasitic (Schistosoma mansoni and Echinococcus granulosus) and free-living (Macrostomun lignano and Schmidtea mediterranea) platyhelminths. Thirty-two putative ABC transporters were identified in the genome of G. salaris, 40 in the genome of P. xenopodis, 46 in the transcriptome of E. nipponicum and 9 in a rather limited ESTs set available for N. melleni. Of the eight ABC subfamilies (A-H) known in metazoans, subfamily H was the only one not found in any monogenean species. In contrast, ABCC was the best represented subfamily. Phylogenetic analyses showed a few cases of one-to-one orthologous relationships, which agree with results from other metazoan species. We found some monogenean ABC members related to subfamilies B, C and G involved in drug resistance in humans. This information may be useful for future functional studies on ABC transporters in monogeneans.

Overfeeding a High-Fat Diet Promotes Sex-Specific Alterations on the Gut Microbiota of the Zebrafish (Danio rerio).

Diet modulates the gut microbiota and is one of the main factors promoting obesity and overweight. In the present study, we investigated the effect of a high-fat diet (HFD) on the gut microbiota of the zebrafish (Danio rerio). Fish were separated into three groups and fed in different regimes: low fat, high fat, and high fat overfed; the experiments were performed on males and females separately. We analyzed more than 2.6 million sequences of variable region V3 of the 16S rRNA gene generated by the Illumina Miniseq platform, clustered to 97% similarity with vsearch and classified with the EzBioCloud database. The weight gain, condition factor (K), and body mass index were calculated as indicators of obesity. Multivariate analysis (PERMANOVA and ANOSIM) and diversity indices (Shannon and Dominance) revealed that overfeeding a HFD disturbs the gut microbiota differently in males and females suggesting that sex is a significant factor (p < 0.05) for the composition of the gut microbiota of zebrafish. The results also indicate that a HFD provided in a basal caloric regime does not promote obesity or alterations in the gut microbiota.

A natural antisense transcript of the fem-1 gene was found expressed in female gonads during the characterization, expression profile, and cellular localization of the fem-1 gene in Pacific white shrimp Penaeus vannamei.

The fem-1 gene in Caenorhabditis elegans is involved in sex differentiation; it is specifically required for all aspects of male development. In this study, the full-length cDNA of the fem-1 (Pvfem-1) gene was isolated from the Pacific whiteleg shrimp Penaeus vannamei. The Pvfem-1 transcript is 3778 nt long and encodes a putative protein (PvFEM-1) of 638 amino acids that presented eight ankyrin repeats. The translated protein showed a significant (P < 0.05) structural similitude by superposition with C. elegans FEM-1 protein. Pvfem-1 expression was evaluated by qPCR and in situ hybridization (ISH) during embryogenesis, larval development, and gonads of both genders in subadult and adult life stages. Pvfem-1 was found expressed in brain, intestine, hepatopancreas, and in the gonads of both genders in subadults and adults when quantified by RT-qPCR. A significant finding was the discovery of a natural antisense transcript (NAT) of Pvfem-1 by ISH. It was present in the oocyte nucleus of subadult female shrimp gonads but was not seen within oocytes from adult females, although it was detected in follicular cells, suggesting a possible post-transcriptional regulation of Pvfem-1 in female gonad. Conversely, in males, no NAT was observed, and Pvfem-1 was found expressed in spermatogonia of both, subadult and adult shrimps indicating a function in male sexual differentiation and gametes generation. This study represents the first step for future functional analysis that is expected to contribute to clarifying the role of Pvfem-1 in sex differentiation and determination.

SNP markers for the genetic characterization of Mexican shrimp broodstocks.

Selective breeding of shrimp has major potential to enhance production traits, including growth and disease resistance. Genetic characterization of broodstock populations is a key element of breeding programs, as it enables decisions on inbreeding restrictions, family structure, and the potential use of genomic selection. Single Nucleotide Polymorphisms (SNPs) are suitable genetic markers for this purpose. A set of SNPs was developed to characterize commercial breeding stocks in Mexico. Individuals from local and imported lines were selected for sequencing using the nextRAD technique, resulting in the identification of 2619 SNPs. Genetic structure analysis showed three to five genetic groups of Ecuadorian and Mexican origins. A subset of 1231 SNPs has potential for stock identification and management. Further, three SNPs were identified as candidate sex-linked markers. The role of SNPs possibly associated with genes related to traits of importance to shrimp farming, such as growth and immune response, should be further investigated.

Molecular variability and genetic structure of white spot syndrome virus strains from northwest Mexico based on the analysis of genomes.

White spot syndrome virus (WSSV) has a ∼300 kb double-stranded DNA genome. It originated in China, spread rapidly through shrimp farms in Asia, and subsequently to America. This study determined complete genome sequences for nine historic WSSV strains isolated from Pacific white shrimp (Litopenaeus vannamei) captured in farm ponds in northwest Mexico (Sinaloa and Nayarit). Genomic DNA was captured by an amplification method using overlapping long-range PCR and sequencing by Ion Torrent-PGM. Complete genome sequences were assembled (length range 255-290 kb) and comparative genome analysis with WSSV strains revealed substantial deletions (3 and 10 kb in two regions) in seven strains, with two strains differing from the rest. Phylogenetic analysis identified that the WSSV strains from the northern area of the state of Sinaloa clustered with strains from China (LC1, LC10, DVI) and Korea (ACF2, ACF4), while those from the southern region of the state of Nayarit (AC1 and JP) differed from both of those and from strains found in Taiwan and Thailand. Our data offer insights into the diversity of the WSSV genome in one country and their divergent origin, suggest that it entered Mexico via multiple routes and that specific genome regions can accommodate substantial deletions without compromising viability.

Methylation status of the putative Pax6 promoter in olive ridley sea turtle embryos with eye defects: An initial approach.

Normal development involves the interplay of genetic and epigenetic regulatory mechanisms. Pax6 is an eye-selector factor responsible for initiating the regulatory cascade for the development of the eyes. For the olive ridley sea turtle (Lepidochelys olivacea), a threatened species, eye malformations have been reported. In order to study the DNA methylation status of the putative promoter of the Pax6 gene in embryos with ocular malformations, an exploratory study was carried out in which DNA was isolated from embryos with anophthalmia, microphthalmia, and cyclopia, as well as from their normal counterparts. The 5'-flanking region from the Pax6 gene was isolated, showing two CpG islands (CGIs). The methylation status of CGIs in malformed embryos was compared with that of normal embryos by bisulfite sequencing. Putative transcription factor binding sites and regulatory features were identified. Methylation patterns were observed in both CpG and non-CpG contexts, and were unique for each malformed embryo; in the CpG context, an embryo with cyclopia showed a methylated cytosine upstream the CGI-1 not present in other embryos, an embryo with left anophthalmia presented two methylated cytosines in the CGI-1, whereas an embryo with left anophthalmia and right microphthalmia showed two methylated cytosines in the CGI-2. Normal embryos did not show methylated cytosines in the CGI-1, but one of them showed one methylcytosine in the CGI-2. Methylated transcription factor-binding sites may affect Pax6 expression associated to the cellular response to environmental compounds and hypoxia, signal transduction, cell cycle, lens physiology and development, as well as the transcription rate. Although preliminary, these results suggest that embryos with ocular malformations present unique DNA methylation patterns in the putative promoter of the Pax6 gene in L. olivacea, and probably those subtle, random changes in the methylation status can cause (at least in part) the aberrant phenotypes observed in these embryos.

Isolation of the sex-determining gene Sex-lethal (Sxl) in Penaeus (Litopenaeus) vannamei (Boone, 1931) and characterization of its embryogenic, gametogenic, and tissue-specific expression.

The Pacific white shrimp Penaeus vannamei is the most cultured shrimp species around the world. Because females grow larger than males, the culture of 'only females' is of great interest, but knowledge on sex determination and differentiation is required for producing only females. In an effort to obtain information associated with reproduction in P. vannamei, transcriptomic data from female gonads was generated, and partial sequences of a transcript were identified as Sex-lethal (Sxl). Its characterization indicated that, differently from other penaeids in which this gene has been isolated, there are six isoforms of the Sxl transcript in P. vannamei (PvanSxl 1-6). These isoforms result from alternative splicing at three splice sites (SS1, SS2, SS3). The first splice-site is unique to P. vannamei, as it has not been reported for other Arthropod species; the second splice-site (SS2) is common among crustaceans, and the third splice-site (SS3) is also unique to P. vannamei and when spliced-out, it is always together with SS2. All isoforms are expressed during embryogenesis as well as gametogenesis of both genders. The two shorter isoforms, PvanSxl-5 and PvanSxl-6, which result from the splicing of SS2 and SS3, were found mostly expressed in adult testis, but PvanSxl-6 was also expressed in oocytes during gametogenesis. During oogenesis, the second largest isoform, PvanSxl-2, which splices-out only SS1, and PvanSxl-4 that splices-out SS1 and SS2 were highly expressed. These two isoforms were also highly expressed during embryonic development. In situ hybridization allowed pinpointing more specifically the cells where the PvanSxl transcripts were expressed. During embryogenesis, hybridization was observed from the one-cell stage embryo to late gastrula. In the female gonad in previtellogenesis, hybridization occurred in the nucleus of oocytes, whereas in secondary vitellogenesis the transcript also hybridized cytoplasmic granules and cortical crypts. Finally, in situ hybridization corroborated the expression of PvanSxl also in the male gonad during spermatogenesis, mostly occurring in the cytoplasm from spermatogonia and spermatocytes.