Penaeid shrimp genome provides insights into benthic adaptation and frequent molting

Identification of putative ingestion-related olfactory receptor genes in the Chinese mitten crab (Eriocheir japonica sinensis)

BACKGROUND

Olfaction plays a central role in mating, spawning, obtaining food and escaping predators, which is essential for survival and reproduction of animals. The nature of the olfactory perception in crabs, which is a major group of crustaceans, has remained elusive.

OBJECTIVE

This project aims to explore the molecular mechanism of olfaction in crabs and further improve our understanding of olfactory perception in crustaceans.

METHODS

The olfactory receptors and ingestion-related gene expression in Eriocheir japonica sinensis were studied by transcriptomic techniques. The de novo assembly, annotation and functional evaluation were performed with bioinformatics tools.

RESULTS

A series of chemosensory receptors associated with olfaction were identified including 33 EsIRs, 24 EsIGluRs, 58 EsVIGluRs, 1 EsOR and 1 EsGC-D. We found IRs were key odorant receptors demonstrating a specific species evolutionary trend in crustaceans. Furthermore, we identified ORs in E. j. sinensis and Litopenaeus vannamei. The incomplete EsOR and LvOR1 structures implied that ORs exist in crustaceans, and may have been degenerated or even lost in the olfactory evolutionary process. In addition, comparative transcriptome analysises demonstrated two possible olfactory transduction pathways of E. j. sinensis: the cGMP-mediated olfactory pathway related to vegetable odor molecules and the cAMP-mediated olfactory pathway related to meat odor molecules. The above results were consistent with its omnivorous ingestion of E. j. sinensis.

CONCLUSIONS

Our study revealed the unique olfactory molecular mechanism of omnivorous crabs and provided valuable information for further functional research on the chemoreception mechanisms in crustaceans.

A chromosome-level genome of the mud crab (Scylla paramamosain estampador) provides insights into the evolution of chemical and light perception in this crustacean

Mud crabs, found throughout the Indo-Pacific region, are coastal species that are important fisheries resources in many tropical and subtropical Asian countries. Here, we present a chromosome-level genome assembly of a mud crab (Scylla paramamosain). The genome is 1.55 Gb (contig N50 191 kb) in length and encodes 17,821 proteins. The heterozygosity of the assembled genome was estimated to be 0.47%. Effective population size analysis suggested that an initial large population size of this species was maintained until 200 thousand years ago. The contraction of cuticle protein and opsin genes compared with Litopenaeus vannamei is assumed to be correlated with shell hardness and light perception ability, respectively. Furthermore, the analysis of three chemoreceptor gene families, the odorant receptor (OR), gustatory receptor (GR) and ionotropic receptor (IR) families, suggested that the mud crab has no OR genes and shows a contraction of GR genes and expansion of IR genes. The numbers of the three gene families were similar to those in three other decapods but different from those in two nondecapods and insects. In addition, IRs were more diversified in decapods than in nondecapod crustaceans, and most of the expanded IRs in the mud crab genome were clustered with the antennal IR clades. These findings suggested that IRs might exhibit more diverse functions in decapods than in nondecapods, which may compensate for the smaller number of GR genes. Decoding the S. paramamosain genome not only provides insight into the genetic changes underpinning ecological traits but also provides valuable information for improving the breeding and aquaculture of this species.

Sex Determination and Differentiation in Decapod and Cladoceran Crustaceans: An Overview of Endocrine Regulation

Mechanisms underlying sex determination and differentiation in animals are known to encompass a diverse array of molecular clues. Recent innovations in high-throughput sequencing and mass spectrometry technologies have been widely applied in non-model organisms without reference genomes. Crustaceans are no exception. They are particularly diverse among the Arthropoda and contain a wide variety of commercially important fishery species such as shrimps, lobsters and crabs (Order Decapoda), and keystone species of aquatic ecosystems such as water fleas (Order Branchiopoda). In terms of decapod sex determination and differentiation, previous approaches have attempted to elucidate their molecular components, to establish mono-sex breeding technology. Here, we overview reports describing the physiological functions of sex hormones regulating masculinization and feminization, and gene discovery by transcriptomics in decapod species. Moreover, this review summarizes the recent progresses of studies on the juvenile hormone-driven sex determination system of the branchiopod genus Daphnia, and then compares sex determination and endocrine systems between decapods and branchiopods. This review provides not only substantial insights for aquaculture research, but also the opportunity to re-organize the current and future trends of this field.

Simple sequence repeats drive genome plasticity and promote adaptive evolution in penaeid shrimp

Simple sequence repeats (SSRs) are rare (approximately 1%) in most genomes and are generally considered to have no function. However, penaeid shrimp genomes have a high proportion of SSRs (>23%), raising the question of whether these SSRs play important functional and evolutionary roles in these SSR-rich species. Here, we show that SSRs drive genome plasticity and adaptive evolution in two penaeid shrimp species, Fenneropenaeus chinensis and Litopenaeus vannamei. Assembly and comparison of genomes of these two shrimp species at the chromosome-level revealed that transposable elements serve as carriers for SSR expansion, which is still occurring. The remarkable genome plasticity identified herein might have been shaped by significant SSR expansions. SSRs were also found to regulate gene expression by multi-omics analyses, and be responsible for driving adaptive evolution, such as the variable osmoregulatory capacities of these shrimp under low-salinity stress. These data provide strong evidence that SSRs are an important driver of the adaptive evolution in penaeid shrimp.

Comparison of growth performance and biochemical components between parent and hybrid offspring in the oriental river prawn, Macrobrachium nipponense

Macrobrachium nipponense, as one of the large-yield farmed shrimp, is facing germplasm degradation. Genetic improvement through hybridization is one of the effective methods to solve this problem. In this study, using a three-line hybrid strategy, two-hybrid F1 populations were obtained using three local populations of M. nipponense as parents for crossbreeding. Five populations were then cultured for 3 months. Growth rate performance was measured by the hepatosomatic index, weight gain, body length growth rate and special growth rate. Biochemical components were also assessed. The results showed that the survival rate and growth performance of the hybrid progeny were better than those of the parents. The levels of triglycerides, total cholesterol, glycogen and lactic acid of the hybrid population were higher than those of the parents. This was consistent with variation in the activity of four digestive enzymes. Compared with the results of the fatty acid and amino acid analysis, it was found that the contents of polyunsaturated fatty acids such as eicosapentaenoic acid, docosapentaenoic acid and docosahexaenoic acid and eight essential amino acids in the hybrid populations were significantly higher than those of their parents, and the contents of flavor amino acids were higher. The expression level of molting genes related to the growth of the parent populations was lower than that of the hybrids. These results show that crossbreeding is effective for the genetic improvement of M. nipponense germplasm. Hybrids showed advantages in growth and nutrition and multigenerational breeding will be required to form a stable germplasm.

The Pontastacus leptodactylus (Astacidae) Repeatome Provides Insight Into Genome Evolution and Reveals Remarkable Diversity of Satellite DNA

Pontastacus leptodactylus is a native European crayfish species found in both freshwater and brackish environments. It has commercial importance for fisheries and aquaculture industries. Up till now, most studies concerning P. leptodactylus have focused onto gaining knowledge about its phylogeny and population genetics. However, little is known about the chromosomal evolution and genome organization of this species. Therefore, we performed clustering analysis of a low coverage genomic dataset to identify and characterize repetitive DNA in the P. leptodactylus genome. In addition, the karyogram of P. leptodactylus (2n = 180) is presented here for the first time consisting of 75 metacentric, 14 submetacentric, and a submetacentric/metacentric heteromorphic chromosome pair. We determined the genome size to be at ~18.7 gigabase pairs. Repetitive DNA represents about 54.85% of the genome. Satellite DNA repeats are the most abundant type of repetitive DNA, making up to ~28% of the total amount of repetitive elements, followed by the Ty3/Gypsy retroelements (~15%). Our study established a surprisingly high diversity of satellite repeats in P. leptodactylus. The genome of P. leptodactylus is by far the most satellite-rich genome discovered to date with 258 satellite families described. Of the five mapped satellite DNA families on chromosomes, PlSAT3-411 co-localizes with the AT-rich DAPI positive probable (peri)centromeric heterochromatin on all chromosomes, while PlSAT14-79 co-localizes with the AT-rich DAPI positive (peri)centromeric heterochromatin on one chromosome and is also located subterminally and intercalary on some chromosomes. PlSAT1-21 is located intercalary in the vicinity of the (peri)centromeric heterochromatin on some chromosomes, while PlSAT6-70 and PlSAT7-134 are located intercalary on some P. leptodactylus chromosomes. The FISH results reveal amplification of interstitial telomeric repeats (ITRs) in P. leptodactylus. The prevalence of repetitive elements, especially the satellite DNA repeats, may have provided a driving force for the evolution of the P. leptodactylus genome.

