Gene set based association analyses for the WSSV resistance of Pacific white shrimp Litopenaeus vannamei

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Circulating Phylotypes of White Spot Syndrome Virus in Bangladesh and Their Virulence

White Spot Syndrome Virus (WSSV) has emerged as one of the most prevalent and lethal viruses globally and infects both shrimps and crabs in the aquatic environment. This study aimed to investigate the occurrence of WSSV in different ghers of Bangladesh and the virulence of the circulating phylotypes. We collected 360 shrimp (Penaeus monodon) and 120 crab (Scylla sp.) samples from the south-east (Cox’s Bazar) and south-west (Satkhira) coastal regions of Bangladesh. The VP28 gene-specific PCR assays and sequencing revealed statistically significant (p < 0.05, Kruskal–Wallis test) differences in the prevalence of WSSV in shrimps and crabs between the study areas (Cox’s Bazar and Satkhira) and over the study periods (2017–2019). The mean Log load of WSSV varied from 8.40 (Cox’s Bazar) to 10.48 (Satkhira) per gram of tissue. The mean values for salinity, dissolved oxygen, temperature and pH were 14.71 ± 0.76 ppt, 3.7 ± 0.1 ppm, 34.11 ± 0.38 °C and 8.23 ± 0.38, respectively, in the WSSV-positive ghers. The VP28 gene-based phylogenetic analysis showed an amino-acid substitution (E→G) at the 167th position in the isolates from Cox’s Bazar (referred to as phylotype BD2) compared to the globally circulating one (BD1). Shrimp PL artificially challenged with BD1 and BD2 phylotypes with filtrates of tissue containing 0.423 × 109 copies of WSSV per mL resulted in a median LT50 value of 73 h and 75 h, respectively. The in vivo trial showed higher mean Log WSSV copies (6.47 ± 2.07 per mg tissue) in BD1-challenged shrimp PL compared to BD2 (4.75 ± 0.35 per mg tissue). Crabs infected with BD1 and BD2 showed 100% mortality within 48 h and 62 h of challenge, respectively, with mean Log WSSV copies of 12.06 ± 0.48 and 9.95 ± 0.37 per gram tissue, respectively. Moreover, shrimp antimicrobial peptides (AMPs), penaeidin and lysozyme expression were lower in the BD1-challenged group compared to BD2 challenged shrimps. These results collectively demonstrated that relative virulence properties of WSSV based on mortality rate, viral load and expression of host immune genes in artificially infected shrimp PL could be affected by single aa substitution in VP28.

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.

Transcriptome analysis of pacific white shrimp (Penaeus vannamei) intestines and hepatopancreas in response to Enterocytozoon hepatopenaei (EHP) infection

Penaeus vannamei is the most economically important species of shrimp cultured worldwide. Enterocytozoon hepatopenaei (EHP) is an emerging pathogen that severely affects the growth and development of shrimps. In this study, the transcriptome differences between EHP-infected and uninfected shrimp were investigated through next-generation sequencing. The unigenes were assembled with the reads from all the four libraries. The differentially expressed genes (DEGs) of intestines and hepatopancreas were analyzed. There were 2,884 DEGs in the intestines and 2,096 DEGs in the hepatopancreas. The GO and KEGG enrichment analysis indicated that DEGs were significantly enriched in signaling pathways associated with nutritional energy metabolism and mobilizing autoimmunity. Moreover, the results suggested the downregulation of key genes in energy synthesis pathways contributed greatly to shrimp growth retardation; the upregulation of immune-related genes enhanced the resistance of shrimp against EHP infection. This study provided identified genes and pathways associated with EHP infection revealing the molecular mechanisms of growth retardation.

Association between ALFPm3 single nucleotide polymorphism and white spot syndrome virus resistance in black tiger shrimp Penaeus monodon

Here single nucleotide polymorphisms (SNPs) were associated with white spot syndrome virus (WSSV) resistance in black tiger shrimp Penaeus monodon. SNPs were identified by single-strand conformation polymorphism (SSCP) screening and DNA sequencing of shrimp sampled from 3 families (100 shrimp per family) challenged with WSSV. Shrimp that died over the 14 d challenge trial were designated susceptible, with those remaining alive on Day 14 designated resistant. To compare SNPs, 10 samples from the susceptible and resistant groups, each comprising DNA pooled from 3 shrimp, were amplified by polymerase chain reaction (PCR) using primers to 12 selected genes and screened by SSCP. SNPs were only identified in the anti-lipopolysaccharide factor 3 (ALFPm3) gene product. Analysis of complete ALFPm3 gene sequences confirmed the existence of 3 SNPs (g.934C>G, g.1186A>G, and g.1898C>G) that were polymorphic between the susceptible and resistant groups. Further analyses using specific tetra-primer amplification refractory mutation system PCR primer sets associated these 3 SNPS, and particularly the g.1186A>G SNP, with WSSV resistance. This SNP thus has potential for use as a DNA marker to select for WSSV resistance in P. monodon breeding programs.

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.

