Penaeid shrimp genome provides insights into benthic adaptation and frequent molting

Improved genomic resources for the black tiger prawn (Penaeus monodon)

World production of farmed crustaceans was 7.8 million tons in 2016. While only making up approximately 10% of world aquaculture production, crustaceans are generally high-value species and can earn significant export income for producing countries. Viet Nam is a major seafood producing country earning USD 7.3 billion in 2016 in export income with shrimp as a major commodity. However, there is a general lack of genomic resources available for shrimp species, which is challenging to obtain due to the need to deal with large repetitive genomes, which characterize many decapod crustaceans. The first tiger prawn (P. monodon) genome assembly was assembled in 2016 using the standard Illumina PCR-based pair-end reads and a computationally-efficient but relatively suboptimal assembler, SOAPdenovo v2. As a result, the current P. monodon draft genome is highly fragmented (> 2 million scaffolds with N₅₀ length of <1000 bp), exhibiting only moderate genome completeness (< 35% BUSCO complete single-copy genes). We sought to improve upon the recently published P. monodon genome assembly and completeness by generating Illumina PCR-free pair-end sequencing reads to eliminate genomic gaps associated with PCR-bias and performing de novo assembly using the updated MaSurCA de novo assembler. Furthermore, we scaffolded the assembly with low coverage Nanopore long reads and several recently published deep Illumina transcriptome paired-end sequencing data, producing a final genome assembly of 1.6 Gbp (1,211,364 scaffolds; N₅₀ length of 1982 bp) with an Arthropod BUSCO completeness of 96.8%. Compared to the previously published P. monodon genome assembly from China (NCBI Accession Code: NIUS01), this represents an almost 20% increase in the overall BUSCO genome completeness that now consists of more than 90% of Arthropod BUSCO single-copy genes. The revised P. monodon genome assembly (NCBI Accession Code: VIGR01) will be a valuable resource to support ongoing functional genomics and molecular-based breeding studies in Vietnam.

CRISPR/Cas9-mediated deletion of one carotenoid isomerooxygenase gene (EcNinaB-X1) from Exopalaemon carinicauda

During the immune defense reaction of invertebrate, a plenty of reactive oxygen species (ROS) could be induced to product. Though ROS can kill foreign invaders, the accumulation of these reactive molecules in animals will cause serious cell damage. Carotenoids could function as scavengers of oxygen radicals. In this research, cDNA and genomic DNA of one carotenoid isomerooxygenase gene (named EcNinaB-X1) were cloned from Exopalaemon carinicauda. EcNinaB-X1 gene was composed of 12 exons and 11 introns. EcNinaB-X1 knock-out (KO) prawns were produced via CRISPR/Cas9 technology and the change of their phenotypes were analyzed. Of the 400 injected one-cell stage embryos with cas9 mRNA and one sgRNA targeting the first exon of EcNinaB-X1 gene, 26 EcNinaB-X1-KO prawns were generated and the mutant rate reached 6.5% after embryo injection. The EcNinaB-X1-KO prawns had significant lower mortality than those in wild-type group when the prawns were challenged with Vibrio parahaemolyticus or Aeromonas hydrophila. In conclusion, we first demonstrate the function of the carotenoid isomerooxygenase gene in immune defense of E. carinicauda by performing directed, heritable gene mutagenesis.

Fine Mapping of the High-pH Tolerance and Growth Trait-Related Quantitative Trait Loci (QTLs) and Identification of the Candidate Genes in Pacific White Shrimp (Litopenaeus vannamei)

High-pH tolerance and growth are important traits for the shrimp culture industry in areas with saline-alkali water. In the present study, an F1 full-sib family of Pacific white shrimp (Litopenaeus vannamei) was generated with a new "semidirectional cross" method, and double-digest restriction site-associated DNA sequencing (ddRAD-Seq) technology was applied to genotype the 2 parents and 148 progenies. A total of 3567 high-quality markers were constructed for the genetic linkage map, and the total map length was 4161.555 centimorgans (cM), showing 48 linkage groups (LGs) with an average interlocus length of 1.167 cM. With a constrained logarithm of odds (LOD) score ≥ 2.50, 12 high-pH tolerance and 2 growth (body weight) QTLs were located. L. vannamei genomic scaffolds were used to assist with the detection of 21 stress- and 5 growth-related scaffold genes. According to the high-pH transcriptome data of our previous study, 6 candidate high-pH response genes were discovered, and 5 of these 6 genes were consistently expressed with the high-pH transcriptome data, validating the locations of the high-pH tolerance trait-related QTLs in this study. This paper is the first report of fine-mapping high-pH tolerance and growth (body weight) trait QTLs in one L. vannamei genetic map. Our results will further benefit marker-assisted selection work and might be useful for promoting genomic research on the shrimp L. vannamei.

A novel ML domain-containing protein (SpMD2) functions as a potential LPS receptor involved in anti-Vibrio immune response