Characterization and Expression Analysis of Insulin Growth Factor Binding Proteins (IGFBPs) in Pacific White Shrimp Litopenaeus vannamei

The insulin signaling (IIS) pathway plays an important role in the metabolism, growth, development, reproduction, and longevity of an organism. As a key member of the IIS pathway, insulin-like growth factor binding proteins (IGFBPs) are widely distributed a family in invertebrates and vertebrates that are critical in various aspects of physiology. As an important mariculture species, the growth of Pacific white shrimp, Litopenaeus vannamei, is one of the most concerning characteristics in this area of study. In this study, we identified three IGFBP genes in the genome of L. vannamei and analyzed their gene structures, phylogenetics, and expression profiles. LvIGFBP1 was found to contain three domains (the insulin growth factor binding (IB) domain, the Kazal-type serine proteinase inhibitor (Kazal) domain, and the immunoglobulin C-2 (IGc2) domain), while LvIGFBP2 and LvIGFBP3 only contained a single IB domain. LvIGFBP1 exhibited high expression in most tissues and different developmental stages, while LvIGFBP2 and LvIGFBP3 were only slightly expressed in hemocytes. The RNA interference of LvIGFBP1 resulted in a significantly smaller increment of body weight than that of control groups. These results will improve our understanding of the conservative structure and function of IGFBPs and show potential applications for the growth of shrimp.

The Neurobiology of Ocean Change - insights from decapod crustaceans

The unprecedented rate of carbon dioxide accumulation in the atmosphere has led to increased warming, acidification and oxygen depletion in the world's oceans, with projected impacts also on ocean salinity. In this perspective article, we highlight potential impacts of these factors on neuronal responses in decapod crustaceans. Decapod crustaceans comprise more than 8,800 marine species which have colonized a wide range of habitats that are particularly affected by global ocean change, including estuarine, intertidal, and coastal areas. Many decapod species have large economic value and high ecological importance because of their global invasive potential and impact on local ecosystems. Global warming has already led to considerable changes in decapod species' behavior and habitat range. Relatively little is known about how the decapod nervous system, which is the ultimate driver of all behaviors, copes with environmental stressors. We use select examples to summarize current findings and evaluate the impact of current and expected environmental changes. While data indicate a surprising robustness against stressors like temperature and pH, we find that only a handful of species have been studied and long-term effects on neuronal activity remain mostly unknown. A further conclusion is that the combined effects of multiple stressors are understudied. We call for greater research efforts towards long-term effects on neuronal physiology and expansion of cross-species comparisons to address these issues.

A Novel Forkhead Box Protein P (FoxP) From Litopenaeus vannamei Plays a Positive Role in Immune Response

The forkhead box protein P (FoxP) family members have been known to be important for regulation of immune responses in vertebrates, but their roles in invertebrate immunity remain unclear. In this study, a novel FoxP gene (LvFoxP) was identified from Pacific white shrimp Litopenaeus vannamei and functionally studied in the context of immune response. Possessing a conserved FoxP coiled-coil domain and a forkhead domain, LvFoxP shared homology to vertebrate FoxP family members, in particular FoxP1. Expression of LvFoxP was detectable in all the examined tissues and could be up-regulated by immune challenge in gill and hemocytes. The LvFoxP protein was present in both the cytoplasm and nucleus of hemocytes and could be nuclear-translocated upon immune stimulation. Silencing of LvFoxP increased the susceptibility of shrimp to infections by Vibrio parahaemolyticus and white spot syndrome virus (WSSV) and down-regulated the expression of multiple components of NF-κB and JAK-STAT pathways and almost all the examined immune effector genes. Moreover, the phagocytic activity of hemocytes from LvFoxP-silenced shrimp against V. parahaemolyticus was decreased. These suggested that LvFoxP could play a positive role in immune response. The current study may provide novel insights into the immunity of invertebrates and the functional evolution of the FoxP family.

The LARK protein is involved in antiviral and antibacterial responses in shrimp by regulating humoral immunity

The LARK proteins containing a C2HC-type zinc finger motif and two RNA recognition motifs are conserved across vertebrates and invertebrates. Previous studies have suggested that invertebrate LARKs and their mammalian counterparts, the RBM4 proteins, regulate gene expression by affecting RNA stability and post-transcriptional processing, participating in multiple life processes. In the current study, the LARK gene from Pacific white shrimp Litopenaeus vannamei was identified and functionally explored in the context of immunity. The LARK protein was mainly present in the nucleus of its expression vector-transfected S2 cells, and the LARK mRNA was detectable in all the tested shrimp tissues. Expression of LARK in gill was up-regulated by immune stimulation with various pathogens. In vivo experiments demonstrated that LARK played positive roles in both antiviral and antibacterial responses and silencing of LARK could make shrimp more susceptible to infection with Vibrio parahaemolyticus and white spot syndrome virus (WSSV). Although silencing of LARK did not affect the phagocytic activity of hemocytes, it regulated expression of many components of the NF-κB and JAK-STAT pathways and a series of immune function proteins. These suggested that LARK could be mainly involved in regulation of humoral immunity. The current study could help reveal the roles of LARK/RBM4 in immunity and further explore the regulatory mechanisms of shrimp immunity.

Segregation distortion: high genetic load suggested by a Chinese shrimp family under high-intensity selection

Segregation distortion is a common phenomenon found in most genetic mapping studies and is an important resource to dissect the mechanism of action in gene loci that cause deviation. Marine animals possess high genetic diversity and genomic heterozygosity, they therefore are ideal model organisms to study segregation distortion induced by selection. In the present study, we constructed a full-sib family of Fenneropenaeus chinensis and exerted high-intensity selection on 10,000 incipient progenies. 2b-RAD method was employed in remaining 273 individuals to develop genome-wide SNPs for segregating analysis and 41,612 SNPs were developed. 50.77% of 32,229 high-quality representative markers deviated from the expected Mendelian ratio. Results showed that most of these distorted markers (91.57%) were influenced at zygotic level. Heterozygote excess (53.07%) and homozygous deletions (41.96%) may both play an important role, sum of which explained 95.03% of distortion after fertilization. However, further results identified highly probable linkage among deleterious alleles, which may account for a considerable portion of heterozygote excess rather than single locus with heterozygote advantage. Results of this study support a major role of deleterious alleles in genetic load, thus in favor of partial dominance hypothesis. It would also offer necessary recommendations for the formulation of breeding strategy in shrimps.

Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei)

Background

Ammonia is one of the most common toxicological environment factors affecting shrimp health. Although ammonia tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available.

Results

In this study, we constructed a high-density genetic map of the Pacific white shrimp (Litopenaeus vannamei) using specific length amplified fragment sequencing (SLAF-seq). The constructed genetic map contained 17,338 polymorphic markers spanning 44 linkage groups, with a total distance of 6360.12 centimorgans (cM) and an average distance of 0.37 cM. Using this genetic map, we identified a quantitative trait locus (QTL) that explained 7.41–8.46% of the phenotypic variance in L. vannamei survival time under acute ammonia stress. We then sequenced the transcriptomes of the most ammonia-tolerant and the most ammonia-sensitive individuals from each of four genetically distinct L. vannamei families. We found that 7546 genes were differentially expressed between the ammonia-tolerant and ammonia-sensitive individuals. Using QTL analysis and the transcriptomes, we identified one candidate gene (annotated as an ATP synthase g subunit) associated with ammonia tolerance.

Conclusions

