Circulatory white spot syndrome virus in South-West region of Bangladesh from 2014 to 2017: molecular characterization and genetic variation

AMB Express, Feb 2018

White Spot Syndrome Virus (WSSV), the etiological agent of White Spot Disease (WSD) is a major impediment for shrimp aquaculture in the worldwide. A critical threshold level of WSSV load in infected shrimp is an important trait for disease manifestation and WSSV transmission in cultured shrimp and subsequently make outbreaks. The present study investigated 120 naturally infected cultured shrimp samples by SYBR Green based qPCR assay for WSD diagnosis and quantification of WSSV load. Among them, 94 samples resulted a variable count of WSSV load ranging from 2.1 × 108 to 2.64 × 1014 copies/g of shrimp tissue. The severity of WSSV infection was assessed based on the established critical threshold load of WSSV in shrimp tissue. Compared to the established critical threshold value of WSSV load in shrimp tissue, our findings showed the horrifying scenario of the severity of WSSV infection in cultured shrimps of Bangladesh that was found to be above the critical limit to initiate an outbreak in the Bangladeshi shrimp aquaculture industry. The latest phylogenetic pattern was altered from the former monophyletic history among WSSVs of Bangladesh due to a variation at 500th nucleotide of VP28 coding gene. Viruses characterized from recent outbreaks in 2015 and 2017 displayed amino acid substitution at position 167 (G→E) on the surface of VP28 protein which has demonstrated the probable replacement of indigenous virus pool. Therefore, it is imperative to take initiative for the management and prevention of WSSV outbreak to sustain shrimp aquaculture in South-West region of Bangladesh.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

Circulatory white spot syndrome virus in South-West region of Bangladesh from 2014 to 2017: molecular characterization and genetic variation

