Genome-wide alternative splicing profile in the posterior kidney of brown trout (Salmo trutta) during proliferative kidney disease.

BACKGROUND: The cnidarian myxozoan parasite Tetracapsuloides bryosalmonae causes chronic proliferative kidney disease (PKD) in salmonids. This parasite is a serious threat to wild and cultured salmonids. T. bryosalmonae undergoes intra-luminal sporogonic development in the kidney of brown trout (Salmo trutta) and the viable spores are released via urine. We investigated the alternative splicing pattern in the posterior kidney of brown trout during PKD. RESULTS: RNA-seq data were generated from the posterior kidney of brown trout collected at 12 weeks post-exposure to T. bryosalmonae. Subsequently, this data was mapped to the brown trout genome. About 153 significant differently expressed alternatively spliced (DEAS) genes, (delta PSI = 5%, FDR P-value < 0.05) were identified from 19,722 alternatively spliced events. Among the DEAS genes, the least and most abundant alternative splicing types were alternative 5' splice site (5.23%) and exon skipping (70.59%), respectively. The DEAS genes were significantly enriched for sodium-potassium transporter activity and ion homeostasis (ahcyl1, atp1a3a, atp1a1a.1, and atp1a1a.5). The protein-protein interaction network analysis enriched two local network clusters namely cation transporting ATPase C-terminus and Sodium/potassium ATPase beta chain cluster, and mixed inclusion of Ion homeostasis and EF-hand domain cluster. Furthermore, the human disease-related salmonella infection pathway was significantly enriched in the protein-protein interaction network. CONCLUSION: This study provides the first baseline information about alternative splicing in brown trout during PKD. The generated data lay a foundation for further functional molecular studies in PKD - brown trout infection model. The information generated from the present study can help to develop therapeutic strategies for PKD in the future.

Disease mortality in domesticated animals is predicted by host evolutionary relationships.

Infectious diseases of domesticated animals impact human well-being via food insecurity, loss of livelihoods, and human infections. While much research has focused on parasites that infect single host species, most parasites of domesticated mammals infect multiple species. The impact of multihost parasites varies across hosts; some rarely result in death, whereas others are nearly always fatal. Despite their high ecological and societal costs, we currently lack theory for predicting the lethality of multihost parasites. Here, using a global dataset of >4,000 case-fatality rates for 65 infectious diseases (caused by microparasites and macroparasites) and 12 domesticated host species, we show that the average evolutionary distance from an infected host to other mammal host species is a strong predictor of disease-induced mortality. We find that as parasites infect species outside of their documented phylogenetic host range, they are more likely to result in lethal infections, with the odds of death doubling for each additional 10 million years of evolutionary distance. Our results for domesticated animal diseases reveal patterns in the evolution of highly lethal parasites that are difficult to observe in the wild and further suggest that the severity of infectious diseases may be predicted from evolutionary relationships among hosts.

Depletion of MHC supertype during domestication can compromise immunocompetence.

The major histocompatibility complex (MHC) plays an important role in infectious disease resistance. The presence of certain MHC alleles and functionally similar groups of MHC alleles (i.e., supertypes) has been associated with resistance to particular parasite species. Farmed and domesticated fish stocks are often depleted in their MHC alleles and supertype diversity, possibly as a consequence of artificial selection for desirable traits, inbreeding (loss of heterozygosity), genetic drift (loss of allelic diversity) and/or reduced parasite biodiversity. Here we quantify the effects of depletion of MHC class II genotype and supertype variation on resistance to the parasite Gyrodactylus turnbulli in guppies (Poecilia reticulata). Compared to the descendants of wild-caught guppies, ornamental fish had a significantly reduced MHC variation (i.e., the numbers of MHC alleles and supertypes per individual, and per population). In addition, ornamental fish were significantly more susceptible to G. turnbulli infections, accumulating peak intensity 10 times higher than that of their wildtype counterparts. Four out of 13 supertypes were associated with a significantly reduced parasite load, and the presence of some supertypes had a dramatic effect on the intensity of infection. Remarkably, the ornamental and wildtype fish differed in the supertypes that were associated with parasite resistance. Analysis with a genetic algorithm showed that resistance-conferring supertypes of the wildtype and ornamental fish shared two unique amino acids in the peptide-binding region of the MHC that were not found in any other alleles. These data show that the supertype demarcation captures some, but not all, of the variation in the immune function of the alleles. This study highlights the importance of managing functional MHC diversity in livestock, and suggests there might be some immunological redundancy among MHC supertypes.

