The myxozoan Tetracapsuloides bryosalmonae is a widely spread endoparasite that causes proliferative kidney disease (PKD) in salmonid fish. We developed an in silico pipeline to separate transcripts of T. bryosalmonae from the kidney tissue of its natural vertebrate host, brown trout (Salmo trutta). After stringent filtering, we constructed a partial transcriptome assembly T. bryosalmonae, comprising 3427 transcripts. Based on homology-restricted searches of the assembled parasite transcriptome and Atlantic salmon (Salmo salar) proteome, we identified four protein targets (Endoglycoceramidase, Legumain-like protease, Carbonic anhydrase 2, Pancreatic lipase-related protein 2) for the development of anti-parasitic drugs against T. bryosalmonae. Earlier work of these proteins on parasitic protists and helminths suggests that the identified anti-parasitic drug targets represent promising chemotherapeutic candidates also against T. bryosalmonae, and strengthen the view that the known inhibitors can be effective in evolutionarily distant organisms. In addition, we identified differentially expressed T. bryosalmonae genes between moderately and severely infected fish, indicating an increased abundance of T. bryosalmonae sporogonic stages in fish with low parasite load. In conclusion, this study paves the way for future genomic research in T. bryosalmonae and represents an important step towards the development of effective drugs against PKD.

The extent of geographic genetic variation is the result of several processes such as mutation, gene flow, selection and drift. Processes that structure the populations of parasite species are often directly linked to the processes that influence the host. Here, we investigate the genetic population structure of the ectoparasite Gyrodactylus thymalli Žitňan, 1960 (Monogenea) collected from grayling (Thymallus thymallus L.) throughout the river Glomma, the largest watercourse in Norway. Parts of the mitochondrial dehydrogenase subunit 5 (NADH 5) and cytochrome oxidase I (COI) genes from 309 G. thymalli were analysed to study the genetic variation and investigated the geographical distribution of parasite haplotypes. Three main clusters of haplotypes dominated the three distinct geographic parts of the river system; one cluster dominated in the western main stem of the river, one in the eastern and one in the lower part. There was a positive correlation between pairwise genetic distance and hydrographic distance. The results indicate restricted gene flow between sub-populations of G. thymalli, most likely due to barriers that limit upstream migration of infected grayling. More than 80% of the populations had private haplotypes, also indicating long-time isolation of sub-populations. According to a molecular clock calibration, much of the haplotype diversity of G. thymalli in the river Glomma has developed after the last glaciation.