Skip to main content
Download PDF
  • Research Article
  • Open access
  • Published:

Moderate heritability of hepatopancreatic parvovirus titre suggests a new option for selection against viral diseases in banana shrimp (Fenneropenaeus merguiensis) and other aquaculture species

Genetics Selection Evolution volume 48, Article number: 64 (2016) Cite this article

Abstract

Background

In shrimp farming, major production losses are caused by viruses. Hepatopancreatic parvovirus (HPV) is one of the viral pathogens that infect banana shrimp (Fenneropenaeus merguiensis). HPV is thought to slow down growth and cause mortality in the juvenile stages of banana shrimp. Genetic improvement through selection of shrimp resistant to viral diseases is one approach to address this issue. This is the first detailed report on an aquaculture species that investigates whether viral titre varies among families and is heritable, and thus whether viral titre per se is a possible candidate trait for selection to produce resistant stock.

Results

HPV titre was measured by quantitative polymerase chain reaction of DNA extracted from 1137 offspring (from 48 full-sib families). Estimated heritability of HPV titre, based on the linear animal mixed model, was moderate (h2 = 0.41). Genetic correlations of HPV with body traits (weight, length and width of body, head and tail) ranged from −0.13 to −0.38. HPV titre was negatively correlated with raw and cooked body colour (−0.33 and −0.43, respectively).

Conclusions

This is the first study based on a large dataset that provides evidence that viral titre may have a genetic component in penaeid shrimp or even in any aquaculture species. The moderate heritability estimated for this trait suggests that resistance to HPV may be achieved by selecting for low HPV titre. With moderate and negative correlations, selection for resistance to HPV should gradually improve body traits and colour of banana shrimp.

Background

Hepatopancreatic parvovirus (HPV) is a viral pathogen that infects penaeid shrimp and also frequently co-infects shrimp with other viruses [1–3]. HPV was reported to cause mortality in the early life stages of shrimp and to stunt growth at juvenile stages in grow-out ponds [4, 5]. In addition, HPV infection was suggested to slow down growth rates because of the observed negative correlation between HPV severity index and shrimp body length [2, 5]. The slowed down growth results in undersized shrimp at harvest time and thus, reduces either quantity or live weight of total production, and consequently contributes to profit loss.

Various biosecurity practices in hatchery and farming systems, e.g., pathogen screening in broodstock, water filtration and disinfection before stocking, have been advocated to control the spread and impact of shrimp viruses [6]. On the one hand, information on the origins and reservoirs of viruses is often lacking and prevents the use of appropriate tools and methods to detect and prevent infections. On the other hand, vaccination against bacteria and viruses has been used in commercial fish farming, especially in salmonids [7, 8]. However, such measures require an adaptive immune system and penaeid shrimp are thought to have only an innate immune system, precluding vaccination as an option to prevent viral infection in shrimp.

In addition to the biosecurity practices, genetic improvement is another option to control the impact of viral diseases by selecting for disease resistance [9]. The main goal of selection for resistance to viral diseases is to develop breeding lines which are less susceptible to viral infection, or at least that tolerate the effects of infection until harvest. Whether such genetic progress can be obtained depends in part on the magnitude of the additive genetic variation (i.e. heritability) of the trait [10, 11].

Heritabilities of resistance against viral diseases have been reported in several aquaculture species. In farmed fish, estimated heritabilities of resistance to viral infection are generally moderate to high and can range from 0.24 to 0.79 [9–12]. In marine shrimp, heritabilities of resistance to viruses have been reported mainly for pacific white shrimp (Litopenaeus vannamei), i.e. against taura syndrome virus (TSV) [13, 14] with moderate heritabilities (h2 = 0.22 to 0.30) and white spot syndrome virus (WSSV) [15, 16] with very low heritabilities (h2 = 0.01 to 0.07).

To date, heritabilities of viral resistance in aquaculture species have been estimated by using survival data after challenge testing. Survival data was considered as a binary variable (i.e. 0 = dead and 1 = alive) that was recorded at the end of the experiment, usually at the time when 50 % of the tested animals are dead or mortality rates flatten out [10, 11, 15]. A classic linear model (a cross-sectional model) is commonly used to handle the binary survival data for genetic estimations; however, this approach does not consider other information in the data such as time until death and other quantitative measures [17]. To better estimate genetic parameters of survival data from challenge testing, a number of more advanced models have been developed, including the threshold model [18], proportional hazard frailty model [19], linear repeatability model [20] and cure model [21]. Such models have been used to predict genetic parameters of survival data in various aquaculture species [9, 11, 16, 22]. For survival from infectious salmon anaemia (ISA) in field outbreaks, the predictive ability of the linear repeatability model was enhanced by 10 to 12 % compared to the classic cross-sectional linear model [23].

Although many comprehensive models have been developed to increase the accuracy of genetic estimates from survival data, setting an experimental infection to obtain survival information still requires a dedicated testing design that considers inoculation approaches, challenge dosage, stock density, testing capacity and operational costs [15]. All intensive attempts to optimize challenge conditions aimed primarily at obtaining records that can represent resistance mechanisms close to those under commercial conditions [16]. Resistance to virus has long been defined as survival from challenge testing, and is commonly applied in farmed fish. However, for shrimp there were concerns on how well survival to acute challenge testing reflects the innate genetic ability to resist to viral diseases because of the lack of a standard protocol for challenge conditions [15, 16]. Lack of standard and repeatable dosages during challenges may contribute to the observed low heritabilities.