In this study, we constructed a high-density genetic map of L. vannamei and identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07254-x.

A double chitin catalytic domain-containing chitinase targeted by c-Jun is involved in immune responses in shrimp

Chitinases are a group of chitin-degrading enzymes widely distributed in organisms. Chitinases containing two chitin catalytic domains have been widely found in arthropods but their functions remain unclear. In this study, a member of these chitinases from Litopenaeus vannamei (dChi) was identified and functionally studied in the context of immunity. The promoter of dChi contained activator protein 1 (AP-1) binding sites and could be regulated by c-Jun. The recombinant dChi protein showed no bacteriostatic activity in vitro but knockdown of dChi in vivo increased the mortality of shrimp and the bacterial load in tissues after Vibrio parahaemolyticus infection, suggesting that dChi could play a positive role in antibacterial responses. However, silencing of dChi expression significantly decreased the mortality of WSSV-infected shrimp and down-regulated the viral load in tissues, indicating that dChi could facilitate WSSV infection. We further demonstrated that dChi was involved in regulation of the bacterial phagocytosis of hemocytes and expression of a series of immune related transcription factors and antimicrobial peptides. These indicated that the roles of dChi in antibacterial responses and anti-WSSV responses in vivo could result from its regulatory effects on the immune system. Taken together, the current study suggests that double chitin catalytic domain-containing chitinases could be important players in immune regulation in crustaceans.

Optimization of high molecular weight DNA extraction methods in shrimp for a long-read sequencing platform

Marine organisms are important to global food security as they are the largest source of animal proteins feeding mankind. Genomics-assisted aquaculture can increase yield while preserving the environment to ensure sufficient and sustainable production for global food security. However, only few high-quality genome sequences of marine organisms, especially shellfish, are available to the public partly because of the difficulty in the sequence assembly due to the complex nature of their genomes. A key step for a successful genome sequencing is the preparation of high-quality high molecular weight (HMW) genomic DNA. This study evaluated the effectiveness of five DNA extraction protocols (CTAB, Genomic-tip, Mollusc DNA, TIANamp Marine Animals DNA, and Sbeadex livestock kits) in obtaining shrimp HMW DNA for a long-read sequencing platform. DNA samples were assessed for quality and quantity using a Qubit fluorometer, NanoDrop spectrophotometer and pulsed-field gel electrophoresis. Among the five extraction methods examined without further optimization, the Genomic-tip kit yielded genomic DNA with the highest quality. However, further modifications of these established protocols might yield even better DNA quality and quantity. To further investigate whether the obtained genomic DNA could be used in a long-read sequencing application, DNA samples from the top three extraction methods (CTAB method, Genomic-tip and Mollusc DNA kits) were used for Pacific Biosciences (PacBio) library construction and sequencing. Genomic DNA obtained from Genomic-tip and Mollusc DNA kits allowed successful library construction, while the DNA obtained from the CTAB method did not. Genomic DNA isolated using the Genomic-tip kit yielded a higher number of long reads (N50 of 14.57 Kb) than those obtained from Mollusc DNA kits (N50 of 9.74 Kb). Thus, this study identified an effective extraction method for high-quality HMW genomic DNA of shrimp that can be applied to other marine organisms for a long-read sequencing platform.

Ecdysis triggering hormone modulates molt behaviour in the redclaw crayfish Cherax quadricarinatus, providing a mechanistic evidence for conserved function in molt regulation across Pancrustacea

Molting enables growth and development across ecdysozoa. The molting process is strictly controlled by hormones - ecdysteroids. Ecdysteroidogenesis occurs in theprothoracic glands and stimulated by prothoracicotropic hormone in insects, while it ensues in the Y-organ and regulated by the molt inhibiting hormone in crustaceans. A peak in ecdysteroids in the hemolymph induces a cascade of multiple neuropeptides including Ecdysis Triggering Hormone (ETH) and Corazonin. The role of ETH is well defined in controlling the molt process in insects, but it is yet to be defined in crustaceans. In this study, we investigated the behavioral response of intermolt crayfish to ETH and Corazonin injections as well as the impact of ETH on the molt period using in vivo assays. Injection of Corazonin and ETH resulted in a clear and immediate eye twitching response to these two neuropeptides. The Corazonin injection induced eye twitching in slow and asynchronous manner, while ETH injection caused eye twitching in a relatively fast and synchronous way. A single injection of ETH to crayfish resulted in a remarkable prolong molt period, at twice the normal molting cycle, suggesting that ETH plays a key role in controlling the molt cycle in decapod crustaceans. Given the key significance of ETH in molt regulation and its plausible application in pest control, we characterized ETH across the pancrustacean orders. Bioinformatic analysis shows the mature ETH sequence is identical in all studied decapod species. ETH can be classified into specific groups based on the associated motif in each insect order and shows an insect motif -KxxPRx to be conserved in crustaceans.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Chemosensory Basis of Feeding Behavior in Pacific White Shrimp, Litopenaeus vannamei

AbstractThe Pacific white shrimp, Litopenaeus vannamei, is important as the principal species in the worldwide aquaculture of shrimp. It has also become a model in the study of crustacean biology, especially because it is one of the first decapod crustaceans to have its genome sequenced. This study examined an aspect of the sensory biology of this shrimp that is important in its aquaculture, by describing its peripheral chemical sensors and how they are used in acquiring and consuming food pellets. We used scanning electron microscopy to describe the diversity of sensilla on the shrimp's major chemosensory organs: antennules, antennae, mouthparts, and legs. Using behavioral studies on animals with selective sensory ablations, we then explored the roles that these chemosensory organs play in the shrimp's search for, and acquisition and ingestion of, food pellets. We found that the antennules mediate odor-activated searching for pellets, with both the lateral and medial antennular flagella contributing to this behavior and thus demonstrating that both aesthetasc (olfactory) and distributed chemosensors on the antennules can mediate this behavior. Once the shrimp finds and grasps the food pellet, the antennular chemoreceptors no longer play a role, and then the chemoreceptors on the mouthparts and legs control ingestion of the pellets. This sequence of chemosensory control of feeding in L. vannamei, a dendrobranchiate crustacean with small antennules and an ability to live and feed in both benthic and pelagic environments, is generally similar to that of the better-studied, large-antennuled, benthic reptantian crustaceans, including spiny lobsters (Achelata), clawed lobsters and crayfish (Astacidea), and crabs (Meirua).

Construction of a High-Density Genetic Map and Identification of Quantitative Trait Loci for Nitrite Tolerance in the Pacific White Shrimp (Litopenaeus vannamei)

Nitrite is a major environmental toxin in aquaculture systems that disrupts multiple physiological functions in aquatic animals. Although nitrite tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available. In this study, we constructed a high-density genetic map of Litopenaeus vannamei with 17,242 single nucleotide polymorphism markers spanning 6,828.06 centimorgans (cM), with an average distance of 0.4 cM between adjacent markers on 44 linkage groups (LGs). Using this genetic map, we identified two markers associated with nitrite tolerance. We then sequenced the transcriptomes of the most nitrite-tolerant and nitrite-sensitive individuals from each of four genetically distinct L. vannamei families (LV-I–4). We found 2,002, 1,983, 1,954, and 1,867 differentially expressed genes in families LV-1, LV-2, LV-3, and LV-4, respectively. By integrating QTL and transcriptomics analyses, we identified a candidate gene associated with nitrite tolerance. This gene was annotated as solute carrier family 26 member 6 (SLC26A6). RNA interference (RNAi) analysis demonstrated that SLC26A6 was critical for nitrite tolerance in L. vannamei. The present study increases our understanding of the molecular mechanisms underlying nitrite tolerance in shrimp and provides a basis for molecular-marker-assisted shrimp breeding.

Immune activation by a multigene family of lectins with variable tandem repeats in oriental river prawn (Macrobrachium nipponense)

