Quantitative genetics of disease resistance in vaccinated and unvaccinated Atlantic salmon (Salmo salar L.)

The genetics of resilience and its relationships with egg production traits and antibody traits in chickens

BACKGROUND

Resilience is the capacity of an animal to be minimally affected by disturbances or to rapidly return to its initial state before exposure to a disturbance. Resilient livestock are desired because of their improved health and increased economic profit. Genetic improvement of resilience may also lead to trade-offs with production traits. Recently, resilience indicators based on longitudinal data have been suggested, but they need further evaluation to determine whether they are indeed predictive of improved resilience, such as disease resilience. This study investigated different resilience indicators based on deviations between expected and observed egg production (EP) by exploring their genetic parameters, their possible trade-offs with production traits, and their relationships with antibody traits in chickens.

METHODS

Egg production in a nucleus breeding herd environment based on 1-week-, 2-week-, or 3-week-intervals of two purebred chicken lines, a white egg-laying (33,825 chickens) and a brown egg-laying line (34,397 chickens), were used to determine deviations between observed EP and expected average batch EP, and between observed EP and expected individual EP. These deviations were used to calculate three types of resilience indicators for two life periods of each individual: natural logarithm-transformed variance (ln(variance)), skewness, and lag-one autocorrelation (autocorrelation) of deviations from 25 to 83 weeks of age and from 83 weeks of age to end of life. Then, we estimated their genetic correlations with EP traits and with two antibody traits.

RESULTS

The most promising resilience indicators were those based on 1-week-intervals, as they had the highest heritability estimates (0.02-0.12) and high genetic correlations (above 0.60) with the same resilience indicators based on longer intervals. The three types of resilience indicators differed genetically from each other, which indicates that they possibly capture different aspects of resilience. Genetic correlations of the resilience indicator traits based on 1-week-intervals with EP traits were favorable or zero, which means that trade-off effects were marginal. The resilience indicator traits based on 1-week-intervals also showed no genetic correlations with the antibody traits, which suggests that they are not informative for improved immunity or vice versa in the nucleus environment.

CONCLUSIONS

This paper gives direction towards the evaluation and implementation of resilience indicators, i.e. to further investigate resilience indicator traits based on 1-week-intervals, in breeding programs for selecting genetically more resilient layer chickens.

Host genetic variation explains reduced protection of commercial vaccines against Piscirickettsia salmonis in Atlantic salmon

Vaccination is a widely used control strategy to prevent Piscirickettsia salmonis causing disease in salmon farming. However, it is not known why all the currently available commercial vaccines generally fail to protect against this pathogenic bacteria. Here, we report, from two different populations, that between-family variation is a strong intrinsic factor that determines vaccine protection for this disease. While in some full-sib families, the protection added by vaccination increased the survival time in 13 days in comparison with their unvaccinated siblings; in other families, there was no added protection by vaccination or even it was slightly negative. Resistance to P. salmonis, measured as days to death, was higher in vaccinated than unvaccinated fish, but only a moderate positive genetic correlation was obtained between these traits. This disputes a previous hypothesis, that stated that both traits were fully controlled by the same genes, and challenges the use of unvaccinated fish as gold standard for evaluating and selecting fish resistant to P. salmonis, particularly if the offspring will be vaccinated. More studies are necessary to evaluate if variation in the host immune response to vaccination could explain the between-family differences in resistance observed in vaccinated fish.

Body Weight Deviations as Indicator for Resilience in Layer Chickens

Resilience is the capacity of an animal to be minimally affected by disturbances or to rapidly return to the state pertained before exposure to a disturbance. Less resilient animals are expected to be more susceptible to environmental perturbations, such as diseases, and will consequently show more and/or greater fluctuations in production than more resilient animals. Natural antibodies (NAb) are antibodies recognizing antigens without previous exposure to these, and are hypothesized to be an indication of general disease resistance. The objective of this research was to investigate genetic parameters of resilience indicators based on standardized body weight (BW) deviations and to investigate its relation with immunity (i.e. NAb) and disease resistance. Keyhole limpet hemocyanin-binding NAb were measured in layer chickens, which were selectively bred for high and low keyhole limpet hemocyanin-binding NAb levels during six generations. In addition, BW data of these layers were collected on a four-weekly interval from 4 weeks of age until 32 weeks of age. Standardized deviations of BW from an individual were compared to lines’ average BW (i.e. across individuals), and these were used to calculate resilience indicators: natural logarithm-transformed variance [ln(variance)], skewness, and lag-one autocorrelation of deviations (i.e. all within an individual). Heritabilities of resilience indicators were between 0.09 and 0.11. Genetic correlations between the three resilience indicators were between -0.20 and 0.40 (with high SE), which might suggest that the resilience indicators capture different aspects of resilience. Genetic correlations between resilience indicators and NAb were close to zero, which suggests that the resilience indicators and NAb capture different aspects of immunity. This might indicate that, in this dataset, environmental perturbations are only to a small extent affected by disease incidence, possibly due to a lack of disease occurrence. However, a lower estimated breeding value for ln(variance) was predictive for lower lesion scores after an avian pathogenic Escherichia coli inoculation and vice versa. In conclusion, this study shows that there is genetic variation in resilience indicators based on BW deviations in layer chickens, which opens up possibilities to improve resilience by means of selective breeding.