An alternative to measuring the phenotypic impact on the host of the virus is to measure the viral titre itself, so that viral titre may be the selection criterion, provided it is heritable. Indeed, Van Koevering et al. [24] and Bosio et al. [25] reported that virus titres are heritable traits, i.e. for spindle streak mosaic virus in wheat and dengue-2 virus in mosquito, respectively. In the current study, we propose, for the first time, a new approach and selection criterion for an aquaculture species i.e. shrimp, which consists in simply measuring viral titre within and among shrimp families. This approach may have several advantages over classical challenge testing and may overcome many of the complications that arise from classical challenge testing. The first possible advantage is that measurements of viral titre can be quantitative, which should add power and resolution to the analyses and selection. Second, viral titre may reflect the innate ability to resist to viruses rather than an ability to survive. There are also a number of other appealing aspects of using viral titre as a measure for viral resistance which have been discussed in our previous study [26]. For example, it does not require a challenge facility to obtain viral titre data; the data can be obtained directly from commercial operations, which makes it highly relevant and applicable to commercial conditions. The high relevance of such viral titre data to commercial conditions can be important especially in shrimp farming systems for which the effectiveness of routine biosecurity practices is hampered by the limited available knowledge about the nature of the viruses (as mentioned above), and thus, often, viral infections cannot be effectively avoided. The lack of need for dedicated experimental challenge facilities reduces costs due to testing operations and facility maintenance. The opportunities to select for viral titre under full commercial conditions apply routinely to shrimp farming in Australia where the future broodstock are taken from commercial grow-out ponds; i.e. broodstock can be selected from the same pond as the animals tested for viral titre based on records of their full-sibs; however such opportunities may not apply so well to other industries that do not take their broodstock from commercial production systems.

One may postulate that using viral titre as a selection criterion irrespective of its heritability cannot yield improvement in survival, due to poor correlations between survival and viral titre and thus that classical challenge testing cannot be replaced by screening for virus titre. However, from a logical point of view, the idea that selection of resistant individuals will reduce or eliminate virus titre and thus, reduce or eliminate the negative impacts of a given virus on a given aquaculture species, is appealing.

Overall, in view of all these potentially desirable aspects of using viral titre as a selection criterion, and notwithstanding the efficacy of viral titre versus survival as a selection criterion, it is remarkable that, until this report and that of Knibb et al. [26], there have been no attempts to estimate the heritability of viral titre in an aquaculture species. Here, we address this knowledge gap by using a large dataset and viral copy number as a quantitative measurement to detect putative genetic resistance to HPV and to estimate the heritability of viral titre and its genetic correlations with body, carcass and flesh quality traits in banana shrimp.

Methods

Animal samples

Hepatopancreas samples of 1137 banana shrimp were collected at harvest time after 140 days of grow-out; shrimp had an average body weight of 16.5 g. These 1137 individuals were the offspring of 48 full-sib families with family size ranging from 9 to 30 offspring. They were cultured within a single grow-out pond in a commercial banana shrimp farm at Cardwell, Australia. Details on rearing conditions, breeding practices, genotyping and pedigree construction are in [27, 28]. The genotypes, based on DNA microsatellites, were used to assign animals to full-sib groups and, thus to establish the pedigree.

Quantitative polymerase chain reaction

Essentially standard curves were produced using known DNA concentrations, and the DNA amounts in the new samples, which were estimated with reference to the standard curve. The genomic DNA of 1137 banana shrimp was extracted using a DNeasy Blood and Tissue Kit (QIAGEN) following the manufacturer’s instructions. DNA concentration was measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific), and normalised to a concentration of 800 ng/μL. HPV titres were quantified based on an absolute quantitative PCR method using the SensiMix HRM kit (Bioline) and validated HPV primers, HPV140F (5′-CTACTCCAATGGAAACTTCTGAGC-3′), HPV140R (5′-GTGGCGTTGGAAGGCACTTC-3′) [29]. Quantitative polymerase chain reaction (qPCR) reactions were performed with a Rotor-Gene 6000 thermal cycler (Corbett Research). Details about standard curve, components of the qPCR reaction mix and cycling conditions are in [26]. Briefly, a standard curve was generated based on tenfold serially-diluted known concentrations of purified HPV PCR products. A final volume of 10 μL reaction mix contained 1 × SensiMix, 0.4 μL Evagreen dye, 3.5 mM MgCl2, 250 nM primers and 4 μL (i.e. 160 ng) normalised DNA. Cycling conditions were as follows: a first step at 95 °C for 10 min, 35 cycles at 95, 60 and 72 °C each for 25 s. The HPV titres were expressed as copy number per μg DNA based on the standard curve.

The hepatopancreas samples were separated into two main groups. The first group consisted of 897 samples for which HPV titres had not been quantified. The second group included 240 samples for which HPV titres had been determined in our previous study [26]. In this study, HPV titres of 135 samples in the second group were quantified for a second time, whereas the remaining 105 samples (that represented five qPCR batches from our previous work) were simply used in the statistical analyses. In the current study, 20 qPCR batches were assayed over a period of five months for HPV detection with an average of 52 samples per batch. In each batch, samples from at least 23 different families were randomly selected. In total, there were 25 qPCR batches represented in the subsequent analyses, 20 batches from the current study and five batches from the previous study [26].

In our previous study [26], each sample was analysed two to three times (replicates) and average HPV titres were used in the statistical analyses. In the current study, we re-analysed 135 samples from the earlier study to verify the consistency of the HPV results. Quantifications of the HPV titre of the 135 samples were conducted once (no replicates). Our statistical analysis showed that there was a very high, close to 1, correlation between the two studies. Statistical data indicated that measuring each sample once was sufficient to determine the HPV titre accurately. Thus, in this study, HPV titres of 897 samples that had not been previously quantified, were obtained from a single HPV titre quantification. As a further ongoing check for consistency of HPV titres, each new batch of tests included in duplicate either three random samples selected from the previous batch, or three known samples having a low [log(HPV) ~ 3], moderate [log(HPV) ~ 5], and high HPV titre [log(HPV) ~ 8], respectively.

Statistical analyses

