Qian Zhang1,2,3†
Yang Yu1,2†
Quanchao Wang1
Fei Liu1Zheng Luo1,3Chengsong Zhang1,2
Xiaojun Zhang1,2Hao Huang4
Jianhai Xiang1,2
Fuhua Li1,2,5*- 1Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
- 2Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao, China
- 3University of Chinese Academy of Sciences, Beijing, China
- 4Hainan Grand Suntop Ocean Breeding Co., Ltd., Wenchang, China
- 5Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao, China
Acute hepatopancreatic necrosis disease (AHPND) is a major bacterial disease in Pacific white shrimp Litopenaeus vannamei farming, which is caused by Vibrio parahaemolyticus. AHPND has led to a significant reduction of shrimp output since its outbreak. Selective breeding of disease-resistant broodstock is regarded as a key strategy in solving the disease problem. Understanding the relationship between genetic variance and AHPND resistance is the basis for marker-assisted selection in shrimp. The purpose of this study was to identify single nucleotide polymorphisms (SNPs) associated with the resistance against AHPND in L. vannamei. In this work, two independent populations were used for V. parahaemolyticus challenge and the resistant or susceptible shrimp were evaluated according to the survival time after Vibrio infection. The above two populations were genotyped separately by a SNP panel designed based on the target sequencing platform using a pooling strategy. The SNP panel contained 508 amplicons from DNA fragments distributed evenly along the genome and some immune-related genes of L. vannamei. By analyzing the allele frequency in the resistant and susceptible groups, 30 SNPs were found to be significantly associated with the resistance of the shrimp against V. parahaemolyticus infection (false discovery rate corrected at P < 0.05). Three SNPs were further validated by individual genotyping in all samples of population 1. Our study illustrated that target sequencing and pooling sequencing were effective in identifying the markers associated with economic traits, and the SNPs identified in this study could be used as molecular markers for breeding disease-resistant shrimp.
Introduction
The Pacific white shrimp Litopenaeus vannamei, which naturally distributed along the Pacific coasts of Central and South America, has become the primary cultivated shrimp species into many regions of the world (Lu et al., 2017). Since it was introduced in China into the 1980s, it has become the most dominant aquaculture shrimp species in China (Lu et al., 2016; Wang et al., 2017). In recent years, due to intensive culture and environmental deterioration, infectious diseases caused by viruses and bacteria have led to serious production loss (Wang et al., 2015). Acute hepatopancreatic necrosis disease (AHPND), also called early mortality syndrome, is a devastating disease that usually occurs in 35 days after stocking in cultivation ponds, and the mortality reached as high as 40% to 100% (Joshi et al., 2014; Lun et al., 2018). During disease outbreak, the shrimp showed clinical signs of atrophied hepatopancreas and empty midgut (Nunan et al., 2014). It was reported that Vibrio parahaemolyticus carrying the Photorhabdus insect-related (Pirvp) binary toxin was the pathogenic agent for AHPND (Lee et al., 2015).
Many immune-related genes were found to be involved in the defense of shrimp against V. parahaemolyticus infection. The Rho signaling pathway was suggested to be helpful for AHPND pathogenesis in shrimp through transcriptome analysis (Ng et al., 2018). Infection of V. parahaemolyticus could strongly activate the genes involved in cell growth and anti-apoptosis (Zheng et al., 2018). In addition, the expression of genes encoding immune effectors, including lectins and antimicrobial peptides (AMPs), was significantly up-regulated after V. parahaemolyticus challenge (Soonthornchai et al., 2016; Qi et al., 2017; Maralit et al., 2018; Qin et al., 2018). Although investigation on the mechanism of AHPND outbreak in shrimp has been conducted, there is still no effective strategy to control this disease.
Breeding disease-resistant broodstock was regarded as an efficient approach (Moss et al., 2012). Compared to the traditional breeding technique, marker-assisted selection (MAS) is more efficient, time-saving, and free from environmental impact (Lu et al., 2018; Wang et al., 2017). Previous studies indicated that gene polymorphism was closely associated with economic traits (Wang et al., 2016a; Li et al., 2017). Immune-related genes were usually considered as optimal candidates for the selection of markers associated with resistance to pathogens (Guo et al., 2013; Hao et al., 2015). Single nucleotide polymorphisms (SNPs) are a type of widely used markers for MAS due to its high polymorphism and abundance in the genome. It has already been used in aquatic species to facilitate selection breeding and speed up the discovery of genes related to the economic trait, such as high growth (Lv et al., 2015; Tsai et al., 2015), disease resistance (Yue et al., 2012), precocity (Xu et al., 2016), and body conformation (Geng et al., 2017).
Over the past years, molecular genetic studies on the penaeid shrimp L. vannamei have made a great progress, such as identification of high-throughput SNPs from the transcriptome (Yu et al., 2014) and construction of high-density linkage map (Yu et al., 2015; Jones et al., 2017). These studies provided fundamental information to identify quantitative trait genes and quantitative trait nucleotides (QTNs) for the AHPND-resistant trait. QTNs have been identified in L. vannamei for several traits, such as growth-related trait (Andriantahina et al., 2013), ammonia-tolerant trait (Lu et al., 2018), and WSSV resistance trait (Liu et al., 2014a). Screening for SNPs or genes associated with AHPND is the basic work to reveal the molecular mechanism of shrimp defense for bacteria infection and it will also be of benefit to MAS or gene-assisted selection in shrimp. However, markers related to the AHPND-resistant trait of shrimp are seldom reported until now. In the present study, genetic association analysis using genotyping by target sequencing method and bulked segregant analysis was performed to screen for disease resistance-related SNPs. The bulks were built by selection of extreme samples from populations (Zou et al., 2016). Consequently, several SNPs or genes were identified to be related to the resistance of the shrimp against V. parahaemolyticus. The results would provide new insights into the molecular basis of disease defense and accelerate the breeding of shrimp with disease resistance.
Materials and Methods
Experimental Shrimp
Healthy shrimp were obtained from Hainan Grand Suntop Ocean Breeding Co., Ltd. (Wenchang, China) and maintained in filtered seawater at a temperature of 26 ± 1°C with continuous aeration. Before the experiment, the shrimp were kept in the aquarium for 3 days to make them acclimate to the laboratory conditions. The shrimp used in this study were the progeny of mating of two breeding lines. The two breeding lines were bred by 7 years of artificial selection, respectively. The progenies produced by the cross of the breeding lines were genetically similar. Although two populations were consisted of multiple families, their genetic backgrounds were very similar.
V. parahaemolyticus Challenge Tests
V. parahaemolyticus was isolated from the diseased shrimp. The Vibrio strain was cultured in tryptic soy broth with 2.0% sodium chloride liquid media. The cultured bacteria were further checked by positive polymerase chain reaction (PCR) amplification of PirAvp and PirBvp genes (Han et al., 2015). Bacterial titer was counted with a hemocytometer under a light microscope. The infection dose was set at 2.5 × 104 colony-forming units (CFU)/g (body weight). A total of 236 shrimp with a body weight of 3.98 ± 1.9 g and a body length of 6.89 ± 1.35 cm were used as the first population. About 10 μl bacterial suspension in phosphate-buffered saline (PBS) with 1.0 × 105 CFU of V. parahaemolyticus was injected into each shrimp at the site between the IV and V abdominal segments. A total of 270 shrimp with a body weight of 5.25 ± 2.9 g and a body length of 7.45 ± 1.56 cm were used as the second population, and each shrimp was injected with 1.3 × 105 CFU of pathogenic bacteria. At the same time, 60 shrimp were set aside as the control group and injected with an equal volume of sterile PBS. The mortality was recorded for 8 days. The dead shrimp were stored at −80°C until DNA extraction. Susceptible and resistant phenotypes were evaluated according to the survival time after challenge (Sawayama et al., 2017). In each population, 60 samples that died at an early time [16 h postinfection (hpi)] after injection of Vibrio were taken as the highly susceptible group, and 60 samples that survived from the infection or died at a later time [6 days postinfection (dpi)] were selected as the resistant group.
Loads of V. parahaemolyticus in Vivo
To explore the dynamic change of V. parahaemolyticus in shrimp, 150 healthy shrimp with a body weight of 5.15 ± 1.6 g were injected with 2.5 × 104 CFU/g (body weight). To determine the pathogenic bacteria count in surviving shrimp, six shrimp were randomly chosen from the tank at 0, 1, 2, 3, 8, and 12 dpi. Various tissues were collected from each shrimp, including hepatopancreas, gill, and muscle. Three repetitions were used each time, and two shrimp were mixed together in each repetition. Tissues were dissected aseptically and then homogenized in PBS. Serial 10-fold dilutions of the supernatant with PBS were produced, and 100-μl solution was inoculated onto thiosulfate citrate bile salts sucrose (TCBS) agar. Plates were incubated at 30°C for 18 h. Colonies were counted to estimate the number of viable bacteria per pool of shrimp tissue (Sotorodriguez et al., 2015).
SNP Genotyping Based on Next-Generation Target Sequencing
DNA was extracted from the muscle of shrimp with the Tiangen Plant Genomic DNA Extraction Kit (Tiangen, Beijing, China) according to the manufacturer’s instructions. DNA concentration was assessed using a NanoDrop 2000 spectrophotometer (Thermo Scientific, USA). A SNP genotyping panel based on the AmpliSeq™ method was designed. This panel contained 508 amplicons including 60 immune-related genes with 148 amplicons, and the other 360 amplicons were selected with even distribution along the shrimp genome (almost 10 cM per marker along the genome). The primer sequences of the amplicons are displayed in Supplemental Table S1.
Each population was divided into two DNA pools from shrimp with susceptibility and two DNA pools with resistance against V. parahaemolyticus. Every DNA pool was generated by combining an equal amount of DNA from 30 susceptible or resistant shrimp. These DNA pools were amplified and next-generation sequencing libraries were constructed using DNA Seq Library Preparation Kit for Amplicon Sequencing-Illumina Compatible (Gnomegen Technologies, San Diego, CA, USA). The amplicons were sequenced using Illumina Hiseq 2500. The raw reads were filtered and then mapped to the target sequence using BWA (version 0.7.12; Li, 2011). SNPs were genotyped using the Genome Analysis Toolkit. The read depth of each allele in the pools was extracted from the VCF file and used to estimate the allele frequency in each pool.
Association Analysis and Candidate Gene Identification
The difference in allele frequency between the susceptible and resistant groups was analyzed using χ2 test. The analysis was performed in the combined population (population 1 + population 2) and P values were corrected using the false discovery rate (FDR) according to the method described by Benjamini and Hochberg (1995). The threshold of adjusted P was set to 0.05 (5% FDR correction). In view of the consistency of favorable alleles in the two populations, SNPs that showed opposite allele frequencies in two populations were filtered out. Statistical tests were performed using the R statistical software (R Core Team, 2018). Screened SNPs were mapped to the genome of L. vannamei (Zhang et al., 2019), and the genes around SNPs were examined for candidate genes according to their locations and functions.
Validation of the Pooling Genotyping Result
To further validate the pooling genotyping result, three SNPs (ALF6-1__22575_510_57, ALF6-1__22575_510__224, and Marker2060_197) were amplified by individual genotyping method in all samples of population 1. Primers were designed based on the flanking sequence of SNPs (Table 1). The PCR program was as follows: 1 cycle of denaturation at 94°C for 5 min, 35 cycles of denaturation at 94°C for 30 s, annealing at 58°C for 30 s, and extension at 72°C for 45 s followed by an extension at 72°C for 7 min. The confirmed PCR products were sequenced by Sanger sequencing in Tsingke Biotech (Qingdao, China). The genotype of each sample was determined based on the sequencing chromatograms.
Table 1 Primer sequences and corresponding annealing temperature of genes.
Transcription Analysis on Two Candidate Genes
Based on the annotation of SNPs located in the genes, two candidate genes, including anti-lipopolysaccharide factor (LvALF6) and phosphatidylinositol 3-kinase (PI3K) regulatory subunit α isoform X2 (LvPI3K), were selected for further analysis. Their expression patterns in the shrimp during V. parahaemolyticus infection were analyzed. Healthy shrimp with a body weight of 8.39 ± 2.2 g and a body length of 8.89 ± 0.78 cm were injected with 1 × 105 CFU of V. parahaemolyticus. The hepatopancreas, hemocytes, and gill of shrimp were collected separately at 1, 6, 24, 48, and 72 hpi. Nine shrimp were collected at each time point and three individuals were put together as one sample. The detailed steps of total RNA extraction and cDNA synthesis were the same as described previously by Wang et al. (2016b). The expression patterns in different tissues were detected by SYBR Green-based quantitative real-time PCR with the primers shown in Table 1. The program was as follows: 95°C for 2 min followed by 40 cycles of 95°C for 15 s, 56°C for 15 s, and 72°C for 20 s, and the melting curve analysis was added to the end of each PCR.
Results
AHPND Challenge Test
The cumulative mortality of two populations was 52.12% and 86.29%, respectively (Figure 1). In population 1, 60 samples that died earlier were selected as the susceptible group, whereas 60 surviving samples were selected as the resistant group. Similarly, 60 shrimp that died first were collected from population 2, whereas 37 surviving shrimp and 23 shrimp that died later were selected as the resistant group. No mortality was observed in the PBS group. The loads of V. parahaemolyticus in shrimp at different infection stages are shown in Table 2. No V. parahaemolyticus was detected in the uninfected shrimp. The number of bacteria in shrimp tissues increased gradually after injection. In hepatopancreas and muscle, the bacterial loads reached 106 and 105 CFU/g at 2 dpi, respectively. The density of bacteria in the gill reached the peak with a bacterial load of 103 CFU/g. At 2 dpi, the bacterial loads decreased slightly. Afterward, it decreased continually in the surviving shrimp and reached 102 to 103 CFU/g in the end.
Figure 1 Accumulated mortality rates of populations 1 and 2 infected with V. parahaemolyticus and PBS. PL1, population 1; PL2, population 2.
Table 2 Loads of V. parahaemolyticus in hepatopancreas, gill, and muscle of shrimp.
SNPs Associated With AHPND Resistance
A total of 1566 SNPs were kept after filtering out SNPs with a minor allele frequency of less than 0.05. The raw data were submitted to the SRA database with accession numbers from SRR9016244 to SRR9016251. Through association analysis on the two populations, 40 significantly different SNPs were identified. After filtering out 10 false-positive sites that showed opposite favorable alleles in two separate populations (Figure 2), 30 SNPs were retained as significant markers (Table 3). The Manhattan plot of -log10(original P) and the Q–Q plot were supplied in the supplementary file (Supplemental Figures S1 and S2). Among these 30 markers, 4 SNPs were located in previously selected immune genes and 26 SNPs were located in the fragment distributed along the genome. As several SNPs were located in the same fragment, a total of 17 independent fragments were screened out to be related with the resistance. Figure 3 shows the -log10(adjusted P) of three SNPs with the most significant difference, including loci Marker15416_294, Marker8720_486, and Marker1077_61 located in the fragment of Marker15416, Marker8720, and Marker1077, respectively.
Figure 2 -log10(adjusted P) of 1556 SNPs (filtering out 10 false-positive sites) in a combined population and the threshold for significance is FDR-corrected P < 0.05.
Table 3 χ2 test and allele frequency distribution of SNP markers with significant difference (FDR-corrected P < 0.05) in two populations.
Figure 3 -log10(adjusted P) of three SNPs, including loci Marker15416_294, Marker8720_486, and Marker1077_61 located in the fragment of Marker15416, (A) Marker8720 (B), and Marker1077 (C), respectively.
To further validate the association result, three SNPs were genotyped by Sanger sequencing in all individuals of population 1. The genotype at the SNP site was determined based on sequencing chromatograms, which is shown in Supplemental Figure S3. As for ALF6-1__22575_510_57, 219 samples of population 1 were successfully genotyped. The number of each genotype was 65, 58, and 96 for the A/A, A/T, and T/T genotypes, respectively (Table 4). The A/A genotype accounted for 37.5% in the susceptible group, which was higher than that in the resistant group (20.2%). The T/T genotype accounted for 39.2% in the susceptible group, which was lower than that in the resistant group (49.5%). χ2 test showed that the shrimp carrying allele T were more resistant than those carrying allele A (P = 0.0204). Another SNP in the same gene was ALF6-1__22575_510_224, which also showed significant difference in 219 individuals. The percentage of G in the susceptible group was lower than that in the resistant group (P = 0.0293). For Marker2060_197, T was the favorable allele (P = 0.031).
Table 4 Distribution of the three markers in the susceptible and resistant groups of population 1.
Genomic Regions Associated With AHPND
The corresponding flanking sequences of the above 30 SNPs with significant difference were obtained by blasting the SNP sequence with the assembled reference genome (Supplemental Table S2). Among these SNPs, there were only 10 markers located in areas with annotated genes. The surrounding candidate genes are listed in Table 5. SNPs ALF6-1__22575_510_57 and ALF6-1__22575_510_224 were both located in LvALF6. The SNP Unigene19157_All__1806_348_223 was located in gene LvPI3K. In addition, Marker66_168 and Marker66_291 were near the ubiquitin carboxyl-terminal hydrolase. Some other genes, such as low molecular weight phosphotyrosine protein phosphatase-like and calpain-B-like, were also found as candidate genes.
Table 5 Information of associated regions and candidate gene identification.
According to the gene annotation and association P value, we focused on two candidate genes. One was LvALF6, where two SNPs with significant difference were located. The other was LvPI3K, which contained one SNP with significant difference designated as Unigene19157_All__1806_348_223 (Table 3). For SNP Unigene19157_All__1806_348_223, the average frequency of C allele was 36% in the susceptible group and 20% in the resistant group. Thus, C allele was identified as a deleterious allele for disease resistance. We further analyzed the expression profiles of the two candidate genes after Vibrio challenge (Figure 4). The expression level of LvALF6 in hepatopancreas and gill was significantly up-regulated after infection (Figure 4A, B), and it was up-regulated at 6 hpi and then down-regulated at the late infection stage in hemocytes (Figure 4C). The expression level of LvPI3K regulatory subunit in hepatopancreas was significantly down-regulated at early infection stage (Figure 4D), whereas it was obviously up-regulated in gill at 48 and 72 hpi (Figure 4E). LvPI3K showed no expression difference in hemocytes after infection (Figure 4F).
Figure 4 Expression profile of LvALF6 and LvPI3K regulatory subunit in hepatopancreas, gill, and hemocytes at different times after V. parahaemolyticus infection (P < 0.01). Relative transcription level of LvALF in hepatopancreas (A), gill (B) and hemocytes (C); Relative transcription level of LvPI3K in hepatopancreas (D), gill (E) and hemocytes (F).
Discussion
In recent years, AHPND has been considered as one of the main problems that hinders the rapid development of the shrimp aquaculture industry. However, no effective method has been reported to control this disease until now. It is widely believed that the selection for a disease-resistant host is one of the valid control measures to prevent disease outbreak in shrimp farming (Cock et al., 2009). A good example is selective breeding of Taura syndrome virus-resistant lines in L. vannamei (Argue et al., 2002; Moss et al., 2012). Recently, Ge et al. (2018) reported a three-generation selective breeding of ridgetail prawn (Exopalaemon carinicauda) to improve the resistance of AHPND infection. However, no research about the selection of AHPND-resistant broodstock was reported for L. vannamei until now. Currently, MAS was known as a method that can greatly accelerate the process of breeding. In red sea bream, it suggested that the RSIVD-resistant trait was controlled by one major quantitative trait loci and could be useful for MAS (Sawayama et al., 2017). Seven SNP markers were identified as markers for the selection of V. parahaemolyticus infection resistance in clam (Meretrix meretrix; Nie et al., 2015). Compared to other species, there are few studies about molecular markers associated with AHPND in shrimp. Thus, screening of SNP markers associated with disease resistance is the first step to MAS.
In this study, we carried out the V. parahaemolyticus challenge experiment to obtain AHPND-resistant and AHPND-susceptible populations. Most of the mortality occurred within the first 48 hpi, whereas there was no mortality in the control shrimp. By investigating the dynamic change of pathogenic bacteria in vivo, healthy shrimp were confirmed to be V. parahaemolyticus free. The amount of V. parahaemolyticus in infected shrimp was higher than that in surviving shrimp, and V. parahaemolyticus density was rapidly increased at 1 or 2 dpi, which was in accordance with the highest mortality rate during this time. In the terminal stage of the disease, a decrease of V. parahaemolyticus counts was observed due to the immune system mechanisms (Khimmakthong and Sukkarun, 2017).
The approaches of discovering SNPs include whole genome sequencing, large-scale amplicon sequencing, transcriptome sequencing, gene-enriched genome sequencing, and so on (Henry et al., 2012). Here, SNP screening and genotyping were performed using a targeting sequencing panel. It contained genome-wide distributed fragments and immune-related genes, which include immune effective factors or immune signaling pathways such as the JAK/STAT, Toll, and IMD pathways. The target sequencing approach based on multiplex PCR is an efficient approach to analyze genetic variation in specific genomic regions. It is cost-effective, flexible, high throughput, and suitable for non-model species. We applied this method for genotyping genomic regions of interest using pooling genotyping strategy in L. vannamei. Subsequently, individual genotyping result proved that pooling sequencing using the designed amplicon panel displayed high sensitivity in detecting allele frequency.
To reduce false-positives, two populations were designed to provide validation for each other. By comparing the allele frequency between the two populations, some SNPs that showed opposite favorable alleles in two separate populations were removed. Among 30 identified associated markers, 14 SNPs were located in LG33. Totally, eight, four, and two SNPs with significant difference were found in the fragment of Marker1077, Marker4976, and Marker66, respectively. We speculate that most of these markers were closely linked. However, these regions have no complete genome annotation. It may result from two reasons: one is that these markers might be closely linked to a candidate gene nearby and the other is that the genome assembly might not be complete; therefore, no candidate genes were identified in this region. With the improvement of the genome assembly, these genomic regions deserve further analysis. Totally, 10 markers were located in the gene region. Several SNPs showing the most significant difference, including loci Marker15416_294, Marker8720_486, and Marker1077_61, were not annotated yet. However, this does not affect the usefulness of SNPs as molecular markers for breeding disease-resistant broodstock.
Among these annotated genes, ALF is a kind of AMP with broad-spectrum activities against bacterial pathogens. A number of ALFs were identified and characterized in E. carinicauda (Lv et al., 2017) and Chinese shrimp Fenneropenaeus chinensis (Li et al., 2015). Different ALFs exhibited diverse antibacterial and antiviral activities in E. carinicauda (Lv et al., 2018) and F. chinensis (Yang et al., 2015). In our previous studies, mutations in LvALF gene have been regarded as the important immune factor in L. vannamei with WSSV resistance (Liu et al., 2014a; Liu et al., 2014b). Several polymorphisms of PtALF in the swimming crab (Portunus trituberculatus) were reported to be associated with resistance and susceptibility to Vibrio alginolyticus (Li et al., 2013). In this study, the expression level of LvALF6 was obviously up-regulated in the shrimp infected by bacteria compared to those in the control group. Therefore, data suggested that LvALF6 played an important role in the disease resistance of shrimp.
In addition, SNPs in LvPI3K were also identified to be associated with V. parahaemolyticus resistance. The PI3K-Akt pathway is an important intracellular signaling pathway involved in various cellular functions, including anti-apoptosis, protein synthesis, glucose metabolism, and cell cycling (Ruan et al., 2014). Several key molecules in this pathway have been reported to be associated with the invasion of some viruses in shrimp (Su et al., 2014; Zhang et al., 2016). Previous studies have reported that V. alginolyticus challenge induced an up-regulation of LvPI3K expression (Kong et al., 2018), which illustrated that LvPI3K might play an important role in Vibrio infection.
In summary, we analyzed 1566 SNPs distributed along genome and candidate genes and identified 30 SNPs associated with shrimp resistance to AHPND. The results proved that the target sequencing approach was a useful method for genotyping interested genomic regions and the pooling genotyping method was more cost-effective than sequencing each individual separately, whereas the results were accurate and convinced. The method established in this study and the identified SNPs could be used for the selective breeding of AHPND-resistant broodstock in L. vannamei. The candidate disease-resistant genes will provide information in dissecting the resistance response of L. vannamei against V. parahaemolyticus.
Data Availability
All datasets generated for this study are included in the manuscript and the supplementary files.
Author Contributions
QZ and YY conducted the experiment and data processing. JX and FHL conceived and supervised the project. QW contributed to statistical analysis. HH and CZ prepared and cultured the experimental animals. ZL participated in the extraction of genomic DNA. FL contributed to prepare V. parahaemolyticus. QZ, YY, and FHL prepared the manuscript. All authors have read and approved the manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (31830100), the National Key R&D Program of China (2018YFD0901301 and 2018YFD0900103), and China Agriculture Research System-48.
Conflict of Interest Statement
Author HH was employed by company Hainan Grand Suntop Ocean Breeding Co.,Ltd.Wenchang, China. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00700/full#supplementary-material
Supplemental Figure S1 | Manhattan plot of -log10(original P).
Supplemental Figure S2 | Q-Q plot for P values.
Supplemental Figure S3 | Chromatograms of different individuals with different allelic variants at the SNP loci (arrows).
References
Andriantahina, F., Liu, X., Huang, H. (2013). Genetic map construction and quantitative trait locus (QTL) detection of growth-related traits in Litopenaeus vannamei for selective breeding applications. PLoS One 8, e75206. doi: 10.1371/journal.pone.0075206
Argue, B. J., Arce, S. M., Lotz, J. M., Moss, S. M. (2002). Selective breeding of Pacific white shrimp (Litopenaeus vannamei) for growth and resistance to Taura syndrome virus. Aquaculture 204, 447–460. doi: 10.1016/S0044-8486(01)00830-4
Benjamini, Y., Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Statist. Soc. B. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Cock, J., Gitterle, T., Salazar, M., Rye, M. (2009). Breeding for disease resistance of Penaeid shrimps. Aquaculture 286, 1–11. doi: 10.1016/j.aquaculture.2008.09.011
Ge, Q., Li, J., Duan, Y., Li, J., Sun, M., Zhao, F. (2016). Isolation of prawn (Exopalaemon carinicauda) lipopolysaccharide and β-1,3-glucan binding protein gene and its expression in responding to bacterial and viral infections. J. Ocean U. China 15, 288–296. doi: 10.1007/s11802-016-2741-2
Ge, Q., Li, J., Li, J., Wang, J., Li, Z. (2018). Immune response of Exopalaemon carinicauda infected with an AHPND-causing strain of Vibrio parahaemolyticus. Fish Shellfish Immun. 74, 223–234. doi: 10.1016/j.fsi.2017.12.042
Geng, X., Liu, S., Yuan, Z., Jiang, Y., Zhi, D., Liu, Z. (2017). A genome-wide association study reveals that genes with functions for bone development are associated with body conformation in catfish. Mar. Biotechnol. 19, 570–578. doi: 10.1007/s10126-017-9775-3
Guo, L., Zhao, X., Zhang, Y., Wang, Z., Zhong, M., Li, S., et al. (2013). Evidences of SNPs in the variable region of hemocyanin Ig-like domain in shrimp Litopenaeus vannamei. Fish Shellfish Immun. 35, 1532–1538. doi: 10.1016/j.fsi.2013.08.024
Han, J. E., Tang, K. F. J., Tran, L. H., Lightner, D. V. (2015). Photorhabdus insect-related (Pir) toxin-like genes in a plasmid of Vibrio parahaemolyticus, the causative agent of acute hepatopancreatic necrosis disease (AHPND) of shrimp. Dis. Aquat. Organ. 113, 33. doi: 10.3354/dao02830
Hao, G., Lin, F., Mu, C., Li, R., Yao, J., Yuan, X., et al. (2015). SNP E4-205 C/T in C-type lectin of Portunus trituberculatus is association with susceptibility/resistance to Vibrio alginolyticus challenge. Aquaculture 442, 125–131. doi: 10.1016/j.aquaculture.2015.02.007
Henry, R. J., Edwards, M., Waters, D. L., Gopala, K. S., Bundock, P., Sexton, T. R., et al. (2012). Application of large-scale sequencing to marker discovery in plants. J. Biosci. 37, 829–841. doi: 10.1007/s12038-012-9253-z
Jones, D. B., Jerry, D. R., Khatkar, M. S., Raadsma, H. W., van der Steen, H., Prochaska, J., et al. (2017). A comparative integrated gene-based linkage and locus ordering by linkage disequilibrium map for the Pacific white shrimp, Litopenaeus vannamei. Sci. Rep. 7, 16. doi: 10.1038/s41598-017-10515-7
Joshi, J., Srisala, J., Truong, V. H., Chen, I. T., Nuangsaeng, B., Suthienkul, O., et al. (2014). Variation in Vibrio parahaemolyticus isolates from a single Thai shrimp farm experiencing an outbreak of acute hepatopancreatic necrosis disease (AHPND). Aquaculture 428, 297–302. doi: 10.1016/j.aquaculture.2014.03.030
Khimmakthong, U., Sukkarun, P. (2017). The spread of Vibrio parahaemolyticus in tissues of the Pacific white shrimp Litopenaeus vannamei analyzed by PCR and histopathology. Microb. Pathog. 113, 107. doi: 10.1016/j.micpath.2017.10.028
Kong, J., Wei, W., Liang, Q., Qiao, X., Kang, H., Liu, Y., et al. (2018). Identifying the function of LvPI3K during the pathogenic infection of Litopenaeus vannamei by Vibrio alginolyticus. Fish Shellfish Immun. 76, 355–367. doi: 10.1016/j.fsi.2018.03.016
Lee, C., Chen, I., Yang, Y., Ko, T., Huang, Y., Huang, J., et al. (2015). The opportunistic marine pathogen Vibrio parahaemolyticus becomes virulent by acquiring a plasmid that expresses a deadly toxin. Proc. Natl. Acad. Sci. U. S. A. 112, 10798–10803. doi: 10.1073/pnas.1503129112
Li, H. (2011). A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27, 2987–2993. doi: 10.1093/bioinformatics/btr509
Li, S., Guo, S., Li, F., Xiang, J. (2015). Functional diversity of anti-lipopolysaccharide factor isoforms in shrimp and their characters related to antiviral activity. Mar. Drugs 13, 2602–2616. doi: 10.3390/md13052602
Li, W., Zhu, Z., Chern, M., Yin, J., Yang, C., Ran, L., et al. (2017). A natural allele of a transcription factor in rice confers broad-spectrum blast resistance. Cell 170, 114–126 e115. doi: 10.1016/j.cell.2017.06.008
Li, X., Cui, Z., Liu, Y., Song, C., Shi, G., Wang, C. (2013). Polymorphisms of anti-lipopolysaccharide factors in the swimming crab Portunus trituberculatus and their association with resistance/susceptibility to Vibrio alginolyticus. Fish Shellfish Immun. 34, 1560–1568. doi: 10.5897/AJB11.2017
Liu, J., Yu, Y., Li, F., Zhang, X., Xiang, J. (2014a). A new anti-lipopolysaccharide factor (ALF) gene with its SNP polymorphisms related to WSSV-resistance of Litopenaeus vannamei. Fish Shellfish Immun. 39, 24–33. doi: 10.1016/j.fsi.2014.04.009
Liu, J., Yu, Y., Li, F., Zhang, X., Xiang, J. (2014b). A new ALF from Litopenaeus vannamei and its SNPs related to WSSV resistance. Chin. J. Oceanol. Limnol. 32, 1232–1247. doi: 10.1007/s00343-015-4010-4
Lu, X., Kong, J., Meng, X., Cao, B., Luo, K., Dai, P., et al. (2018). Identification of SNP markers associated with tolerance to ammonia toxicity by selective genotyping from de novo assembled transcriptome in Litopenaeus vannamei. Fish Shellfish Immun. 73, 158–166. doi: 10.1016/j.fsi.2017.12.005
Lu, X., Luan, S., Cao, B., Meng, X., Sui, J., Dai, P., et al. (2017). Estimation of genetic parameters and genotype-by-environment interactions related to acute ammonia stress in Pacific white shrimp (Litopenaeus vannamei) juveniles at two different salinity levels. PLoS One 12, e0173835. doi: 10.1371/journal.pone.0173835
Lu, X., Luan, S., Luo, K., Meng, X., Li, W., Sui, J., et al. (2016). Genetic analysis of the Pacific white shrimp (Litopenaeus vannamei): heterosis and heritability for harvest body weight. Aquac. Res. 47, 3365–3375. doi: 10.1111/are.12820
Lun, J., Liu, D., Liu, T., Zhang, S., Dong, Y., Li, C., et al. (2018). Evaluation of outer membrane protein U (OmpU) as a novel capture target of Vibrio parahaemolyticus and rapid detection of acute hepatopancreatic necrosis disease (AHPND) using PCR combined with immunomagnetic separation. Aquaculture 485, 225–232. doi: 10.1016/j.aquaculture.2017.11.046
Lv, J., Zhang, D., Gao, B., Liu, P., Li, J. (2015). Transcriptome and MassARRAY analysis for identification of transcripts and SNPs for growth traits of the swimming crab Portunus trituberculatus. Gene 566, 229–235. doi: 10.1016/j.gene.2015.04.050
Lv, X., Li, S., Liu, F., Li, F., Xiang, J. (2017). Identification and function analysis of an anti-lipopolysaccharide factor from the ridgetail prawn Exopalaemon carinicauda. Dev. Comp. Immunol. 70, 128–134. doi: 10.1016/j.dci.2017.01.010
Lv, X., Li, S., Zhang, C., Xiang, J., Li, F. (2018). Multiple isoforms of anti-lipopolysaccharide factors and their antimicrobial functions in the ridgetail prawn Exopalaemon carinicauda. Mar. Drugs 16, 13. doi: 10.3390/md16050145
Maralit, B. A., Jaree, P., Boonchuen, P., Tassanakajon, A., Somboonwiwat, K. (2018). Differentially expressed genes in hemocytes of Litopenaeus vannamei challenged with Vibrio parahaemolyticus AHPND (VPAHPND) and VPAHPND toxin. Fish Shellfish Immun. 81, 284–296. doi: 10.1016/j.fsi.2018.06.054
Moss, S. M., Moss, D. R., Arce, S. M., Lightner, D. V., Lotz, J. M. (2012). The role of selective breeding and biosecurity in the prevention of disease in penaeid shrimp aquaculture. J. Invertebr. Pathol. 110, 247–250. doi: 10.1016/j.jip.2012.01.013
Ng, T. H., Lu, C., Lin, S., Chang, C., Tran, L., Chang, W., et al. (2018). The Rho signalling pathway mediates the pathogenicity of AHPND-causing V. parahaemolyticus in shrimp. Cell Microbiol. 20, 1–12. doi: 10.1111/cmi.12849
Nie, Q., Yue, X., Liu, B. (2015). Development of Vibrio spp. infection resistance related SNP markers using multiplex SNaPshot genotyping method in the clam Meretrix meretrix. Fish Shellfish Immun. 43, 469–476. doi: 10.1016/j.fsi.2015.01.030
Nunan, L., Lightner, D., Pantoja, C., Gomez-Jimenez, S. (2014). Detection of acute hepatopancreatic necrosis disease (AHPND) in Mexico. Dis. Aquat. Organ. 111, 81–86. doi: 10.3354/dao02776
Qi, C., Wang, L., Liu, M., Jiang, K., Wang, M., Zhao, W., et al. (2017). Transcriptomic and morphological analyses of Litopenaeus vannamei intestinal barrier in response to Vibrio paraheamolyticus infection reveals immune response signatures and structural disruption. Fish Shellfish Immun. 70, 437–450. doi: 10.1016/j.fsi.2017.09.004
Qin, Z., Babu, V. S., Wan, Q., Zhou, M., Liang, R., Muhammad, A., et al. (2018). Transcriptome analysis of Pacific white shrimp (Litopenaeus vannamei) challenged by Vibrio parahaemolyticus reveals unique immune-related genes. Fish Shellfish Immun. 77, 164–174. doi: 10.1016/j.fsi.2018.03.030
R Core Team, (2018). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. URL https://www.R-project.org/.
Ruan, L., Liu, R., Xu, X., Shi, H. (2014). Molecular cloning and characterization of a threonine/serine protein kinase lvakt from Litopenaeus Vannamei. Chin. J. Oceanol. Limnol. 32, 792–798. doi: 10.1007/s00343-014-3295-z
Sawayama, E., Tanizawa, S., Kitamura, S. I., Nakayama, K., Ohta, K., Ozaki, A., et al. (2017). Identification of quantitative trait loci for resistance to RSIVD in red sea bream (Pagrus major). Mar. Biotechnol. 19, 1–13. doi: 10.1007/s10126-017-9779-z
Soonthornchai, W., Chaiyapechara, S., Klinbunga, S., Thongda, W., Tangphatsornruang, S., Yoocha, T., et al. (2016). Differentially expressed transcripts in stomach of Penaeus monodon in response to AHPND infection. Dev. Comp. Immunol. 65, 53–63. doi: 10.1016/j.dci.2016.06.013
Sotorodriguez, S. A., Gomezgil, B., Lozanoolvera, R., Betancourtlozano, M., Moralescovarrubias, M. S. (2015). Field and experimental evidence of Vibrio parahaemolyticus as the causative agent of acute hepatopancreatic necrosis disease of cultured shrimp (Litopenaeus vannamei) in northwestern Mexico. Appl. Environ. Microb. 81, 1689–1699. doi: 10.1128/AEM.03610-14
Su, M. A., Huang, Y. T., Chen, I., Lee, D. Y., Hsieh, Y. C., Li, C. Y., et al. (2014). A An invertebrate Warburg effect: a shrimp virus achieves successful replication by altering the host metabolome via the PI3K-Akt-mTOR pathway. PLoS Pathog. 10 (6), e1004196. doi: 10.1371/journal.ppat.1004196
Tsai, H. Y., Hamilton, A., Tinch, A. E., Guy, D. R., Gharbi, K., Stear, M. J., et al. (2015). Genome wide association and genomic prediction for growth traits in juvenile farmed Atlantic salmon using a high density SNP array. Bmc Genomics 16, 969. doi: 10.1186/s12864-015-2117-9
Wang, L., Wang, X., Liu, J., Chen, C., Liu, Y., Wang, W. (2015). Rab from the white shrimp Litopenaeus vannamei: characterization and its regulation upon environmental stress. Ecotoxicology 24, 1765. doi: 10.1007/s10646-015-1481-1
Wang, Q., Zhang, X., Yu, Y., Li, F., Xiang, J. (2017). Predictive ability of genomic selection models for breeding value estimation on growth traits of Pacific white shrimp Litopenaeus vannamei. Chin. J. Oceanol. Limnol. 35, 1–9. doi: 10.1007/s00343-017-6038-0
Wang, X., Wang, H., Liu, S., Ferjani, A., Li, J., Yan, J., et al. (2016a). Genetic variation in ZmVPP1 contributes to drought tolerance in maize seedlings. Nat. Genet. 48, 1233–1241. doi: 10.1038/ng.3636
Wang, Z., Li, S., Li, F., Xie, S., Xiang, J. (2016b). Identification and function analysis of a novel vascular endothelial growth factor, LvVEGF3, in the Pacific whiteleg shrimp Litopenaeus vannamei. Dev. Comp. Immunol. 63, 111–120. doi: 10.1016/j.dci.2016.05.020
Xu, Z., Tang, L., Li, Y., Ge, J., Pan, J. (2016). Identification of SNPs in the 5′-flanking region and 3′-UTR of the MIH gene and their association with precocity of the Chinese mitten crab Eriocheir sinensis. Aquac. Res. 47, 992–1000. doi: 10.1111/are.12559
Yang, H., Li, S., Li, F., Lv, X., Xiang, J. (2015). Recombinant expression and functional analysis of an isoform of anti-lipopolysaccharide factors (FcALF5) from Chinese shrimp Fenneropenaeus chinensis. Dev. Comp. Immunol. 53, 47–54. doi: 10.1016/j.dci.2015.06.015
Yu, Y., Wei, J., Zhang, X., Liu, J., Liu, C., Li, F., et al. (2014). SNP discovery in the transcriptome of white Pacific shrimp Litopenaeus vannamei by next generation sequencing. PLoS One 9, e87218. doi: 10.1371/journal.pone.0087218
Yu, Y., Zhang, X., Yuan, J., Li, F., Chen, X., Zhao, Y., et al. (2015). Genome survey and high-density genetic map construction provide genomic and genetic resources for the Pacific white shrimp Litopenaeus vannamei. Sci. Rep. 5, 15612. doi: 10.1038/srep15612
Yue, X., Wang, H., Huang, X., Wang, C., Chai, X., Wang, C., et al. (2012). Single nucleotide polymorphisms in i-type lysozyme gene and their correlation with Vibrio-resistance and growth of clam Meretrix meretrix based on the selected resistance stocks. Fish Shellfish Immun. 33, 559–568. doi: 10.1016/j.fsi.2012.06.007
Zhang, X., Song, X., Huang, J. (2016). Impact of Vibrio parahaemolyticus and white spot syndrome virus (WSSV) co-infection on survival of penaeid shrimp Litopenaeus vannamei. Chin. J. Oceanol. Limnol. 34, 1278–1286. doi: 10.1007/s00343-016-5165-3
Zhang, X., Yuan, J., Sun, Y., Li, S., Gao, Y., Yu, Y., et al. (2019). Penaeid shrimp genome provides insights into benthic adaptation and frequent molting. Nat. Commun. 10, 356. doi: 10.1038/s41467-018-08197-4
Zheng, Z., Wang, F., Aweya, J. J., Li, R., Yao, D., Zhong, M., et al. (2018). Comparative transcriptomic analysis of shrimp hemocytes in response to acute hepatopancreas necrosis disease (AHPND) causing Vibrio parahemolyticus infection. Fish Shellfish Immun. 74, 10–18. doi: 10.1016/j.fsi.2017.12.032
Keywords: single nucleotide polymorphism, target sequencing, association analysis, disease resistance, penaeid shrimp
Citation: Zhang Q, Yu Y, Wang Q, Liu F, Luo Z, Zhang C, Zhang X, Huang H, Xiang J and Li F (2019) Identification of Single Nucleotide Polymorphisms Related to the Resistance Against Acute Hepatopancreatic Necrosis Disease in the Pacific White Shrimp Litopenaeus vannamei by Target Sequencing Approach. Front. Genet. 10:700. doi: 10.3389/fgene.2019.00700
Received: 28 February 2019; Accepted: 03 July 2019;
Published: 02 August 2019.
Edited by:
Gen Hua Yue, Temasek Life Sciences Laboratory, SingaporeReviewed by:
Jian Xu, Chinese Academy of Fishery Sciences (CAFS), ChinaBaoqing Ye, Temasek Life Sciences Laboratory, Singapore
Copyright © 2019 Zhang, Yu, Wang, Liu, Luo, Zhang, Zhang, Huang, Xiang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fuhua Li, fhli@qdio.ac.cn
†These authors have contributed equally to this work.