Identification of variants associated with hard clam, Mercenaria mercenaria, resistance to Quahog Parasite Unknown disease.

Severe losses in aquacultured and wild hard clam (Mercenaria mercenaria) stocks have been previously reported in the northeastern United States due to a protistan parasite called QPX (Quahog Parasite Unknown). Previous work demonstrated that clam resistance to QPX is under genetic control. This study identifies single nucleotide polymorphism (SNP) associated with clam survivorship from two geographically segregated populations, both deployed in an enzootic site. The analysis contrasted samples collected before and after undergoing QPX-related mortalities and relied on a robust draft clam genome assembly. ~200 genes displayed significant variant enrichment at each sampling point in both populations, including 18 genes shared between both populations. Markers from both populations were identified in genes related to apoptosis pathways, protein-protein interaction, receptors, and signaling. This research begins to identify genetic markers associated with clam resistance to QPX disease, leading the way for the development of resistant clam stocks through marker-assisted selection.

Modulation of local and systemic immune responses in brown trout (Salmo trutta) following exposure to Myxobolus cerebralis.

Myxobolus cerebralis, the etiological agent of Whirling Disease (WD), is a freshwater myxozoan parasite with considerable economic and ecological relevance for salmonids. There are differences in disease susceptibility between species and strains of salmonids. Recently, we have reported that the suppressor of cytokine signaling SOCS1 and SOCS3 are key in modulating rainbow trout (Oncorhynchus mykiss) immune responses and that resistant fish apparently exhibit effective Th17 cell response after exposure to M. cerebralis. It is unclear whether such molecules and pathways are also involved in the immune response of M. cerebralis infected brown trout (Salmo trutta). Hence, this study aimed to explore their role during immune modulation in infected brown trout, which is considered resistant to this parasite. Fish were exposed to the triactinomyxon (TAM) stages of M. cerebralis and quantitative real-time PCR (RT-qPCR) was carried out to examine local (caudal fin) and systemic (head kidney, spleen) immune transcriptional changes associated with WD over time in infected and control fish. All of the immune genes in the three tissues studied were differentially expressed in infected fish at multiple time points. Brown trout reduced the parasite load and demonstrated effective immune responses, likely by keeping pro-inflammatory and anti-inflammatory cytokines in balance whilst stimulating efficient Th17-mediated immunity. This study increases knowledge on the brown trout immune response to M. cerebralis and helps us to understand the underlying mechanisms of WD resistance.

Effects of Enteromyxum spp. (Myxozoa) infection in the regulation of intestinal E-cadherin: Turbot against gilthead sea bream.

Enteromyxoses are relevant diseases for turbot and gilthead sea bream aquaculture. The myxozoan parasites invade the intestinal mucosa, causing a cachectic syndrome associated with intestinal barrier alteration; nonetheless, their pathological impact is different. Turbot infected by Enteromyxum scophthalmi develop more severe intestinal lesions, reaching mortality rates of 100%, whereas in E. leei-infected gilthead sea bream, the disease progresses slowly, and mortality rates are lower. The mechanisms underlying the different pathogenesis are still unclear. We studied the distribution and expression changes of E-cadherin, a highly conserved protein of the adherens junctions, in the intestine of both species by immunohistochemistry and quantitative PCR, using the same immunohistochemical protocol and common primers. The regular immunostaining pattern observed in control fish turned into markedly irregular in parasitized turbot, showing an intense immunoreaction at the host-parasite interface. Nevertheless, E-cadherin gene expression was not significantly modulated in this species. On the contrary, no evident changes in the protein distribution were noticed in gilthead sea bream, whereas a significant gene downregulation occurred in advanced infection. The results contribute to the understanding of the different host-parasite interactions in enteromyxoses. Host and parasite cells appear to establish diverse relationships in these species, which could underlie the different pathological picture.