The shrimp HPV titres that were quantified by qPCR were transformed to the base 10 logarithm termed log(HPV). All preliminary statistics were analysed using the IBM SPSS Statistics 22 software [30]. The log(HPV) values for shrimp samples that were higher than 3 standard deviations of the mean (\({\bar{\text{x}}} \, \pm \, 3\sigma\)) were considered as outliers, and excluded from the statistical analyses. The observed \({\bar{\text{x}}} \, \pm \, 3\sigma\) of log(HPV) was 4.74 ± 3.81, ranging from the minimum value of 0.93 to the maximum of 8.55. Log(HPV) values for shrimp samples were all within the permitted range, which indicates that there were no outliers in the data, and thus no HPV titres for the tested shrimp were excluded from statistical analyses.

We used histograms to graphically represent the distribution of the HPV titre data and assess its normality. The empirical frequencies of log(HPV) followed the theoretical curve drawn from a normal distribution, which indicated that the HPV titre data was approximately normally-distributed. ANOVA analysis was conducted to examine statistical differences of log(HPV) means between the 48 families studied. The general linear model was used to investigate the relationships between fixed factors and log(HPV). The fixed factors included sampling locations around the pond, sampling time, the person who sampled the hepatopancreas, sex, DNA quality, qPCR batch and qPCR run date. Among these factors, only the effect of qPCR batch was statistically significant (P < 0.05) and was included in the model to estimate heritability of log(HPV). The proportion of log(HPV) variation explained by the qPCR batch in the statistical model was only 3.8 %.

Mixed model analyses

Heritability (h2)

The heritability of log(HPV) was estimated using a linear animal mixed model in ASReml 4.1 [31]. The model was defined as:

$${\text{Y}}_{\text{mn}} = \upmu + {\text{Q}}_{\text{m}} + {\text{a}}_{\text{n}} + {\text{e}}_{\text{mn}} ,$$
(1)

where Ymn = log(HPV), μ = estimated population mean, Qm = qPCR batch; m = 1–25, an = additive genetic effect of individual shrimp in the pedigree; n = 1 to 1137, emn = residual.

Phenotypic (rp) and genotypic (rg) correlations

The phenotypic and genotypic correlations of log(HPV) with commercial traits such as body weight, body length, head length, body width, abdominal tail weight, meat yield, raw and cooked colour, and ‘flesh streaks’, were estimated using a series of two-trait animal mixed models in ASReml 4.1 [31]. ASReml enabled the specification of different effects for different traits in the statistical models. For body and flesh quality traits, the model included the fixed effects of sampling batch, sampling time, operator and sex [28]. The effect fitted for log(HPV) was only qPCR batch as described in Eq. 1. The general form of the two-trait model was defined as:

$${\text{Y}}_{\text{ijklmn}} = \upmu + {\text{S}}_{\text{i}} + {\text{T}}_{\text{j}} + {\text{O}}_{\text{k}} + {\text{S}}_{\text{l}} + {\text{Q}}_{\text{m}} + {\text{a}}_{\text{n}} + {\text{e}}_{\text{ijklmn}} ,$$
(2)

where Yijklmn = traits under study, μ = estimated population mean, Si = sampling batch within a single grow-out pond; i = 1–8, Tj = sampling time; j = AM and PM, Ok = operator; k = 1 and 2, S1 = Sex; l = M and F, Qm = qPCR batch; m = 1–25, an = additive genetic effect of individual shrimp in the pedigree; n = 1 to 1137, eijklmn = residual.

Testing different analysis models for estimating heritability (h2)

In this study, heritability of log(HPV) was estimated by using the model that included the significant effect i.e. qPCR batch (Eq. 1). In addition, we fitted a different model that included other fixed effects as described in Eq. 2 to estimate the heritability of log(HPV). The statistical difference in heritability estimates between the two models was examined by using the z-score [32]:

$${\text{z}} = \frac{{{\text{x}}_{\text{i}} - {\text{x}}_{\text{j}} }}{{\left( {\upsigma_{\text{i}}^{2} + \upsigma_{\text{j}}^{2} } \right)^{0.5} }},$$

where xi and xj are the estimated heritabilities based on Eqs. 1 and 2, and σi and σj are their respective standard errors. The observed z-score was tested against a large-sample normal distribution (α = 0.05).

Results

The log(HPV) values averaged for each of the 48 families ranged from 3.52 ± 0.16 to 6.29 ± 0.20, and showed significant differences across families (ANOVA F47,1089 = 7.334, P < 0.001; Fig. 1). These significant differences in log(HPV) across families were also confirmed by the log-likelihood ratio test (P < 0.001). The log(HPV) of the family with the highest HPV titre was more than 2.5 fold higher than that of the family with the lowest HPV titre. There were also large differences, up to sixfold, in log(HPV) among individual shrimp within families. An example of log(HPV) variation within the family with the lowest (i.e. highest HPV resistance) and that with the highest HPV titre is in Fig. 2.

Fig. 1
figure 1

Average log(HPV) showing significant differences (P < 0.001) across the 48 full-sib shrimp families

Full size image
Fig. 2
figure 2

Variation in log(HPV) among individual shrimp within a family. Low: family with the lowest average HPV titre (i.e. the highest HPV resistance); high: family with the highest average HPV titre

Full size image

Reproducibility and correlation of qPCR results

A subset of 135 of the 240 samples that were previously assayed was re-analysed in this study. The correlation of HPV titres between the two studies was highly significant (P < 0.001) with a correlation coefficient (r) close to 1, r = 0.98, (Fig. 3). This consistency between both studies indicated a high level of accuracy in the quantification of HPV titres. In spite of the high reproducibility of HPV titres between both studies, we still expected some differences due to systematic factors such as qPCR being done by different persons, on different dates or times, and qPCR-reaction-mix preparation. In this study, we considered the combination of these factors as the effect of qPCR batch.

Fig. 3
figure 3

Scatter plot of the linear relationship of log(HPV) between the previous [26] (N = 135) and present study (N = 135). X-axis (Jan–Feb 2014) represent qPCR results from the previous study and Y-axis (Apr–Jun 2015) qPCR results for the same samples obtained in the current study

Full size image

Heritability (h2) and correlations (r) between traits