Genomic regions with repeated sequences are unstable and prone to rapid DNA diversification. However, the role of tandem repeats within the coding region is not fully characterized. Here, we have identified a new hypervariable C-type lectin gene family with different numbers of tandem repeats (Rlecs; R means repeat) in oriental river prawn (Macrobrachium nipponense). Two types of repeat units (33 or 30 bp) are identified in the second exon, and the number of repeat units vary from 1 to 9. Rlecs can be classified into 15 types through phylogenetic analysis. The amino acid sequences in the same type of Rlec are highly conservative outside the repeat regions. The main differences among the Rlec types are evident in exon 5. A variable number of tandem repeats in Rlecs may be produced by slip mispairing during gene replication. Alternative splicing contributes to the multiplicity of forms in this lectin gene family, and different types of Rlecs vary in terms of tissue distribution, expression quantity and response to bacterial challenge. These variations suggest that Rlecs have functional diversity. The results of experiments on sugar binding, microbial inhibition and clearance, regulation of antimicrobial peptide gene expression and prophenoloxidase activation indicate that the function of Rlecs with the motif of YRSKDD in innate immunity is enhanced when the number of tandem repeats increases. Our results suggest that Rlecs undergo gene expansion through gene duplication and alternative splicing, which ultimately leads to functional diversity.

Gonadulins, the fourth type of insulin-related peptides in decapods

Insulin and related peptides play important roles in the regulation of growth and reproduction. Until recently three different types of insulin-related peptides had been identified from decapod crustaceans. The identification of two novel insulin-related peptides from Sagmariasus verreauxi and Cherax quadricarinatus suggested that there might a fourth type. Publicly available short read archives show that orthologs of these peptides are commonly present in these animals. Most decapods have two genes coding such peptides, but Penaeus species have likely only one and some palaemonids have three. Interestingly, expression levels can vary more than thousand-fold in the gonads of Portunus trituberculatus, where gonadulin 1 is expressed by the testis and gonadulin 2 by the ovary. Although these peptides are also expressed in other tissues, the occasionally very high expression in the gonads led to them being called gonadulins.

The "IAG-Switch"-A Key Controlling Element in Decapod Crustacean Sex Differentiation

The androgenic gland (AG)-a unique crustacean endocrine organ that secretes factors such as the insulin-like androgenic gland (IAG) hormone-is a key player in crustacean sex differentiation processes. IAG expression induces masculinization, while the absence of the AG or a deficiency in IAG expression results in feminization. Therefore, by virtue of its universal role as a master regulator of crustacean sexual development, the IAG hormone may be regarded as the sexual "IAG-switch." The switch functions within an endocrine axis governed by neuropeptides secreted from the eyestalks, and interacts downstream with specific insulin receptors at its target organs. In recent years, IAG hormones have been found-and sequenced-in dozens of decapod crustacean species, including crabs, prawns, crayfish and shrimps, bearing different types of reproductive strategies-from gonochorism, through hermaphroditism and intersexuality, to parthenogenesis. The IAG-switch has thus been the focus of efforts to manipulate sex developmental processes in crustaceans. Most sex manipulations were performed using AG ablation or knock-down of the IAG gene in males in order to sex reverse them into "neo-females," or using AG implantation/injecting AG extracts or cells into females to produce "neo-males." These manipulations have highlighted the striking crustacean sexual plasticity in different species and have permitted the manifestation of either maleness or femaleness without altering the genotype of the animals. Furthermore, these sex manipulations have not only facilitated fundamental studies of crustacean sexual mechanisms, but have also enabled the development of the first IAG-switch-based monosex population biotechnologies, primarily for aquaculture but also for pest control. Here, we review the crustacean IAG-switch, a unique crustacean endocrine mechanism, from the early discoveries of the AG and the IAG hormone to recent IAG-switch-based manipulations. Moreover, we discuss this unique early pancrustacean insulin-based sexual differentiation control mechanism in contrast to the extensively studied mechanisms in vertebrates, which are based on sex steroids.

The ins and outs of host-microsporidia interactions during invasion, proliferation and exit

Microsporidia are a large group of fungal-related obligate intracellular parasites. They are responsible for infections in humans as well as in agriculturally and environmentally important animals. Although microsporidia are abundant in nature, many of the molecular mechanisms employed during infection have remained enigmatic. In this review, we highlight recent work showing how microsporidia invade, proliferate and exit from host cells. During invasion, microsporidia use spore wall and polar tube proteins to interact with host receptors and adhere to the host cell surface. In turn, the host has multiple defence mechanisms to prevent and eliminate these infections. Microsporidia encode numerous transporters and steal host nutrients to facilitate proliferation within host cells. They also encode many secreted proteins which may modulate host metabolism and inhibit host cell defence mechanisms. Spores exit the host in a non-lytic manner that is dependent on host actin and endocytic recycling proteins. Together, this work provides a fuller picture of the mechanisms that these fascinating organisms use to infect their hosts.

Transcriptomic data on the transgenerational exposure of the keystone amphipod Gammarus locusta to simvastatin

The use of transcriptomics data brings new insights and works as a powerful tool to explore the molecular mode of action (MoA) of transgenerational inheritance effects of contaminants of emerging concern. Therefore, in this dataset, we present the transcriptomic data of the transgenerational effects of environmentally relevant simvastatin levels, one of the most prescribed human pharmaceuticals, in the keystone amphipod species Gammarus locusta. In summary, G. locusta juveniles were maintained under simvastatin exposure up to adulthood (exposed group - F0E) and the offspring of F0E were transferred to control water for the three subsequent generations (transgenerational group - F1T, F2T and F3T). To gain insights into the biological functions and canonical pathways transgenerationally disrupted by simvastatin, a G. locusta de novo transcriptome assembly was produced and the transcriptomic profiles of three individual G. locusta females, per group, over the four generations (F0 to F3) - solvent control groups (F0.C, F1.C, F2.C and F3.C), F0 320 ng/L simvastatin exposed group (F0.320E) and F1 to F3 320 transgenerational group (F1.320T; F2.320T and F3.320T) - were analyzed. Briefly, Illumina HiSeq™ 2500 platform was used to perform RNA sequencing, and due to the unavailability of G. locusta genome, the RNA-seq datasets were assembled de novo using Trinity and annotated with Trinotate software. After assembly and post-processing steps, 106093 transcripts with N50 of 2371 bp and mean sequence length of 1343.98 bp was produced. BUSCO analyses showed a transcriptome with gene completeness of 97.5 % Arthropoda library profile. The Bowtie2, RSEM and edgeR tools were used for the differential gene expression (DEGs) analyses that allowed the identification of a high quantity of genes differentially expressed in all generations. Finally, to identify the main metabolic pathways affected by the transgenerational effects of SIM across all generations, the DGEs genes were blasted onto KEGG pathways database using the KAAS webserver. The data furnished in this article allows a better molecular understanding of the transgenerational effects produced by simvastatin in the keystone amphipod G. locusta and has major implications for hazard and risk assessment of pharmaceuticals and other emerging contaminants. This article is related to the research article entitled "Transgenerational inheritance of chemical-induced signature: a case study with simvastatin [1].

Illuminating the impact of diel vertical migration on visual gene expression in deep-sea shrimp

Diel vertical migration (DVM) of marine animals represents one of the largest migrations on our planet. Migrating fauna are subjected to a variety of light fields and environmental conditions that can have notable impacts on sensory mechanisms, including an organism's visual capabilities. Among deep-sea migrators are oplophorid shrimp that vertically migrate hundreds of metres to feed in shallow waters at night. These species also have bioluminescent light organs that emit light during migrations to aid in camouflage. The organs have recently been shown to contain visual proteins (opsins) and genes that infer light sensitivity. Knowledge regarding the impacts of vertical migratory behaviour, and fluctuating environmental conditions, on sensory system evolution is unknown. In this study, the oplophorid Systellaspis debilis was either collected during the day from deep waters or at night from relatively shallow waters to ensure sampling across the vertical distributional range. De novo transcriptomes of light-sensitive tissues (eyes/photophores) from the day/night specimens were sequenced and analysed to characterize opsin diversity and visual/light interaction genes. Gene expression analyses were also conducted to quantify expression differences associated with DVM. Our results revealed an expanded opsin repertoire among the shrimp and differential opsin expression that may be linked to spectral tuning during the migratory process. This study sheds light on the sensory systems of a bioluminescent invertebrate and provides additional evidence for extraocular light sensitivity. Our findings further suggest opsin co-expression and subsequent fluctuations in opsin expression may play an important role in diversifying the visual responses of vertical migrators.