Identification and Expression Profiling of Toll-Like Receptors of Brown Trout (Salmo trutta) during Proliferative Kidney Disease.

Proliferative kidney disease is an emerging disease among salmonids in Europe and North America caused by the myxozoan parasite Tetracapsuloides bryosalmonae. The decline of endemic brown trout (Salmo trutta) in the Alpine streams of Europe is fostered by T. bryosalmonae infection. Toll-like receptors (TLRs) are a family of pattern recognition receptors that acts as sentinels of the immune system against the invading pathogens. However, little is known about the TLRs' response in salmonids against the myxozoan infection. In the present study, we identified and evaluated TLR1, TLR19, and TLR13-like genes of brown trout using data-mining and phylogenetic analysis. The expression pattern of TLRs was examined in the posterior kidney of brown trout infected with T. bryosalmonae at various time points. Typical Toll/interleukin-1 receptor protein domain was found in all tested TLRs. However, TLR13-like chr2 had a short amino acid sequence with no LRR domain. Phylogenetic analysis illustrated that TLR orthologs are conserved across vertebrates. Similarly, a conserved synteny gene block arrangement was observed in the case of TLR1 and TLR19 across fish species. Interestingly, all tested TLRs showed their maximal relative expression from 6 to 10 weeks post-exposure to the parasite. Our results suggest that these TLRs may play an important role in the innate defense mechanism of brown trout against the invading T. bryosalmonae.

Survival of Whirling-Disease-Resistant Rainbow Trout Fry in the Wild: A Comparison of Two Strains.

Introduced pathogens can affect fish populations, and three main factors affect disease occurrence: the environment, host, and pathogen. Manipulating at least one of these factors is necessary for controlling disease. Myxobolus cerebralis, the parasite responsible for salmonid whirling disease, became established in Colorado during the 1990s and caused significant declines in wild Rainbow Trout Oncorhynchus mykiss populations. Attempts to re-establish Rainbow Trout have focused on manipulating salmonid host resistance. A Rainbow Trout strain known as GR × CRR was developed for stocking in Colorado by crossing a whirling-disease-resistant strain known as the German Rainbow Trout (GR) with the Colorado River Rainbow Trout (CRR). The GR × CRR fish exhibit resistance similar to that shown by GR, and survival and reproduction were expected to be similar to those of CRR. One disadvantage of stocking GR × CRR is that outcrossing and backcrossing could decrease resistance, and laboratory studies have indicated that this can occur. A potential disadvantage of stocking pure GR is lower survival due to domestication. To compare fry survival between the strains, a field experiment was conducted in 1.6-km reaches of nine Colorado streams. Each stream was stocked in August 2014 with 5,000 GR × CRR and 5,000 GR individuals. In October 2014, April 2015, and August 2015, apparent survival was assessed. Two laboratory predation experiments were also conducted. The field experiment revealed that short-term apparent survival was influenced by stream, and growth rate was influenced by strain and stream. However, after 12 months, there was no difference in apparent survival or growth rate between the GR and GR × CRR strains. Laboratory experiments showed that survival did not differ between the strains when confronted with Brown Trout Salmo trutta predation. Our results indicate that the GR strain is a viable option for stocking in streams where M. cerebralis is enzootic. Further evaluation is needed to determine whether GR fish will survive to maturity and reproduce.