The heritability of log(HPV) was moderate, i.e. 0.41 ± 0.08 (Table 1), which indicated that this trait was heritable and could be potentially improved by selective breeding. Phenotypic (r p) and genotypic (r g) correlations of log(HPV) with body and carcass traits are in Table 2. Phenotypic correlations between log(HPV) and the other traits studied were low (−0.01 to −0.14), whereas genetic correlations were generally moderate and ranged from −0.13 to 0.5. Phenotypic and genetic correlations of log(HPV) with body and carcass traits were negative except those with meat yield which were positive. The genetic correlation between log(HPV) and yellow hepatopancreas was moderate and positive. All estimated genetic correlations had greater magnitudes than phenotypic correlations.

Table 1 Number of samples (N), mean, standard deviation (SD) and heritability (h2) for HPV titre
Full size table
Table 2 Phenotypic (r P) and genetic (r g) correlations of log(HPV) with commercial traits
Full size table

Effect of using different models on heritability (h2)

The additional analysis that included other systematic effects in the model (Eq. 2) yielded heritability estimates for log(HPV) that were very similar to those based on the model that included only the effect of qPCR batch (Eq. 1) (0.403 ± 0.080 versus 0.410 ± 0.080, z = 0.097, P > 0.05).

Discussion

Reproducibility and feasibility

Overall, we found that qPCR for HPV detection in banana shrimp was highly reproducible when the same sample was tested even months apart. Such high reproducibility of qPCR has been reported elsewhere [33, 34] and this is of practical importance since the estimation of genetic parameters requires a large number of samples, and many months of analyses.

The current work was completed over several months and faster turnaround is probably needed for applied selection programs. However, recent technologies and qPCR machines that are capable of analysing 96 samples, even up to hundreds of samples at a time, have become available, at relatively low costs; our calculations are that the same work can now be completed in weeks instead of months. DNA assays such as immunofluorescence titration and virus particle counts are fast and inexpensive methods to measure viral titres, although perhaps at the cost of accuracy [24, 25]. Moreover, should interest shift from challenge testing to viral titre as selection criterion, it is also anticipated that the identification of single nucleotide polymorphisms (SNPs) specific to full-sib families with low or high HPV titres from differential gene expression analyses [35] would be a useful alternative since SNPs represent an even lower-cost selection criterion; however, to date, quantitative trait loci (QTL) for HPV resistance have not been identified.

Substantial variation among families and moderate heritability

This is the first report, apart from our preliminary study [26], which provides an estimated heritability for viral titre in an aquaculture species. In contrast to the previously reported low estimated heritabilities for survival to virus challenge in shrimp, our estimated heritability for viral titre was moderate. It is a little premature to suggest that a moderate heritability would be the general expectation from future similar studies in different species and with different viruses. However, there are several discussable aspects of our approach, which suggest that it may yield higher heritabilities than those based on survival data from challenge testing in shrimp. For example, on the one hand, survival data may be influenced by substantial uncontrolled environmental variation to the challenge testing conditions (e.g. stock density and infection dosages due to variable feeding) [15]. We found a high level of reproducibility in the quantification of viral titre for given animals and the between-batch effects contributed to less than 5 % of the variation in viral titre. On the other hand, there is the potential for much uncontrolled environmental variation in viral titre among shrimp in the ponds, and thus, natural outbreaks may have more variable environmental effects than controlled challenge tests, but the prima facie evidence of a moderate heritability indicates that this uncontrolled variation was not sufficient to obscure genetic effects. Another advantage of using viral titre from qPCR analysis to assess resistance to virus is that shrimp samples can be directly collected at harvest times under farm condition, and transferred to the laboratory for viral quantification. This procedure of testing the animals from commercial conditions does not require challenge facilities, thus investments related to challenge testing (e.g. challenge facilities, testing and maintenance) are reduced or even completely avoided. This approach also guaranties that the assessed parameters are relevant to commercial operations, whereas there is always concern about whether dedicated challenge testing is fully relevant to commercial production. The moderate heritability of viral titre suggests that selective breeding for low viral titre is possible, and may result in shrimp that are resistant to the virus.

Not only was there substantial HPV titre variation among families, there was also substantial variation among individuals within each family, which suggests that selection should consider both within- and between-family variation.

Genetic associations between HPV titre and commercial traits

The genetic correlation between HPV titre and body weight was negative and the magnitude of the correlation was moderate (−0.33), which suggest that selection to increase HPV resistance (i.e. by reducing HPV titre) would improve growth performances (e.g. increase body weight) in banana shrimp and vice versa. Our results contrast with those of a previous study on pacific white shrimp for which highly negative genetic correlations between body weight and WSSV resistance were found for two breeding lines (−0.55 and −0.64) [15]. In marine fish breeding, genetic correlations between body weight and resistance to pathogens vary according to the infective agent. Low genetic correlations (r g ~ 0) were reported between body weight and resistance to viruses and bacteria, e.g., viral nervous necrosis (VNN), vibriosis and furunculosis in various cultured marine fishes [36–38]. However, a moderately negative genetic correlation (−0.32) was reported in Atlantic salmon between body weight and resistance to parasitic disease, C. rogercresseyi [39], while adverse genetic correlations of body weight with resistance to viral haemorrhagic septicaemia (VHS) were low to moderate (−0.14 to −0.23) in rainbow trout (Oncorhynchus mykiss) [40]. To some degree, the genetic correlations may depend on how the data are collected: growth measured on infected animals may be a different trait to growth measured on uninfected animals, and the former may rather be a measure of disease resistance than growth per se.