Genomic selection for white spot syndrome virus resistance in whiteleg shrimp boosts survival under an experimental challenge test

White spot syndrome virus (WSSV) causes major worldwide losses in shrimp aquaculture. The development of resistant shrimp populations is an attractive option for management of the disease. However, heritability for WSSV resistance is generally low and genetic improvement by conventional selection has been slow. This study was designed to determine the power and accuracy of genomic selection to improve WSSV resistance in Litopenaeus vannamei. Shrimp were experimentally challenged with WSSV and resistance was evaluated as dead or alive (DOA) 23 days after infestation. All shrimp in the challenge test were genotyped for 18,643 single nucleotide polymorphisms. Breeding candidates (G₀) were ranked on genomic breeding values for WSSV resistance. Two G₁ populations were produced, one from G₀ breeders with high and the other with low estimated breeding values. A third population was produced from “random” mating of parent stock. The average survival was 25% in the low, 38% in the random and 51% in the high-genomic breeding value groups. Genomic heritability for DOA (0.41 in G₁) was high for this type of trait. The realised genetic gain and high heritability clearly demonstrates large potential for further genetic improvement of WSSV resistance in the evaluated L. vannamei population using genomic selection.

A SNP of the hemocyanin gene (LvHc) is a marker for high growth and ammonia-tolerance in Pacific white shrimp Litopenaeus vannamei

Expression levels of hemocyanin (LvHc), activating transcription factor 4 (LvAtf4), glutathione S-transferase (LvGst), caspase 2 (LvCasp2) and anti-lipopolysaccharide factor (LvAlf) were examined in the hepatopancreas of Pacific white shrimp Litopenaeus vannamei juveniles exposed to a lethal concentration of ammonia-N (32.15 mg/l). The expression levels of all transcripts except LvAlf were significantly greater (P < 0.05) in tolerant shrimp (Lv-AT; N = 30) that survived up to 72 h post treatment (hpt) than in susceptible shrimp (Lv-AS24 and Lv-AS72; N = 45 and 15), that died within 24 h or between 24 and 72 hpt, respectively. Subsequently, effects of non-lethal concentrations of ammonia-N (control, 10 and 20 mg/l) on the expression of LvHc in juvenile shrimp were examined. Compared to the control, expression levels of LvHc transcripts in hemocytes and the hepatopancreas of tested shrimp changed after exposure to ammonia-N. One SNP (C > T₅₄₅) was found in the LvHc₃₂₂ gene segment. Real-time PCR amplification of specific alleles (real-time PASA) was developed for detection of C > T₅₄₅ genotypes. Juveniles in the lethal exposure test that carried a C/T₅₄₅ genotype showed a greater average body weight and total length (8.46 ± 0.36 g and 10.05 ± 0.16 cm) than those with a C/C₅₄₅ genotype (7.48 ± 0.31 g and 9.60 ± 0.13 cm) (P < 0.05). Similar results were found in the second generation (G2) of a growth-improved stock (3 and 4 families of BIOTEC-G2-L1 and BIOTEC-G2-L2) and in commercially farmed shrimp (2 groups). Accordingly, expression levels and SNP of LvHc can serve as markers for selection high growth performance in ammonia-tolerant L. vannamei.

The krüppel-like factor of Penaeus vannamei negatively regulates transcription of the small subunit hemocyanin gene as part of shrimp immune response

Hemocyanin is a multifunctional respiratory glycoprotein, which has also been implicated in other biological functions in shrimp. Moreover, recent studies have revealed that hemocyanin is also involved in a broad range of immune-related activities in shrimp. However, in spite of the considerable interest in unraveling the reasons behind the multiple immune-related functions of hemocyanin, little is known about its transcriptional regulation. Here, DNA pull-down and Liquid Chromatography - Tandem Mass Spectrometry (LC-MS/MS) analyses were used to isolate and identify the putative transcription factor(s) that are involved in the transcriptional regulation of the small subunit hemocyanin gene of Penaeus vannamei (PvHMCs). Krüppel-like factor (designated PvKruppel), a zinc finger transcription factor homolog in P. vannamei, was identified among the putative transcription factors, while bioinformatics analysis revealed the presence of Krüppel-like factor binding site (KLF motif) on the core promoter region of PvHMCs. Mutational analysis and electrophoretic mobility shift assay (EMSA) confirmed that PvKruppel could bind to the KLF motif on the core promoter region of PvHMCs. Moreover, in response to lipopolysaccharide (LPS), Vibrio parahaemolyticus and white spot syndrome virus (WSSV) challenge, transcript levels of PvKruppel and PvHMCs were negatively correlated. Furthermore, overexpression of PvKruppel significantly reduced the promoter activity of PvHMCs, while PvKruppel knockdown by RNA interference or lipopolysaccharides (LPS) stimulation resulted in a significant increase in the transcript level of PvHMCs. Taken together, our present study provides mechanistic insights into the transcriptional regulation of PvHMCs by PvKruppel during shrimp immune response to pathogens.