The myeloid differentiation protein 2 (MD2)-related lipid-recognition (ML) proteins display diverse biological functions in host immunity and lipid metabolism by interacting with different lipids. Human MD2, an indispensable accessory protein in TLR4 signaling pathway, specifically recognizes lipopolysaccharides (LPS), thereby leading to the activation of TLR4 signaling pathway to produce many effectors that participate in inflammatory and immuneresponses against Gram-negative bacteria. Toll and immune deficiency (IMD) pathways are first characterized in Drosophila and are reportedly present in crustaceans, but the recognition and activation mechanism of these signaling pathways in crustaceans remains unclear. In the present study, a novel ML protein was characterized in mud crab (Scylla paramamosain) and designated as SpMD2. The complete SpMD2 cDNA sequence is 1114 bp long with a 465 bp open reading frame; it encodes a protein that contains 154 amino acids (aa). In the deduced protein, a signal peptide (1-21 aa residues) and a ML domain (43-151 aa residues) were predicted. SpMD2 shared a similar three-dimensional structure and a close evolutionary relationship with human MD2. SpMD2 was highly expressed in gills, hemocytes, intestine, and hepatopancreas and was upregulated in gills and hemocytes after challenges with bacteria, thereby suggesting its involvement in antibacterial defense. Western blot assay showed that SpMD2 possesses strong binding activities to different bacteria and two fungi. ELISA demonstrated that SpMD2 exhibits binding abilities to LPS, lipid A, peptidoglycan (PGN), and lipoteichoic acid (LTA). Its binding ability to LPS and lipid A were stronger than to PGN or LTA, implying that SpMD2 was an important LPS-binding protein in mud crab. Bacterial clearance assay revealed that the pre-incubation of Vibrio parahemolyticus with SpMD2 facilitates bacterial clearance in vivo and that knockdown of SpMD2 dramatically suppresses the bacterial clearance and decreases the expression of several antimicrobial peptides (AMPs). Furthermore, SpMD2 overexpression could enhance the promoter activity of SpALF2. These results revealed that SpMD2 affects bacterial clearance by regulating AMPs. Thus, by binding to LPS and by regulating AMPs, SpMD2 may function as a potential receptor, which is involved in the recognition and activation of a certain immune signaling pathway against Gram-negative bacteria. This study provides new insights into the diverse functions of ML proteins and into the antibacterial mechanisms of crustaceans.

Doing More with Less: A Comparison of 16S Hypervariable Regions in Search of Defining the Shrimp Microbiota

The shrimp has become the most valuable traded marine product in the world, and its microbiota plays an essential role in its development and overall health status. Massive high-throughput sequencing techniques using several hypervariable regions of the 16S rRNA gene are broadly applied in shrimp microbiota studies. However, it is essential to consider that the use of different hypervariable regions can influence the obtained data and the interpretation of the results. The present study compares the shrimp microbiota structure and composition obtained by three types of amplicons: one spanning both the V3 and V4 hypervariable regions (V3V4), one for the V3 region only (V3), and one for the V4 region only (V4) using the same experimental and bioinformatics protocols. Twenty-four samples from hepatopancreas and intestine were sequenced and evaluated using the GreenGenes and silva reference databases for clustering and taxonomic classification. In general, the V3V4 regions resulted in higher richness and diversity, followed by V3 and V4. All three regions establish an apparent clustering effect that discriminates between the two analyzed organs and describe a higher richness for the intestine and a higher diversity for the hepatopancreas samples. Proteobacteria was the most abundant phyla overall, and Cyanobacteria was more common in the intestine, whereas Firmicutes and Actinobacteria were more prevalent in hepatopancreas samples. Also, the genus Vibrio was significantly abundant in the intestine, as well as Acinetobacter and Pseudomonas in the hepatopancreas suggesting these taxa as markers for their respective organs independently of the sequenced region. The use of a single hypervariable region such as V3 may be a low-cost alternative that enables an adequate description of the shrimp microbiota, allowing for the development of strategies to continually monitor the microbial communities and detect changes that could indicate susceptibility to pathogens under real aquaculture conditions while the use of the full V3V4 regions can contribute to a more in-depth characterization of the microbial composition.

The Complete Genome of an Endogenous Nimavirus (Nimav-1_LVa) From the Pacific Whiteleg Shrimp Penaeus (Litopenaeus) Vannamei

White spot syndrome virus (WSSV), the lone virus of the genus Whispovirus under the family Nimaviridae, is one of the most devastating viruses affecting the shrimp farming industry. Knowledge about this virus, in particular, its evolution history, has been limited, partly due to its large genome and the lack of other closely related free-living viruses for comparative studies. In this study, we reconstructed a full-length endogenous nimavirus consensus genome, Nimav-1_LVa (279,905 bp), in the genome sequence of Penaeus (Litopenaeus) vannamei breed Kehai No. 1 (ASM378908v1). This endogenous virus seemed to insert exclusively into the telomeric pentanucleotide microsatellite (TAACC/GGTTA)n. It encoded 117 putative genes, with some containing introns, such as g012 (inhibitor of apoptosis, IAP), g046 (crustacean hyperglycemic hormone, CHH), g155 (innexin), g158 (Bax inhibitor 1 like). More than a dozen Nimav-1_LVa genes are involved in the pathogen-host interactions. We hypothesized that g046, g155, g158, and g227 (semaphorin 1A like) were recruited host genes for their roles in immune regulation. Sequence analysis indicated that a total of 43 WSSV genes belonged to the ancestral/core nimavirus gene set, including four genes reported in this study: wsv112 (dUTPase), wsv206, wsv226, and wsv308 (nucleocapsid protein). The availability of the Nimav-1_LVa sequence would help understand the genetic diversity, epidemiology, evolution, and virulence of WSSV.

A crustacean annotated transcriptome (CAT) database

BACKGROUND

Decapods are an order of crustaceans which includes shrimps, crabs, lobsters and crayfish. They occur worldwide and are of great scientific interest as well as being of ecological and economic importance in fisheries and aquaculture. However, our knowledge of their biology mainly comes from the group which is most closely related to crustaceans - insects. Here we produce a de novo transcriptome database, crustacean annotated transcriptome (CAT) database, spanning multiple tissues and the life stages of seven crustaceans.

DESCRIPTION

A total of 71 transcriptome assemblies from six decapod species and a stomatopod species, including the coral shrimp Stenopus hispidus, the cherry shrimp Neocaridina davidi, the redclaw crayfish Cherax quadricarinatus, the spiny lobster Panulirus ornatus, the red king crab Paralithodes camtschaticus, the coconut crab Birgus latro, and the zebra mantis shrimp Lysiosquillina maculata, were generated. Differential gene expression analyses within species were generated as a reference and included in a graphical user interface database at http://cat.sls.cuhk.edu.hk/. Users can carry out gene name searches and also access gene sequences based on a sequence query using the BLAST search function.

