Litopenaeus vannamei Transcriptome Profile of Populations Evaluated for Growth Performance and Exposed to White Spot Syndrome Virus (WSSV)

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

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

Long noncoding RNA profiling in hepatopancreas of Pacific white shrimp and its role in response to white spot syndrome virus infection

Long noncoding RNA (lncRNA) is a potential regulator of biological processes, including immunity, reproduction, and development. Although several transcriptome studies have focused on responses of viral infections in several organisms, the role of lncRNAs in viral responses in shrimp is still unclear. Therefore, this work aimed to identify putative lncRNAs and study their role in white spot syndrome virus (WSSV) infection in white shrimp. The hepatopancreas transcriptome from WSSV infected shrimp was analyzed in silico to identify putative lncRNAs. Among 221,347 unigenes of the de novo assembled transcriptome, 44,539 putative lncRNAs were identified, 32 of which were differentially expressed between WSSV-infected and control shrimp. Five candidate lncRNAs were validated for their expressions in shrimp tissues and in response to WSSV infection. Lnc164 was chosen for further investigation of its role in WSSV infection. Knockdown of lnc164 prolonged survival of shrimp when challenged with WSSV, suggesting a role in shrimp immunity. In addition, lnc164 was not directly involved in the control of total hemocytes and viral loads in hemolymph of WSSV-infected shrimp. A set of lnc164-regulated genes was obtained by RNA sequencing among which 251 transcripts were differentially expressed between lnc164 knockdown and control shrimp. Six immune-related genes were validated for their expression profiles. Our work sheds light on lncRNA profiles in L. vannamei in response to WSSV infection and paves the way to a functional study of lnc164 in host antiviral response.

The pyruvate kinase of the whiteleg shrimp Litopenaeus vannamei: Gene structure and responses to short term hypoxia

The shrimp Litopenaeus vannamei is the main farmed crustaceans worldwide. This crustacean suffers environmental changes in oxygen availability that affect its energy metabolism. Pyruvate kinase (PK) catalyzes the last reaction of glycolysis and is key for the regulation of glycolysis and gluconeogenesis. There is ample knowledge about mammalian PK, but in crustaceans, the information is very scarce. In this study, we analyzed in silico the structures of the PK gene and protein. Also, the effects of hypoxia on gene expression, enzymatic activity, glucose, and lactate in hepatopancreas and muscle were analyzed. The PK gene is 15,103 bp and contains 11 exons and 10 introns, producing four mRNA variants by alternative splicing and named PK1, PK2, PK3 and PK4, and two proteins with longer C-terminus and two with a 12 bp insertion. The promoter contains putative binding sites for transcription factors (TF) that are typically involved in stress responses. The deduced amino acid sequences contain the classic domains, binding sites for allosteric effectors and potential reversible phosphorylation residues. Protein modeling indicates a homotetramer with highly conserved structure. The effect of hypoxia for 6 and 12 h showed tissue-specific patterns, with higher expression, enzyme activity and lactate in muscle, but higher glucose in hepatopancreas. Changes in response to hypoxia were detected at 12 h in expression with induction in muscle and reduction in hepatopancreas, while enzyme activity was maintained, and glucose and lactate decreased. These results show rapid changes in expression and metabolites, while enzyme activity was maintained to cope with short-term hypoxia.

Comparison of Gene Expression Between Resistant and Susceptible Families Against VPAHPND and Identification of Biomarkers Used for Resistance Evaluation in Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) has caused a heavy loss to shrimp aquaculture since its outbreak. Vibrio parahaemolyticus (VP_AHPND) is regarded as one of the main pathogens that caused AHPND in the Pacific white shrimp Litopenaeus vannamei. In order to learn more about the mechanism of resistance to AHPND, the resistant and susceptible shrimp families were obtained through genetic breeding, and comparative transcriptome approach was used to analyze the gene expression patterns between resistant and susceptible families. A total of 95 families were subjected to VP_AHPND challenge test, and significant variations in the resistance of these families were observed. Three pairs of resistant and susceptible families were selected for transcriptome sequencing. A total of 489 differentially expressed genes (DEGs) that presented in at least two pairwise comparisons were screened, including 196 DEGs highly expressed in the susceptible families and 293 DEGs in the resistant families. Among these DEGs, 16 genes demonstrated significant difference in all three pairwise comparisons. Gene set enrichment analysis (GSEA) of all 27,331 expressed genes indicated that some energy metabolism processes were enriched in the resistant families, while signal transduction and immune system were enriched in the susceptible families. A total of 32 DEGs were further confirmed in the offspring of the detected families, among which 19 genes were successfully verified. The identified genes in this study will be useful for clarifying the genetic mechanism of shrimp resistance against Vibrio and will further provide molecular markers for evaluating the disease resistance of shrimp in the breeding program.