Improving HPV resistance in shrimp is also expected to increase some carcass characteristics (i.e. abdominal tail weight) as indicated by the negative (favourable) genetic correlations between the two traits. Besides growth performances and carcass characteristics, raw and cooked colours of banana shrimp are expected to be improved (more pronounced dark or red colour) by selecting for HPV resistance. Another finding of our study was the negative genetic correlation between HPV titre and condition of the hepatopancreas (yellow hepatopancreas). This could be a desirable trait from a commercial perspective because selection for reduction in HPV titre (i.e. increase in HPV resistance) will likely decrease the incidence of shrimp with yellow (presumably undesirable) hepatopancreas. Unfortunately, to our knowledge, genetic correlations between viral resistance and condition of the hepatopancreas have never been reported in any other aquaculture species for comparison with our findings.

To date, the nature of the defence mechanisms against viruses in shrimp is still unclear. Various studies [1–3, 41] have reported that HPV was co-expressed with other viral diseases, frequently with monodon-type baculovirus (MBV), WSSV and infectious hypodermal and hematopoietic necrosis virus (IHHNV), but there is no report on genetic correlations between HPV and other viral diseases. This is an important knowledge gap, since selection against one virus may confer resistance to others.

The favourable heritability of HPV titre and genetic correlations that we obtained in this study indicate promising opportunities to improve HPV resistance in banana shrimp without compromising other important traits.

Association of HPV titre with survival

In this study, we do not report any hard data that specifically relate HPV titre to survival. What data that are available are only circumstantial. For example, as reported above, we observed negative genetic and phenotypic correlations of HPV titre with measures of growth; this suggests that HPV titre has an impact on production traits but may not particularly relate to survival. There are data on counts of offspring per family but they likely reflect problems of differential fertility and other genetic or stochastic events in the hatchery and prior. Indeed these counts were not significantly correlated with family means of log(HPV).

One unresolved matter is the degree to which HPV titre and survival are correlated, or not, and whether such a correlation is evident only at a certain threshold titre, i.e. perhaps below this threshold there is no correlation; shrimp may have HPV but are not “sick” or “diseased” or do not show any specific clinical signs. Indeed, various studies [1, 4, 5] have reported that HPV-infected shrimp do not show gross signs of disease and act normally; only impacts on stunt growth in juveniles and mortality in severely-infected post-larvae were observed. In addition, Flegel et al. [2] reported that most black tiger shrimp from commercial grow-out ponds were positive for HPV, while no ponds showed an unusual level of mortality. Therefore, it is likely that moderately HPV-infected shrimp cannot be readily identified just by external appearances; only severely-infected shrimp may present symptoms. Manivannan et al. [1] reported that the symptoms of severe infections in juvenile shrimp were a whitish hepatopancreas, poor growth rate, anorexia, gill fouling and susceptibility to surface as a result of reduced preening activity. Based on literature data, it is hypothesized that health consequences may vary according to HPV titres, which is consistent with some industry operating views and procedures. The industry considers that ponds with high HPV titres risk catastrophic collapse when they reach the threshold at which HPV impacts health, but ponds are harvested well before the onset of collapse, as would be the case in our experiment. Thus, one can argue about whether HPV titres had reached a sufficient level for major impacts to occur on either individuals or the whole pond. Another complication in associating HPV titre with survival per se, is that HPV is often detected with other viruses [3, 41]; mortality may not reflect specifically that due to HPV since the disease may be essentially due to other viruses.

To what degree many of these uncertainties and unknowns can be best resolved by returning to more classical challenge testing is debatable. Perhaps we could test predefined HPV titres against ensuing survival; this approach presumably would identify any effects of HPV titre and relate the effects specifically to HPV rather than to other viral infections. However, we would then be returning to the swath of issues and problems with challenge testing referred to in the discussion such as costs and potential inadequacy with industry practices, etc.

In conclusion of this section, it is possible that selection for viral titre at the levels recorded here may not achieve desirable outcomes for situations of extreme viral titres and catastrophic collapse, i.e. correlation between HPV titres reported here and survival to HPV is poor as the pond evolves towards collapse. Notwithstanding these uncertainties, what we can say ipso facto, is that if we can eliminate the virus by selection, as suggested by this study, then concordantly, we will eliminate the negative effects of the virus.

Mode of transmission (vertical or horizontal)

If transmission is horizontal

There are two major types of transmission, horizontal and vertical, although the two are not mutually exclusive. Horizontal transmission refers to all possible situations in rearing environments that diseases can be transmitted from infected animals to healthy ones in the same facility. There is circumstantial evidence for horizontal transmission of HPV: Catap et al. [42] detected infections in healthy black tiger shrimp (Penaeus monodon) after feeding animals with infected black tiger shrimp postlarvae, while Sivakumar et al. [43] reported that experimentally-infected live artemia given to healthy postlarvae of black tiger shrimp apparently caused infections. Other viruses (e.g. yellow-head virus (YHV) and WSSV) were also reported to be able to transmit and infect cultured shrimp horizontally [44, 45]. All these studies suggest that viruses can be transmitted either through infected live feed or by moribund shrimp being eaten due to the cannibalistic behaviour of shrimp.

In our study, shrimp were cultured communally in the same grow-out pond, thus presumably, they had equal chances to be exposed to diseases excepting prior infections due to vertical transmission. Variation in this type of (horizontal) transmission in the pond is more like an environmental factor, and by itself should not contribute to the additive genetic variation of HPV titre and the repeatable differences among families.

If transmission is vertical

There is a possibility that vertical transmission, from parents to offspring, may have contributed to the differences of HPV titres among families. Several previous studies investigated the potential for vertical transmission, from infected parents to offspring, of different farmed shrimp viral diseases, e.g., MBV and IHHNV. All these studies confirmed the existence of a parental contribution to viral infection in the offspring, however the proportion of virus transmitted to the offspring depended on the severity of the infection in the parents [46–48].