Variance and covariance estimates for resistance to bacterial cold water disease and columnaris disease in two rainbow trout breeding populations1

Family-based selective breeding can be an effective strategy for controlling diseases in aquaculture. This study aimed to estimate (co)variance components for resistance to bacterial cold water disease (BCWD) and columnaris disease (CD) in two unrelated rainbow trout nucleus breeding populations: the USDA, ARS, National Center for Cool and Cold Water Aquaculture odd-year line (ARS-Fp-R), which has been subjected to five generations of selection for improved resistance to BCWD, and the Troutlodge, Inc., May-spawning odd-year line (TLUM), which has been selected for improved growth performance but not for disease resistance. A total of 46,805 and 27,821 pedigree records were available from both populations, respectively. Between 44 and 138 families per generation and population were evaluated under controlled BCWD and CD challenges, providing 32,311 and 17,861 phenotypic records for BCWD resistance, and 13,603 and 9,413 for CD resistance, in the ARS-Fp-R and TLUM populations, respectively. A two-trait animal threshold model assuming an underlying normal distribution for the binary survival phenotypes was used to estimate (co)variance components separately for each population. Resistance to BCWD (h2 = 0.27 ± 0.04 and 0.43 ± 0.08) and CD (h2 = 0.23 ± 0.07 and 0.34 ± 0.09) was moderately heritable in the ARS-Fp-R and TLUM populations, respectively. The genetic correlation between the resistance to BCWD and CD was favorably positive in the ARS-Fp-R (0.40 ± 0.17) and TLUM (0.39 ± 0.18) populations. These findings suggest that both disease resistance traits can be improved simultaneously even if genetic selection pressure is applied to only one of the two traits.

Genetically influenced resistance to stress and disease in salmonids in relation to present-day breeding practice - a short review

While intensive fish production has many advantages, it also has a number of drawbacks as regards disease and stress. To date, there has been no conclusive review of disease resistance at Czech fish farms. The aim of the study was to describe briefly the existing salmonid breeding practice in the Czech Republic and to point out the trends and new possibilities gaining ground around Europe. However, the present situation in the Czech stocks is not rare at all and therefore it is used here as a model example representing numerous breeding practices in Europe. Stress and disease resistance in fish is polygenic and quantitative, making selection for such traits difficult. In recent years, however, fish breeding methods have developed rapidly, with the use of genetic analysis tools, for example, now allowing much greater selection accuracy. Gradual progress in understanding the importance of individual genetic markers offers many new options that can be utilised in breeding practice. New selection methods, such as quantitative trait loci (QTLs) and genomic selection, are increasingly employed in European aquaculture. Next generation sequencing techniques now help in the finding of new and promising QTLs that can be used in assisted selection. This review maps the current progress in improving salmonid resistance to stress and disease in aquaculture and at the same time provides the breeders with a short overview of the latest tools of genetically controlled breeding and of the newest products available at the European market.

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

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 (h(2) = 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.

Rainbow trout (Oncorhynchus mykiss) resistance to columnaris disease is heritable and favorably correlated with bacterial cold water disease resistance

Columnaris disease (CD), caused by Flavobacterium columnare, is an emerging disease affecting rainbow trout aquaculture. Objectives of this study were to 1) estimate heritability of CD resistance in a rainbow trout line (ARS-Fp-R) previously selected 4 generations for improved bacterial cold water disease (BCWD) resistance; 2) estimate genetic correlations among CD resistance, BCWD resistance, and growth to market BW; and 3) compare CD resistance among the ARS-Fp-R, ARS-Fp-S (selected 1 generation for increased BCWD susceptibility), and ARS-Fp-C (selection control) lines. Heritability of CD resistance was estimated using data from a waterborne challenge of 44 full-sib ARS-Fp-R families produced using a paternal half-sib mating design, and genetic correlations were estimated using these data and 5 generations of BCWD resistance, 9-mo BW (approximately 0.5 kg), and 12-mo BW (approximately 1.0 kg) data from 405 ARS-Fp-R full-sib families. The CD and BCWD challenges were initiated at approximately 52 and 84 d posthatch, or approximately 650 and 1,050 degree days (°C × d), respectively. Survival of ARS-Fp-R families ranged from 0 to 48% following CD challenge and heritability estimates were similar between CD (0.17 ± 0.09) and BCWD (0.18 ± 0.03) resistance, and the genetic correlation between these 2 traits was favorable (0.35 ± 0.25). Genetic correlations were small and antagonistic (-0.15 ± 0.08 to -0.19 ± 0.24) between the 2 resistance traits and 9- and 12-mo BW. Two challenges were conducted in consecutive years to compare CD resistance among ARS-Fp-R, ARS-Fp-C, and ARS-Fp-S families. In the first challenge, ARS-Fp-R families (83% survival) had greater CD resistance than ARS-Fp-C (73.5%; P = 0.02) and ARS-Fp-S (68%; P < 0.001) families, which did not differ (P = 0.16). In the second challenge, using an approximately 2.5-fold greater challenge dose, ARS-Fp-R families exhibited greater CD resistance (56% survival) than ARS-Fp-S (38% survival; P = 0.02) families. The favorable genetic correlation between CD and BCWD resistance is supported by greater CD resistance of the ARS-Fp-R line compared to the ARS-Fp-C and ARS-Fp-S lines and suggests that both traits will be improved simultaneously when selection is practiced on only 1 trait. In summary, these data indicate the feasibility of further selective breeding of the BCWD-resistant ARS-Fp-R line for increased CD resistance to produce a double pathogen-resistant line of rainbow trout.