Characterization of Penaeus vannamei mitogenome focusing on genetic diversity

The diversity of the Penaeus vannamei mitochondrial genome has still been poorly characterized, there are no validated mitochondrial markers available for populational studies, and the heteroplasmy has not yet been investigated in this species. In this study, metagenomic reads extracted from the muscle of a single individual were used to assemble the mitochondrial genome (mtDNA). These data associated with mitochondrial genomes previously described allowed to evaluate the inter-individual variability and heteroplasmy. Comparison among 45 mtDNA control regions led to the detection of conserved and variable segments and the characterization of two hypervariable regions. The analysis of diversity revealed mostly low frequency polymorphisms, and heteroplasmy was found in practically all mitochondrial genes, with a high occurrence of indels. These results indicate that the design of mitochondrial markers for P. vannamei must be done with caution. The mapping of conserved and variable regions and the characterization of heteroplasmy presented here will contribute to increasing the efficiency of mitochondrial markers for population or individual studies.

A chromosome-level assembly of the black tiger shrimp (Penaeus monodon) genome facilitates the identification of growth-associated genes

To salvage marine ecosystems from fishery overexploitation, sustainable and efficient aquaculture must be emphasized. The knowledge obtained from available genome sequence of marine organisms has accelerated marine aquaculture in many cases. The black tiger shrimp (Penaeus monodon) is one of the most prominent cultured penaeid shrimps (Crustacean) with an average annual global production of half a million tons in the last decade. However, its currently available genome assemblies lack the contiguity and completeness required for accurate genome annotation due to the highly repetitive nature of the genome and technical difficulty in extracting high-quality, high-molecular weight DNA. Here, we report the first chromosome-level whole-genome assembly of P. monodon. The combination of long-read Pacific Biosciences (PacBio) and long-range Chicago and Hi-C technologies enabled a successful assembly of this first high-quality genome sequence. The final assembly covered 2.39 Gb (92.3% of the estimated genome size) and contained 44 pseudomolecules, corresponding to the haploid chromosome number. Repetitive elements occupied a substantial portion of the assembly (62.5%), the highest of the figures reported among crustacean species. The availability of this high-quality genome assembly enabled the identification of genes associated with rapid growth in the black tiger shrimp through the comparison of hepatopancreas transcriptome of slow-growing and fast-growing shrimps. The results highlighted several growth-associated genes. Our high-quality genome assembly provides an invaluable resource for genetic improvement and breeding penaeid shrimp in aquaculture. The availability of P. monodon genome enables analyses of ecological impact, environment adaptation and evolution, as well as the role of the genome to protect the ecological resources by promoting sustainable shrimp farming.

Global mRNA and miRNA Analysis Reveal Key Processes in the Initial Response to Infection with WSSV in the Pacific Whiteleg Shrimp

White Spot Disease (WSD) presents a major barrier to penaeid shrimp production. Mechanisms underlying White Spot Syndrome Virus (WSSV) susceptibility in penaeids are poorly understood due to limited information related to early infection. We investigated mRNA and miRNA transcription in Penaeus vannamei over 36 h following infection. Over this time course, 6192 transcripts and 27 miRNAs were differentially expressed—with limited differential expression from 3–12 h post injection (hpi) and a more significant transcriptional response associated with the onset of disease symptoms (24 hpi). During early infection, regulated processes included cytoskeletal remodelling and alterations in phagocytic activity that may assist WSSV entry and translocation, novel miRNA-induced metabolic shifts, and the downregulation of ATP-dependent proton transporter subunits that may impair cellular recycling. During later infection, uncoupling of the electron transport chain may drive cellular dysfunction and lead to high mortalities in infected penaeids. We propose that post-transcriptional silencing of the immune priming gene Dscam (downregulated following infections) by a novel shrimp miRNA (Pva-pmiR-78; upregulated) as a potential mechanism preventing future recognition of WSSV that may be suppressed in surviving shrimp. Our findings improve our understanding of WSD pathogenesis in P. vannamei and provide potential avenues for future development of prophylactics and treatments.