CONCLUSIONS

The data generated and deposited in this database offers a valuable resource for the further study of these crustaceans, as well as being of use in aquaculture development.

Insights on Molecular Mechanisms of Ovarian Development in Decapod Crustacea: Focus on Vitellogenesis-Stimulating Factors and Pathways

Vitellogenesis in crustaceans is an energy-consuming process. Though the underlying mechanisms of ovarian maturation in decapod Crustacea are still unclear, evidence indicates the process to be regulated by antagonistically-acting inhibitory and stimulating factors specifically originating from X-organ/sinus gland (XO/SG) complex. Among the reported neuromediators, neuropeptides belonging to the crustacean hyperglycemic hormone (CHH)-family have been studied extensively. The structure and dynamics of inhibitory action of vitellogenesis-inhibiting hormone (VIH) on vitellogenesis have been demonstrated in several species. Similarly, the stimulatory effects of other neuropeptides of the CHH-family on crustacean vitellogenesis have also been validated. Advancement in transcriptomic sequencing and comparative genome analysis has led to the discovery of a large number of neuromediators, peptides, and putative peptide receptors having pleiotropic and novel functions in decapod reproduction. Furthermore, differing research strategies have indicated that neurotransmitters and steroid hormones play an integrative role by stimulating neuropeptide secretion, thus demonstrating the complex intertwining of regulatory factors in reproduction. However, the molecular mechanisms by which the combinatorial effect of eyestalk hormones, neuromediators and other factors coordinate to regulate ovarian maturation remain elusive. These multifunctional substances are speculated to control ovarian maturation possibly via the autocrine/paracrine pathway by acting directly on the gonads or by indirectly exerting their stimulatory effects by triggering the release of a putative gonad stimulating factor from the thoracic ganglion. Acting through receptors, they possibly affect levels of cyclic nucleotides (cAMP and cGMP) and Ca²⁺ in target tissues leading to the regulation of vitellogenesis. The "stimulatory paradox" effect of eyestalk ablation on ovarian maturation continues to be exploited in commercial aquaculture operations, and is outweighed by the detrimental physiological effects of this procedure. In this regard, the development of efficient alternatives to eyestalk ablation based on scientific knowledge is a necessity. In this article, we focus principally on the signaling pathways of positive neuromediators and other factors regulating crustacean reproduction, providing an overview of their proposed receptor-mediated stimulatory mechanisms, intracellular signaling, and probable interaction with other hormonal signals. Finally, we provide insight into future research directions on crustacean reproduction as well as potential applications of such research to aquaculture technology development.

Investigation of Gene Sequence Divergence, Expression Dynamics, and Endocrine Regulation of the Vitellogenin Gene Family in the Whiteleg Shrimp Litopenaeus vannamei

In this report, we studied the vitellogenin gene family in the whiteleg shrimp Litopenaeus vannamei by transcriptomics, bioinformatics, and molecular biology methods. At least three moderately homologous vitellogenin (Vg) genes (i.e. LvVg1, LvVg2, and LvVg3) were identified in the genome. The deduced LvVg proteins consisted of a vitellogenin_N domain, a DUF1943 domain, and a VWD domain typical of most vitellogenins from oviparous animals. LvVg1 was the most abundant Vg expressed in the hepatopancreas and ovary of maturing females. Furthermore, multiple isoforms of LvVg1 were evolved presumably due to the need for rapid Vg production during the rapid phase of vitellogenesis. LvVg transcripts were detected in different larval stages, juveniles, and subadults. During the non-reproductive cycle, LvVg expression in the hepatopancreas peaked at the intermolt stages. During the female vitellogenesis cycle, a two-phase expression pattern of LvVg1 gene was observed in the hepatopancreas and ovary. Moreover, the eyestalk optic nerve, brain, and thoracic ganglion consisted of factors that differentially regulated the expression of the three Vg genes. In addition to their reproduction-related roles, Vg may also be involved in growth and molt-related processes. Phylogenetic analysis revealed the early expansion and separation of these Vg genes, and it is most likely correlated with the expansion of Vg's function. In conclusion, the evolution of multiple LvVg1 isoforms and the acquisition of different Vg genes (i.e. LvVg2 and LvVg3) may occur universally in most decapods. Full information on the total number of Vg genes and precise knowledge on the expression pattern and endocrine regulation of each Vg during all life cycle stages are crucial for us to understand the roles of this emerging gene family in the control of shrimp reproduction and other non-reproductive processes.

Sex-Biased CHHs and Their Putative Receptor Regulate the Expression of IAG Gene in the Shrimp Litopenaeus vannamei

The “eyestalk-androgenic gland (AG)-testis” endocrine axis is involved in male sexual differentiation of crustaceans. The insulin-like androgenic gland hormone (IAG), secreted from the AG, plays a central role in this axis, however key factors upstream the IAG are still poorly understood. Here, two crustacean hyperglycemic hormone (CHH) genes (LvCHH1 and LvCHH2) and their putative receptor guanylate cyclase (LvGC) were identified in Litopenaeus vannamei. LvCHH1 and LvCHH2 belonged to CHH subfamily I members and LvGC was a membrane-bound guanylate cyclase. They were all differentially expressed in eyestalks and gonads of males and females. RNA interference (RNAi) of either LvCHH1 or LvCHH2 increased LvIAG expression, while injection of their recombinant protein decreased LvIAG expression, indicating that LvCHH1 and LvCHH2 are inhibitory factors of LvIAG expression. Yeast two-hybrid assay showed that both LvCHH1 and LvCHH2 interacted with LvGC and their RNAi and recombinant protein injection exerted opposite regulatory effects on the transcriptional expression of LvGC. Meanwhile, knockdown of LvGC increased LvIAG expression. These results suggest that LvGC is the receptor of LvCHH1 and LvCHH2 and they are all involved in male sexual development by regulating LvIAG expression. The present study unveils missing upstream elements in the “eyestalk-AG-testis” endocrine axis in crustacean.