From a scientific perspective, vertical transmission if not swamped by horizontal transmission, could lead to inflated or completely artefactual genetic parameter estimates. Unfortunately, estimation of the maternal effect was not possible with our data structure and we cannot discount that the maternal effect from vertical transmission may have merged with, and inflated, the additive genetic estimates. In other words, there is considerable risk that the heritability is biased but that our design cannot avoid this pitfall. The only publications so far for farmed fish that are communally reared at early stages show that maternal effects are negligible, at least for traits such as growth and condition index [49, 50]. Such results contrast with those from terrestrial livestock for which the offspring after birth require maternal care and, thus, maternal effects are common [51–53]. Nevertheless, until further work is done, we cannot say if the pattern of a negligible maternal effect for farmed fish that are raised communally at early stages will also apply for viral infection in shrimp. This matter can be addressed in future studies that consider the use of maternal half-sib families, and by assessing the viral titres in the female broodstock.

While we have no direct evidence for vertical transmission of HPV in shrimp, if it occurs and contributes to the differences in HPV titres among families, then on-farm selection for families with low HPV titres will still result in lower HPV titres in the next generation through a reduction in the viral titre of the parents.

Conclusions

The additive genetic variation of HPV titre suggests that resistance to HPV in banana shrimp can be improved by selective breeding. Selection for HPV titre may also improve other traits of commercial importance such as body weight, meat yield, condition of the hepatopancreas, and shrimp colour. This is the first major report on viral resistance as a quantitative trait and the genetic parameters of HPV titre, for an aquaculture species.