Siddique et al. AMB Expr Circulatory white spot syndrome virus in South-West region of Bangladesh from 2014 to 2017: molecular characterization and genetic variation Mohammad Anwar Siddique 0 1 Md. Inja‑Mamun Haque 0 3 Santonu Kumar Sanyal 1 2 Anwar Hossain 3 Shuvro Prokash Nandi 1 2 A. S. M. Rubayet Ul Alam 1 Munawar Sultana 1 Mahmud Hasan 3 M. Anwar Hossain 1 0 Mohammad Anwar Siddique and Md. Inja‐ Mamun Haque have 1 Department of Microbiology, University of Dhaka , Dhaka 1000 , Bangladesh 2 Pre‐ sent Address: Department of Microbiology, Jagannath University , Dhaka 1100 , Bangladesh 3 Department of Fisheries, University of Dhaka , Dhaka 1000 , Bangladesh White Spot Syndrome Virus (WSSV), the etiological agent of White Spot Disease (WSD) is a major impediment for shrimp aquaculture in the worldwide. A critical threshold level of WSSV load in infected shrimp is an important trait for disease manifestation and WSSV transmission in cultured shrimp and subsequently make outbreaks. The present study investigated 120 naturally infected cultured shrimp samples by SYBR Green based qPCR assay for WSD diagnosis and quantification of WSSV load. Among them, 94 samples resulted a variable count of WSSV load ranging from 2.1 × 108 to 2.64 × 1014 copies/g of shrimp tissue. The severity of WSSV infection was assessed based on the established critical threshold load of WSSV in shrimp tissue. Compared to the established critical threshold value of WSSV load in shrimp tissue, our findings showed the horrifying scenario of the severity of WSSV infection in cultured shrimps of Bangladesh that was found to be above the critical limit to initiate an outbreak in the Bangladeshi shrimp aquaculture industry. The latest phylogenetic pattern was altered from the former monophyletic history among WSSVs of Bangladesh due to a variation at 500th nucleotide of VP28 coding gene. Viruses characterized from recent outbreaks in 2015 and 2017 displayed amino acid substitution at position 167 (G→E) on the surface of VP28 protein which has demonstrated the probable replacement of indigenous virus pool. Therefore, it is imperative to take initiative for the management and prevention of WSSV outbreak to sustain shrimp aquaculture in South‑ West region of Bangladesh. Cultured shrimp; VP28; Real‑ time qPCR; Severity; Mutation; Phylogeny Introduction White spot syndrome virus (WSSV) is a large dsDNA virus and the causative agent of White spot disease (WSD), has been reported to cause severe mortality in farmed shrimp especially black tiger shrimp, Penaeus monodon in many countries (Balasubramanian et  al. 2008) . Black tiger shrimp is the main coastal aquaculture species and second highest foreign currency earning source in Bangladesh (DoF 2013) . Khulna division is the largest farmed shrimp producer (accounting for 86% of national shrimp aquaculture production) and Satkhira district (Southwestern parts of Bangladesh) accounted for 31% of Khulna shrimp farm production and 26% of national shrimp farm production in 2014 (BBS 2017) . However, shrimp aquaculture industry in Bangladesh has been greatly suffered from viral diseases particularly WSD and losses hundreds of million dollars every year (Lafferty et  al. 2015) . Since its first report in East Asia in 1992, WSSV spread quickly to shrimp-farming areas in Southeast Asia, Central and Latin America, causing major economic damage to shrimp aquaculture (Oidtmann and Stentiford 2011) . As a result, WSD has been listed a devastating viral disease of crustaceans by the Office of International Epizootic (OIE) (Nunan and Lightner 2011) . Cultured penaeid shrimp from post-larvae (PL) to brood stock are highly susceptible to WSSV infection (Vidal et  al. 2001) . Cultured shrimp life span is maximum of four to five  months. During the culture period, shrimp have to tolerate more variable conditions, less controlled temperature and salinity fluctuation (Jeswin et  al. 2015) . Extreme salinities is identified as one of the key environmental stress factor that might intensify the susceptibility of shrimp to viral infection (Ramos-Carreño et al. 2014) . Therefore, recurrent variation of environmental conditions such as temperature, pH and salinity make shrimp more susceptible to WSSV infection and in such conditions viral load has been found at a very high titer. The viral load in infected shrimp is known to be one of the most important traits for transmission  of virus. WSD outbreak is triggered when the viral load in the infected shrimp reached a critical concentration (Durand et  al. 2003). Several studies on WSSV genetic structure and functions have shown that a few envelope proteins are involved in initiation of infection. WSSV genome contains a single copy of VP28 gene (615  bp, 204 amino acids) which encodes for a key structural protein for virion envelope. This envelope protein plays a significant role in progression of systemic infection as it is often involved with viral entry, assembly and budding (Syed Musthaq et  al. 2009) . It infers that, quantification of copies of this gene in extracted DNA samples from WSSV infected shrimp might be indicative of the number of virions and naked genome prior to assembly present in the infected or carrier host. In Bangladesh, it has been observed from field experience that some farmers in South-West region are able to obtain a reasonably good harvest despite the presence of WSSV infection characterized by severe white spot on carapace, while similar observations have also been reported from China (Lan et  al. 2002) . These observations suggested the existence of possible disease resistant shrimp and genetic variants of WSSV (Chakrabarty et  al. 2014a, b; Mallik et  al. 2016; Morse 1994) . Genetic variants arisen due to mutations could be a probable reason behind molecular diagnostic failure mainly owing to inhibition of binding PCR primers; that might lead to false negative PCR results or non-specific PCR products which in return, limits the use of PCR based diagnostic kits (Mendoza-cano and Sanchez-paz 2013; Kwok et  al. 1990) . Thereupon, the comparison of the nucleotide variation in VP28 gene could be considered a possible marker to differentiate WSSV isolates. Furthermore, our earlier studies were confined to detection, not WSSV load quantification which did not depict the complete scenario of circulating WSSV load in cultured shrimp of Bangladesh (Hossain et  al. 2015) . Therefore, this study focused for the first time on the quantification of the circulating WSSV load in cultured shrimp samples in this region. Hence, the work reports the first quantitative assay of WSSV infection in naturally infected cultured shrimp to infer the severity of WSD and describes duly the genetic variation of WSSV circulating in the South-West region of Bangladesh. Materials and methods Shrimp samples collection The moribund shrimp samples (n =  120) were collected randomly from WSSV prevalent zone of Satkhira district, Bangladesh from January 2014 to June 2017. The sample collection techniques and transportation procedure was done by following the protocol described by Hossain et al. (2015). DNA extraction Total DNA was extracted from tissue samples of WSD suspected shrimps by Maxwell® 16 automated nucleic acid extraction system (Promega, Madison, WI, USA) using Maxwell® 16 tissue DNA purification kit (Promega, Madison, WI, USA; Catalog-AS1030) according to manufacturer’s instruction. Concentration and optical purity of the extracted DNA was evaluated using NanoDrop™ 2000. Quantitative real‑time PCR (qPCR) assay Our method was optimized with newly designed primer (WSSV-q28F 5′-TGTGACCAAGACCATCGAAA-3′ and WSSV-q28R 5′-CTTGATTTTGCCCAAGGTGT-3′) following the method developed by Mendoza-Cano and Sánchez-Paz (2013) with some modification (recombinant plasmid based standard instead of purified PCR product based standard). The recombinant plasmid containing VP28 gene (TOPO TA Vector having complete CDS of VP28 gene as an insert) was gel purified using Wizard® SV Gel and PCR Clean-Up System (Promega, Madison, WI, USA). The standard was prepared from serial dilution of purified recombinant plasmid in a linear logarithmic scale of 1.0 × 109 to 102 copies per reaction. All qPCR reactions were run at a final volume of 25 µL in the Applied Biosystems® 7500 Real-Time PCR system (Foster City, CA, USA) using 2× SYBR® Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA), 100  nM of each of the primers and 2.5  µL (10% of total reaction volume) of each template DNA. Thermal cycling parameters were set for an initial denaturation step at 95 °C for 10 min followed by 40 cycles at 95 °C for 15 s for DNA denaturation with subsequent annealing and extension at 53  °C for 30  s. Melt curve analysis was also performed to differentiate the specific amplicon from primer dimer formation or amplification of other non-specific product. Moreover, the qPCR products were also electrophoresed in agarose gel to nullify the presence of spurious amplicon. The present experiment was conducted with duplicates replication to quantify the viral load. WSSV load per gram of tissue sample was calculated according to the following equation: Viral load per gram of tissue  =  [viral load per reaction  ×  (Final Elution volume/volume of template DNA per reaction)  ×  dilution factor]  ±  Standard Deviation (SD). To assess the reproducibility of the standard curve, standard reactions were performed three times independently including duplications of each reaction. The realtime PCR data were analyzed by using 7500 software, version 2.0.6 (Applied Biosystems, Foster City, CA, USA). The data were analyzed by using the statistical program Microsoft Excel 2010 and presented as mean ± SD. Polymerase chain reaction amplification and sequencing of VP28 A conventional PCR was performed to amplify the VP28 gene using GoTaq 2× Hot Start Colorless Master Mix (Promega, USA) with forward and reverse primer (Rout et  al. 2007) . The PCR-amplified products were analyzed by 1% agarose gels containing ethidium bromide. PCR amplicons were purified using the Wizard® SV Gel and PCR Clean-Up System (Promega, USA) and purified PCR products were subjected to an automated cycle sequencing reaction using BigDye® Terminator v3.1 cycle sequencing kit (Applied Biosystems®, USA) according to manufacturer’s instructions with the same primers used in the PCR reaction and analysis of the data was performed in ABI Genetic Analyzer (Applied Biosystems®, USA). Sequence analysis and submission to NCBI GenBank The raw sequence data were assembled using SeqMan version 7.0 (DNASTAR, Inc., Madison, WI, USA) and the assembled sequences were compared with other entries from NCBI GenBank (Benson et  al. 2005) using BLAST (Altschul et al. 1990) search to reveal the identification and matching with VP28 gene of WSSV. Afterwards, assembled sequences of the VP28 gene of WSSV were submitted to NCBI GenBank (Accession No. MF489075–MF489082). Assortment of WSSV sequence dataset and rationale behind the dataset generation A total of 74 sequences were compared in this investigation, among which 18 isolates had origin in Bangladesh (Additional file 1: Table S2). In this study, we have sequenced eight of these isolates from Bangladesh. 56 reference sequences were collected from NCBI GenBank database to represent other countries where WSSV has been identified. We have screened the sequences based on cluster formation in phylogenetic tree to cover each country. After the thorough screening, 42 sequences were selected for final phylogenetic tree reconstruction. In case of mutational analysis, 66 sequences were taken into consideration due to the identical amino acid composition of ten sequences in the dataset. These sequences were generated from the isolates of same location and collection year (Chidambaram/M(3–8)/ India/2009, Kadalur-NM(1–4)/India/2009) as presented in Additional file 1: Table S2. The rationale was taken into account to produce a non-redundant dataset by screening which was necessary to leverage the final analytical processes and for better view of the generating results. Recombination detection We first analyzed the sequences of WSSV with RDP v4.56 program (Martin et  al. 2015) that incorporates RDP (Martin and Rybicki 2000) , GENECONV (Padidam et  al. 1999) , Bootscan (Martin et  al. 2005) , Chimaera (Hall 1999) , MaxChi (Garcia-Boronat et al. 2008) , SiScan (Kabat et al. 1977) and 3Seq (Boni et al. 2007) methods as well as GARD (Kosakovsky Pond et  al. 2006) method of Datamonkey webserver (Delport et  al. 2010) to find the evidence for recombination events to prevent potential biases during phylogenetic  reconstruction. We considered significant events executed by four or more methods using the Bonferroni (Bland and Altman 1995) correction to nullify the false positive results (Faye et al. 2014) . Phylogenetic tree reconstruction jModelTest (version 2.1.10) package (Darriba et al. 2012) was used for substitution model selection by computing likelihood score out of 56 models. The best fitted model was selected using lowest corrected Akaike information criterion (cAIC) (Hurvich and Tsai 1993) and Bayesian information criterion (BIC) values (Schwarz 1978) . MEGA7 (Kumar et  al. 2016) was also implied to check out the best substitution model for the dataset based on AIC (Akaike 1973) and BIC values. Maximum Likelihood (ML) method (Felsenstein 1981) was implied to reconstruct the phylogenetic tree. Fewer than 5% alignment gaps, missing data, and ambiguous bases were allowed at any position as all the sequences were not completely aligned on full range. To observe the statistical significance, 1000 times replication of the branches of the tree was performed. Mutational analysis Variability scores in both nucleotide and amino acid level were calculated using MEGA7 and Seaview v4.6 (Gouy et  al. 2009) in both overall dataset including all the sequences and only Bangladesh originated isolates. BioEdit Sequence Alignment Editor v7.2.5 (Hall 1999) was used for a better visualization of the mutations with the position keeping all the sequences altogether in a single frame. Protein Variability server (PVS) (GarciaBoronat et al. 2008) was used for finding out protein variability index using the Wu-Kabat variability coefficient of the aligned protein sequences (Kabat et  al. 1977) . Using the PVS, a consensus sequence was also obtained from this alignment that also search for the variable site with a threshold level of less than or equal to 1. Then, Skylign (Wheeler et al. 2014) was used to generate sequence logo after inputting the aligned file. In case of Skylign, fragmented sequences were taken into account and the probability of amino acid residue for each position was detected. Protein modeling and alignment On the basis of major mutation in VP28 isolated from Bangladesh, the protein modeling was performed taken PDB id 2ED6 (Tang et  al. 2007) and BAN_SH_BG-1C_2014 as template and target, respectively. SWISS-MODEL server (Biasini et  al. 2014) was used for generating the protein model and then PyMol (DeLano 2002) was used for visualization of the PDB file. Using PyMol, alignment of the Bangladeshi isolate and the reference VP28 (2ED6) focusing on the region (160–170). Using PDB files of target, Ramachandran plot (Ramachandran et  al. 1963) was produced to verify whether the residues were in one of the three regionsfavored, allowed and outlier using RAMPAGE (Chen et  al. 2010) . To assess both the quality and energy criteria of 3D structure, PROSA (Wiederstein and Sippl 2007; Sippl 1993) was used. The 3D structures of the reference template sequence (PDB id—2ED6) and Bangladeshi isolate sequence (BAN_SH_BG-1C_2014) were superimposed on each other to visualize the mutational change in a clear view by implying PyMol. Results Quantification of WSSV in cultured shrimps To quantify viral loads, 120 shrimp samples from 120 different farms were tested (from January 2014 to June 2017). Positive amplification pattern was observed in 94 samples with the presence of specific PCR product of 148  bp confirmed by significant fluorescence signal having a signature asymptotic curve (Fig.  1). No significant fluorescence signal was observed for the duplicated negative control (NTC) and 26 negative samples (compared to standard positive samples). The cycle threshold value (CT) for the negative samples and negative control of amplification was beyond determination index (three times of independent real-time PCR assay with duplication in each sample). On the other hand, cycle threshold value (CT) in all 94 positive samples were between 18 and 33, and all positive samples showed quantitative index in real-time qPCR assay. After amplification, all samples were subjected to melt curve analysis. Melting temperatures of amplicons (Tm) were presented by plotting the values of negative derivation of the fluorescence signal (-Rn) against temperature in degree Celsius. Melt curve analysis of standard (recombinant plasmid based) and WSSV positive field samples (crude DNA) yielded single expected dissociation peak of Tm at 81.11 ± 0.17 °C, indicating specific amplification of target gene (Fig.  2). Another Tm peak (ranged from 65 to 67° C) was evident beyond the expected one after melt curve analysis especially of the negative samples which suggests the presence of possible primer dimer. The agarose gel electrophoresis of the qPCR products nullified the presence of spurious amplicon other than the primer dimer (in case of negative samples). WSSV load and severity of infection in cultured shrimps of Bangladesh Quantitative study revealed a wide range of WSSV load (8.7 × 103 to 1.10 × 1010 copies; mean value) per 2.5 µL of template DNA which is equivalent to 2.1  ×  108 to 2.64  ×  1014  copies/g of tissue (Additional file  1: Table S1). Among 94 positive samples, 87 samples were found to contain WSSV load between the logarithmic scales of 9 and 11. However, WSSV  load in two samples were on the logarithmic scale of 12. In addition, one sample contained WSSV load at 14 logarithmic scale  (phenotypic observation:  highly spotted on entire exoskeleton) and WSSV load of 8 logarithmic scale was observed for another four samples (phenotypic observation:  spot not appeared on exoskeleton of shrimp). The degree of WSSV infection was categorized based on the established critical threshold limit (Table  1) described by Walker et  al. (2011) and Tang and Lightner (2000) . According to Walker et  al. (2011), out of 120 shrimp samples, 94 were WSSV positive in which 73 (77.66%) were graded as moderate to heavy infection and 21 (22.34%) were grouped as light or very light infections (Table 1). In another context based on the report of Tang and Lightner (2000) , the number of WSSV genomes increased over 100-fold between severity levels (Table 1). The severity of infection level G3 (moderate to severe) and G4 (severe) was found in 87 (72.5%) samples; altogether 78.33% samples were infected which stipulates the alarming situation of WSD in the South-West region of Bangladesh. Compiled dataset for phylogenetic and mutational analyses Altogether 74 sequences were taken into account in this study where 18 of these were isolated from Bangladesh considering the location and year of isolation (Additional file 1: Table S2). Among the isolates, 66 of the sequences were chosen for mutational analysis. Observation of phylogenetic history RDP and GARD based recombination detection methods implemented in this study showed no breakpoint event and hence recombination in VP28 sequences along with out-group [Monodon Baculovirus (MBV)] included for further phylogenetic studies. Tamura-Nei substitution model with uniform rate variation and pattern was chosen as best model based on AIC, cAIC and BIC values generated from jModelTest and MEGA7 softwares. The isolates from Bangladesh could be divided into different clusters as shown from the tree (Fig. 3). Two clusters (group 1 and 2) contained 17 of the sequences isolated from Bangladesh that mostly related to isolates from India and Vietnam. Another cluster contained only one isolate from Bangladesh (BAN_SH_KU-2_2017) along with isolates from other countries. The other countries included South Korea, China, Thailand, Iran, Mexico, Netherlands, India, Brazil, Japan, USA, Indonesia and Vietnam. The isolate BAN_SH_SU-3_2015 clustered with other VP28 sequences that was previously reported (Hossain et  al. 2015) to be circulated in Bangladesh. BAN_SH_KU-2_2017 with the isolates of China and South Korea formed a slightly different branch within group 3. VP28 gene sequences of other local samples Analysis of mutation pattern in VP28 In total, there were found 13 mutations (Additional file 2: Figure S1) covering all the isolates from different countries specially focusing on and clustering the Bangladeshi isolates serially  for better visualization of the dissimilarity with other sequences of VP28 (Table  2). White spot syndrome virus strain CN01 (Accession No. NC_003225) was used as reference sequence for mutational analysis unambiguity. Amino acid residues at positions 42 and 167 showed mutations in case of Bangladeshi isolates among which position 167 showed glycine instead of glutamic acid for the isolates circulating in 2014 (Hossain et  al. 2015) and only one isolate from 2015 (BAN_SH_ SU-3_2015) (Additional file 3: Figure S2). Prediction and mapping of exclusive mutation found in Bangladesh After producing PDB file, the quality of the structure was checked using Ramachandran plot (Additional file 4: Figure S3) showing that 96.4% of the residues fell within the favored region. The PROSA analysis of the model showed maximum residues to have negative interaction energy, and the overall Z-score is − 6.4. In the superimposed structure, the region flanking (10  bp) the unique mutation (E→G) was focused (Fig.  4) and the mutated amino acid side chain protruded out of the structure. This position is in between β-barrel no. seven and eight (β7–β8). Discussion Higher severity indicates WSD outbreak in SouthW‑est region of Bangladesh Satkhira is the top-tier shrimp producing district (BBS 2017) , therefore, it is imperative to conduct regular surveillance of WSSV in this area. Our optimized SYBR Green qPCR method was applied to diagnose and quantify WSSV load in diseased cultured shrimp samples. In our previous studies, conventional and realtime PCR based detection of WSSV in cultured shrimp was reported from Bangladesh (Hossain et  al. 2015, 2016) . In this study, quantitative real-time PCR assay was performed on field shrimp samples; to reduce the (See figure on next page) Fig. 3 Molecular phylogenetic analysis (codon based) of WSSV isolated from Bangladesh between 2014 to 2017. Nucleotide sequences of VP28 coding region were used to construct the tree based on the Tamura‑Nei model by Maximum Likelihood method. The analysis included 42 nucleotide sequences among which 18 were from Bangladesh. The sequences generated from isolates of Bangladesh are shown in bold letter, which are placed in three distinct clusters among four. Sequences generated exclusively for this study and submitted from previous study (Hossain et al. 2015) are marked with (filled triangle) and (filled diamond) symbols, respectively. Monodon Baculovirus (Accession No. HQ222840) was taken as outgroup in the tree The country name and year of isolation of the isolates that showed the mutation from the Reference Sequence CN01/China/2000 (Accession No.—NC_003225) Italicized data indicate the specific amino acid variations (frequency) with respective positions and regions in VP28 of WSSV isolates collected from Bangladesh and related countries possible PCR inhibitors, all samples were diluted one logarithmic scale which is supported by a former study (Leal et al. 2013) . The resulted standard deviation value of Tm (± 0.17) °C revealed significant agreement in the amplified PCR products (Mendoza-Cano and SánchezPaz 2013) . Infections at low or even moderate  to  high intensities do not necessarily result in disease, however, once the virus reaches a critical threshold level in some shrimp, disease might occur (Walker et  al. 2011) . Once some shrimp in a pond get infected with WSSV, other shrimp in that pond will be exposed to high level of virus through ingestion or immersion and quickly develop severe disease resulting in rapid production failure. Presence of minimum 105 WSSV copies is necessary to transmit WSSV in shrimp by immersion in cultured condition (Meng et  al. 2010) . A recent report demonstrated that 107.1–107.5 virus per gram of tissue is needed to infect shrimp by peroral inoculation and via feeding (Thuong et  al. 2016) . A study by Jeswin et  al. (2015) has demonstrated that a dose of 6 × 107 WSSVg−1 was sufficient to kill shrimp effectively within 96 h. Analysis of WSSV titer in diseased shrimp under current study, showed that all 94 WSSV positive shrimp samples had viral load above 108  copies/g of tissue. This high viral loads might pave a way to transmit viruses to non-infected cultured shrimp. With poor pond management, the virus causing the disease outbreak can also spread to adjacent ponds in a radiating outbreak, primarily through the transfer of WSSV contaminated water (Walker et al. 2011) . In addition, our personal experiences, farmer’s information and observed field scenario also confirmed WSD outbreak in the present study areas. Sahoo et  al. (2010 ) have reported that moderate or heavy level of WSSV presence in ponds is significantly associated with production outcomes in terms of survival, yield and length of the production cycle. Furthermore, the greater the level of severity of WSSV infection, the faster the rate of mortality. In the study area, most farms practice improved extensive or semi-intensive culture technology and do not follow scientific management. Tidal water is filled in, and drained out, by using the same structure as inlet and outlet of the farm. The water enters into the pond was found to mix with the discharge of the same farm or other farms that would contaminated the water. These factors, collectively, might be attributed as one of the most key sources of WSSV infections in cultured shrimps in Bangladesh. Till date, we have no sustainable treatment technique available for extermination of WSSV. WSSV infected shrimp cannot be cured, only preventive measures could be taken to get rid of this burden. The whole culture environments might be improvised by implementing good aquaculture practices (GAqP) and regular surveillance of critical WSSV load (a load that implies disease outbreak) using quantitative real-time PCR technique that might resist mass mortality due to WSSV infection and huge economic destruction from unexpected outbreak of this disease in Bangladesh. In our previous report (Hossain et  al. 2015) , we presented the high prevalence of WSSV in shrimp farming areas of Bangladesh from where the field samples were collected. Quantitative index of high WSSV titer observed in the current study might be a key underlying cause of the WSD outbreaks in South-West region of Bangladesh. One single nucleotide change has altered the phylogenetic pattern among Bangladeshi isolates In DNA virus, rates of evolutionary change in genome is slow and the nucleotide substitution rate is far lower than RNA virus (Duffy et al. 2008) . Being a DNA virus, WSSV must show the properties of being conserved in both nucleotide and protein level. In case of VP28 of WSSV, the major structural protein coding gene displays conservation of the DNA sequence, even at a level of 99–100% in most cases (You et al. 2004; Molina‐Garza et al. 2008) . In Bangladesh, the isolated viruses showed genetic divergence falling under different clusters. These different clusters consisted of WSSV samples from other countries. From the phylogenetic tree, it can be assumed that all of the isolates collected in 2014 from Bangladesh were quite diverged from the samples collected in the following years and the result correlates with the previous work (Hossain et  al. 2015) . Along with the 2014 isolates, two isolates from 2015 clustered in the same group- BAN_ SH_SU-3_2015 and BAN_SH_AL-2_2015. On the other hand, isolates from 2017 got divided into two groups wherein BAN_SH_KU-2_2017 fell under separate clade from the other four isolates. The sequences from India, Vietnam and Bangladesh fell into all of the major clades, that is not found in case of other countries as shown in the tree. A more conservative approach to phylogeny might give a different result from the one shown in this study. However, the divergence was not that significant in the reconstructed phylogenetic tree as there have no significant variations present among the sequences of VP28. In case of Bangladeshi isolates, only one nucleotide got changed that means 0.3% of the total nucleotide changed. Interestingly, this single change has a major impact on the protein variability of VP28. We have found by calculating the diversity both in nucleotide and protein level among sequences from different geographical locations that 5% and 7.7% of the nucleotides and amino acids were variable. The conservation of sequences of the coding region of structural proteins of different WSSV isolates provides a strong basis for generating highly specific and sensitive nucleic acid and immuno-based molecular detection methodologies for WSSV (You et al. 2004) . As an example, we have used PCR primers designed based on the sequence of VP28 gene for detection of WSSV. Based on the sequence homology or similarity, complete VP28 protein can be applied as potential recombinant antigen to produce vaccine that triggers the adaptive quasi-immune response in shrimp (Namikoshi et al. 2004, Venegas et al. 2000) . However, all WSSV proteins involved in viral infection of shrimp need to be assessed in molecular level to uncover the detailed genetic make-up. The first time report of mutation pattern of VP28: an exclusive mutation was spotted The shrimp (Penaeus monodon) protein PmRab7 is located in endosome that acts as receptor for the WSSV envelop protein VP28 (Sritunyalucksana et  al. 2006; van Hulten et  al. 2001) . VP28 is a major structural envelope protein of WSSV responsible for systemic infection in shrimp forming an important part of “infectome” that is found to be crucial in cell recognition, attaching and penetration into the shrimp cell (Chang et  al. 2010; van Hulten et al. 2001) . VP28 is located on the outer surface of the virus and the transmembrane N-terminal portion of the protein embedded in the envelope followed by the N–terminal α–helix (α1) and then C–terminal nine stranded β–barrel structure-a feature uncommon in other viral proteins (Tang et al. 2007) . From the result of mutational analysis (summarized in Table 2), some unique mutations at different geographical locations were found. Joseph et  al. (2015 ) reported the nucleotide and amino acid variations found in VP28 of WSSV. For example, two unique mutations were captured from the isolates of Thailand in 2006. The positions of the mutated amino acids are in 38 and 145, respectively among which P→L substitution might be crucial as the amino acid situated in the outside of the envelope in a beta-strand (Tang et  al. 2007) . In case of Bangladesh, the unique mutation (E→G) at position 167 that falls between two beta strands of protein that are thought to be involved in receptor recognition (Verma et  al. 2013) . This exclusive mutation was found in the isolates of Bangladesh collected in 2014 and only in one of the isolates of 2015. Another mutation at position 42 (D→G) was found common to all of the isolates of WSSV in Bangladesh. The significance of these mutations in the particular positions of VP28 was not determined as none of the amino acid was found to interact with PmRab7 of shrimp (Verma et al. 2013) . For finding out the exact position and possible role of the amino acids at the uniquely mutated position in a 3D structure of VP28, we performed the alignment of the designed and representative VP28 protein from Bangladesh over the reference protein from PDB consisting of different amino acid other than that was found for the first time in Bangladesh. The region spanning (160–170) amino acid residues positioned on the surface area rather than groove area of VP28 protein (Fig.  4). The change of amino acid (glutamic acid to glycine, E→G) at position 167 may not give any positive effect on the infectivity of VP28 as it does not give much opportunity to interact with other proteins or ions rather than glutamic acid. Further investigation into the interaction of the mutated amino acid with different related proteins or ions will shed light on the possible role. To reduce the risks of WSD  outbreak in shrimp farming, we recommend to stock virus-free PL produced from PCR-negative broods and proper maintenance of poststocking to reduce risks of WSSV infection. Therefore, quantitative Real-time PCR methods can be applied to detect the infected brood shrimp or PL prior to infestation of the virus to its school that results an outbreak. From our observation, it can be suggested that lack of proper quarantine system on import of fish/shrimp and fish/shrimp brood from the WSSV endemic country could be responsible for a WSD outbreak in Bangladesh. The higher WSSV load in our studied samples which is found above the critical limit for dissemination of the virus can subsequently lead to initiate an outbreak. This study describes a unique mutation pattern of VP28 gene of circulating WSSV in South-West region of Bangladesh. Overall, the findings of this research work could be a beacon for rapid screening of WSSV (including the new mutant) in brood and post-larvae of shrimp and may be a way out to mitigate the huge economic loss every year in shrimp farming of Bangladesh. Additional files Additional file 1: Table S1. qPCR profile of the studied samples by our designed and standardized method. Table S2. Accession no., sequence IDs, sources and collection time of the WSSVs sequences of dataset. Additional file 2: Figure S1. Wu‑Kabat protein variability index showed 13 unique mutations as observed from the protruded peaks’ positioned along with the amino acid residue of VP28 protein. Additional file 3: Figure S2. The aligned protein sequence data of VP28 taking into consideration 66 sequences generated by BioEdit Sequence Alignment Editor. The left portion of the frame contained the isolate ids of the sequences. The digit in the upper portion of the figure represented the residue number of protein. The red bar indicates the unique mutation position along with the mutated amino acid. Additional file 4: Figure S3. Ramachandran plot of modeled VP28 protein of WSSV. The number of residues in favored, allowed and outlier region are 96.4, 3.6 and 0%, respectively. On the X and Y axes, φ (phi) and ψ (psi) represent the torsion angles around Cα (alpha carbon) to amine and carboxyl group of different amino acids. Abbreviations WSSV: white spot syndrome virus; WSD: white spot disease; OIE: office of inter‑ national epizootic; NCBI: national center for biotechnology information; BLAST: basic local alignment search tool; cAIC: corrected akaike information criterion; BIC: Bayesian information criterion; AIC: Akaike information criterion; ML: maximum likelihood; PDB: protein data bank; 3D: three dimensional; CT: cycle threshold; Tm: melting temperature; MBV: Penaeus monodon‑type baculovirus; P: proline; L: lysine; E: glutamic acid; G: glycine; D: aspartic acid. Authors’ contributions MAS, MIH, AH and SKS carried out the studies (sampling, molecular and data analysis) and participating in drafted manuscript. SPN and MS partially car‑ ried out work on bench and helped to draft the manuscript. ARA, MIH and MAS carried out work the bioinformatics analysis. MAH and MH developed hypothesis, supervised the whole work and helped to prepare the manuscript. All authors read and approved the final manuscript. Authors’ information MAS is a Lecturer in the Department of Microbiology, University of Dhaka, Bangladesh and has a keen interest in the field of molecular microbiology, virology, genomics and bioinformatics. MIH is a Lecturer in the Department of Fisheries, University of Dhaka, Bangladesh, with an excellent academic performance. His key research inter‑ est is in the field of molecular and computational biology. SKS is a Lecturer in the Department of Microbiology, Jagannath Univer‑ sity, Bangladesh and his research interest is in the field of environmental microbiology. AH is an Assistant professor in the Department of Fisheries, University of Dhaka, Bangladesh. His key research interest is in the field of molecular biology. SPN is an Assistant Professor in the Department of Microbiology, Jagan‑ nath University, Bangladesh and primarily interested in the field of molecular microbiology, virology, cancer biology and computational biology. ARA has completed his B.S. (Honours) and M.S. from the Department of Microbiology, University of Dhaka, Bangladesh. His key research interest is in the field of molecular and computational biology. MS is an Associate Professor in the Department of Microbiology, University of Dhaka, Bangladesh and her main area of research is antibiotic resistance and arsenic microbiome. MH is a Professor in the Department of Fisheries, University of Dhaka, Bangladesh. He is currently working on commercial culture of Tubificid worms. MAH is the leader of the Microbial Genetics and Bioinformatics Laboratory and a Professor in the Department of Microbiology, University of Dhaka, Bang‑ ladesh. He is currently working on antibiotic resistance in bacteria, develop‑ ment of Foot‑and‑Mouth Disease Vaccine, biofilm formation, and heavy metal resistance in bacteria. Acknowledgements The corresponding author thanks UGC and Ministry of Education, Govt. of Bangladesh for providing the grant for the establishment of the facilities used in this work. This research was partially supported by ACI Ltd. Bangladesh. Competing interest The authors declare that they have no competing interests. Availability of data and materials The dataset(s) supporting the conclusions of this article is(are) included within the article and its additional file(s). Consent for publication Not applicable. Ethics approval and consent to participate Not applicable. Funding Not applicable. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in pub‑ lished maps and institutional affiliations. Akaike H ( 1973 ) Information theory and an extension of the maximum likelihood principle . In: Petrov BN , Csaki F (eds) Proceedings of the 2nd international symposium on information theory. Akademiai Kiado , Budapest, pp 267 - 281 Altschul SF , Gish W , Miller W , Myers EW , Lipman DJ ( 1990 ) Basic local alignment search tool . J Mol Biol 215 : 403 - 410 Balasubramanian G , Sarathi M , Venkatesan C , Thomas J , Sahul Hameed A ( 2008 ) Oral administration of antiviral plant extract of Cynodon dactylon on a large scale production against white spot syndrome virus (WSSV) in Penaeus monodon . Aquaculture 279 : 2 - 5 BBS , Yearbook of Agricultural Statistics‑ 2017 ( 28thSeries ) ( 2017 ) Bangladesh Bureau of Statistics, Statistics and Informatics Division (SID), Ministry of Planning, Government of the People's Republic of Bangladesh . BBS, Dhaka Benson DA , Karsch‑Mizrachi I , Lipman DJ , Ostell J , Wheeler DL ( 2005 ) GenBank . Nucleic Acids Res 33 : D34 - D38 Biasini M , Bienert S , Waterhouse A , Arnold K , Studer G , Schmidt T , Kiefer F , Cassarino TG , Bertoni M , Bordoli L , Schwede T ( 2014 ) SWISS‑MODEL: modelling protein tertiary and quaternary structure using evolutionary information . Nucleic Acids Res 42 : W252 - W258 Bland JM , Altman DG ( 1995 ) Multiple significance tests: the Bonferroni method . BMJ 310:170 Boni MF , Posada D , Feldman MW ( 2007 ) An exact nonparametric method for inferring mosaic structure in sequence triplets . Genetics 176 : 1035 - 1047 Chakrabarty U , Dutta S , Mallik A , Mandal N ( 2014a ) White spot syndrome virus (WSSV) and prevalence of disease resistance in a commercially cultured population of Penaeus monodon Fabricius, 1798 (Decapoda , Dendrobranchiata). Crustaceana 87 : 1593 - 1605 Chakrabarty U , Mallik A , Mondal D , Dutta S , Mandal N ( 2014b ) Assessment of WSSV prevalence and distribution of disease‑resistant shrimp among the wild population of Penaeus monodon along the west coast of India . J Inv Path 119 : 12 - 18 . Trans Emer Dis 63 : 1 - 10 Chang YS , Liu WJ , Lee CC , Chou TL , Lee YT , Wu TS , Huang JY , Huang WT , Lee TL , Kou GH , Wang AH ( 2010 ) A 3D model of the membrane protein complex formed by the white spot syndrome virus structural proteins . PLoS ONE 5:e10718 Chen VB , Arendall WB , Headd JJ , Keedy DA , Immormino RM , Kapral GJ , Murray LW , Richardson JS , Richardson DC ( 2010 ) MolProbity: all‑atom structure validation for macromolecular crystallography . Acta Crystallographica Sect D Biol Crystal 66 : 12 - 21 Darriba D , Taboada GL , Doallo R , Posada D ( 2012 ) jModelTest 2: more models, new heuristics and parallel computing . Nat Method 9 : 772 DeLano WL ( 2002 ) The PyMOL molecular graphics system http://pymol .org Delport W , Poon AF , Frost SD , Kosakovsky Pond SL ( 2010 ) Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology . Bioinformatics 26 : 2455 - 2457 DoF ( 2013 ) National fish week 2013 compendium, DoF, ministry of fisheries and livestock, People's Republic of Bangladesh . DoF, Dhaka Duffy S , Shackelton LA , Holmes EC ( 2008 ) Rates of evolutionary change in viruses: patterns and determinants . Nat Rev Genet 9 : 267 - 276 Durand S , Redman R , Mohney L , Tang‑Nelson K , Bonami JR , Lightner DV ( 2003 ) Qualitative and quantitative studies on the relative virus load of tails and heads of shrimp acutely infected with WSSV . Aquaculture 216 : 9 - 18 Faye O , Freire CC , Iamarino A , Faye O , de Oliveira JV , Diallo M , Zanotto PM ( 2014 ) Molecular evolution of Zika virus during its emergence in the 20th century . PLoS Negl Trop Dis 8 : e2636 Felsenstein J ( 1981 ) Evolutionary trees from DNA sequences: a maximum likelihood approach . J Mol Evol 17 : 368 - 376 Garcia‑Boronat M , Diez‑Rivero CM , Reinherz EL , Reche PA ( 2008 ) PVS: a web server for protein sequence variability analysis tuned to facilitate conserved epitope discovery . Nucleic Acids Res 36 : W35 - W41 Gouy M , Guindon S , Gascuel O ( 2009 ) SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building . Mol Biol Evol 27 : 221 - 224 Hall TA ( 1999 ) BioEdit a user‑friendly biological sequence alignment editor and analysis program for Windows 95 /98/NT. Information Retrieval Ltd., London, pp 95 - 98 Hossain A , Nandi SP , Siddique MA , Sanyal SK , Sultana M , Hossain MA ( 2015 ) Prevalence and distribution of white spot syndrome virus in cultured shrimp . Lett Appl Microbiol 60 : 128 - 134 Hossain A , Sanyal SK , Siddique MA , Biswas RK , Sultana M , Hossain MA ( 2016 ) Envelope protein gene VP466‑a target for PCR detection of white spot syndrome virus in shrimp . Bangladesh J Microbiol 31 ( 1 ): 65 - 68 Hurvich CM , Tsai CL ( 1993 ) A corrected akaike information criterion for vector autoregressive model selection . J Time Series Anal 14 : 271 - 279 Jeswin J , Anju A , Thomas PC , Paulton MP , Vijayan KK ( 2015 ) Analysis of viral load between different tissues and rate of progression of white spot syndrome virus (WSSV) in Penaeus monodon . Aqua Res 46 : 2003 - 2012 Joseph TC , Rajan LA , James R , Lalitha K , Surendran P ( 2015 ) Variations of structural protein sequences among geographical isolates of white spot syndrome virus . Int Aquat Res 7 : 85 - 91 Kabat E , Wu T , Bilofsky H ( 1977 ) Unusual distributions of amino acids in complementarity‑ determining (hypervariable) segments of heavy and light chains of immunoglobulins and their possible roles in specificity of antibody‑ combining sites . J Biol Chem 252 : 6609 - 6616 Kosakovsky Pond SL , Posada D , Gravenor MB , Woelk CH , Frost SD ( 2006 ) Automated phylogenetic detection of recombination using a genetic algorithm . Mol Biol Evol 23 : 1891 - 1901 Kumar S , Stecher G , Tamura K ( 2016 ) MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets . Mol Biol Evol 33 : 1870 - 1874 Kwok S , Kellogg D , McKinney N , Spasic D , Goda L , Levenson C ( 1990 ) Effects of primer‑template mismatches on the polymerase chain reaction: human immunodeficiency virus type 1 model studies . Nucleic Acid Res 18 : 999 - 1005 Lafferty KD , Harvell CD , Conrad JM , Friedman CS , Kent ML , Kuris AM , Powell EN , Rondeau D , Saksida SM ( 2015 ) Infectious diseases affect marine fisheries and aquaculture economics . Ann Rev Mar Sci 7 : 471 - 496 Lan Y , Lu W , Xu X ( 2002 ) Genomic instability of prawn white spot bacilliform virus (WSBV) and its association to virus virulence . Virus Res 90 : 269 - 274 Leal C , Carvalho‑ Castro G , Cottorello A , Leite R , Figueiredo H ( 2013 ) Comparative analysis of conventional PCR and real‑time PCR to diagnose shrimp WSD . Braz J Microbiol 44 : 901 - 904 Mallik A , Chakrabarty U , Dutta S , Mondal D , Mandal N ( 2016 ) Study on the distribution of disease‑resistant shrimp identified by DNA markers in respect to WSSV infection in different seasons along the entire east coast of India aiming to prevent white spot disease in penaeus monodon . Trans Emer Dis 63 : 1 - 10 Martin D , Rybicki E ( 2000 ) RDP: detection of recombination amongst aligned sequences . Bioinformatics 16 : 562 - 563 Martin D , Posada D , Crandall K , Williamson C ( 2005 ) A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints . AIDS Res Human Retrov 21 : 98 - 102 Martin DP , Murrell B , Golden M , Khoosal A , Muhire B ( 2015 ) RDP4: detection and analysis of recombination patterns in virus genomes . Virus Evol 1 . Mendoza‑ Cano F , Sánchez‑Paz A ( 2013 ) Development and validation of a quantitative real‑time polymerase chain assay for universal detection of the white spot syndrome virus in marine crustaceans . Virol J 10:186 Meng XH , Jang IK , Seo HC , Cho YR ( 2010 ) A TaqMan real‑time PCR assay for survey of white spot syndrome virus (WSSV) infections in Litopenaeus vannamei postlarvae and shrimp of farms in different grow‑ out seasons . Aquaculture 310 : 32 - 37 Molina‑ Garza ZJ , Galaviz‑Silva L , Rosales‑Encinas J , Alcocer ‑ González J ( 2008 ) Nucleotide sequence variations of the major structural proteins (VP15, VP19, VP26 and VP28) of white spot syndrome virus (WSSV), a pathogen of cultured Litopenaeus vannamei in Mexico . J Fish Dis 31 : 197 - 203 Morse SS ( 1994 ) Towards an evolutionary biology of viruses. The evolutionary biology of viruses . Raven press, New York, pp 1 - 28 Namikoshi A , Wu JL , Yamashita T , Nishizawa T , Nishioka T , Arimoto M ( 2004 ) Vaccination trials with Penaeus japonicus to induce resistance to white spot syndrome virus . Aquaculture 229 : 25 - 35 Nunan LM , Lightner DV ( 2011 ) Optimized PCR assay for detection of white spot syndrome virus (WSSV) . J Virol Methods 171 : 318 - 321 Oidtmann B , Stentiford GD ( 2011 ) White spot syndrome virus (WSSV) concentrations in crustacean tissues-a review of data relevant to assess the risk associated with commodity trade . Transbound Emer Dis 58 : 469 - 482 Padidam M , Sawyer S , Fauquet CM ( 1999 ) Possible emergence of new geminiviruses by frequent recombination . Virology 265 : 218 - 225 Ramachandran GN , Ramakrishnan C , Sasisekharan V ( 1963 ) Stereochemistry of polypeptide chain configurations . J Mol Biol 7 : 95 - 99 Ramos‑ Carreño S , Valencia‑ Yáñez R , Correa‑Sandoval F , Ruíz‑ García N , DíazHerrera F , Giffard‑Mena I ( 2014 ) White spot syndrome virus (WSSV) infection in shrimp (Litopenaeus vannamei) exposed to low and high salinity . Archiv Virol 159 ( 9 ): 2213 - 2222 Rout N , Kumar S , Jaganmohan S , Murugan V ( 2007 ) DNA vaccines encoding viral envelope proteins confer protective immunity against WSSV in black tiger shrimp . Vaccine 25 : 2778 - 2786 Sahoo AK , Mohan C , Shankar K , Corsin F , Turnbull JF , Thakur PC , Hao NV , Morgan KL , Padiyar AP ( 2010 ) Clinical white spot disease status in Penaeus monodon during the middle of the culture period‑its epidemiological significance . J Fis Dis 33 : 609 - 615 Schwarz G ( 1978 ) Estimating the dimension of a model . Ann Stat 6 : 461 - 464 Sippl MJ ( 1993 ) Recognition of errors in three‑ dimensional structures of proteins . Proteins Struct Funct Bioinform 17 : 355 - 362 Sritunyalucksana K , Wannapapho W , Lo CF , Flegel TW ( 2006 ) PmRab7 is a VP28‑binding protein involved in white spot syndrome virus infection in shrimp . J Virol 80 : 10734 - 10742 Syed Musthaq S , Madhan S , Sahul Hameed AS , Kwang J ( 2009 ) Localization of VP28 on the baculovirus envelope and its immunogenicity against white spot syndrome virus in Penaeus monodon . Virology 391 : 315 - 324 Tang KF , Lightner DV ( 2000 ) Quantification of white spot syndrome virus DNA through a competitive polymerase chain reaction . Aquaculture 189 : 11 - 21 Tang X , Wu J , Sivaraman J , Hew CL ( 2007 ) Crystal structures of major envelope proteins VP26 and VP28 from white spot syndrome virus shed light on their evolutionary relationship . J Virol 81 : 6709 - 6717 Thuong K , Tuan V , Li W , Sorgeloos P , Bossier P , Nauwynck H ( 2016 ) Per os infectivity of white spot syndrome virus (WSSV) in white‑legged shrimp (Litopenaeus vannamei) and role of peritrophic membrane . Vet Res 47 : 1 van Hulten MC , Witteveldt J , Snippe M , Vlak JM ( 2001 ) White spot syndrome virus envelope protein VP28 is involved in the systemic infection of shrimp . Virology 285 : 228 - 233 Venegas CA , Nonaka L , Mushiake K , Nishizawa T , Murog K ( 2000 ) Quasiimmune response of Penaeus japonicus to penaeid rod‑shaped DNA virus (PRDV) . Dis Aqua Org 42 : 83 - 90 Verma AK , Gupta S , Verma S , Mishra A , Nagpure NS , Singh SP , Pathak AK , Sarkar UK , Singh SP , Singh M , Seth PK ( 2013 ) Interaction between shrimp and white spot syndrome virus through PmRab7‑ VP28 complex: an insight using simulation and docking studies . J Mol Model 19 : 1285 - 1294 Vidal OM , Granja CB , Aranguren F , Brock JA , Salazar M ( 2001 ) A profound effect of hyperthermia on survival of Litopenaeus vannamei juveniles infected with white spot syndrome virus . J World Aqua Soc 32 : 364 - 372 Walker PJ , Gudkovs N , Mohan C , Raj VS , Pradeep B , Sergeant E , Mohan AC , Ravibabu G , Karunasagar I , Santiago TC ( 2011 ) Longitudinal disease studies in small‑holder black tiger shrimp (Penaeus monodon) ponds in Andhra Pradesh, India. II. Multiple WSSV genotypes associated with disease outbreaks in ponds seeded with uninfected postlarvae . Aquaculture 319 : 18 - 24 Wheeler TJ , Clements J , Finn RD ( 2014 ) Skylign: a tool for creating informative, interactive logos representing sequence alignments and profile hidden Markov models . BMC Bioinf 15 : 7 Wiederstein M , Sippl MJ ( 2007 ) ProSA‑ web: interactive web service for the recognition of errors in three‑ dimensional structures of proteins . Nucleic Acids Res 35 : W407 - W410 You Z , Nadala EC Jr, Yang J , Loh PC ( 2004 ) Conservation of the DNA sequences encoding the major structural viral proteins of WSSV . Dis Aquat Organ 61 : 159 - 163

This is a preview of a remote PDF:

Mohammad Anwar Siddique, Md. Inja-Mamun Haque, Santonu Kumar Sanyal, Anwar Hossain, Shuvro Prokash Nandi, A. S. M. Rubayet Ul Alam, Munawar Sultana, Mahmud Hasan, M. Anwar Hossain. Circulatory white spot syndrome virus in South-West region of Bangladesh from 2014 to 2017: molecular characterization and genetic variation, AMB Express, 2018, 25, DOI: 10.1186/s13568-018-0553-z