Feed competition reduces heritable variation for body weight in Litopenaeus vannamei

Background

Competition is a common social interaction among shrimp and depending on its intensity, it can affect heritable variation and response to selection. Little is known about the variance of indirect genetic effects (IGE) under competitive and non-competitive conditions in shrimp. In this study, we used extended mixed linear models to estimate genetic parameters for the direct genetic effect (DGE) and IGE on body weight in Litopenaeus vannamei raised under ad libitum (AF, non-competitive environment) and restricted (RF, competitive environment) feeding regimes.

Results

Estimates of heritabilities for body weight obtained with a traditional animal model (i.e. without accounting for IGE) were 0.11 ± 0.09 under AF and 0.25 ± 0.11 under RF. With extended animal models that accounted for IGE, the corresponding estimates for body weight were 0.07 ± 0.08 and 0.34 ± 0.11. Thus, heritabilities were higher under the RF regime than under the AF regime, regardless of whether IGE was accounted for or not. The log-likelihood ratio test revealed significant IGE under the RF regime. Although estimates of indirect genetic variance were low (0.0023 ± 0.0013 for AF and 0.0028 ± 0.0012 for RF), they contributed substantially to the total heritable variance: 66.8% for AF and 692.2% for RF. The total heritable variance was smaller under the RF regime (0.7 ± 1.3) than under the AF regime (5.8 ± 2.6) because of the high contribution of the negative covariance between DGE and IGE (− 7.03). Estimates of the correlation between DGE and IGE were 0.32 ± 0.47 under AF and − 0.93 ± 0.15 under RF, those of DGE and IGE for body weight between both regimes were 0.94 ± 0.07 and 0.67 ± 0.20, respectively, and those of IGE for body weight with DGE for survival were − 0.12 ± 0.22 under AF and − 0.58 ± 0.20 under RF.

Conclusions

These results indicate that strong competitive interactions occurred under the RF regime in L. vannamei. Significant reranking and variation in IGE of individuals were observed between the two feeding regimes. Strong competitive interactions reduced the total heritable variation for body weight when food was restricted. These results indicate that the extent of competition among L. vannamei depends on the feeding regime applied and that this competition affects the genetic basis of body weight.

The functional relevance of shrimp C-type lectins in host-pathogen interactions

C-type lectins (CTLs) are key recognition proteins in shrimp immunity. A few years ago we reviewed sequence information, ligand specificity, expression profiles and specific functions of the shrimp CTLs. Since then, multiple integrated studies that implemented biochemical approaches using both the native and recombinant proteins, functional genetic approaches using RNA interference, and mechanistic studies by analyzing protein-protein interactions were carried out. Results from these rigorous studies revealed the functions and mechanisms of action of selected members of the shrimp CTL family. This review focuses on this new knowledge, that includes unique structural aspects, functions, and mechanisms in host-pathogen interactions, the functional relevance of regions other than the C-type lectin domain, and the regulation of transcription of shrimp CTLs. Thus, this review aims to provide a detailed update of recent studies that have contributed to our better understanding of the shrimp immune events that involve CTL functions.

The Polymorphism of LvMMD2 and Its Association with Growth Traits in Litopenaeus vannamei

The Pacific white shrimp Litopenaeus vannamei is one of the major economic aquaculture species. The growth trait is considered as the most important trait in L. vannamei aquaculture. Identification of the genetic components underlying growth-related traits in L. vannamei could be useful for the selective breeding of growth trait. Our previous work identified several growth-related SNPs by genome-wide association study (GWAS). Based on the assembled genome, we identified a new candidate gene (LvMMD2) beside the associated marker. This gene encodes the progestin and AdipoQ receptor 10 (PAQR10) protein. We further investigate the polymorphisms of LvMMD2 and their association with body weight of L. vannamei. By resequencing the coding region of LvMMD2, a total of 8 SNPs were identified, including 6 synonymous mutations and 2 nonsynonymous mutations. Association analyses based on a population of 322 individuals revealed that several SNPs located in the coding region of LvMMD2 were significantly associated with the body weight, especially the nonsynonymous mutation named as MMD_5 contributed the most association to the trait and it could explain 10.5% of phenotypic variance. In addition, several genes involved in growth and development have been identified as LvMMD2-interacting genes. These findings strongly suggested that LvMMD2 might be an important gene regulating the shrimp growth. More importantly, the MMD_5 could be a promising candidate locus for marker-assisted selection (MAS) of the body weight in L. vannamei.

High accuracy of pooled DNA genotyping by 2b-RAD sequencing in the Pacific white shrimp, Litopenaeus vannamei

Using pooled DNA genotyping to estimate the proportional contributions from multiple families in a pooled sample is of particular interest for selective breeding in aquaculture. We compared different pooled libraries with separate 2b-RAD sequencing of Litopenaeus vannamei individuals to assess the effect of different population structures (different numbers of individuals and families) on pooled DNA sequencing, the accuracy of parent sequencing of the DNA pools and the effect of SNP numbers on pooled DNA sequencing. We demonstrated that small pooled DNA genotyping of up to 53 individuals by 2b-RAD sequencing could provide a highly accurate assessment of population allele frequencies. The accuracy increased as the number of individuals and families increased. The allele frequencies of the parents from each pool were highly correlated with those of the pools or the corresponding individuals in the pool. We chose 500-28,000 SNPs to test the effect of SNP number on the accuracy of pooled sequencing, and no linear relationship was found between them. When the SNP number was fixed, increasing the number of individuals in the mixed pool resulted in higher accuracy of each pooled genotyping. Our data confirmed that pooled DNA genotyping by 2b-RAD sequencing could achieve higher accuracy than that of individual-based genotyping. The results will provide important information for shrimp breeding programs.

Genomic Characterization and Expression of Juvenile Hormone Esterase-Like Carboxylesterase Genes in Pacific White Shrimp, Litopenaeus vannamei

The sesquiterpenoid methyl farnesoate (MF), a juvenile hormone (JH) analog, plays important roles in many physiological processes of crustaceans, such as morphogenesis, molting and reproduction. Juvenile hormone esterase-like (JHE-like) carboxylesterase (CXE) is a key enzyme in MF degradation, playing a significant role in regulating MF titer. However, its function is barely known in shrimp. In this study, a total of 21 JHE-like CXEs (LvCXEs) were characterized in Pacific white shrimp Litopenaeus vannamei, based on the full genome and multi-transcriptomic data. LvCXE has a conserved triplet catalytic site (Ser-Glu-His) and a characteristic GxSxG motif. Most LvCXEs were highly expressed in the hepatopancreas, which was the main site for MF degradation. LvCXEs containing a GESAG motif showed a specific expansion in the L. vannamei genome. Those GESAG-containing LvCXEs presented differential expressions at different larvae stages and different molting stages of L. vannamei, which suggested their potential functions in development and molting. Additionally, when the transcription level of CXEs was inhibited, it could lead to failed molt and death of L. vannamei. When we further detected the expression levels of the key ecdysone responsive transcription factors including LvE75, LvBr-C, LvHr3 and LvFtz-f1 after the CXE inhibitor was injected into L. vannamei, they all showed apparent down-regulation. These results suggested that the expansion of LvCXEs in the L. vannamei genome should contribute to the regulation of metamorphosis at larvae stages and frequent molting during the growth of L. vannamei.

Impacts of Methyl Farnesoate and 20-Hydroxyecdysone on Larval Mortality and Metamorphosis in the Kuruma Prawn Marsupenaeus japonicus