Maturation of an Antimicrobial Peptide Inhibits Aeromonas hydrophila Infection in Crayfish

Rapid synthesis and release of active antimicrobial peptides (AMPs) is an important strategy in innate immune. Processing of the precursor into the active form is a common posttranslational modification of AMPs in mammals. However, in invertebrates, the mechanism of AMP maturation is largely unknown. In the current study, to our knowledge, a novel potential AMP, designated as PcnAMP, was identified because of its significant induction by bacterial infection in the red swamp crayfish (Procambarus clarkii). PcnAMP was cleaved into a short fragment postinfection. Using the purified native peptide, this cleavage was found to be mediated by trypsin after synthesis. Proteolysis produced an N-terminal peptide that exerted the antibacterial function. Although the N-terminal peptide did not show significant similarity to any other sequences, it was predicted to have an overall helical structure and high amphipathicity, both of which are typical features of many AMPs. The N-terminal active peptide exhibited a wide spectrum of antimicrobial activity. Atomic force microscope imaging and flow cytometry analysis showed that treatment with the active form of PcnAMP led to the collapse of the bacterial cell wall and permeabilization of the bacterial cell membrane. Thus, this study provided a new candidate for therapeutic agent development, and revealed new insights into the maturation of AMPs in invertebrates.

Ecdysteroid-mimicking compounds act as both agonists and antagonists to the crustacean ecdysone receptor

To characterize the potential endocrine-disrupting chemicals (EDCs) in the environment that interact with the crustacean ecdysone receptor (EcR), we established a method involving in silico modeling/molecular docking and in vitro reporter gene assay. Cherry shrimp (Neocaridina davidi) EcR (NdEcR) and retinoid X receptor (NdRxR) were identified and cloned for use in this method. A theoretical 3D model of NdEcR ligand-binding domain (LBD) was built in silico based on sequence homology with the established X-ray structure of insect EcR. The interaction of the NdEcR LBD with ecdysteroids, diacylhydrazine (DAH) pesticides, and other potential EDCs was evaluated using molecular docking programs. The results revealed that the ligand-binding pocket in the NdEcR LBD was flexible and adaptive for accommodating ligands of different shapes. The agonistic and antagonistic activities of the candidate compounds were further assessed by in vitro reporter gene assay using human cell lines transiently transfected with NdEcR and NdRxR expression plasmids and a reporter plasmid containing synthesized ecdysone response element. The assay was validated by the dose-dependent responses of EcR-mediated gene transcription after treating the transfected cell lines with ecdysteroids, 20-hydroxyecdysone, and ponasterone A. Examination of the candidate compounds using the reporter gene assay revealed restricted functional specificity to ecdysteroids and DAHs. Three of the tested DAH pesticides originally targeting the insect EcR were found to be weak agonists and strong antagonists of NdEcR. These results suggest that DAHs are potential EDCs for crustaceans that disrupt their ecdysteroid signals by functioning as EcR agonists or antagonists.

Transcriptome-wide analysis of immune responses in Eriocheir sinensis hemocytes after challenge with different microbial derivatives

Knowledge about how Eriocheir sinensis interacts with microorganisms in its ambient environment is still lacking. Using RNA-Seq, we determined the most conserved genes and pathways compared with other animals and detected highly-induced immune genes in E. sinensis hemocytes post-in vivo challenge with different microbial derivatives. In total, 33.2 million high-quality reads were generated and assembled into 177,679 contigs. Completeness assessment and functional annotation were performed. Lipopolysaccharide, peptidoglycan, and β-1, 3-glucan stimulation induced 373, 173, and 108 differentially expressed (DE) transcripts, respectively. GO terms such as 'G-protein-coupled receptor binding', 'negative regulation of mitogen-activated protein kinase activity', and 'positive regulation of blood circulation' were enriched in the DE transcripts. Quantitative real-time PCR validated the data for selected genes. Our data contribute to understanding the immune defense mechanism in E. sinensis and the development of the innate immune system, thereby providing insights into disease control and prevention in aquaculture.

Chromosome-level reference genome of X12, a highly virulent race of the soybean cyst nematode Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean that is spreading across major soybean production regions worldwide. Increased SCN virulence has recently been observed in both the United States and China. However, no study has reported a genome assembly for H. glycines at the chromosome scale. Herein, the first chromosome-level reference genome of X12, an unusual SCN race with high infection ability, is presented. Using whole-genome shotgun (WGS) sequencing, Pacific Biosciences (PacBio) sequencing, Illumina paired-end sequencing, 10X Genomics linked reads and high-throughput chromatin conformation capture (Hi-C) genome scaffolding techniques, a 141.01-megabase (Mb) assembled genome was obtained with scaffold and contig N50 sizes of 16.27 Mb and 330.54 kilobases (kb), respectively. The assembly showed high integrity and quality, with over 90% of Illumina reads mapped to the genome. The assembly quality was evaluated using Core Eukaryotic Genes Mapping Approach and Benchmarking Universal Single-Copy Orthologs. A total of 11,882 genes were predicted using de novo, homolog and RNAseq data generated from eggs, second-stage juveniles (J2), third-stage juveniles (J3) and fourth-stage juveniles (J4) of X12, and 79.0% of homologous sequences were annotated in the genome. These high-quality X12 genome data will provide valuable resources for research in a broad range of areas, including fundamental nematode biology, SCN-plant interactions and co-evolution, and also contribute to the development of technology for overall SCN management.