Analysis of transcriptome difference between rapid-growing and slow-growing in Penaeus vannamei

Penaeus vannamei is the principle cultured shrimp species in China. However, with the increase of culture density, the growth difference between individuals is also expanding. Here, we make use of RNA-seq to study the growth mechanisms of P. vannamei. After 120 days, we examined the transcriptomes of rapid-growing individuals (RG) and slow-growing individuals (SG). A total of 2116 and 176 differentially expressed genes (DEGs) were found in SG and RG, respectively. Moreover, the main DEGs are opsin, heat shock protein (HSP), actin, myosin, superoxide dismutase (SOD), cuticle protein, and chitinase. GO analysis further revealed that the DEGs were enriched in biological processes significantly, such as "sensory perception," "sensory perception of light stimulus," "response to stimulus," and "response to stress." Additionally, KEGG enrichment analysis showed that the DEGs were mainly enriched in "pentose and glucuronate interconversions," "amino sugar and nucleotide sugar metabolism," "glycophospholipid biosynthesis," and "glutathione metabolism." Interestingly, the upstream genes in the ecdysone signaling pathway, including molting inhibition hormone (MIH) and crustacean hyperglycemic hormone (CHH), did not differ significantly between RG and SG, which suggests that the cause for the inconsistent growth performance is due to the stress levels rather than the ecdysone signal pathway. In summary, this work provides data that will be useful for future studies on shrimp growth and development.

Multiple transcriptome mining coupled with tissue specific molecular cloning and mass spectrometry provide insights into agatoxin-like peptide conservation in decapod crustaceans

Over the past decade, in silico genome and transcriptome mining has led to the identification of many new crustacean peptide families, including the agatoxin-like peptides (ALPs), a group named for their structural similarity to agatoxin, a spider venom component. Here, analysis of publicly accessible transcriptomes was used to expand our understanding of crustacean ALPs. Specifically, transcriptome mining was used to investigate the phylogenetic/structural conservation, tissue localization, and putative functions of ALPs in decapod species. Transcripts encoding putative ALP precursors were identified from one or more members of the Penaeoidea (penaeid shrimp), Sergestoidea (sergestid shrimps), Caridea (caridean shrimp), Astacidea (clawed lobsters and freshwater crayfish), Achelata (spiny/slipper lobsters), and Brachyura (true crabs), suggesting a broad, and perhaps ubiquitous, conservation of ALPs in decapods. Comparison of the predicted mature structures of decapod ALPs revealed high levels of amino acid conservation, including eight identically conserved cysteine residues that presumably allow for the formation of four identically positioned disulfide bridges. All decapod ALPs are predicted to have amidated carboxyl-terminals. Two isoforms of ALP appear to be present in most decapod species, one 44 amino acids long and the other 42 amino acids in length, both likely generated by alternative splicing of a single gene. In carideans, a gene or terminal exon duplication appears to have occurred, with alternative splicing producing four ALPs, two 44 and two 42 amino acid isoforms. The identification of ALP precursor-encoding transcripts in nervous system-specific transcriptomes (e.g., Homarus americanus brain, eyestalk ganglia, and cardiac ganglion assemblies, finding confirmed using RT-PCR) suggests that members of this peptide family may serve as locally-released and/or hormonally-delivered neuromodulators in decapods. Their detection in testis- and hepatopancreas-specific transcriptomes suggests that members of the ALP family may also play roles in male reproduction and innate immunity/detoxification.

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.

A novel cuticle protein involved in WSSV infection to the Pacific white shrimp Litopenaeus vannamei