Mapping and validation of a major QTL affecting resistance to pancreas disease (salmonid alphavirus) in Atlantic salmon (Salmo salar)

Pancreas disease (PD), caused by a salmonid alphavirus (SAV), has a large negative economic and animal welfare impact on Atlantic salmon aquaculture. Evidence for genetic variation in host resistance to this disease has been reported, suggesting that selective breeding may potentially form an important component of disease control. The aim of this study was to explore the genetic architecture of resistance to PD, using survival data collected from two unrelated populations of Atlantic salmon; one challenged with SAV as fry in freshwater (POP 1) and one challenged with SAV as post-smolts in sea water (POP 2). Analyses of the binary survival data revealed a moderate-to-high heritability for host resistance to PD in both populations (fry POP 1 h²~0.5; post-smolt POP 2 h²~0.4). Subsets of both populations were genotyped for single nucleotide polymorphism markers, and six putative resistance quantitative trait loci (QTL) were identified. One of these QTL was mapped to the same location on chromosome 3 in both populations, reaching chromosome-wide significance in both the sire- and dam-based analyses in POP 1, and genome-wide significance in a combined analysis in POP 2. This independently verified QTL explains a significant proportion of host genetic variation in resistance to PD in both populations, suggesting a common underlying mechanism for genetic resistance across lifecycle stages. Markers associated with this QTL are being incorporated into selective breeding programs to improve PD resistance.

Genetic variation in bacterial kidney disease (BKD) susceptibility in Lake Michigan Chinook Salmon and its progenitor population from the Puget Sound

Mass mortality events in wild fish due to infectious diseases are troubling, especially given the potential for long-term, population-level consequences. Evolutionary theory predicts that populations with sufficient genetic variation will adapt in response to pathogen pressure. Chinook Salmon Oncorhynchus tshawytscha were introduced into Lake Michigan in the late 1960s from a Washington State hatchery population. In the late 1980s, collapse of the forage base and nutritional stress in Lake Michigan were thought to contribute to die-offs of Chinook Salmon due to bacterial kidney disease (BKD). Previously, we demonstrated that Lake Michigan Chinook Salmon from a Wisconsin hatchery have greater survival following BKD challenge relative to their progenitor population. Here, we evaluated whether the phenotypic divergence of these populations in BKD susceptibility was due to selection rather than genetic drift. Comparison of the overall magnitude of quantitative trait to neutral marker divergence between the populations suggested selection had occurred but a direct test of quantitative trait divergence was not significant, preventing the rejection of the null hypothesis of differentiation through genetic drift. Estimates of phenotypic variation (VP ), additive genetic variation (VA ) and narrow-sense heritability (h (2)) were consistently higher in the Wisconsin relative to the Washington population. If selection had acted on the Wisconsin population there was no evidence of a concomitant loss of genetic variation in BKD susceptibility. The Renibacterium salmoninarum exposures were conducted at both 14°C and 9°C; the warmer temperature accelerated time to death in both populations and there was no evidence of phenotypic plasticity or a genotype-by-environment (G × E) interaction. High h (2) estimates for BKD susceptibility in the Wisconsin population, combined with a lack of phenotypic plasticity, predicts that future adaptive gains in BKD resistance are still possible and that these adaptive gains would be stable under the temperature range evaluated here.