Biochemical Characterization of a Novel α/β-Hydrolase/FSH from the White Shrimp Litopenaeus vannamei

(1) Background: Lipases and esterases are important enzymes that share the α/β hydrolase fold. The activity and cellular localization are important characteristics to understand the role of such enzymes in an organism. (2) Methods: Bioinformatic and biochemical tools were used to describe a new α/β hydrolase from a Litopenaeus vannamei transcriptome (LvFHS for Family Serine Hydrolase). (3) Results: The enzyme was obtained by heterologous overexpression in Escherichia coli and showed hydrolytic activity towards short-chain lipid substrates and high affinity to long-chain lipid substrates. Anti-LvFHS antibodies were produced in rabbit that immunodetected the LvFSH enzyme in several shrimp tissues. (4) Conclusions: The protein obtained and analyzed was an α/β hydrolase with esterase and lipase-type activity towards long-chain substrates up to 12 carbons; its immunodetection in shrimp tissues suggests that it has an intracellular localization, and predicted roles in energy mobilization and signal transduction.

What can transcriptomics reveal about the phylogenetic/structural conservation, tissue localization, and possible functions of CNMamide peptides in decapod crustaceans?

Over the past several years, in silico analyses of arthropod genomes/transcriptomes have led to the identification of several previously unknown peptide families. The CNMamides are one such peptide group, having been discovered via computational analyses of the fruit fly, Drosophila melanogaster, genome; both a CNMamide precursor and receptor were identified. Recently, a CNMamide family member, VMCHFKICNLamide (disulfide bridging between the cysteine residues), was predicted via in silico mining of a crayfish, Procambarus clarkii, transcriptome, suggesting the presence of this peptide group in members of the Decapoda. Here, using publically accessible transcriptomic data, the phylogenetic/structural conservation, tissue localization, and possible functions of the CNMamide family in decapods were explored. Evidence for CNMamide precursors was found for members of each decapod infraorder for which significant sequence data are available, suggesting a ubiquitous conservation of the CNMamide family in the Decapoda. For the Penaeoidea, Caridea, Astacidea and Achelata, the isoform of CNMamide originally identified from P. clarkii appears to be ubiquitously conserved; in members of the Brachyura, VMCHFKICNMamide (disulfide bridging between the cysteine residues) is the native isoform. Interestingly, the decapod CNMamide gene appears to also have a splice variant in which the carboxy-terminal portion of the preprohormone containing the CNMamide peptide is replaced by one containing a different disulfide bridged peptide that is structurally unrelated to it; this second peptide shows considerable conservation within, but variation among, decapod infraorders. A highly conserved putative CNMamide receptor was identified from members of the Penaeoidea, Astacidea and Brachyura. Phylogenetic analyses support the annotation of the decapod receptor as a true member of the CNMamide receptor family. The presence of precursor and receptor transcripts in both nervous system- and reproductive tissue-specific transcriptomes suggests CNMamides serve as modulators of decapod neural and reproductive control systems.

Temporal changes in transcriptome profile provide insights of White Spot Syndrome Virus infection in Litopenaeus vannamei

Shrimp aquaculture is severely affected by WSSV. Despite an increasing effort to understand host/virus interaction by characterizing changes in gene expression (GE) following WSSV infection, the majority of published studies have focussed on a single time-point, providing limited insight on the development of host-pathogen interaction over the infection cycle. Using RNA-seq, we contrasted GE in gills of Litopenaeus vannamei at 1.5, 18 and 56 hours-post-infection (hpi), between WSSV-challenged and control shrimps. Time course analysis revealed 5097 differentially expressed genes: 63 DEGs were viral genes and their expression in WSSV group either peaked at 18 hpi (and decreased at 56 hpi) or increased linearly up to 56 hpi, suggesting a different role played by these genes during the course of infection. The remaining DEGs showed that WSSV altered the expression of metabolic, immune, apoptotic and cytoskeletal genes and was able to inhibit NF-κB and JAK/STAT pathways. Interestingly, GE changes were not consistent through the course of infection but were dynamic with time, suggesting the complexity of host-pathogen interaction. These data offer novel insights into the cellular functions that are affected during the course of infection and ultimately provide a valuable resource towards our understanding of the host-pathogen dynamics and its variation with time.

Genome-Wide Analysis of Alternative Splicing Provides Insights Into Stress Response of the Pacific White Shrimp Litopenaeus vanname

Alternative splicing (AS) can enhance transcript diversity dramatically and play an important role in stress adaptation. Limited researches of AS have been reported in the Pacific white shrimp (Litopenaeus vannamei), which is an important aquaculture species in the world. Here, we performed a genome-wide identification of AS events in L. vannamei based on eight transcriptomes. We identified 38,781 AS events in the shrimp genome, and some of them were validated by polymerase chain reaction experiments. These AS events correspond to 9,209 genes, accounting for 36% of protein-coding genes in the shrimp genome. The number of AS events increased after virus or bacteria infection and low salinity stress. Type 1 AS genes (AS was initially activated) were mainly enriched in substance and energy metabolism, such as carbon metabolism and amino metabolism. However, type 2 AS genes (AS events changed) displayed specific enrichment under different stress challenges. Specifically, type 2 AS genes under biotic stresses were mainly enriched in the pathogenic pathway and immune network, and the AS genes under low salinity stress were significantly enriched for betalain biosynthesis. In summary, our study indicates that AS events are complex in shrimp and may be related to stress adaptation. These results will provide valuable resource for functional genomic studies on crustaceans.

Molecular characterization of a novel white spot syndrome virus response protein (dubbed LvWRP) from Litopenaeus vannamei

