Gene expression response to a nematode parasite in novel and native eel hosts

Subtle transcriptomic response of Eurasian perch (Perca fluviatilis) associated with Triaenophorus nodulosus plerocercoid infection

Determining the physiological effects of parasites and characterizing genes involved in host responses to infections are essential to improving our understanding of host-parasite interactions and their ecological and evolutionary consequences. This task, however, is complicated by high diversity and complex life histories of many parasite species. The use of transcriptomics in the context of wild-caught specimens can help ameliorate this by providing both qualitative and quantitative information on gene expression patterns in response to parasites in specific host organs and tissues. Here, we evaluated the physiological impact of the widespread parasite, the pike tapeworm (Triaenophorus nodulosus), on its second intermediate host, the Eurasian perch (Perca fluviatilis). We used an RNAseq approach to analyse gene expression in the liver, the target organ of T. nodulosus plerocercoids, and spleen which is one of the main immune organs in teleost fishes. We compared perch collected from multiple lakes consisting of individuals with (n = 8) and without (n = 6) T. nodulosus plerocercoids in the liver. Results revealed a small number of differentially expressed genes (DEGs, adjusted p-value ≤0.05) in both spleen (n = 22) and liver (n = 10). DEGs in spleen consisted of mostly upregulated immune related genes (e.g., JUN, SIK1, THSB1), while those in the liver were often linked to metabolic functions (e.g., FABP1, CADM4, CDAB). However, Gene Ontology (GO) analysis showed lack of functional enrichment among DEGs. This study demonstrates that Eurasian perch displays a subtle response at a gene expression level to T. nodulosus plerocercoid infection. Given that plerocercoids are low-metabolic activity transmission stages, our results suggest that moderate T. nodulosus plerocercoid infection most likely does not provoke an extensive host immune response and have relatively low physiological costs for the host. Our findings illustrate that not all conspicuous infections have severe effects on host gene regulation.

Effects of the invasive swim bladder parasite Anguillicola crassus on health and condition indicators in the European eel

Parasites negatively affect biological processes within their hosts, which may alter for example health, growth, and reproductive ability. Non-native invasive parasites, in particular, may have large effects on the endemic hosts, given that the hosts lack evolved specific defences against such parasites. The swim bladder nematode Anguillicola crassus, an invasive parasite originating from Asia, is found in the European eel (Anguilla anguilla, L. 1758), since the 1980s. We investigated whether A. crassus affected several indicators related to health of the European eel (spleen- and liver size, body fat content and relative condition). Our results indicate that during the continental residency of the eels, infection by A. crassus had no major negative impacts on the investigated health indicators at the generally low infection intensities present in this study (median 2-3 visible parasites). Given that many of the adult eels were found to have swim bladder damage, concerns about their spawning migration through deeper oceanic environments can still be raised. To allow further investigations, we suggest that quantification of swim bladder damage should be implemented in eel-monitoring programs. Compared to other parasite pressure parameters, swim bladder damage provides additional information about past infections and future problems.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Rat and fish peripheral blood leukocytes respond distinctively to Anisakis pegreffii (Nematoda, Anisakidae) crude extract

Infective third-stage larvae (L3) of the marine nematode Anisakis pegreffii cause inflammation and clinical symptoms in humans, their accidental host, that subside and self-resolve in a couple of weeks after L3 die. To characterise the differences in an early immune response of a marine vs. terrestrial host, we stimulated peripheral blood leukocytes (PBLs) of fish (paratenic host) and rat (accidental, human-model host) with A. pegreffii crude extract and analysed PBL transcriptomes 1 and 12 h post-stimulation. Fish and rat PBLs differentially expressed 712 and 493 transcripts, respectively, between 1 and 12 h post-stimulation (false discovery rate, FDR <0.001, logFC >2). While there was a difference in the highest upregulated transcripts between two time-points, the same Gene Ontologies, biological processes (intracellular signal transduction, DNA-dependent transcription, and DNA-regulated regulation of transcription), and molecular functions (ATP and metal ion binding) were enriched in the two hosts, showing an incrementing dynamic between 1 and 12 h. This suggests that the two distinct hosts employ qualitatively different transcript cascades only to achieve the same effect, at least during an early innate immunity response. Activation of later immunity elements and/or a combination of other host's intrinsic conditions may contribute to the death of L3 in the terrestrial host.

Genome-wide methylation in the panmictic European eel (Anguilla anguilla)