Comparison of Sampling and Detection Methods for Chinook Salmon and Steelhead Naturally Infected with Myxobolus cerebralis.

Myxobolus cerebralis (Mc) is a myxozoan parasite causing whirling disease in hatchery- and natural-origin salmonids. To minimize spread of this parasite and the incidence of its associated disease, fish health professionals routinely screen fish for Mc before stocking or moving the fish to Mc-free waters. Sample collection for Mc traditionally entails lethal sampling of cranial tissue followed by pepsin-trypsin digest (PTD) and screening of the sample for mature myxobolid myxospores (PTD method); however, nonlethal sampling methods would be advantageous in some circumstances, such as when dealing with rare or otherwise valuable fish. Accordingly, we compared Mc detections in cranial cartilage by using the PTD method with PCR assays of fin biopsies collected from juvenile Chinook Salmon Oncorhynchus tshawytscha and adult steelhead O. mykiss. Cranial samples were also analyzed using PCR methods for comparative purposes. Results indicated that Mc could be detected by PCR in fin clips, but the results generated by this approach differed significantly from those associated with PTD- and/or PCR-based analysis of cranial cartilage samples. Polymerase chain reaction-based analysis-of individual head samples and head digest pools in both species as well as fins in steelhead-yielded more positive detections than PTD analysis alone. The PCR-based analysis of head and fin tissues yielded different Mc detection rates in both species, but the nature of the detection disparity varied depending on the species and/or life stage of the fish. We conclude that for lethal cranial samples, neither PTD nor PCR should be used alone, but using these techniques in concert may provide the most complete and accurate estimation of Mc presence in a group of salmonids. If imperiled or highly valuable fish are in question, nonlethal fin samples may be used to generate some information regarding Mc status, with the understanding that parasite DNA detections do not necessarily signify mature infections or disease.

Telomere length and antioxidant defense associate with parasite-induced retarded growth in wild brown trout.

Early growth conditions can have profound impacts on individuals' development, growth and physiology, with subsequent long-term consequences for individuals' fitness and life expectancy. Telomere length (TL) has been suggested to indicate both individual fitness and life expectancy in wide range of species, as the telomere attrition rate at early age can be accelerated due to exposure to various stressors, including parasites and inflammatory diseases, which increase production of reactive oxygen species (ROS) and influence antioxidant (AO) levels. We investigated impacts of Tetracapsuloides bryosalmonae infection, a causative agent of proliferative kidney disease (PKD), on AO status and TL in a natural population of juvenile brown trout (Salmo trutta). The fish with higher parasite load showed more severe kidney hyperplasia, anemia and smaller body size compared to less parasitized fish. Furthermore, fish with severe PKD symptoms had lower SOD-, CAT- and GST activity than fish with milder kidney hyperplasia. However, parasite load was not directly correlated either with AOs or with TL. Smaller fish showed shorter TLs, potentially reflecting lower individual quality. The fish, which were less sensitive to parasite-induced impaired growth, quantified as parasite load-adjusted fork length, showed also longer TLs, lower GR- and GST activity and less GSHtot compared to more sensitive fish. These results provide novel knowledge about the impacts of the PKD in brown trout at the molecular level and support the idea that TL may reflect individual quality and ability to cope with parasitic infections.

Does the parasite-mediated selection drive the MHC class IIB diversity in wild populations of European chub (Squalius cephalus)?