Physiological functions of juvenile hormone (JH) and molting hormone have been demonstrated in insects. JH, molting hormone and their mimics (insect growth regulators, IGRs) show endocrine-disrupting effects not only on target pest insects but also on other arthropod species such as crustaceans. However, little is known about the endocrine-disrupting effects of IGRs on benthic crustaceans. In this study, laboratory experiments were conducted to investigate effects of representative innate JH in crustaceans (methyl farnesoate, MF) and molting hormone (20-hydroxyecdysone, 20E, active form of ecdysteroid) on larval stages of the kuruma prawn Marsupenaeus japonicus, which is a decapod crustacean living in warm seawater. Larval development of kuruma prawn progresses in the order of nauplius, zoea, mysis, and then post-larvae with molting and metamorphosis, but it is unknown whether both MF and 20E have crucial roles in metamorphosis and molting of this species. Treatments of either MF or 20E on shrimp larvae were attempted at each developmental stage and those effects were validated. In terms of EC50 values between mortality and metamorphosis, there were apparent differences in the transition from nauplius to zoea (MF: 7.67 and 0.12 μM; 20E: 3.84 and 0.06 μM in survival and metamorphic rates, respectively). In contrast, EC50 values in MF and 20E treatments showed high consistency in the transitions between zoea to mysis (EC50 values for survival; MF: 1.25 and 20E: 0.22 μM), and mysis to post-larvae (EC50 values for survival; MF: 0.65 and 20E: 0.46 μM). These data suggest that nauplius has strong resistance against exposure to MF and 20E. Moreover, both chemicals induced high mortality triggered by the disruption of molting associated with metamorphosis. To our knowledge, this is the first experimental evidence that investigates in vivo physiological functions of MF and 20E in the larval stages of kuruma prawn, shedding light on not only ecotoxicological impacts of IGRs released into nature, but also endocrine mechanisms underlying larval development with metamorphosis in benthic decapod crustaceans.

Shedding the Light on Litopenaeus vannamei Differential Muscle and Hepatopancreas Immune Responses in White Spot Syndrome Virus (WSSV) Exposure

White Spot Syndrome Virus (WSSV) is one of the main threats to farming Litopenaeus vannamei, the most important crustacean commercialized in aquaculture worldwide. Here, we performed RNA-seq analyses in hepatopancreas and muscle from WSSV-negative (healthy) and WSSV-positive (unhealthy) L. vannamei, previously exposed to the virus, to obtain new insights about the molecular basis of resistance to WSSV. We detected 71% of our reads mapped against the recently described L. vannamei genome. This is the first report mapping RNA-seq transcripts from shrimps exposed to WSSV against the species reference genome. Differentially expressed gene (DEG) analyses were performed for four independent comparisons, and 13,338 DEGs were identified. When the redundancies and isoforms were disregarded, we observed 8351 and 6514 DEGs, respectively. Interestingly, after crossing the data, we detected a common set of DEGs for hepatopancreas and healthy shrimps, as well as another one for muscle and unhealthy shrimps. Our findings indicate that genes related to apoptosis, melanization, and the Imd pathway are likely to be involved in response to WSSV, offering knowledge about WSSV defense in shrimps exposed to the virus but not infected. These data present potential to be applied in further genetic studies in penaeids and other farmed shrimp species.

Molecular and Functional Diversity of Crustin-Like Genes in the Shrimp Litopenaeus vannamei

Crustins are crustacean cationic cysteine-rich antimicrobial peptides that contain one or two whey acidic protein (WAP) domain(s) at the carboxyl terminus and mainly show antimicrobial and/or proteinase inhibitory activities. Here, we performed genome and transcriptome screening and identified 34 full-length crustin-like encoding genes in Litopenaeus vannamei. Multiple sequence analysis of the deduced mature peptides revealed that these putative crustins included 10 type Ia, two type Ib, one type Ic, 11 type IIa, three type IIb, four type III, one type IV, one type VI, and one type VII. These putative crustins were clustered into different groups. Phylogenetic analysis, considering their domain composition, showed that different types of crustin-like genes in crustaceans might be originated from the WAP core region, along with sequence insertion, duplication, deletion, and amino acid substitution. Tissue distribution analysis suggested that most crustin-like genes were mainly detected in immune-related tissues while several crustin-like genes exhibited tissue-specific expression patterns. Quantitative PCR analysis on 15 selected crustin-like genes showed that most of them were apparently upregulated after Vibrio parahaemolyticus or white spot syndrome virus (WSSV) infection. One type Ib crustin-like gene, mainly expressed in the ovary, showed the highest expression levels before the gastrula stage and was hardly detected after the limb bud stage, suggesting that it was a maternal immune effector. Collectively, the present data revealed the molecular and functional diversity of crustins and their potential evolutionary routes in crustaceans.

The Seminavis robusta genome provides insights into the evolutionary adaptations of benthic diatoms

Benthic diatoms are the main primary producers in shallow freshwater and coastal environments, fulfilling important ecological functions such as nutrient cycling and sediment stabilization. However, little is known about their evolutionary adaptations to these highly structured but heterogeneous environments. Here, we report a reference genome for the marine biofilm-forming diatom Seminavis robusta, showing that gene family expansions are responsible for a quarter of all 36,254 protein-coding genes. Tandem duplications play a key role in extending the repertoire of specific gene functions, including light and oxygen sensing, which are probably central for its adaptation to benthic habitats. Genes differentially expressed during interactions with bacteria are strongly conserved in other benthic diatoms while many species-specific genes are strongly upregulated during sexual reproduction. Combined with re-sequencing data from 48 strains, our results offer insights into the genetic diversity and gene functions in benthic diatoms.

The sterol regulatory element binding protein homolog of Penaeus vannamei modulates fatty acid metabolism and immune response

The sterol regulatory element binding proteins (SREBPs) transcription factors family, which regulate the expression of genes involved in cellular lipid metabolism and homeostasis, have recently been implicated in various physiological and pathophysiological processes such as immune regulation and inflammation in vertebrates. Consistent with other invertebrates, we identified a single SREBP ortholog in Penaeus vannamei (designated PvSREBP) with transcripts ubiquitously expressed in tissues and induced by lipopolysaccharide (LPS), Vibrio parahaemolyticus and Streptococcus iniae. In vivo RNA interference (RNAi) of PvSREBP attenuated the expression of several fatty acid metabolism-related genes (i.e., cyclooxygenase (PvCOX), lipoxygenase (PvLOX), fatty acid binding protein (PvFABP) and fatty acid synthase (PvFASN)), which consequently decreased the levels of total polyunsaturated fatty acids (ΣPUFAs). In addition, PvSREBP silencing decreased transcript levels of several immune-related genes such as hemocyanin (PvHMC) and trypsin (PvTrypsin), as well as genes encoding for heat-shock proteins (i.e., PvHSP60, PvHSP70 and PvHSP90). Moreover, in silico analysis revealed the presence of SREBP binding motifs on the promoters of most of the dysregulated genes, while shrimp depleted of PvSREBP were more susceptible to V. parahaemolyticus infection. Collectively, we demonstrated the involvement of shrimp SREBP in fatty acids metabolism and immune response, and propose that PvSREBP and PvHMC modulate each other through a feedback mechanism to establish homeostasis. The current study is the first to show the dual role of SREBP in fatty acid metabolism and immune response in invertebrates and crustaceans.

Sex-Specific Transcriptome Sequencing of Zoea I Larvae and Identification of Sex-Linked Genes Using Bulked Segregant Analysis in Pacific White Shrimp Litopenaeus vannamei

In order to learn the sex determination and differentiation process in shrimp, we performed sex-specific transcriptome sequencing of the Zoea I larvae in Litopenaeus vannamei (L. vannamei) using a DNA/RNA co-extraction method. The sex-specific genes and the sex-linked SNPs were identified. De novo assembly of all the clean reads generated 41,270 unigenes with an average length of 1026 bp and an N50 of 1922 bp. A total of 60 differentially expressed genes (DEGs) between females and males at Zoea I stage were obtained, in which 41 unigenes had annotations. DEGs that might be related to sex development in L. vannamei were analyzed. Besides, paralogous genes of some known sex-related genes in model organisms such as sex-lethal (Sxl) showed no sex-biased difference, which indicated they may not play roles in sex development at Zoea I stage. Candidate sex-linked SNPs were validated in six populations of L. vannamei, and SNPs in Unigene0020898 and Unigene0020336 were fully linked to genders, suggesting the two genes located in the sex-determining region of L. vannamei. Our report provided the sex-specific gene expression in early developmental stage of L. vannamei, and found some sex-biased genes that may participate in the sex development process. We also detected two sex-linked genes which may be located in the sex determination region of L. vannamei through bulked segregant analysis, and the SNPs on these genes also provided a method for genetic sex identification at transcriptional level. Our findings will facilitate further researches on molecular mechanism of sex determination and differentiation in shrimps.