The role of methylation in adaptive, developmental and speciation processes has attracted considerable interest, but interpretation of results is complicated by diffuse boundaries between genetic and non-genetic variation. We studied whole genome genetic and methylation variation in the European eel, distributed from subarctic to subtropical environments, but with panmixia precluding genetically based local adaptation beyond single-generation responses. Overall methylation was 70.9%, with hypomethylation predominantly found in promoters and first exons. Redundancy analyses involving juvenile glass eels showed 0.06% and 0.03% of the variance at SNPs to be explained by localities and environmental variables, respectively, with GO terms of genes associated with outliers primarily involving neural system functioning. For CpGs 2.98% and 1.36% of variance was explained by localities and environmental variables. Differentially methylated regions particularly included genes involved in developmental processes, with hox clusters featuring prominently. Life stage (adult versus glass eels) was the most important source of inter-individual variation in methylation, likely reflecting both ageing and developmental processes. Demethylation of transposable elements relative to pure European eel was observed in European X American eel hybrids, possibly representing postzygotic barriers in this system characterized by prolonged speciation and ongoing gene flow. Whereas the genetic data are consistent with a role of single-generation selective responses, the methylation results underpin the importance of epigenetics in the life cycle of eels and suggests interactions between local environments, development and phenotypic variation mediated by methylation variation. Eels are remarkable by having retained eight hox clusters, and the results suggest important roles of methylation at hox genes for adaptive processes.

Assembly-free rapid differential gene expression analysis in non-model organisms using DNA-protein alignment

BACKGROUND

RNA-seq is being increasingly adopted for gene expression studies in a panoply of non-model organisms, with applications spanning the fields of agriculture, aquaculture, ecology, and environment. For organisms that lack a well-annotated reference genome or transcriptome, a conventional RNA-seq data analysis workflow requires constructing a de-novo transcriptome assembly and annotating it against a high-confidence protein database. The assembly serves as a reference for read mapping, and the annotation is necessary for functional analysis of genes found to be differentially expressed. However, assembly is computationally expensive. It is also prone to errors that impact expression analysis, especially since sequencing depth is typically much lower for expression studies than for transcript discovery.

RESULTS

We propose a shortcut, in which we obtain counts for differential expression analysis by directly aligning RNA-seq reads to the high-confidence proteome that would have been otherwise used for annotation. By avoiding assembly, we drastically cut down computational costs - the running time on a typical dataset improves from the order of tens of hours to under half an hour, and the memory requirement is reduced from the order of tens of Gbytes to tens of Mbytes. We show through experiments on simulated and real data that our pipeline not only reduces computational costs, but has higher sensitivity and precision than a typical assembly-based pipeline. A Snakemake implementation of our workflow is available at: https://bitbucket.org/project_samar/samar .

CONCLUSIONS

The flip side of RNA-seq becoming accessible to even modestly resourced labs has been that the time, labor, and infrastructure cost of bioinformatics analysis has become a bottleneck. Assembly is one such resource-hungry process, and we show here that it can be avoided for quick and easy, yet more sensitive and precise, differential gene expression analysis in non-model organisms.

Host gene expression in wildlife disease: making sense of species-level responses

Emerging infectious diseases are significant threats to wildlife conservation, yet the impacts of pathogen exposure and infection can vary widely among host species. As such, conservation biologists and disease ecologists have increasingly aimed to understand species-specific host susceptibility using molecular methods. In particular, comparative gene expression assays have been used to contrast the transcriptomic responses of disease-resistant and disease-susceptible hosts to pathogen exposure. This work usually assumes that the gene expression responses of disease-resistant species will reveal the activation of molecular pathways contributing to host defence. However, results often show that disease-resistant hosts undergo little gene expression change following pathogen challenge. Here, we discuss the mechanistic implications of these "null" findings and offer methodological suggestions for future molecular studies of wildlife disease. First, we highlight that muted transcriptomic responses with minimal immune system recruitment may indeed be protective for nonsusceptible hosts if they limit immunopathology and promote pathogen tolerance in systems where susceptible hosts suffer from genetic dysregulation. Second, we argue that overly narrow investigation of responses to pathogen exposure may overlook important, constitutively active molecular pathways that underlie species-specific defences. Finally, we outline alternative study designs and approaches that complement interspecific transcriptomic comparisons, including intraspecific gene expression studies and genomic methods to detect signatures of selection. Collectively, these insights will help ecologists extract maximal information from conservation-relevant transcriptomic data sets, leading to a deeper understanding of host defences and, ultimately, the implementation of successful conservation interventions.

Encapsulation of Anguillicola crassus reduces the abundance of adult parasite stages in the European eel (Anguilla anguilla)

Encapsulation of the parasitic nematode Anguillicola crassus Kuwahara, Niimi & Hagaki is commonly observed in its native host, the Japanese eel (Anguilla japonica Temminck & Schlegel). Encapsulation has also been described in a novel host, the European eel (A. anguilla L.), and there is evidence that encapsulation frequency has increased since the introduction of A. crassus. We examined whether encapsulation of A. crassus provides an advantage to its novel host in Lake Müggelsee, NE Germany. We provide the first evidence that encapsulation was associated with reduced abundance of adult A. crassus. This pattern was consistent in samples taken 3 months apart. There was no influence of infection on the expression of the two metabolic genes studied, but the number of capsules was negatively correlated with the expression of two mhc II genes of the adaptive immune response, suggesting a reduced activation. Interestingly, eels that encapsulated A. crassus had higher abundances of two native parasites compared with non-encapsulating eels. We propose that the response of A. anguilla to infection by A. crassus may interfere with its reaction to other co-occurring parasites.

Corrigendum

[This corrects the article DOI: 10.1002/ece3.5728.].