The genes of major histocompatibility complex (MHC) provide an excellent opportunity to study host-parasite relationships because they are expected to evolve in response to parasites and variation in parasite communities. In this study, we investigated the potential role of parasite-mediated selection acting on MHC class IIB (DAB) genes in European chub (Squalius cephalus) natural populations. We found significant differences between populations in metazoan parasites, neutral and adaptive genetic diversities. The analyses based on pairwise data revealed that populations with dissimilar MHC allelic profiles were geographically distant populations with significantly different diversity in microsatellites and a dissimilar composition of parasite communities. The results from the generalized estimating equations method (GEE) on the level of individuals revealed that metazoan parasite load in European chub was influenced by the diversity of DAB alleles as well as by the diversity of neutral genetic markers and host traits reflecting condition and immunocompetence. The multivariate co-inertia analysis showed specific associations between DAB alleles and parasite species. DAB1-like alleles were more involved in associations with ectoparasites, while DAB3-like alleles were positively associated with endoparasites which could suggest potential differences between DAB genes caused by different selection pressure. Our study revealed that parasite-mediated selection is not the only variable affecting MHC diversity in European chub; however, we strongly support the role of neutral processes as the main driver of DAB diversity across populations. In addition, our study contributes to the understanding of the evolution of MHC genes in wild living fish.

Immunohistochemical detection and gene expression of TNFα in turbot (Scophthalmus maximus) enteromyxosis.

Enteromyxum scophthalmi (Myxozoa) constitutes one of the most devastating pathogens for turbot (Scophthalmus maximus, L.) aquaculture. This parasite causes a severe intestinal parasitosis that leads to a cachectic syndrome with high morbidity and mortality rates for which no therapeutic options are available. Presence of inflammatory infiltrates, increased apoptotic rates and epithelial detaching have been described at intestinal level, as well as leukocyte depletion in lymphohaematopoietic organs. Previous investigations on enteromyxosis in turbot showed the high susceptibility of this species to the parasite and reported the existence of a dysregulated immune response against the parasite. The pleiotropic cytokine tumour necrosis factor alpha (TNFα) plays a major role in immune response and is involved in a wide range of biological activities. In teleost, the gene expression of this cytokine has been found regulated under several pathological conditions. Teleost TNFα shows some analogous functions with its mammalian counterparts, but the extent of its activities is still poorly understood. Cytokines are generally considered as a double-edge sword and TNFα has been implicated in the pathogenesis of different inflammatory diseases as well as in wasting syndromes described in mammals. The aim of this work was to analyse the expression of TNFα during enteromyxosis with molecular (Q-PCR) and morphological (immunohistochemistry) tools. Kidney, spleen and pyloric caeca from turbot with moderate and severe infections were analysed and compared to healthy naïve fish. TNFα expression was increased in both spleen and kidney in the earlier stages of the disease, whereas in severely infected fish, the expression decreased, especially in kidney. At the intestinal level, an increase in the number of TNFα-positive cells was noticed, which was proportional to the infiltration of inflammatory cells. The results demonstrate the involvement of TNFα in the immune response to E. scophthalmi in turbot, which could be related to the development of the clinic signs and lesions.

RNA-seq analysis reveals significant transcriptome changes in turbot (Scophthalmus maximus) suffering severe enteromyxosis.

BACKGROUND: Enteromyxosis caused by the intestinal myxozoan parasite Enteromyxum scophthalmi is a serious threat for turbot (Scophthalmus maximus, L.) aquaculture, causing severe catarrhal enteritis leading to a cachectic syndrome, with no therapeutic options available. There are still many aspects of host-parasite interaction and disease pathogenesis that are yet to be elucidated, and to date, no analysis of the transcriptomic changes induced by E. scophthalmi in turbot organs has been conducted. In this study, RNA-seq technology was applied to head kidney, spleen and pyloric caeca of severely infected turbot with the aim of furthering our understanding of the pathogenetic mechanisms and turbot immune response against enteromyxosis. RESULTS: A huge amount of information was generated with more than 23,000 identified genes in the three organs, amongst which 4,762 were differently expressed (DE) between infected and control fish. Associate gene functions were studied based on gene ontology terms and available literature, and the most interesting DE genes were classified into five categories: 1) immune and defence response; 2) apoptosis and cell proliferation; 3) iron metabolism and erythropoiesis; 4) cytoskeleton and extracellular matrix and 5) metabolism and digestive function. The analysis of down-regulated genes of the first category revealed evidences of a connexion failure between innate and adaptive immune response, especially represented by a high number of DE interferon-related genes in the three organs. Furthermore, we found an intense activation of local immune response at intestinal level that appeared exacerbated, whereas in kidney and spleen genes involved in adaptive immune response were mainly down-regulated. The apoptotic machinery was only clearly activated in pyloric caeca, while kidney and spleen showed a marked depression of genes related to erythropoiesis, probably related to disorders in iron homeostasis. The genetic signature of the causes and consequences of cachexia was also demonstrated by the down-regulation of the genes encoding structural proteins and those involved in the digestive metabolism. CONCLUSIONS: This transcriptomic study has enabled us to gain a better understanding of the pathogenesis of enteromyxosis and identify a large number of DE target genes that bring us closer to the development of strategies designed to effectively combat this pathogen.