White spot syndrome, which is caused by white spot syndrome virus (WSSV), is a highly contagious disease of penaeid shrimp. However, there is currently incomplete understanding of the infection mechanism and pathogenesis of WSSV. In this study, a novel gene of a previously uncharacterized WSSV response protein (LvWRP) in Litopenaeus vannamei was identified and characterized. The LvWRP gene has an open reading frame (ORF) of 879 bp encoding a putative protein of 292 amino acids. Sequence analysis revealed that LvWRP shared 24.9% identity with an uncharacterized protein of Penaeus monodon nudivirus. Real-time qPCR analysis showed that LvWRP was ubiquitously expressed in shrimp tissues, with transcript levels induced in hemocytes upon immune challenge with Vibrio parahaemolyticus, Streptoccocus iniae, lipopolysaccharide (LPS), and WSSV. In addition, RNA interference-mediated knockdown of LvWRP followed by WSSV challenge revealed significant decrease in the transcript levels of WSSV IE1 and VP28 genes coupled with a reduction in WSSV copies in shrimp hemocytes. Moreover, depletion of LvWRP followed by WSSV challenge significantly increased the transcript levels of Vago4 and Vago5 as well as increased the phosphorylation of STAT, while hemocytes apoptosis in terms of caspase 3/7 activity was decreased. These results suggest that LvWRP is important for WSSV replication in shrimp, and therefore one of the vital host factors in WSSV infection.

Proteomic characterization of the hepatopancreas in the Pacific white shrimp Litopenaeus vannamei under cold stress: Revealing the organism homeostasis mechanism

To understand the homeostasis mechanism of crustacean hepatopancreas to cold stress, iTRAQ proteomics based on the genome database of Litopenaeus vannamei (L. vannamei) was applied to investigate proteins changes and variety of the hepatopancreas during cold stress stage in this study. A total of 4062 distinct proteins were identified, 137 differentially expressed proteins (DEPs) including 62 differentially up-regulated proteins (DUPs) and 75 differentially down-regulated proteins (DDPs) were identified in G1 (18 °C) compared with CK (28 °C), 359 DEPs including 131 DUPs and 228 DDPs were identified in G2 (13 °C for 24 h) compared with CK. Based on bioinformatics analysis, the cold tolerance of L. vannamei might be related to energy metabolism such as amino acid, carbohydrate, lipid, and oxidative phosphorylation. Moreover, shrimp immunity was declined during cold stress stage. However, L. vannamei could cope with cold stress by enhancing the production of ATP and UFA. Notably, arginine kinase, heat shock proteins, and histones may act as positive regulators in L. vannamei under cold stress. Ten randomly selected proteins were used for validation using qRT-PCR and the expressions on the transcription level for most of the genes were similar to the results of iTRAQ. These results indicated that L. vannamei can maintain the organism homeostasis by a series of orderly regulatory process during cold stress. Furthermore, the results can provide guidance for shrimp farming.

Production of WW males lacking the masculine Z chromosome and mining the Macrobrachium rosenbergii genome for sex-chromosomes

The cultivation of monosex populations is common in animal husbandry. However, preselecting the desired gender remains a major biotechnological and ethical challenge. To achieve an efficient biotechnology for all-female aquaculture in the economically important prawn (Macrobrachium rosenbergii), we achieved – for the first time – WW males using androgenic gland cells transplantation which caused full sex-reversal of WW females to functional males. Crossing the WW males with WW females yielded all-female progeny lacking the Z chromosome. We now have the ability to manipulate – by non-genomic means – all possible genotype combinations (ZZ, WZ and WW) to retain either male or female phenotypes and hence to produce monosex populations of either gender. This calls for a study of the genomic basis underlying this striking sexual plasticity, questioning the content of the W and Z chromosomes. Here, we report on the sequencing of a high-quality genome exhibiting distinguishable paternal and maternal sequences. This assembly covers ~ 87.5% of the genome and yielded a remarkable N50 value of ~ 20 × 10⁶ bp. Genomic sex markers were used to initiate the identification and validation of parts of the W and Z chromosomes for the first time in arthropods.

Isolation and characterisation of 14 novel microsatellite markers through Next Generation Sequencing for the commercial Atlantic seabob shrimp Xiphopenaeus kroyeri

Assessing population genetic structure is a crucial step to support fisheries and conservation management. DNA microsatellite molecular markers are a widely used tool in population genotyping. In the present study, we characterised and developed 14 novel polymorphic microsatellite markers for a decapod crustacean, the Atlantic seabob shrimp Xiphopenaeus kroyeri (Heller, 1862), through rapid and cost-effective Illumina shotgun sequencing and a Galaxy-based bioinformatic pipeline. We genotyped 60 individuals from 2 populations with the newly developed microsatellites, resulting in the detection of 3 to 29 alleles per locus. Four loci deviated from Hardy-Weinberg equilibrium. Cross-amplification in a cryptic congeneric species was successful for eight loci (57%). The microsatellite loci developed in this study will be highly relevant for genetic and evolutionary studies of X. kroyeri, and for the stock management of this commercially exploited species.

Identification of Single Nucleotide Polymorphisms Related to the Resistance Against Acute Hepatopancreatic Necrosis Disease in the Pacific White Shrimp Litopenaeus vannamei by Target Sequencing Approach

Acute hepatopancreatic necrosis disease (AHPND) is a major bacterial disease in Pacific white shrimp Litopenaeus vannamei farming, which is caused by Vibrio parahaemolyticus. AHPND has led to a significant reduction of shrimp output since its outbreak. Selective breeding of disease-resistant broodstock is regarded as a key strategy in solving the disease problem. Understanding the relationship between genetic variance and AHPND resistance is the basis for marker-assisted selection in shrimp. The purpose of this study was to identify single nucleotide polymorphisms (SNPs) associated with the resistance against AHPND in L. vannamei. In this work, two independent populations were used for V. parahaemolyticus challenge and the resistant or susceptible shrimp were evaluated according to the survival time after Vibrio infection. The above two populations were genotyped separately by a SNP panel designed based on the target sequencing platform using a pooling strategy. The SNP panel contained 508 amplicons from DNA fragments distributed evenly along the genome and some immune-related genes of L. vannamei. By analyzing the allele frequency in the resistant and susceptible groups, 30 SNPs were found to be significantly associated with the resistance of the shrimp against V. parahaemolyticus infection (false discovery rate corrected at P < 0.05). Three SNPs were further validated by individual genotyping in all samples of population 1. Our study illustrated that target sequencing and pooling sequencing were effective in identifying the markers associated with economic traits, and the SNPs identified in this study could be used as molecular markers for breeding disease-resistant shrimp.

