- Open Access
Parental genetic diversity of brown trout (Salmo trutta m. fario) brood stock affects offspring susceptibility to whirling disease
Parasites & Vectorsvolume 8, Article number: 141 (2015)
Whirling disease, caused by the myxozoan parasite Myxobolus cerebralis, has high economical and ecological importance worldwide. Susceptibility to the disease varies considerably among salmonid species. In brown trout (Salmo trutta) the infection is usually subclinical with low mortality, which increases the risk of parasite dissemination, especially when farm fish are used for stocking natural habitats. The influence of intraspecific genetic differences (especially the level of homozygosity) on susceptibility is unknown. Therefore, we examined the possible correlations between parental genetic diversity and offspring susceptibility of brown trout stocks to whirling disease.
Two brown trout brood stocks from a German and a Hungarian fish farm were genetically characterized using microsatellite and lineage-specific genetic markers. The individual inbreeding coefficient f and pairwise relatedness factor r were estimated based on eight microsatellite markers. Brood stock populations were divided into groups according to low and high f and r value estimates and subjected to selective fertilization. The offspring from these separate groups were exposed to M. cerebralis actinospores, and the infection prevalence and intensity was measured and statistically analysed.
The analysis of phylogeographic lineage heritage revealed high heterogeneity in the Hungarian brood stock since > 50% of individuals were Atlantic-Danubian hybrids, while only pure Atlantic-descending specimens were detected in the German population. Based on f msat and r msat estimations, classified non-inbred (NIB), inbred (IB) and a group of closely related fish (REL) were created. The susceptibility of their offspring varied considerably. Although there was no significant difference in the prevalence of M. cerebralis infection, the mean intensity of infection differed significantly between NIB and IB groups. In REL and IB groups, a high variability was observed in infection intensity. No external clinical signs were observed in the exposed brown trout groups.
Our findings indicate that the allelic diversity of brown trout brood stock may constitute a significant factor in disease susceptibility, i.e. the intensity of parasite infection in the subsequent generation.
The effects of inbreeding on salmonid reproduction, growth and survival have an extensive record in scientific literature. Numerous studies support the thesis of fitness-reducing negative effects, which may occur in fish with an increase in inbreeding [1-3]. The reasons for elevated inbreeding are usually the selection of less than optimal (and thus more closely related) individuals in artificial reproduction, the selection of certain desired phenotypes, the selection of spawning fish only during short sections of the spawning period, the use of only few female spawners due to sufficient amounts of eggs [4-6], the mixing of gametes of different brood stock animals [7,8], as well as differences in survival rates or sole re-use of certain parts of the brood stock population [9-12].
Besides these, selective breeding is often utilized for controlling economically relevant diseases. Through propagation of disease-survivors as brood stock, or the genetic manipulation of fish stocks, the disease resistance of fish populations can be increased [13-16]. In salmonids, resistant strains were developed against common diseases such as furunculosis  and viral haemorrhagic septicaemia (VHS) . Ceratonova (syn Ceratomyxa) shasta was the first myxozoan parasite to which resistant strains were developed and evaluated [19,20]. In genetic studies on various rainbow trout (Oncorhynchus mykiss) strains with different susceptibility, several genome regions were found to be involved in resistance to C. shasta. The complex, multigenic trait of disease resistance has been confirmed for another myxozoan parasite, Myxobolus cerebralis [21-23]. The parasite responsible for whirling disease causes serious declines in wild and farmed salmonid populations worldwide. M. cerebralis infects a range of salmonid hosts that vary in susceptibility [24-27]. Fetherman et al.  showed that the introduction of a “M. cerebralis-resistant” rainbow trout strain to a natural habitat contributed to the reduction of infection prevalence and disease severity despite the low survival rate of stocked fish. Fish age and size also affect parasite development. The susceptibility of salmonid fry to the parasite decreases with age and growth [24,29-31]. Surprisingly, resistance is not associated with the level of skeletal ossification, but rather with other age- and size-related factors, such as the stage of development of the central nervous system . Brown trout (Salmo trutta) is assumed to be the original host, which evolved with the parasite , and is considered one of the least susceptible salmonid host species. Infection in brown trout usually proceeds without clinical signs and modest to absent mortality. However, subclinical infections may increase the risk of the dissemination of the parasite into natural environments, especially in Europe, where hatchery-reared brown trout are used for stocking natural habitats, as a major part of conservation efforts . Therefore, the role of reservoir hosts needs to be considered carefully and factored into management plans and epidemics as suggested for the monogenean ectoparasite Gyrodactylus salaris .
In the present study, we examined the correlations between parental genetic diversity and susceptibility traits of brown trout to the whirling disease parasite, M. cerebralis, for the first time. Marker heterozygosity-based grouping of brood stock was conducted on the basis of microsatellite and lineage-derived genetic markers, and the susceptibility of the offspring from assorted parental fish was evaluated in experimental infection trials in order to measure the effect of different homozygosity levels on the susceptibility of the resulting offspring.
Brood stock tagging and sampling
Fish were sampled in a Hungarian and a German trout hatchery in October, 2010, and in June, 2011. All examined fish were individually tagged with PIT tags (Loligo Systems). From the brown trout (BT) brood stock of Lillafüred trout hatchery in Hungary (48°6′59.22"N, 20°34′ 46.21"E), 167 specimens of 3–4 year-old fish were examined. From the brood stock of the trout hatchery in Aufseß, Germany (49°52′47.01"N, 11°13′41.43"E), 195 individuals (approx. equal number of males and females) from the same age group were sampled. During sampling, fin clips of approximately 1 cm2 were cut from the caudal fin of all specimens, photos were taken of every fish, and the total length was also recorded. Fin clips were fixed in 70% ethanol and stored at +4°C for molecular analysis.
DNA extraction and lineage analysis
For DNA extraction, an approximately 5 mm2 fin clip cut was dried in a Speed Vac Concentrator (Savant), and homogenized in a 1.5 ml microcentrifuge tube with a sterile pestle (Eppendorf) in ultrapure MilliQ water. Extraction was carried out as described by Estoup et al. . Briefly, 500 μl 10% Chelex 100 resin solution (BioRad) containing 0.3 mg/ml proteinase K was added to the homogenate and incubated under constant shaking at 55°C for 1 h. Digestion was terminated by heating to 100°C for 15 min. DNA extracts were stored at −20°C until further use. An approximately 1088 bp fragment of the mitochondrial DNA control region (mtDNA CR) was PCR-amplified using the primer pair RiBa (5’-CAC CCT TAA CTC CCA AAG CTA AG-3’)  and HN20 (5’-GTG TTA TGC TTT AGT TAA GC-3’) . From nuclear DNA, a 428 bp fragment of the lactate-dehydrogenase gene C1 region (LDH-C1) was amplified with the primer pair Ldhxon3F (5’-GGC AGC CTC TTC CTC AAA ACG CCC AA-3’) and Ldhxon4R (5’-CAA CCT GCT CTC TCC CTC CTG CTG ACG AA-3’) . The total volume of PCR reactions was 25 μl, which contained 15–40 ng DNA (0.5 μl), 1× Taq PCR reaction buffer (Fermentas, Thermo Scientific), 1.5 mM MgCl2, 0.2 mM dNTP mix (Sigma), 0.5 μM of each primer and 1.25 units of recombinant Taq DNA Polymerase (Fermentas, Thermo Scientific). For PCR conditions, the protocol described by Marić et al.  was followed.
The PCR-restriction fragment length polymorphism (PCR-RFLP) was conducted as per Marić et al. . The amplified PCR products of mtDNA CR were digested with the restriction enzyme SatI (Fermentas, Thermo Scientific), while BseLI (Fermentas, Thermo Scientific) was used for LDH-C1 PCR products. The RFLP pattern was visualized on a 1.5% agarose gel in 0.5× TAE buffer stained with GelRed nucleic acid stain (Biotium). Photos were taken using a Gel Logic 212 Imaging System (KODAK).
Eight microsatellite (msat) loci were examined on the tagged brood stock fish populations (Table 1). The DNA of msat Str-15, Str-60, Ssa-85 and Ssa-197 were amplified using a tetraplex PCR assay developed and optimized in the present study. The 25 μl total volume of PCR reaction contained approx. 50 ng template DNA, 1× AmpliTaq Gold PCR buffer (Life Technologies), 3 mM MgCl2, 40 μM dNTP (Sigma), 0.5 μM (for Str-15)/0.25 μM (for Ssa-197)/0.125 μM (for Ssa-85 and Str-60) of each primer and 1 u AmpliTaq Gold Taq DNA polymerase (Life Technologies). A duplex PCR to amplify msat Sso-197 and Str-543 was performed under similar conditions with 0.5 μM of each primer. Microsatellites SsoSL-438 and OKI-10 were amplified with simplex PCR assays containing half of the amount of MgCl2 and dNTP as mentioned above (1.5 mM and 20 μM, respectively). The PCR program started with an initial denaturation at 94°C for 5 min, followed by 35 cycles of 94°C/30 s, 55 or 60°C as indicated in Table 1 for 30 s and 72°C/30 s. The final elongation step was performed at 72°C for 3 min.
The size of the msat loci were estimated with DNA fragment analysis on ABI 3100 Genetic Analyzer on POP6 polymer at 60°C using a GeneScan500 ROX size standard. Then, msat allele size detection was conducted with the software PeakScanner v1.0 (Life Technologies). Raw size values were normalized based on the length of the nucleotide (nt) repeats in the msat loci; e.g. for 2 nt repeats “CT” in Str-15, the size 199 was corrected to 200 as 1 nt difference is shorter than one 2 nt-repeat, therefore the two size values were treated identically (as they belong to the same number of 2-nt repeats).
Pairwise relatedness (r msat) and individual inbreeding coefficients (f msat) were estimated using the triadic likelihood estimator (TrioML) in the software COANCESTRY by J. Wang .
Grouping of fish for selective fertilization
On the basis of the obtained data regarding the genetic status of BT brood stock, four groups were selected. Classified non-inbred (NIB) and inbred (IB) groups of parent fish in Aufseß were distinguished based on f msat values. One group composed of a dozen NIB individuals of genetically related male and female parent fish (REL) was identified for individual fertilization based on r msat values. Four male and female individuals from the Atlantic-Danubian hybrid (LBT) brood stock group were used for selective fertilization in the Lillafüred trout hatchery.
Infection trial with M. cerebralis
Offspring of BT brood stock groups IB, NIB and REL were obtained from the trout hatchery in Aufseß, Germany, whereas LBT and Steelhead strain rainbow trout (Oncorhynchus mykiss; LRBT) fingerlings originated from the Lillafüred trout hatchery, Hungary (Table 2). Fish were kept in a parasite-free environment in both hatcheries before transport. Two month old fry were transported to the laboratory. They were kept in flow-through aquaria at 15 ± 2°C water temperature, and were fed with commercial trout food. The fish were exposed at the age of 3 months.
Infective Myxobolus cerebralis triactinomyxon spores (TAMs) were obtained from our in vivo laboratory oligochaete cultures. TAMs were harvested by filtering the water from the culture containers through 20 μm mesh. TAMs used for the infection trials were less than 48 h old [47,48]. Fish were exposed individually to 3000 freshly filtered TAMs/fish in 500 ml dechlorinated tap water at 15°C for 3 h. Fifty-three individuals per group were exposed. Non-exposed control fish were kept under the same conditions with the exception of TAMs. After individual exposure, fish groups were kept separately in flow-through aquaria and fed daily. The infection experiment was terminated 4 months post exposure. Fish were euthanized with 200 mg.L−1 tricaine-methanesulfonate (MS222, Sigma), and kept frozen at −20°C until further examination.
Fish were decapitated and whole heads were divided horizontally. The cartilage and minced skeletal elements of a half head (including gill arches) were transferred into a 2 ml microcentrifuge tube containing a metal ball, and after the addition of 800 μl distillated water, the sample was homogenized with a TissueLyser LT (Qiagen) at 50 Hz for 10 min. Instead of quantifying parasite DNA of all developmental stages using species-specific PCR, morphologically intact M. cerebralis myxospores were counted which are transmission stages and play a key role in the propagation of the parasite. Spores were counted in 20 μl tissue homogenate on a microscopic slide with 20 × 20 mm coverslip under a light microscope (Zeiss Axiostar Plus) in three replicates.
Statistical analyses were performed using the R program for Windows (R commander 2.15.1). A chi-square test was used for analyzing the difference in infection prevalence among BT groups. Difference of mean infection intensities among groups were tested by one-way ANOVA. Pairwise comparisons were made by Dunnett’s test with IB as the reference group. Non-infected individuals were excluded from the statistical analysis of infection intensity as suggested by Rózsa et al. .
We declare that the treatment of fish complied with the relevant Hungarian legislation (Section 49 of Act No. XXVIII/1998 on the protection and preservation of animals, and Executive decree No. 40/2013) and the Directive 2010/63/EU on the protection of animals used for scientific purposes and the guidelines and recommendations of the Federation of Laboratory Animal Science Associations.
As a result of PCR-RFLP, one genotype was detected using mtDNA CR, and three genotypes with LDH-C1. For all examined samples, the PCR-RFLP pattern characteristic of the Atlantic (At) lineage was detected, as the amplified fragment of mtDNA CR was cut into two fragments, a 654 and a 434 bp long one. The pure Danubian (Da) genotype, whose amplified mtDNA CR fragment is not cut by the restriction endonuclease SatI, was not present in the collected fish samples. In the German brown trout population, only the Atlantic allele was detected with the PCR-RFLP of LDH-C1. In every examined fish sample, the endonuclease BseLI cut the 428 bp LDH-C1 fragment at one restriction site, resulting in 353 and 75 bp fragments. Thus, 100% of the examined fish from Germany belonged to the LDH-C1_At-At (Atlantic) genotype.
In the Hungarian brood stock population, both Atlantic (LDH-C1_At) and Danubian (LDH-C1_Da) LDH-C1 alleles were present, and the PCR-RFLP patterns of At (LDH-C1_At-At), Da (LDH-C1_Da-Da) and At-Da hybrid (Hyb) genotypes (LDH-C1_At-Da) were distinguishable (Table 3). Half of the Hungarian fish population showed the hybrid LDH-C1 pattern, while the genotype Da was detected in 25 cases of 144 (17%). All of the fish having LDH-C1_Da-Da (Da) showed the At pattern using mtDNA, which might be explained by back-crossing of hybrid parents. The gender difference was recognizable since more females than males (23% and 12.5% respectively) belonged to the Da genotype, which was the opposite for the hybrids (Table 3).
The size range of the examined msat loci is displayed in Table 1. The highest allele size variability was detected for OKI-10, for which 23 different msat length values have been observed in the examined brood stock. The slightest length variation was detected for Str-15 and Str-60 (4 different msat lengths each).
The estimated f msat values varied between 0 and 0.2091 for females, and 0 and 0.6197 for males. The examined brood stock specimens were divided into groups with relatively low and high individual inbreeding coefficients with thresholds resulting from the f msat estimation. These were f msat > 0.05 for more homozygous (IB) and f msat < 0.02 for more heterozygous (NIB) females and f msat > 0.1 for more homozygous (IB) and f msat < 0.02 for more heterozygous (NIB) males. The latter threshold resulted from significantly higher individual f msat values for NIB males, while fish with intermediate values were excluded from the experiments. The frequency distribution of f msat values showed that while more than 40% of males could be classified as inbred based on our markers, only about 30% of females had f msat values that made them candidates for the IB group (Figure 1).
The assignment of closely related males to selected females from the NIB group (all individuals’ f msat < 0.02) for experimental mating of related and non-related parents (paired samples from single female egg batches) was based on r msat estimation. Non-related parents had an r msat < 0.0002, close-sib (highly related) pairs had a r msat > 0.2. The r msat values of male–female pairs considered related (REL) varied between 0.2208 and 0.4481.
The inbreeding coefficient f msat of randomly selected NIB and IB offspring individuals (ten of each group) has also been estimated and showed high variance within the groups as expected (NIB: 0.0002 - 0.3414; IB: 0.0004 - 0.1353).
Infection trial with M. cerebralis
The exposed BT groups showed no clinical signs of infection. The infection prevalence was highest (91%) in the positive control group (LRBT) with moderate clinical signs in the affected fish specimens (blackened tail, growth retardation and slight spinal deformations). In the non-exposed, negative control fish, M. cerebralis myxospores were not detected. The prevalence of infection varied between 44 and 69% among BT groups, but the difference was not significant between any of the BT groups (Figure 2). However, the intensity of M. cerebralis infection showed a more remarkable difference. The observed mean spore number was highest for the IB group (274 ± 344.54; mean spore number ± standard deviance, SD) and lowest for LBT (52 ± 85.41). The group NIB and REL were mid-ranged regarding mean spore numbers (91 ± 126.98 and 174 ± 274.72, respectively), while in the LRBT group an extreme variance in spore numbers was observed (8517 ± 25793.39). When comparing all values to that of the group IB as a reference, the infection intensity of the NIB and the At-Da hybrid group, LBT was significantly lower than that of the reference group (p < 0.005). The related group, REL did not differ significantly from the reference IB group (Figure 3).
Inbreeding reduces the rates of reproduction and survival in most regularly outbred species . But inbreeding depression can not only lead to a strong reduction of growth in salmonids . According to other studies, inbreeding causes increased embryo mortality, and affects spawning age and egg counts of female rainbow trout [4,52-54]. Interestingly, inbred males had no effect on hatching success, whereas inbreeding in female spawners had fatal consequences for developing eggs and caused poor hatching rates . Our findings support these observations, as more IB females failed to ovulate in the spawning period than in the NIB group.
Until now, heterozygosity-fitness correlations remain in criticism. The subject of controversy is the question whether relationships between degree of homozygosity and observed fitness parameters are detectable and comprehensible at all. According to Balloux et al.  increased heterozygosity is a rather poor indicator for genome-wide inbreeding and also a large number of markers can reflect the inbreeding coefficient f only under extreme conditions of kinship mating. Even pedigree data, as preferred by most authors, computationally tend to underestimate the real degree of inbreeding. The decisive advantage of the marker-based approach is the possibility to investigate the existing inbreeding (−depression) in animals with identical origin (and thus same expected f) on the basis of individual values. To circumvent the mentioned problems, especially when using microsatellite data, Wang  proposed a new probability estimator, TrioML which is better correlated with fitness parameters and weighted for allelic diversity of a given population. Therefore, in the present study, the TrioML estimator was chosen for the calculation of f msat to give a more reliable estimation of the homozygosity level of the examined brood stock.
The signs of whirling disease and their severity vary between fish host species. Among other species such as Atlantic salmon (Salmo salar), brown trout usually is an asymptomatic carrier of M. cerebralis and can act as a vector for the parasite . Approximately 20 d post exposure, the parasite locates in the cartilage (primarily in the opercula, cranium and spine) causing lesions and inflammation, and potentially damaging the skeletal structure of the fish . Locations of lesions vary among fish species and can explain variability in disease severity and survival among different species. Cartilage lesions, observed by histology in rainbow trout, tend to concentrate in the cranium around the brain, but can be found throughout the cartilage in the body in young fry . However, in brown trout, lesions concentrate in fin rays and gill arches [58,59]. Under experimental conditions, Hedrick et al.  compared the relative susceptibility of age-matched brown and rainbow trout fry to M. cerebralis. Using different TAM doses (from 10 to 10000 TAMs/fish), they compared the effect of exposure dose on the prevalence of clinical signs, the severity of microscopic lesions, the presence and number of spores in the cartilage of infected fish, and the relative risk of infection for each fish species. In our study, a dose of 3000 TAMs/fish was used for exposure and a significant difference was detected both in prevalence and intensity of infection, similarly to the findings of Hedrick et al.  when using > 1000 TAMs/fish.
Although Nichols et al.  and Fetherman et al.  presumed that the geographic distribution of fish strains might have an influence on the genetic basis of resistance to both C. shasta and M. cerebralis, the correlation between the phylogeographic composition of salmonid populations and their resistance to the diseases has not been studied so far. Our study is also the first to compare the susceptibility of different genetic lineages of brown trout. The examined German brood stock represents a homogenous population, while the Hungarian population shows high heterogeneity concerning its phylogeographic lineages. Genetic analysis determined that the German population is Atlantic-descendant, while the Hungarian brood stock merely contains a higher ratio of At markers than the autochthonous Da one, and thereby > 50% of brood stock is comprised of At-Da hybrids. And the heterogenous group LBT composed of At and At-Da hybrid parent fish was the least susceptible to M. cerebralis of all examined fish groups.
As a key finding, the fully heterozygous group NIB was significantly less susceptible to the parasite than the more homozygous IB group. Thus, our findings suggest that there is a positive correlation between degree of parental heterozygosity and disease resistance. Higher variability in infection intensity was observed in the groups IB and REL compared to NIB and LBT, albeit the prevalence of infection did not differ significantly among the examined brown trout groups. The variations in groups IB and REL was mainly caused by a small fraction of specimens (<5%) that showed extremely heavy infections compared to their conspecifics. This was the case in the LRBT group as well, of course with even higher variance in spore numbers. Individual differences in resistance to certain myxozoan-caused diseases including whirling disease are well-known . In most cases, the individual differences in immune competence or the development of the innate immune system are the causes of variations in disease resistance. But in the present study, individual variations in the susceptibility of offspring might be also caused by the allele recombination of parent fish pairs/groups that may result higher variations in the homozygosity level of the offspring population than in that of the parental one. The findings of the genetic analysis of randomly selected offspring individuals support this, as we showed that IB brood fish also had NIB offspring and vice versa. The brood stock fish in REL group were heterozygous, but genetically closely related individuals were paired, therefore the relatively high susceptibility of REL offspring can be explained by the “forced” inbreeding caused by the selective fertilization. The opposite overall tendency could be observed in the IB group, in which only the parents were inbred and the selective fertilization probably outbred the offspring, at least to a certain extent. In the offspring of IB and REL, not only the elevated homozygosity may be a reason for the increased parasite susceptibility, but the resulting allele depletion that lowers immune diversity and versatility. The variations in infection intensity were less remarkable in the offspring of the heterozygous groups NIB and LBT, probably because these brood stock fish were non-related and the heterozygosity level did not decrease considerably in offspring.
Our findings indicate that the homozygosity level and thus the genetic diversity of brood stock has a significant influence on the intensity of parasite infection in the subsequent generations. As there is no effective treatment against whirling disease as yet, preventive measures against the parasite should be accompanied by regular genetic diversification of brood stock under controlled conditions to avoid increased levels of homozygosity in brood stock populations, even if barely measurable.
Charlesworth D, Charlesworth B. Inbreeding depression and its evolutionary consequences. Ann Rev Ecol Syst. 1987;18:237–68.
Falconer DS, Mackay TFC. Introduction to Quantitative Genetics. Essex: Longman Group Ltd; 1996. p. 464.
Lynch M, Walsh B. Genetics and Analysis of Quantitative Traits. Sunderland, MA: Sinauer Associates Inc; 1998. p. 980.
Kincaid HL. Inbreeding in fish populations used for aquaculture. Aquaculture. 1983;33:215–27.
Flagg TA, Waknitz FW, Maynard DJ, Milner GB, Mahnken CVW. The effect of hatcheries on native coho salmon populations in the lower Columbia River. Am Fish Soc Symp. 1995;15:366–75.
Wang SZ, Hard J, Utter F. Salmonid inbreeding: a review. Rev Fish Biol Fish. 2002;11:301–19.
Tave D. Genetics for Fish Hatchery Managers. Westport, Connecticut: AVI Publishing Company Inc; 1986. p. 299.
Withler RE, Beacham TD. Genetic consequences of the simultaneous or sequential addition of semen from multiple males during hatchery spawning of chinook salmon (Oncorhynchus tshawytscha). Aquaculture. 1994;12:11–23.
Simon RC. Management techniques to minimize the loss of genetic variability in hatchery fish populations. Am Fish Soc Symp. 1991;10:487–94.
Kincaid HL. An evaluation of inbreeding and effective population size in salmonid broodstocks in federal and state hatcheries. Am Fish Soc Symp. 1995;15:193–204.
Geiger HA, Smoker WW, Zhivotovsky LA, Gharrett AJ. Variability of family size and marine survival in pink salmon (Oncorhynchus gorbuscha) has implications for conservation biology and human use. Can J Fish Aquat Sc. 1997;54:2684–90.
Hard JJ, Connell L, Hershberger WK, Harrell LW. Genetic variation in mortality of chinook salmon (Oncorhynchus tshawytscha) during a bloom of the marine alga Heterosigma akashiwo. J Fish Biol. 2000;56:1387–97.
Herman RL. Prevention and control of fish diseases in hatcheries. In: Snieszko SF, editor. Symposium on diseases of fisheries and shellfishes. Special Publication 5. Bethesda, MD: American Fisheries Society; 1970. p. 3–15.
Price DJ. Genetics of susceptibility and resistance to disease in fishes. J Fish Biol. 1985;26:509–19.
Chevassus B, Dorson M. Genetics of resistance to disease in fishes. Aquaculture. 1990;85:83–107.
Wiegertjes GF, Stet RJM, Parmentier HK, van Muiswinkel WB. Immunogenetics of disease resistance in fish: a comparative approach. Dev Comp Immunol. 1996;20:365–82.
Ehlinger NF. Selective breeding of trout for resistance to furunculosis. N Y Fish Game J. 1977;24:25–36.
Slierendrecht WJ, Olesen NJ, Juul-Madsen HR, Lorenzen N, Henryon M, Berg P, et al. Rainbow trout offspring with different resistance to viral haemorrhagic septicaemia. Fish Shellfish Immunol. 2001;11:155–67.
Ibarra AM, Hedrick RP, Gall GAE. Genetic analysis of rainbow trout susceptibility to the myxosporean Ceratomyxa shasta. Aquaculture. 1994;120:239–62.
Nichols KM, Bartholomew J, Thorgaard GH. Mapping multiple genetic loci associated with Ceratomyxa shasta resistance in Oncorhynchus mykiss. Dis Aquat Org. 2003;56:145–54.
Schisler GJ, Myklebust KA, Hedrick RP. Inheritance of Myxobolus cerebralis resistance among F1-generation crosses of whirling disease resistant and susceptible rainbow trout strains. J Aquat Anim Health. 2006;18:109–15.
Baerwald MR, Petersen JL, Hedrick RP, Schisler GJ, May B. A major effect quantitative trait locus for whirling disease resistance identified in rainbow trout (Oncorhynchus mykiss). Heredity. 2011;106:920–6.
Fetherman ER, Winkelman DL, Schisler GJ, Antolin MF. Genetic basis of differences in myxospore count between whirling disease-resistant and susceptible strains of rainbow trout. Dis Aquat Org. 2012;102:97–106.
Halliday MM. The biology of Myxosoma cerebralis: the causative organism of whirling disease of salmonids. J Fish Biol. 1976;9:339–57.
Hoffman GL. Myxobolus cerebralis, a worldwide cause of salmonid whirling disease. J Aquat Anim Health. 1990;2:30–7.
Hedrick RP, El-Matbouli M, Adkison MA, MacConnell E. Whirling disease: re-emergence among wild trout. Immunol Rev. 1998;166:365–76.
MacConnell E, Vincent ER. The effects of Myxobolus cerebralis on the salmonid host. In: Bartholomew JL, Wilson JC, editors. Whirling disease: Reviews and current topics. American Fisheries Society, Symposium 29. MD: Bethesda; 2002. p. 95–107.
Fetherman ER, Winkelman DL, Baerwald MR, Schisler GJ. Survival and Reproduction of Myxobolus cerebralis-Resistant Rainbow Trout Introduced to the Colorado River and Increased Resistance of Age-0 Progeny. PLoS One. 2014;9:e96954.
Markiw ME. Whirling disease: Earliest susceptible age of rainbow trout to the triactinomyxid of Myxobolus cerebralis. Aquaculture. 1991;92:1–6.
Markiw ME. Experimentally induced whirling disease. I. Dose response of fry and adults of rainbow trout exposed to the triactinomyxon stage of Myxobolus cerebralis. J Aquat Anim Health. 1992;4:40–3.
El-Matbouli M, Fischer-Scherl T, Hoffmann RW. Present knowledge of the life cycle, taxonomy, pathology, and therapy of some Myxosporea spp. important for freshwater fish. Ann Rev Fish Dis. 1992;3:367–402.
Ryce EKN, Zale AV, MacConnell E, Nelson M. Effects of fish age versus size on the development of whirling disease in rainbow trout. Dis Aquat Org. 2005;63:69–76.
Hoffman GL. Intercontinental and transcontinental dissemination and transfaunation of fish parasites with emphasis on whirling disease (Myxosoma cerebralis). In: Snieszko SF, editor. Symposium on diseases of fisheries and shellfishes. Special Publication 5. Bethesda, MD: American Fisheries Society; 1970. p. 69–81.
Laikre L, Antunes A, Apostolidis A, Berrebi P, Duguid A, Ferguson A, et al. Conservation genetic management of brown trout (Salmo trutta) in Europe. Report by the concerted action on identification, management and exploitation of genetic resources in the brown trout (Salmo trutta). Silkeborg: TROUTCONCERT (EU FAIR CT97–3882). Danmarks fiskeriundersrgelser; 1999. p. 1–91.
Paladini G, Hansen H, Williams CF, Taylor NGH, Rubio-Mejía OL, Denholm SJ, et al. Reservoir hosts for Gyrodactylus salaris may play a more significant role in epidemics than previously thought. Parasites Vectors. 2014;7:576.
Estoup A, Largiader CR, Perrot E, Chourrout D. Rapid one-tube DNA extraction for reliable PCR detection of fish polymorphic markers and transgenes. Mol Marine Biol Biotech. 1996;5:295–8.
Snoj A, Jug T, Melkic E, Susnik S, Pohar J, Dovc P, et al. Mitochondrial and microsatellite DNA analysis of marble trout in Slovenia. J Freshwater Biol (Quaderni ETP). 2000;29:5–11.
Bernatchez L. The evolutionary history of brown trout Salmo trutta L. inferred from phylogeographic, nested clade, and mismatch analyses of mitochondrial DNA variation. Evolution. 2001;55:351–79.
McMeel OM, Hoey EM, Ferguson A. Partial nucleotide sequences, and routine typing by polymerase chain reaction-restriction fragment length polymorphism, of the brown trout (Salmo trutta) lactate dehydrogenase, LDH-C1-C1*90 and *100 alleles. Mol Ecol. 2001;10:29–34.
Marić S, Simonović P, Razpet A. Genetic characterization of brood stock brown trout from Bled fish-farm, Slovenia. Period Biologorum. 2010;112:145–8.
Estoup A, Presa P, Krieg F, Vaiman D, Guyomard R. CT)n and (GT)n microsatellites: a new class of genetic markers for Salmo trutta L. (brown trout. Heredity. 1993;71:488–96.
Presa P, Guyomard R. Conservation of microsatellites in three species of salmonids. J Fish Biol. 1996;49:1326–9.
Slettan A, Olsaker I, Lie Ø. Atlantic salmon, Salmo salar, microsatellites at the SSOSL25, SSOSL85, SSOSL311, SSOSL417 loci. Anim Gen. 1995;26:281–2.
O’Reilly PT, Hamilton LC, McConnell SK, Wright JM. Rapid analysis of genetic variation in Atlantic salmon (Salmo salar) by PCR multiplexing of dinucleotide and tetranucleotide microsatellites. Can J Fisheries Aquat Sci. 1996;53:2292–8.
Smith CT, Koop BF, Nelson RJ. Isolation and characterization of coho salmon (Oncorhynchus kisutch) microsatellites and their use in other salmonids. Mol Ecol. 1998;7:1614–6.
Wang J. COANCESTRY: A program for simulating, estimating and analysing relatedness and inbreeding coefficients. Mol Ecol Resources. 2011;11:141–5.
Kallert DM, El-Matbouli M. Differences in viability and reactivity of actinospores of three myxozoan species upon ageing. Folia Parasitol. 2008;55:105–10.
Eszterbauer E, Kallert DM, Grabner D, El-Matbouli M. Differentially expressed parasite genes involved in host recognition and invasion of the triactinomyxon stage of Myxobolus cerebralis (Myxozoa). Parasitology. 2009;136:367–77.
Rózsa L, Reiczigel J, Majoros G. Quantifying parasites in samples of hosts. J Parasitol. 2000;86:228–32.
Keller LF, Waller DM. Inbreeding effects in wild populations. Trends Ecol Evol. 2002;17:230–41.
McKay LR, McMillan I, Sadler SE, Moccia RD. Effects of mating system on inbreeding levels and selection response in Salmonid aquaculture. Aquaculture. 1992;100:100–1.
Aulstad D, Gjedrem T, Skjervold H. Genetic and environmental sources of variation in length and weight of rainbow trout (Salmo gairdneri). J Fish Res Board Canada. 1972;29:237–41.
Gjerde B, Gunnes K, Gjedrem T. Effect of inbreeding on survival and growth in rainbow trout. Aquaculture. 1983;34:327–32.
Su GS, Liljedahl LE, Gall GAE. Effects of inbreeding on growth and reproductive traits in rainbow trout (Oncorhynchus mykiss). Aquaculture. 1996;142:139–48.
Balloux F, Amos W, Coulson T. Does heterozygosity estimate inbreeding in real populations? Mol Ecol. 2004;13:3021–31.
Steinbach Elwell LC, Stromberg KE, Ryce EKN, Bartholomew JL. Whirling Disease in the United States. A summary of progress in research and management 2009. Bozeman, Montana: Trout Unlimited and The Whirling Disease Foundation; 2009. p. 61.
El-Matbouli M, Hoffmann RW, Mandok C. Light and electron microscopic observations on the route of the triactinomyxon-sporoplasm of Myxobolus cerebralis from epidermis into rainbow trout cartilage. J Fish Biol. 1995;46:919–35.
Hedrick RP, McDowell TS, Gay M, Marty GD, Georgiadis MP, MacConnell E. Comparative susceptibility of rainbow trout Onchorhynchus mykiss and brown trout Salmo trutta to Myxobolus cerebralis, the cause of salmonid whirling disease. Dis Aquat Org. 1999;37:173–83.
Baldwin TJ, Vincent ER, Silflow RM, Stanek D. Myxobolus cerebralis infection in rainbow trout (Oncorhynchus mykiss) and brown trout (Salmo trutta) exposed under natural stream conditions. J Vet Diag Invest. 2000;12:312–21.
Gómez D, Bartholomew J, Sunyer JO. Biology and mucosal immunity to myxozoans. Dev Comp Immunol. 2014;43:243–56.
We thank the János Bolyai Research Scholarship of the Hungarian Academy of Sciences, Hungary for EE, and the Alexander von Humboldt-Foundation for a Feodor-Lynen Stipend for DMK. Thanks are due to J. Reiczigel for his helpful advice for the statistical analysis. The study was supported by the Hungarian Scientific Research Fund (OTKA grant No. K112301).
Freeware software: R version 2.15.1 (2012-06-22) -- “Roasted Marshmallows”. The R Foundation for Statistical Computing. ISBN 3-900051-07-0. Platform: i386-pc-mingw32/i386 (32-bit).
The authors declare that they have no competing interests.
EE and DMK conceived and designed the study. DMK, BF, GH executed the selective fish fertilization. ZT performed the lineage analysis, and EE carried out the microsatellite analysis. BF executed the infection trial. CFG and GZ prepared the samples for parasite spore counting. EE and BF counted the parasites and evaluated the data. EE drafted the manuscript and all authors read and approved its final version.