Interleukin gene expression is strongly modulated at the local level in a fish-parasite model.

The goal of this work was to identify interleukin (IL)-related genes in the gilthead sea bream (GSB) (Sparus aurata L.) and how they are modulated by the parasite Enteromyxum leei, a myxozoan that causes severe enteritis with a strong inflammatory response. A Blast-X search of our transcriptomic GSB database (www.nutrigroup-iats.org/seabreamdb) identified 16 new sequences encompassing seven ILs (IL-7, IL-8, IL-10, IL-12ß, IL-15, IL-18, and IL-34), the interleukin enhancer-binding factor 2 (ILF2), and eight IL receptors (IL-R); IL-R1, IL-6RA, IL-6RB, IL-8RA, IL-10RA, IL-10RB, IL-18R1, and IL-22R. Except for ILF2, their expression, plus that of IL-1ß, IL-1R2, IL-6, and TNF-α (from public repositories), were analysed by 96-well PCR array of samples of blood, spleen, head kidney, and intestine of GSB that were anally intubated with E. leei (recipient group, RCPT). Only the expression profile of the intestine of RCPT fish showed significant difference as compared to samples from PBS-inoculated fish. At 17 days post intubation (dpi), the expression of key pro-inflammatory ILs, such as IL-8, IL-8R, IL-12ß, and TNFα was significantly up-regulated, whereas at 64 dpi, anti-inflammatory IL expression (IL-6, IL-6RB, IL-7, IL-10, IL-10RA, and IL-15) was predominant. These results indicate a modification of the IL expression at late times post infection, probably to protect the fish intestine from the parasite and damage inflicted by an excessive inflammatory response. Furthermore, the response is mainly mediated at the local level as no significant changes were detected in blood, spleen and head kidney.

Ecological genomics of host behavior manipulation by parasites.

Among the vast array of niche exploitation strategies exhibited by millions of different species on Earth, parasitic lifestyles are characterized by extremely successful evolutionary outcomes. Some parasites even seem to have the ability to 'control' their host's behavior to fulfill their own vital needs. Research efforts in the past decades have focused on surveying the phylogenetic diversity and ecological nature of these host-parasite interactions, and trying to understand their evolutionary significance. However, to understand the proximal and ultimate causes of these behavioral alterations triggered by parasitic infections, the underlying molecular mechanisms governing them must be uncovered. Studies using ecological genomics approaches have identified key candidate molecules involved in host-parasite molecular cross-talk, but also molecules not expected to alter behavior. These studies have shown the importance of following up with functional analyses, using a comparative approach and including a time-series analysis. High-throughput methods surveying different levels of biological information, such as the transcriptome and the epigenome, suggest that specific biologically-relevant processes are affected by infection, that sex-specific effects at the level of behavior are recapitulated at the level of transcription, and that epigenetic control represents a key factor in managing life cycle stages of the parasite through temporal regulation of gene expression. Post-translational processes, such as protein-protein interactions (interactome) and post translational modifications (e.g. protein phosphorylation, phosphorylome), and processes modifying gene expression and translation, such as interactions with microRNAs (microRNAome), are examples of promising avenues to explore to obtain crucial insights into the proximal and ultimate causes of these fascinating and complex inter-specific interactions.