The Two NF-κB Pathways Regulating Bacterial and WSSV Infection of Shrimp

The outbreak of diseases ordinarily results from the disruption of the balance and harmony between hosts and pathogens. Devoid of adaptive immunity, shrimp rely largely on the innate immune system to protect themselves from pathogenic infection. Two nuclear factor-κB (NF-κB) pathways, the Toll and immune deficiency (IMD) pathways, are generally regarded as the major regulators of the immune response in shrimp, which have been extensively studied over the years. Bacterial infection can be recognized by Toll and IMD pathways, which activate two NF-κB transcription factors, Dorsal and Relish, respectively, to eventually lead to boosting the expression of various antimicrobial peptides (AMPs). In response to white-spot-syndrome-virus (WSSV) infection, these two pathways appear to be subverted and hijacked to favor viral survival. In this review, the recent progress in elucidating microbial recognition, signal transduction, and effector regulation within both shrimp Toll and IMD pathways will be discussed. We will also highlight and discuss the similarities and differences between shrimps and their Drosophila or mammalian counterparts. Understanding the interplay between pathogens and shrimp NF-κB pathways may provide new opportunities for disease-prevention strategies in the future.

Molecular cloning of the Rab7 effector RILP (Rab-interacting lysosomal protein) in Litopenaeus vannamei and preliminary analysis of its role in white spot syndrome virus infection

To investigate the role of the Rab7 effector RILP (Rab-interacting lysosomal protein) in white spot syndrome virus (WSSV) infection, the full-length cDNA of RILP (LvRILP) was cloned in Litopenaeus vannamei, which consists of 1595 bp and encodes a polypeptide of 411 amino acids. Sequence analysis and multiple sequence alignment displayed that LvRILP contained a conserved RILP region from 277 amino acid to 325 amino acid. Both the LvRILP and Rab7 mRNA were most highly expressed in stomach and most lowly expressed in hemocyte, which were significantly up-regulated and exhibited similar kinetics post WSSV infection. The interaction of Rab7 with LvRILP was verified by both GST Pull-down and ELISA. Meanwhile, the results of Pull-down assays showed that the GST-tagged VP28 (GST-VP28), His-tagged Rab7 (His-Rab7) and His-RILP formed a tripartite complex. After silencing by specific LvRILP dsRNA, the LvRILP mRNA level exhibited a significant reduction, and the expression levels of three WSSV genes ie1, wsv477 and vp28 all exhibited decreases at 24, 36 and 48 h post WSSV infection. These results suggested that the Rab7 effector RILP was involved in WSSV infection.

Comparative physiology and aquaculture: Toward Omics-enabled improvement of aquatic animal health and sustainable production

Omics-technologies have revolutionized biomedical research over the past two decades, and are now poised to play a transformative role in aquaculture. This article serves as an introduction to a Virtual Special Issue of Comparative Biochemistry and Physiology - Part D: Genomics and Proteomics (CBPD), with the objective to showcase the state-of-the-science for Omics in aquaculture. In this editorial, we describe the role that Omics can play in aquaculture, and provide a synopsis for each of the Special Issue articles that use these technologies to improve aquaculture practices. Current genomic resources available for some aquaculture species are also described. The number of datasets is impressive for species such as Atlantic salmon and rainbow trout, totaling in the thousands (NCBI Gene Expression Omnibus and Sequence Read Archive). We present a conceptual framework that describes how Omics can be leveraged to understand complex responses of aquatic animals in culture for relevant physiological outcomes, such as fecundity, growth, and immunity. Lastly, knowledge gaps and new directions are identified to address current obstacles in aquaculture. Articles in this Special Issue on aquaculture in CBPD highlight the diversity and scope of Omics in aquaculture. As the technology becomes more cost-effective, it is anticipated that genomics, transcriptomics, proteomics, metabolomics and lipidomics will play increasingly important roles in stock diagnostics (e.g. genetics, health, performance). The timing is right, as global concerns are reaching critical levels over food availability/security and water restrictions for humankind.

Genome Scan for Genomic Regions and Genes Associated with Growth Trait in Pacific White Shrimp Litopeneaus vannamei

The Pacific white shrimp Litopeneaus vannmei (L. vannmei) is a predominant aquaculture shrimp species worldwide, and it is considered as the aquaculture species with the highest single output value. Advances in selective breeding have accelerated the development of L. vannmei aquaculture. Recently, the genome-wide association studies (GWAS) have been applied in aquaculture animals and markers associated with economic traits were identified. In this study, we focused on the growth trait of L. vannamei and performed GWAS to identify SNPs or genes associated with growth. Genomic regions in linkage group 7, 27, 33, and 38 were identified to be associated with body weight and body length of the shrimp. Further, candidate gene association analysis was performed in two independent populations and the result demonstrated that the SNPs in the genes protein kinase C delta type and ras-related protein Rap-2a were significantly associated with the growth trait of L. vannamei. This study showed that GWAS analysis is an efficient approach for screening trait-related markers or genes. The genomic regions and genes identified in this study are essential for further fine mapping of growth-related genes. The identified markers will provide useful information for marker-assisted selection in L. vannamei.