References

  1. Manivannan S, Otta SK, Karunasagar I, Karunasagar I. Multiple viral infection in Penaeus monodon shrimp postlarvae in an Indian hatchery. Dis Aquat Organ. 2002;48:233–6.

    Article  PubMed  Google Scholar 

  2. Flegel TW, Nielsen L, Thamavit V, Kongtim S, Pasharawipas T. Presence of multiple viruses in non-diseased, cultivated shrimp at harvest. Aquaculture. 2004;240:55–68.

    Article  Google Scholar 

  3. Joseph TC, James R, Anbu Rajan L, Surendran PK, Lalitha KV. Occurrence of viral pathogens in Penaeus monodon post-larvae from aquaculture hatcheries. Data Brief. 2015;4:170–6.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Spann KM, Adlard RD, Hudson DA, Pyecroft SB, Jones TC, Voigt MOC. Hepatopancreatic parvo-like virus (HPV) of Penaeus japonicus cultured in Australia. Dis Aquat Organ. 1997;31:239–41.

    Article  Google Scholar 

  5. Flegel TW, Thamavit V, Pasharawipas T, Alday-Sanz V. Statistical correlation between severity of hepatopancreatic parvovirus infection and stunting of farmed black tiger shrimp (Penaeus monodon). Aquaculture. 1999;174:197–206.

    Article  Google Scholar 

  6. Walker PJ, Mohan CV. Viral disease emergence in shrimp aquaculture: origins, impact and the effectiveness of health management strategies. Rev Aquacult. 2009;1:125–54.

    Article  Google Scholar 

  7. Salonius K, Simard N, Harland R, Ulmer JB. The road to licensure of a DNA vaccine. Curr Opin Investig Drugs. 2007;8:635–41.

    CAS  PubMed  Google Scholar 

  8. Drangsholt TMK, Gjerde B, Ødegård J, Finne-Fridell F, Evensen Ø, Bentsen HB. Quantitative genetics of disease resistance in vaccinated and unvaccinated Atlantic salmon (Salmo salar L.). Heredity (Edinb). 2011;107:471–7.

    Article  CAS  Google Scholar 

  9. Henryon M, Berg P, Olesen NJ, Kjær TE, Slierendrecht WJ, Jokumsen A, et al. Selective breeding provides an approach to increase resistance of rainbow trout (Onchorhynchus mykiss) to the diseases, enteric redmouth disease, rainbow trout fry syndrome, and viral haemorrhagic septicaemia. Aquaculture. 2005;250:621–36.

    Article  Google Scholar 

  10. Kjøglum S, Henryon M, Aasmundstad T, Korsgaard I. Selective breeding can increase resistance of Atlantic salmon to furunculosis, infectious salmon anaemia and infectious pancreatic necrosis. Aquac Res. 2008;39:498–505.

    Article  Google Scholar 

  11. Ødegård J, Olesen I, Dixon P, Jeney Z, Nielsen HM, Way K, et al. Genetic analysis of common carp (Cyprinus carpio) strains. II: resistance to koi herpesvirus and Aeromonas hydrophila and their relationship with pond survival. Aquaculture. 2010;304:7–13.

    Article  Google Scholar 

  12. Olesen I, Hung D, Ødegård J. Genetic analysis of survival in challenge tests of furunculosis and ISA in Atlantic salmon. Genetic parameter estimates and model comparisons. Aquaculture. 2007;272:S297–8.

    Article  Google Scholar 

  13. Fjalestad KT, Carr WC, Lotz J, Sweeney JN, Gjedrem T. Genetic variation and selection response in body weight and disease resistance in Pacific white shrimp (Penaeus vannamei). Aquaculture. 1999;173:10.

    Google Scholar 

  14. Argue BJ, Arce SM, Lotz JM, Moss SM. Selective breeding of Pacific white shrimp (Litopenaeus vannamei) for growth and resistance to Taura Syndrome Virus. Aquaculture. 2002;204:447–60.

    Article  Google Scholar 

  15. Gitterle T, Salte R, Gjerde B, Cock J, Johansen H, Salazar M, et al. Genetic (co)variation in resistance to White Spot Syndrome Virus (WSSV) and harvest weight in Penaeus (Litopenaeus) vannamei. Aquaculture. 2005;246:139–49.

    Article  Google Scholar 

  16. Gitterle T, Gjerde B, Cock J, Salazar M, Rye M, Vidal O, et al. Optimization of experimental infection protocols for the estimation of genetic parameters of resistance to White Spot Syndrome Virus (WSSV) in Penaeus (Litopenaeus) vannamei. Aquaculture. 2006;261:501–9.

    Article  Google Scholar 

  17. Ødegård J, Olesen I, Gjerde B, Klemetsdal G. Evaluation of statistical models for genetic analysis of challenge test data on furunculosis resistance in Atlantic salmon (Salmo salar): prediction of field survival. Aquaculture. 2006;259:116–23.

    Article  Google Scholar 

  18. Gianola D, Foulley J. Sire evaluation for ordered categorical data with a threshold model. Genet Sel Evol. 1983;15:201–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ducrocq V, Casella G. A Bayesian analysis of mixed survival models. Genet Sel Evol. 1996;28:505–29.

    Article  PubMed Central  Google Scholar 

  20. Meuwissen THE, Veerkamp RF, Engel B, Brotherstone S. Single and multitrait estimates of breeding value for survival using sire and animal models. Anim Sci. 2002;75:15–24.

    Article  Google Scholar 

  21. Farewell VT. The use of mixture models for the analysis of survival data with long-term survivors. Biometrics. 1982;38:1041–6.

    Article  CAS  PubMed  Google Scholar 

  22. Ødegård J, Gitterle T, Madsen P, Meuwissen THE, Yazdi MH, Gjerde B, et al. Quantitative genetics of taura syndrome resistance in pacific white shrimp (Penaeus vannamei): a cure model approach. Genet Sel Evol. 2011;43:14.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Ødegård J, Olesen I, Gjerde B, Klemetsdal G. Evaluation of statistical models for genetic analysis of challenge-test data on ISA resistance in Atlantic salmon (Salmo salar): prediction of progeny survival. Aquaculture. 2007;266:70–6.

    Article  Google Scholar 

  24. Van Koevering M, Haufler Z, Fulbright D, Isleib T, Everson E. Heritability of resistance in winter wheat to wheat spindle streak mosaic virus. Phytopathology. 1987;77:742–4.

    Article  Google Scholar 

  25. Bosio CF, Beaty BJ, Black WC 4th. Quantitative genetics of vector competence for dengue-2 virus in Aedes aegypti. Am J Trop Med Hyg. 1998;59:965–70.

    CAS  PubMed  Google Scholar 

  26. Knibb W, Quinn J, Kuballa A, Powell D, Remilton C, Nguyen NH. Yearly, pond, lineage and family variation of hepatopancreatic parvo-like virus (HPV) copy number in banana shrimp Fenneropenaeus merguiensis. J Invertebr Pathol. 2015;128:73–9.

    Article  PubMed  Google Scholar 

  27. Knibb W, Whatmore P, Lamont R, Quinn J, Powell D, Elizur A, et al. Can genetic diversity be maintained in long term mass selected populations without pedigree information? A case study using banana shrimp Fenneropenaeus merguiensis. Aquaculture. 2014;428–9:71–8.

    Article  Google Scholar 

  28. Nguyen NH, Quinn J, Powell D, Elizur A, Thoa NP, Nocillado J, et al. Heritability for body colour and its genetic association with morphometric traits in Banana shrimp (Fenneropenaeus merguiensis). BMC Genet. 2014;15:132.

    Article  PubMed  PubMed Central  Google Scholar 

  29. La Fauce KA, Layton R, Owens L. TaqMan real-time PCR for detection of hepatopancreatic parvovirus from Australia. J Virol Methods. 2007;140:10–6.

    Article  PubMed  Google Scholar 

  30. IBM Corp. IBM SPSS statistics for windows, version 22.0. Armonk: IBM Corp; 2013.

    Google Scholar 

  31. Gilmour AR, Gogel BJ, Cullis BR, Welham SJ, Thompson R. ASReml user guide release 4.1 structural specification. Hemel Hempstead: VSN International Ltd; 2015.

    Google Scholar 

  32. Nguyen NH, Khaw HL, Ponzoni RW, Hamzah A, Kamaruzzaman N. Can sexual dimorphism and body shape be altered in Nile tilapia (Oreochromis niloticus) by genetic means? Aquaculture. 2007;272:S38–46.

    Article  Google Scholar 

  33. Heid CA, Stevens J, Livak KJ, Williams PM. Real time quantitative PCR. Genome Res. 1996;6:986–94.

    Article  CAS  PubMed  Google Scholar 

  34. Hein I, Klein D, Lehner A, Bubert A, Brandl E, Wagner M. Detection and quantification of the iap gene of Listeria monocytogenes and Listeria innocua by a new real-time quantitative PCR assay. Res Microbiol. 2001;152:37–46.

    Article  CAS  PubMed  Google Scholar 

  35. Powell D, Knibb W, Nguyen NH, Elizur A. Transcriptional profiling of Banana Shrimp Fenneropenaeus merguiensis with differing levels of viral load. Integr Comp Biol. 2016. doi:10.1093/icb/icw029.

    PubMed  Google Scholar 

  36. Perry GML, Tarte P, Croisetière S, Belhumeur P, Bernatchez L. Genetic variance and covariance for 0 + brook charr (Salvelinus fontinalis) weight and survival time of furunculosis (Aeromonas salmonicida) exposure. Aquaculture. 2004;235:263–71.

    Article  Google Scholar 

  37. Silverstein JT, Vallejo RL, Palti Y, Leeds TD, Rexroad CE 3rd, Welch TJ, et al. Rainbow trout resistance to bacterial cold-water disease is moderately heritable and is not adversely correlated with growth. J Anim Sci. 2009;87:860–7.

    Article  CAS  PubMed  Google Scholar 

  38. Bangera R, Ødegård J, Præbel AK, Mortensen A, Nielsen HM. Genetic correlations between growth rate and resistance to vibriosis and viral nervous necrosis in Atlantic cod (Gadus morhua L.). Aquaculture. 2011;317:67–73.

    Article  Google Scholar 

  39. Yáñez JM, Lhorente JP, Bassini LN, Oyarzún M, Neira R, Newman S. Genetic co-variation between resistance against both Caligus rogercresseyi and Piscirickettsia salmonis, and body weight in Atlantic salmon (Salmo salar). Aquaculture. 2014;433:295–8.

    Article  Google Scholar 

  40. Henryon M, Jokumsen A, Berg P, Lund I, Pedersen PB, Olesen NJ, et al. Genetic variation for growth rate, feed conversion efficiency, and disease resistance exists within a farmed population of rainbow trout. Aquaculture. 2002;209:59–76.

    Article  Google Scholar 

  41. Umesha RK, Uma A, Otta SK, Karunasagar I, Karunasagar I. Detection by PCR of hepatopancreatic parvovirus (HPV) and other viruses in hatchery-reared Penaeus monodon postlarvae. Dis Aquat Organ. 2003;57:141–5.

    Article  CAS  PubMed  Google Scholar 

  42. Catap ES, Lavilla-Pitogo CR, Maeno Y, Travina RD. Occurrence, histopathology and experimental transmission of hepatopancreatic parvovirus infection in Penaeus monodon postlarvae. Dis Aquat Organ. 2003;57:11–7.

    Article  CAS  PubMed  Google Scholar 

  43. Sivakumar VK, Sarathi M, Venkatesan C, Sivaraj A, Hameed ASS. Experimental exposure of Artemia to hepatopancreatic parvo-like virus and subsequent transmission to post-larvae of Penaeus monodon. J Invertebr Pathol. 2009;102:191–5.

    Article  CAS  PubMed  Google Scholar 

  44. Hamano K, Miyoshi T, Aue-umneoy D, Srisapoome P, Maeno Y, Tsutsui I. Waterborne and cannibalism-mediated transmission of the Yellow head virus in Penaeus monodon. Aquaculture. 2015;437:161–6.

    Article  Google Scholar 

  45. Haryadi D, Verreth JAJ, Verdegem MCJ, Vlak JM. Transmission of white spot syndrome virus (WSSV) from Dendronereis spp. (Peters) (Nereididae) to penaeid shrimp. J Fish Dis. 2015;38:419–28.

    Article  CAS  PubMed  Google Scholar 

  46. Cowley JA, Hall MR, Cadogan LC, Spann KM, Walker PJ. Vertical transmission of gill-associated virus (GAV) in the black tiger prawn Penaeus monodon. Dis Aquat Organ. 2002;50:95–104.

    Article  PubMed  Google Scholar 

  47. Motte E, Yugcha E, Luzardo J, Castro F, Leclercq G, Rodrı́guez J, et al. Prevention of IHHNV vertical transmission in the white shrimp Litopenaeus vannamei. Aquaculture. 2003;219:57–70.

    Article  Google Scholar 

  48. Kanjanasopa D, Pongtippatee P, Vanichviriyakit R, Wongtripop S, Pradeep PJ, Withyachumnarnkul B. An evidence on trans-ovarian transmission of Monodon baculovirus (MBV) infection in Penaeus monodon. Aquaculture. 2015;443:5–10.

    Article  Google Scholar 

  49. Ninh NH, Ponzoni RW, Nguyen NH, Woolliams JA, Taggart JB, McAndrew BJ, et al. A comparison of communal and separate rearing of families in selective breeding of common carp (Cyprinus carpio): estimation of genetic parameters. Aquaculture. 2011;322–3:39–46.

    Article  Google Scholar 

  50. Whatmore P, Nguyen NH, Miller A, Lamont R, Powell D, D’Antignana T, et al. Genetic parameters for economically important traits in yellowtail kingfish Seriola lalandi. Aquaculture. 2013;400–1:77–84.

    Article  Google Scholar 

  51. Mandal A, Neser FWC, Rout PK, Roy R, Notter DR. Estimation of direct and maternal (co)variance components for pre-weaning growth traits in Muzaffarnagari sheep. Livest Sci. 2006;99:79–89.

    Article  Google Scholar 

  52. van Arendonk JAM, van Rosmeulen C, Janss LLG, Knol EF. Estimation of direct and maternal genetic (co) variances for survival within litters of piglets. Livest Prod Sci. 1996;46:163–71.

    Article  Google Scholar 

  53. Ahunu BK, Arthur PF, Kissiedu HWA. Genetic and phenotypic parameters for birth and weaning weights of purebred and crossbred Ndama and West African Shorthorn cattle. Livest Prod Sci. 1997;51:165–71.

    Article  Google Scholar 