Vitellogenin Receptor (VgR) Mediates Oocyte Maturation and Ovarian Development in the Pacific White Shrimp (Litopenaeus vannamei)

Oocyte maturation and ovarian development are sequentially coordinated events critical to reproduction. In the ovaries of adult oviparous animals such as birds, bony fish, insects, and crustaceans, vitellogenin receptor (VgR) is a plasma membrane receptor that specifically mediates vitellogenin (Vg) transport into oocytes. Accumulation of Vg drives sexual maturation of the female crustaceans by acting as a pivotal regulator of nutritional accumulation within oocytes, a process known as vitellogenesis. However, the mechanisms by which VgR mediates vitellogenesis are still not fully understood. In this study, we first identified a unique VgR (Lv-VgR) and characterized its genomic organization and protein structural domains in Litopenaeus vannamei, a predominant cultured shrimp species worldwide. This newly identified Lv-VgR phylogenetically forms a group with VgRs from other crustacean species within the arthropod cluster. Duplicated LBD/EGFD regions are found exclusively among arthropod VgRs but not in paralogs from vertebrates and nematodes. In terms of expression patterns, Lv-VgR transcripts are specifically expressed in ovaries of female shrimps, which increases progressively during ovarian development, and rapidly declines toward embryonic development. The cellular and subcellular locations were For analyzed by in situ hybridization and immunofluorescence, respectively. The Lv-VgR mRNA was found to be expressed in the oocytes of ovaries, and Lv-VgR protein was found to localize in the cell membrane of maturing oocytes while accumulation of the ligand Vg protein assumed an even cytoplasmic distribution. Silencing of VgR transcript expression by RNAi was effective for stunting ovarian development. This present study has thus provided new insights into the regulatory roles of VgR in crustacean ovarian development.

Crystal structure of a C-type lysozyme from Litopenaeus vanamei exhibiting a high binding constant to its chitotriose inhibitor

Although information about invertebrate lysozymes is scarce, these enzymes have been described as components of the innate immune system, functioning as antibacterial proteins. Here we describe the first thermodynamic and structural study of a new C-type lysozyme from a Pacific white shrimp Litopenaeus vannamei (LvL), which has shown high activity against both Gram (+) and Gram (-) bacteria including Vibrio sp. that is one of the most severe pathogens in penaeid shrimp aquaculture. Compared with hen egg-white lysozyme, its sequence harbors a seven-residue insertion from amino acid 97 to 103, and a nine-residue extension at the C-terminus only found in penaeid crustaceans, making this enzyme one of the longest lysozyme reported to date. LvL was crystallized in the presence and absence of chitotriose. The former crystallized as a monomer in space group P6₁ and the latter in P2₁2₁2₁ with two monomers in the asymmetric unit. Since the enzyme crystallized at a pH where lysozyme activity is deficient, the ligand could not be observed in the P6₁ structure; therefore, we performed a docking simulation with chitotriose to compare with the hen egg lysozyme crystallized in the presence of the ligand. Remarkably, additional amino acids in LvL caused an increase in the length of α-helix H4 (residues 97-103) that is directly related to ligand recognition. The K_a for chitotriose (4.1 × 10⁵ M^-1), as determined by Isothermal Titration Calorimetry, was one order of magnitude higher than those for lysozymes from hen and duck eggs. Our results revealed new interactions of chitiotriose with residues in helix H4.

Harnessing genomics to fast-track genetic improvement in aquaculture

Aquaculture is the fastest-growing farmed food sector and will soon become the primary source of fish and shellfish for human diets. In contrast to crop and livestock production, aquaculture production is derived from numerous, exceptionally diverse species that are typically in the early stages of domestication. Genetic improvement of production traits via well-designed, managed breeding programmes has great potential to help meet the rising seafood demand driven by human population growth. Supported by continuous advances in sequencing and bioinformatics, genomics is increasingly being applied across the broad range of aquaculture species and at all stages of the domestication process to optimize selective breeding. In the future, combining genomic selection with biotechnological innovations, such as genome editing and surrogate broodstock technologies, may further expedite genetic improvement in aquaculture.

Exceptional diversity of opsin expression patterns in Neogonodactylus oerstedii (Stomatopoda) retinas

Stomatopod crustaceans possess some of the most complex animal visual systems, including at least 16 spectrally distinct types of photoreceptive units (e.g., assemblages of photoreceptor cells). Here we fully characterize the set of opsin genes expressed in retinal tissues and determine expression patterns of each in the stomatopod Neogonodactylus oerstedii Using a combination of transcriptome and RACE sequencing, we identified 33 opsin transcripts expressed in each N. oerstedii eye, which are predicted to form 20 long-wavelength-sensitive, 10 middle-wavelength-sensitive, and three UV-sensitive visual pigments. Observed expression patterns of these 33 transcripts were highly unusual in five respects: 1) All long-wavelength and short/middle-wavelength photoreceptive units expressed multiple opsins, while UV photoreceptor cells expressed single opsins; 2) most of the long-wavelength photoreceptive units expressed at least one middle-wavelength-sensitive opsin transcript; 3) the photoreceptors involved in spatial, motion, and polarization vision expressed more transcripts than those involved in color vision; 4) there is a unique opsin transcript that is expressed in all eight of the photoreceptive units devoted to color vision; and 5) expression patterns in the peripheral hemispheres of the eyes suggest visual specializations not previously recognized in stomatopods. Elucidating the expression patterns of all opsin transcripts expressed in the N. oerstedii retina reveals the potential for previously undocumented functional diversity in the already complex stomatopod eye and is a first step toward understanding the functional significance of the unusual abundance of opsins found in many arthropod species' visual systems.

The First Genome Survey of the Antarctic Krill (Euphausia superba) Provides a Valuable Genetic Resource for Polar Biomedical Research

The world-famous Antarctic krill (Euphausia superba) plays a fundamental role in the Antarctic food chain. It resides in cold environments with the most abundant biomass to support the Antarctic ecology and fisheries. Here, we performed the first genome survey of the Antarctic krill, with genomic evidence for its estimated genome size of 42.1 gigabases (Gb). Such a large genome, however, is beyond our present capability to obtain a good assembly, although our sequencing data are a valuable genetic resource for subsequent polar biomedical research. We extracted 13 typical protein-coding gene sequences of the mitochondrial genome and analyzed simple sequence repeats (SSRs), which are useful for species identification and origin determination. Meanwhile, we conducted a high-throughput comparative identification of putative antimicrobial peptides (AMPs) and antihypertensive peptides (AHTPs) from whole-body transcriptomes of the Antarctic krill and its well-known counterpart, the whiteleg shrimp (Penaeus vannamei; resident in warm waters). Related data revealed that AMPs/AMP precursors and AHTPs were generally conserved, with interesting variations between the two crustacean species. In summary, as the first report of estimated genome size of the Antarctic krill, our present genome survey data provide a foundation for further biological research into this polar species. Our preliminary investigations on bioactive peptides will bring a new perspective for the in-depth development of novel marine drugs.

Comparative transcriptomic analysis of deep- and shallow-water barnacle species (Cirripedia, Poecilasmatidae) provides insights into deep-sea adaptation of sessile crustaceans

BACKGROUND

Barnacles are specialized marine organisms that differ from other crustaceans in possession of a calcareous shell, which is attached to submerged surfaces. Barnacles have a wide distribution, mostly in the intertidal zone and shallow waters, but a few species inhabit the deep-sea floor. It is of interest to investigate how such sessile crustaceans became adapted to extreme deep-sea environments. We sequenced the transcriptomes of a deep-sea barnacle, Glyptelasma gigas collected at a depth of 731 m from the northern area of the Zhongjiannan Basin, and a shallow-water coordinal relative, Octolasmis warwicki. The purpose of this study was to provide genetic resources for investigating adaptation mechanisms of deep-sea barnacles.

