First Insights into Population Structure and Genetic Diversity Versus Host Specificity in Trypanorhynch Tapeworms Using Multiplexed Shotgun Genotyping

Abstract

Theory predicts relaxed host specificity and high host vagility should contribute to reduced genetic structure in parasites while strict host specificity and low host vagility should increase genetic structure. Though these predictions are intuitive, they have never been explicitly tested in a population genomic framework. Trypanorhynch tapeworms, which parasitize sharks and rays (elasmobranchs) as definitive hosts, are the only order of elasmobranch tapeworms that exhibit considerable variability in their definitive host specificity. This allows for unique combinations of host use and geographic range, making trypanorhynchs ideal candidates for studying how these traits influence population-level structure and genetic diversity. Multiplexed shotgun genotyping (MSG) data sets were generated to characterize component population structure and infrapopulation diversity for a representative of each trypanorhynch suborder: the ray-hosted Rhinoptericola megacantha (Trypanobatoida) and the shark-hosted Callitetrarhynchus gracilis (Trypanoselachoida). Adults of R. megacantha are more host-specific and less broadly distributed than adults of C. gracilis, allowing correlation between these factors and genetic structure. Replicate tapeworm specimens were sequenced from the same host individual, from multiple conspecific hosts within and across geographic regions, and from multiple definitive host species. For R. megacantha, population structure coincided with geography rather than host species. For C. gracilis, limited population structure was found, suggesting a potential link between degree of host specificity and structure. Conspecific trypanorhynchs from the same host individual were found to be as, or more, genetically divergent from one another as from conspecifics from different host individuals. For both species, high levels of homozygosity and positive FIS values were documented.

Keywords: Callitetrarhynchus gracilis; Rhinoptericola megacantha; RAD-seq; Trypanorhyncha; elasmobranch tapeworms; population genomics.