Dietary vegetable oils do not alter the intestine transcriptome of gilthead sea bream (Sparus aurata), but modulate the transcriptomic response to infection with Enteromyxum leei.

BACKGROUND: Studies conducted with gilthead sea bream (Sparus aurata L.) have determined the maximum dietary replacement of fish meal and oil without compromising growth or product quality. The present study aimed to analyze the effect of the nutritional background on fish health and fish fed plant protein-based diets with fish oil (FO diet) or a blend of vegetable oils (66VO diet) were exposed for 102 days to the intestinal myxosporean parasite Enteromyxum leei, and the intestine transcriptome was analyzed with a customized oligo-microarray of 7,500 annotated genes. RESULTS: Infection prevalence was high and similar in the two diet groups, but the outcome of the disease was more pronounced in fish fed the 66VO diet. No differences were found in the transcriptome of both diet control groups, whereas the number of differentially expressed genes in infected groups was considerable. K-means clustering of these differentially expressed genes identified four expression patterns that reflected the progression of the disease with the magnitude of the fold-change being higher in infected 66VO fish. A positive correlation was found between the time of infection and the magnitude of the transcriptional change within the 66VO group, being higher in early infected animals. Within this diet group, a strong up-regulation of many components of the immune specific response was evidenced, whereas other genes related to complement response and xenobiotic metabolism were down-regulated. CONCLUSIONS: The high replacement of fish oil by vegetable oils in practical fish feeds did not modify the intestine transcriptome of gilthead sea bream, but important changes were apparent when fish were exposed to the myxosporean E. leei. The detected changes were mostly a consequence rather than a cause of the different disease progression in the two diet groups. Hence, the developed microarray constitutes an excellent diagnostic tool to address changes associated with the action of intestinal pathogens, but lacks a prognostic value to predict in advance the different susceptibility of growing fish to the current pathogen.

Genomic regions showing copy number variations associate with resistance or susceptibility to gastrointestinal nematodes in Angus cattle.

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. We previously reported an initial analysis of copy number variations (CNVs) in Angus cattle selected for resistance or susceptibility to gastrointestinal nematodes. In this study, we performed a large-scale analysis of CNVs using SNP genotyping data from 472 animals of the same population. We detected 811 candidate CNV regions, which represent 141.8 Mb (~4.7%) of the genome. To investigate the functional impacts of CNVs, we created 2 groups of 100 individual animals with extremely low or high estimated breeding values of eggs per gram of feces and referred to these groups as parasite resistant (PR) or parasite susceptible (PS), respectively. We identified 297 (~51 Mb) and 282 (~48 Mb) CNV regions from PR and PS groups, respectively. Approximately 60% of the CNV regions were specific to the PS group or PR group of animals. Selected PR- or PS-specific CNVs were further experimentally validated by quantitative PCR. A total of 297 PR CNV regions overlapped with 437 Ensembl genes enriched in immunity and defense, like WC1 gene which uniquely expresses on gamma/delta T cells in cattle. Network analyses indicated that the PR-specific genes were predominantly involved in gastrointestinal disease, immunological disease, inflammatory response, cell-to-cell signaling and interaction, lymphoid tissue development, and cell death. By contrast, the 282 PS CNV regions contained 473 Ensembl genes which are overrepresented in environmental interactions. Network analyses indicated that the PS-specific genes were particularly enriched for inflammatory response, immune cell trafficking, metabolic disease, cell cycle, and cellular organization and movement.

Genetic basis of differences in myxospore count between whirling disease-resistant and susceptible strains of rainbow trout.