Genetic heterogeneity of within-family variance of body weight in Atlantic salmon (Salmo salar)

BACKGROUND

Canalization is defined as the stability of a genotype against minor variations in both environment and genetics. Genetic variation in degree of canalization causes heterogeneity of within-family variance. The aims of this study are twofold: (1) quantify genetic heterogeneity of (within-family) residual variance in Atlantic salmon and (2) test whether the observed heterogeneity of (within-family) residual variance can be explained by simple scaling effects.

RESULTS

Analysis of body weight in Atlantic salmon using a double hierarchical generalized linear model (DHGLM) revealed substantial heterogeneity of within-family variance. The 95% prediction interval for within-family variance ranged from ~0.4 to 1.2 kg2, implying that the within-family variance of the most extreme high families is expected to be approximately three times larger than the extreme low families. For cross-sectional data, DHGLM with an animal mean sub-model resulted in severe bias, while a corresponding sire-dam model was appropriate. Heterogeneity of variance was not sensitive to Box-Cox transformations of phenotypes, which implies that heterogeneity of variance exists beyond what would be expected from simple scaling effects.

CONCLUSIONS

Substantial heterogeneity of within-family variance was found for body weight in Atlantic salmon. A tendency towards higher variance with higher means (scaling effects) was observed, but heterogeneity of within-family variance existed beyond what could be explained by simple scaling effects. For cross-sectional data, using the animal mean sub-model in the DHGLM resulted in biased estimates of variance components, which differed substantially both from a standard linear mean animal model and a sire-dam DHGLM model. Although genetic differences in canalization were observed, selection for increased canalization is difficult, because there is limited individual information for the variance sub-model, especially when based on cross-sectional data. Furthermore, potential macro-environmental changes (diet, climatic region, etc.) may make genetic heterogeneity of variance a less stable trait over time and space.

Survival, growth and sexual maturation in Atlantic salmon exposed to infectious pancreatic necrosis: a multi-variate mixture model approach

Background

Outbreaks of infectious pancreatic necrosis (IPN) in Atlantic salmon can result in reduced growth rates in a fraction of the surviving fish (runts). Genetic and environmental variation also affects growth rates within different categories of healthy animals and runts, which complicates identification of runts. Mixture models are commonly used to identify the underlying structures in such data, and the aim of this study was to develop Bayesian mixture models for the genetic analysis of health status (runt/healthy) of surviving fish from an IPN outbreak.

Methods

Five statistical models were tested on data consisting of 10 972 fish that died and 3959 survivors with recorded growth data. The most complex models (4 and 5) were multivariate normal-binary mixture models including growth, sexual maturity and field survival traits. Growth rate and liability of sexual maturation were treated as two-component normal mixtures, assuming phenotypes originated from two potentially overlapping distributions, (runt/normal). Runt status was an unobserved binary trait. These models were compared to mixture models with fewer traits (Models 2 and 3) and a classical linear animal model for growth (Model 1).

Results

Assuming growth as a mixture trait improved the predictive ability of the statistical model considerably (Model 2 vs. 1). The final models (4 and 5) yielded the following results: estimated (underlying) heritabilities were moderate for growth in healthy fish (0.32 ± 0.04 and 0.35 ± 0.05), runt status (0.39 ± 0.07 and 0.36 ± 0.08) and sexual maturation (0.33 ± 0.05), and high for field survival (0.47 ± 0.03 and 0.48 ± 0.03). Growth in healthy animals, runt status and survival showed consistent favourable genetic associations. Sexual maturation showed an unfavourable non-significant genetic correlation with runt status, but favourable genetic correlations with other traits. The estimated fraction of healthy fish was 81-85%. The estimated breeding values for runt status and (normal) growth were consistent for the most complex models (4 and 5), but showed imperfect correlations with estimated breeding values from the simpler models.

Conclusions

Modelling growth in IPN survivors as a mixture trait improved the predictive ability of the model compared with a classical linear model. The results indicated considerable genetic variation in health status among survivors. Mixture modelling may be useful for the genetic analysis of diseases detected mainly through indicator traits.