A Novel Candidate Gene Associated With Body Weight in the Pacific White Shrimp Litopenaeus vannamei

Improvements of growth traits are always the focus in selective breeding programs for the Pacific white shrimp Litopenaeus vannamei (L. vannamei). Identification of growth-related genes or markers can contribute to the application of modern breeding technologies, and thus accelerate the genetic improvement of growth traits. The aim of this study was to identify the genes and molecular markers associated with the growth traits of L. vannamei. A population of 200 individuals was genotyped using 2b-RAD techniques for genome-wide linkage disequilibrium (LD) analysis and genome-wide association study (GWAS). The results showed that the LD decayed fast in the studied population, which suggest that it is feasible to fine map the growth-related genes with GWAS in L. vannamei. One gene designated as LvSRC, encoding the class C scavenger receptor (SRC), was identified as a growth-related candidate gene by GWAS. Further targeted sequencing of the candidate gene in another population of 322 shrimps revealed that several non-synonymous mutations within LvSRC were significantly associated with the body weight (P < 0.01), and the most significant marker (SRC_24) located in the candidate gene could explain 13% of phenotypic variance. The current results provide not only molecular markers for genetic improvement in L. vannamei, but also new insights for understanding the growth regulation mechanism in penaeid shrimp.

Genome-Wide Identification and Expression Profiles of Myosin Genes in the Pacific White Shrimp, Litopenaeus vannamei

As the main structural protein of muscle fiber, myosin is essential for multiple cellular processes or functions, especially for muscle composition and development. Although the shrimp possess a well-developed muscular system, the knowledge about the myosin family in shrimp is far from understood. In this study, we performed comprehensive analysis on the myosin genes in the genome of the Pacific white shrimp, Litopenaeus vannamei. A total of 29 myosin genes were identified, which were classified into 14 subfamilies. Among them, Myo2 subfamily was significantly expanded in the penaeid shrimp genome. Most of the Myo2 subfamily genes were primarily expressed in abdominal muscle, which suggested that Myo2 subfamily genes might be responsible for the well-developed muscular system of the penaeid shrimp. In situ hybridization detection showed that the slow-type muscle myosin gene was mainly localized in pleopod muscle and superficial ventral muscle of the shrimp. This study provides valuable insights into the evolutionary and functional characterization of myosin genes in shrimps, which provides clues for us to understand the well-developed muscular system of shrimp.

Characterization of a Lymphoid Organ Specific Anti-lipopolysaccharide Factor From Shrimp Reveals Structure-Activity Relationship of the LPS-Binding Domain

Anti-lipopolysaccharide factor (ALF) is a kind of important antimicrobial peptides with broad-spectrum antimicrobial activities. The LPS-binding domain (LBD) contributes to the major antimicrobial activity of ALF. However, LBDs from different ALFs share low sequence similarity. The general character of LBDs needs to be elucidated to understand the molecular mechanism of their function and facilitate LBD-original drug design. Here we identified a lymphoid organ specifically expressed ALF, designated as FcALF8, from the Chinese shrimp Fenneropenaeus chinensis. The synthetic LBD peptide of FcALF8 (LBD8) showed strong antibacterial activities to the pathogenic Vibrio, such as Vibrio alginolyticus, Vibrio harveyi, and Photobacterium damselae with a MIC value of 0.5–1, 1–2, and 1–2 μM, respectively. FcALF8 knock-down using dsRNA led to significant increase of the viable bacteria in the lymphoid organ and hepatopancreas of shrimp upon V. harveyi infection. On the contrary, the proliferation of V. harveyi in the shrimp lymphoid organ and hepatopancreas significantly decreased after infected by LBD8 pre-incubated V. harveyi. Sequence alignments showed that the LBDs from 39 ALFs shared only two identical cysteine residues. However, 17 of the total 22 LBD residues showed high similarity when the amino acids were classified into hydrophobic and hydrophilic ones. A further activity analysis on modified LBD8 peptides showed that the antibacterial activity of LBD8 was lost after linearization and apparently weakened after changing the amino acid property at certain positions. The data indicated that the disulfide bond and amino acid property contributed to the conservation of the functional domain. To the best of our knowledge, this is the first identified ALFs specifically expressed in the lymphoid organ of shrimp with strong antibacterial activity. The present data will give creative instructions for the design of LBD-originated antimicrobial agents.

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.

High-Density Genetic Linkage Maps Provide Novel Insights Into ZW/ZZ Sex Determination System and Growth Performance in Mud Crab (Scylla paramamosain)

Mud crab, Scylla paramamosain is one of the most important crustacean species in global aquaculture. To determine the genetic basis of sex and growth-related traits in S. paramamosain, a high-density genetic linkage map with 16,701 single nucleotide polymorphisms (SNPs) was constructed using SLAF-seq and a full-sib family. The consensus map has 49 linkage groups, spanning 5,996.66 cM with an average marker-interval of 0.81 cM. A total of 516 SNP markers, including 8 female-specific SNPs segregated in two quantitative trait loci (QTLs) for phenotypic sex were located on LG32. The presence of female-specific SNP markers only on female linkage map, their segregation patterns and lower female: male recombination rate strongly suggest the conformation of a ZW/ZZ sex determination system in S. paramamosain. The QTLs of most (90%) growth-related traits were found within a small interval (25.18–33.74 cM) on LG46, highlighting the potential involvement of LG46 in growth. Four markers on LG46 were significantly associated with 10–16 growth-related traits. BW was only associated with marker 3846. Based on the annotation of transcriptome data, 11 and 2 candidate genes were identified within the QTL regions of sex and growth-related traits, respectively. The newly constructed high-density genetic linkage map with sex-specific SNPs, and the identified QTLs of sex- and growth-related traits serve as a valuable genetic resource and solid foundation for marker-assisted selection and genetic improvement of crustaceans.