We used a quantitative genetics approach and estimated broad sense heritability (h2b) of myxospore count and the number of genes involved in myxospore formation to gain a better understanding of how resistance to Myxobolus cerebralis, the parasite responsible for whirling disease, is inherited in rainbow trout Oncorhynchus mykiss. An M. cerebralis-resistant strain of rainbow trout, the German Rainbow (GR), and a wild, susceptible strain of rainbow trout, the Colorado River Rainbow (CRR), were spawned to create 3 intermediate crossed populations (an F1 cross, F2 intercross, and a B2 backcross between the F1 and the CRR). Within each strain or cross, h2b was estimated from the between-family variance of myxospore counts using full-sibling families. Estimates of h2b and average myxospore counts were lowest in the GR strain, F1 cross, and F2 intercross (h2b = 0.34, 0.42, and 0.34; myxospores fish-1 = 275, 9566, and 45780, respectively), and highest in the B2 backcross and CRR strain (h2b = 0.93 and 0.89; myxospores fish-1 = 97865 and 187595, respectively). Comparison of means and a joint-scaling test suggest that resistance alleles arising from the GR strain are dominant to susceptible alleles from the CRR strain. Resistance was retained in the intermediate crosses but decreased as filial generation number increased (F2) or backcrossing occurred (B2). The estimated number of segregating loci responsible for differences in myxospore count in the parental strains was 9 ± 5. Our results indicate that resistance to M. cerebralis is a heritable trait within these populations and would respond to either artificial selection in hatcheries or natural selection in the wild.

A major effect quantitative trait locus for whirling disease resistance identified in rainbow trout (Oncorhynchus mykiss).

Whirling disease, caused by the pathogen Myxobolus cerebralis, leads to skeletal deformation, neurological impairment and under certain conditions, mortality of juvenile salmonid fishes. The disease has impacted the propagation and survival of many salmonid species over six continents, with particularly negative consequences for rainbow trout. To assess the genetic basis of whirling disease resistance in rainbow trout, genome-wide mapping was initiated using a large outbred F(2) rainbow trout family (n=480) and results were confirmed in three additional outbred F(2) families (n=96 per family). A single quantitative trait locus (QTL) region on chromosome Omy9 was identified in the large mapping family and confirmed in all additional families. This region explains 50-86% of the phenotypic variance across families. Therefore, these data establish that a single QTL region is capable of explaining a large percentage of the phenotypic variance contributing to whirling disease resistance. This is the first genetic region discovered that contributes directly to the whirling disease phenotype and the finding moves the field closer to a mechanistic understanding of resistance to this important disease of salmonid fish.

Genome-Wide Association Study Using Whole-Genome Sequence Data for Fertility, Health Indicator, and Endoparasite Infection Traits in German Black Pied Cattle.

This genome-wide association study (GWAS) aimed to identify sequence variants (SVs) and candidate genes associated with fertility and health in endangered German Black Pied cattle (DSN) based on whole-genome sequence (WGS) data. We used 304 sequenced DSN cattle for the imputation of 1797 genotyped DSN to WGS. The final dataset included 11,413,456 SVs of 1886 cows. Cow traits were calving-to-first service interval (CTFS), non-return after 56 days (NR56), somatic cell score (SCS), fat-to-protein ratio (FPR), and three pre-corrected endoparasite infection traits. We identified 40 SVs above the genome-wide significance and suggestive threshold associated with CTFS and NR56, and three important potential candidate genes (ARHGAP21, MARCH11, and ZNF462). For SCS, most associations were observed on BTA 25. The GWAS revealed 61 SVs, a cluster of 10 candidate genes on BTA 13, and 7 pathways for FPR, including key mediators involved in milk fat synthesis. The strongest associations for gastrointestinal nematode and Dictyocaulus viviparus infections were detected on BTA 8 and 24, respectively. For Fasciola hepatica infections, the strongest associated SVs were located on BTA 4 and 7. We detected 200 genes for endoparasite infection traits, related to 16 pathways involved in host immune response during infection.