Download references

Authors’ contributions

CP conducted the qPCR experiments, performed the statistical analyses and drafted the manuscript. NHN participated in the design and data collection and contributed to the manuscript, JQ carried out the experiment in our previous study and contributed to the dataset in the current study, WK conceived the study, participated in the design and data collection and contributed to the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We acknowledge the support from Cooperative Research Centre grant (Australian Seafood CRC: Project No. 2009/724) to Wayne Knibb. We would like to thank Prof Abigail Elizur, Dr. Josephine Nocillado, and Ngo Phu Thoa for their contributions to sample and data collection, Daniel Powell for data collection and his helpful advices during laboratory works, and Seafarm at Cardwell.

Competing interests

The authors declare that they have no competing interests.

Availability of data and material

The datasets supporting the conclusions of this article are available upon request.

Ethics approval and consent to participate

Under Australian and Queensland legislation crustaceans are exempted from requiring ethics approval.

Author information

Authors and Affiliations

  1. University of the Sunshine Coast, Locked Bag 4, Maroochydore DC, QLD, 4558, Australia

    Chontida Phuthaworn, Nguyen Hong Nguyen, Jane Quinn & Wayne Knibb

Authors
  1. Chontida Phuthaworn
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Nguyen Hong Nguyen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jane Quinn
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Wayne Knibb
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Wayne Knibb.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Phuthaworn, C., Nguyen, N.H., Quinn, J. et al. Moderate heritability of hepatopancreatic parvovirus titre suggests a new option for selection against viral diseases in banana shrimp (Fenneropenaeus merguiensis) and other aquaculture species. Genet Sel Evol 48, 64 (2016). https://doi.org/10.1186/s12711-016-0243-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-016-0243-8

Keywords

Genetics Selection Evolution

ISSN: 1297-9686