Identification of Functional Gene Modules Associated With STAT-Mediated Antiviral Responses to White Spot Syndrome Virus in Shrimp

White spot syndrome virus (WSSV) is one of the major threats to shrimp aquaculture. It has been found that the signal transducer and activator of transcription (STAT) protein plays an important role in the antiviral immunity of shrimp with a WSSV infection. However, the mechanism that underlies the STAT-mediated antiviral responses in shrimp, against WSSV infection, remains unclear. In this work, based on the gene expression profiles of shrimp with an injection of WSSV and STAT double strand RNA (dsRNA), we constructed a gene co-expression network for shrimp and identified the gene modules that are possibly responsible for STAT-mediated antiviral responses. These gene modules are found enriched in the regulation of the viral process, JAK-STAT cascade and the regulation of immune effector process pathways. The gene modules identified here provide insights into the molecular mechanism that underlies the STAT-mediated antiviral response of shrimp, against WSSV.

WSV181 inhibits JAK/STAT signaling and promotes viral replication in Drosophila

The Janus kinase/signal transducers and activators of transcription (JAK/STAT) pathway plays a critical role in host defense against viral infections. Here, we report the use of the Drosophila model system to investigate the modulation of the JAK/STAT pathway by the white spot syndrome virus (WSSV) protein WSV181. WSV181 overexpression in transgenic flies resulted in the downregulation of STAT92E and STAT92E-targeted genes. This result indicates that WSV181 can suppress JAK/STAT signaling by controlling STAT92E expression. An infection experiment was carried out on transgenic Drosophila infected with Drosophila C virus and on Litopenaeus vannamei injected with recombinant WSV181 and WSSV. The increased viral load and suppressed transcript levels of JAK/STAT pathway components indicate that WSV181 can promote viral proliferation by inhibiting the JAK/STAT pathway. This study provided evidence for the role of WSV181 in viral replication and revealed a new mechanism through which WSSV evades host immunity to maintain persistent infection.

Identification and immune-related analysis of SNPs in Litopenaeus vannamei Toll3 receptor

Tolls and Toll-like receptors (TLRs), as innate immune-recognition receptors that recognize molecular patterns associated with microbial pathogens, play a critical role in antimicrobial immune responses. Here, we report on single nucleotide polymorphisms (SNPs) of Litopenaeus vannamei Toll3 (LvToll3). Multiple sequence alignment of the L. vannamei Toll3 Leucine rich repeat C-terminal domain (LvToll3-LRR-CT) with other L. vannamei Tolls LRR-CT domains showed 39.23% - 43.96% homology at the nucleic acid level and 20.31% - 30.00% identity at the amino acid level. Analysis of different shrimp tissues by polymerase chain reaction denaturing gradient gel electrophoresis (PCR-DGGE) revealed that LvToll3-LRR-CT had genetic polymorphisms at both the genomic deoxyribonucleic acid (gDNA) and complementary deoxyribonucleic acid (cDNA) levels. Further, high-throughput sequencing analysis confirmed the presence of 8 non-synonymous SNP (nsSNP) and 1 nsSNPs with frequency greater than 1% at the gDNA level, while 13 nsSNPs and 2 nsSNPs with frequency greater than 1% at the cDNA level. In silico analysis revealed that the α-helix secondary structure and tertiary structure of LvToll3 changed when 3 SNPs (C2039T, T2041C, T2228C) were mutated. Interestingly, 2 novel bands on PCR-DGGE, which were identified as 2 nsSNPs (C2140A, T2186A) were observed following challenge with Streptococcus iniae but not with Vibrio parahaemolyticus or White spot syndrome virus (WSSV). Moreover, the secondary and tertiary structures of LvToll3 changed when the nsSNP T2186A was mutated. The present findings therefore provide novel insight into the molecular basis of shrimp innate immune response to pathogens through the generation of specific SNPs.

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.