As the most productive crustacean species in aquaculture, Litopenaeus vannamei is seriously threatened by white spot syndrome virus (WSSV), which has caused huge economic damage in the past decades. Shrimp cuticle proteins are the important components in the frontier target tissues, including cuticle and the chitinous lining of the digestive tract. In present study, a novel cuticle protein gene, named LvCPAP1, was isolated and demonstrated to play an important role in WSSV infection. The deduced amino acid sequence of LvCPAP1 contained a signal peptide and a conserved chitin-binding domain type 2 (ChBD2). Tissue distribution analysis revealed that LvCPAP1 was predominantly expressed in epidermis and stomach. The transcription levels of LvCPAP1 in epidermis and stomach were significantly regulated upon WSSV challenge. DsRNA silencing of LvCPAP1 decreased the in vivo WSSV copy numbers and the death rate of shrimp after WSSV infection, indicating that LvCPAP1 might facilitate WSSV invasion. In addition, the interaction between LvCPAP1 and the major envelop protein VP24 of WSSV was revealed by yeast two-hybrid system and further confirmed by dot blot and pull-down assays. The present study implied that cuticle protein LvCPAP1 might favor the entry process of WSSV, which provided new clues for understanding the role of cuticle proteins during virus infection.

A new insight to biomarkers related to resistance in survived-white spot syndrome virus challenged giant tiger shrimp, Penaeus monodon

The emergence of diseases such as white spot disease has become a threat to Penaeus monodon cultivation. Although there have been a few studies utilizing RNA-Seq, the cellular processes of host-virus interaction in this species remain mostly anonymous. In the present study, P. monodon was challenged with WSSV by intramuscular injection and survived for 12 days. The effect of the host gene expression by WSSV infection in the haemocytes, hepatopancreas and muscle of P. monodon was studied using Illumina HiSeq 2000. The RNA-Seq of cDNA libraries was developed from surviving WSSV-challenged shrimp as well as from normal healthy shrimp as control. A comparison of the transcriptome data of the two groups showed 2,644 host genes to be significantly up-regulated and 2,194 genes significantly down-regulated as a result of the infection with WSSV. Among the differentially expressed genes, our study discovered HMGB, TNFSF and c-Jun in P. monodon as new potential candidate genes for further investigation for the development of potential disease resistance markers. Our study also provided significant data on the differential expression of genes in the survived WSSV infected P. monodon that will help to improve understanding of host-virus interactions in this species.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

Identification of SNPs potentially related to immune responses and growth performance in Litopenaeus vannamei by RNA-seq analyses

Litopenaeus vannamei is one of the most important shrimp species for worldwide aquaculture. Despite this, little genomic information is available for this penaeid and other closely related taxonomic crustaceans. Consequently, genes, proteins and their respective polymorphisms are poorly known for these species. In this work, we used the RNA sequencing technology (RNA-seq) in L. vannamei shrimp evaluated for growth performance, and exposed to the White Spot Syndrome Virus (WSSV), in order to investigate the presence of Single Nucleotide Polymorphisms (SNPs) within genes related to innate immunity and growth, both features of great interest for aquaculture activity. We analyzed individuals with higher and lower growth rates; and infected (unhealthy) and non-infected (healthy), after exposure to WSSV. Approximately 7,000 SNPs were detected in the samples evaluated for growth, being 3,186 and 3,978 exclusive for individuals with higher and lower growth rates, respectively. In the animals exposed to WSSV we found about 16,300 unique SNPs, in which 9,338 were specific to non-infected shrimp, and 7,008 were exclusive to individuals infected with WSSV and symptomatic. In total, we describe 4,312 unigenes containing SNPs. About 60% of these unigenes returned GO blastX hits for Biological Process, Molecular Function and Cellular Component ontologies. We identified 512 KEGG unique KOs distributed among 275 pathways, elucidating the majority of metabolism roles related to high protein metabolism, growth and immunity. These polymorphisms are all located in coding regions, and certainly can be applied in further studies involving phenotype expression of complex traits, such as growth and immunity. Overall, the set of variants raised herein enriches the genomic databases available for shrimp, given that SNPs originated from nextgen are still rare for this relevant crustacean group, despite their huge potential of use in genomic selection approaches.