RESULTS

Totals of 62,470 and 51,585 unigenes were assembled for G. gigas and O. warwicki, respectively, and functional annotation of these unigenes was made using public databases. Comparison of the protein-coding genes between the deep- and shallow-water barnacles, and with those of four other shallow-water crustaceans, revealed 26 gene families that had experienced significant expansion in G. gigas. Functional annotation showed that these expanded genes were predominately related to DNA repair, signal transduction and carbohydrate metabolism. Base substitution analysis on the 11,611 single-copy orthologs between G. gigas and O. warwicki indicated that 25 of them were distinctly positive selected in the deep-sea barnacle, including genes related to transcription, DNA repair, ligand binding, ion channels and energy metabolism, potentially indicating their importance for survival of G. gigas in the deep-sea environment.

CONCLUSIONS

The barnacle G. gigas has adopted strategies of expansion of specific gene families and of positive selection of key genes to counteract the negative effects of high hydrostatic pressure, hypoxia, low temperature and food limitation on the deep-sea floor. These expanded gene families and genes under positive selection would tend to enhance the capacities of G. gigas for signal transduction, genetic information processing and energy metabolism, and facilitate networks for perceiving and responding physiologically to the environmental conditions in deep-sea habitats. In short, our results provide genomic evidence relating to deep-sea adaptation of G. gigas, which provide a basis for further biological studies of sessile crustaceans in the deep sea.

Air Exposure Affects Physiological Responses, Innate Immunity, Apoptosis and DNA Methylation of Kuruma Shrimp, Marsupenaeus japonicus

Air exposure stress is a common phenomenon for commercial crustacean species in aquaculture and during waterless transportation. However, the antioxidant responses to air exposure discussed in previous studies may be insufficient to present the complexities involved in this process. The comprehensive immune responses, especially considering the immune genes, cell apoptosis, and epigenetic changes, are still unknown. Accordingly, we investigated the multifaceted responses of Marsupenaeus japonicus to air exposure. The results showed that the expression profiles of the apoptosis genes (e.g., IAP, TXNIP, caspase, and caspase-3) and the hypoxia-related genes (e.g., hsp70, hif-1α, and HcY) were all dramatically induced in the hepatopancreas and gills of M. japonicus. Heart rates, T-AOC (total antioxidant capacity) and lactate contents showed time-dependent changes upon air exposure. Air exposure significantly induced apoptosis in hepatopancreas and gills. Compared with the control group, the apoptosis index (AI) of the 12.5 h experimental group increased significantly (p < 0.05) in the hepatopancreas and gills. Most individuals in the experimental group (EG, 12.5 h) had lower methylation ratios than the control group (CG). Air exposure markedly reduced the full-methylation and total-methylation ratios (31.39% for the CG and 26.46% for the EG). This study provided a comprehensive understanding of the antioxidant responses of M. japonicus considering its physiology, innate immunity, apoptosis, and DNA methylation levels, and provided theoretical guidance for waterless transportation.

Microecological Koch's postulates reveal that intestinal microbiota dysbiosis contributes to shrimp white feces syndrome

BACKGROUND

Recently, increasing evidence supports that some complex diseases are not attributed to a given pathogen, but dysbiosis in the host intestinal microbiota (IM). The full intestinal ecosystem alterations, rather than a single pathogen, are associated with white feces syndrome (WFS), a globally severe non-infectious shrimp disease, while no experimental evidence to explore the causality. Herein, we conducted comprehensive metagenomic and metabolomic analysis, and intestinal microbiota transplantation (IMT) to investigate the causal relationship between IM dysbiosis and WFS.

RESULTS

Compared to the Control shrimp, we found dramatically decreased microbial richness and diversity in WFS shrimp. Ten genera, such as Vibrio, Candidatus Bacilloplasma, Photobacterium, and Aeromonas, were overrepresented in WFS, whereas 11 genera, including Shewanella, Chitinibacter, and Rhodobacter were enriched in control. The divergent changes in these populations might contribute the observation that a decline of pathways conferring lipoic acid metabolism and mineral absorption in WFS. Meanwhile, some sorts of metabolites, especially lipids and organic acids, were found to be related to the IM alteration in WFS. Integrated with multiomics and IMT, we demonstrated that significant alterations in the community composition, functional potentials, and metabolites of IM were closely linked to shrimp WFS. The distinguished metabolites which were attributed to the IM dysbiosis were validated by feed-supplementary challenge. Both homogenous selection and heterogeneous selection process were less pronounced in WFS microbial community assembly. Notably, IMT shrimp from WFS donors eventually developed WFS clinical signs, while the dysbiotic IM can be recharacterized in recipient shrimp.

CONCLUSIONS

Collectively, our findings offer solid evidence of the causality between IM dysbiosis and shrimp WFS, which exemplify the 'microecological Koch's postulates' (an intestinal microbiota dysbiosis, a disease) in disease etiology, and inspire our cogitation on etiology from an ecological perspective. Video abstract.

Litopenaeus vannamei sirtuin 6 homolog (LvSIRT6) is involved in immune response by modulating hemocytes ROS production and apoptosis

The histone deacetylase, sirtuin 6 (SIRT6), plays an essential role in the regulation of oxidative stress, mitochondrial function and inflammation in mammals. However, the specific role of SIRT6 in invertebrate immunity has not been reported. Here, we characterized for the first time, a sirtuin 6 homolog in Litopenaeus vannamei (LvSIRT6), with full-length cDNA of 2919 bp and 1536 bp open reading frame (ORF) encoding a putative protein of 511 amino acids, which contains a typical SIR2 domain. Sequence and phylogenetic analysis revealed that LvSIRT6 shares a close evolutionary relationship with SIRT6 from invertebrates. Real-time quantitative PCR analysis of LvSIRT6 transcripts revealed that they were ubiquitously expressed in shrimp and induced in hepatopancreas and hemocytes upon challenge with Vibrio parahaemolyticus, Streptococcus iniae, lipopolysaccharide (LPS), and white spot syndrome virus (WSSV), suggesting the involvement of LvSIRT6 in shrimp immune response. Moreover, knockdown of LvSIRT6 decreased mitochondrial membrane potential and increased total ROS level in hemocytes, especially upon V. parahaemolyticus challenge. Depletion of LvSIRT6 also increased hemocytes apoptosis in terms of decreased expression of pro-survival LvBcl-2, but increased expression of pro-apoptotic LvBax and LvCytochrome C, coupled with high LvCaspase3/7 activity. Shrimp were rendered more susceptible to V. parahaemolyticus infection upon LvSIRT6 knockdown. Taken together, our present data suggest that LvSIRT6 plays an important role in shrimp immune response by modulating hemocytes ROS production and apoptosis during pathogen challenge.

A novel WSSV responsive plasma protein from Litopenaeus vannamei contributes to hemocytes anti-apoptosis

Litopenaeus vannamei, as an important marine aquaculture species, has attracted more and more attentions in past several years. More recently people got its genome fine mapping, which unveiled a gene treasure. In this study, we have identified a novel trypsin-like protein which came from previous WSSV-infected shrimp plasma iTRAQ data. This protein is a 39 kDa protein with 363 amino acids. It contains a conserved trypsin-domain and could be strongly induced with WSSV infection. Interestingly, knockdown of this protein made shrimps vulnerable to WSSV infection. Further exploration unveiled that this fragility was probably due to the fact that knockdown of this protein could cause shrimp hemocytes apoptosis, which indicated that this protein played key roles in preventing shrimp hemocytes from apoptosis. To further explore how LvTLAP protected shrimp hemocytes from apoptosis, GST pull down assay was applied to screen LvTLAP interacting protein in shrimp plasma. L. vannamei growth and transformation-dependent-like protein (LvGTD-like protein) was identified as a LvTLAP interacting protein, which played proapoptotic roles in cells. Thus, a possible explanation for LvTLAP anti-apoptosis activity was that this protein could block LvGTD-like protein proapoptotic activity to protect shrimp hemocytes from death. In general, our study has uncovered a novel WSSV responsive shrimp plasma protein, which played key roles in shrimp hemocytes anti-apoptosis and shrimp against WSSV infection.