When parasites persist: tapeworms survive host extinction and reveal waves of dispersal across Beringia

Phylogenomics of American pika (Ochotona princeps) lineage diversification

Quaternary climate oscillations have profoundly influenced current species distributions. For many montane species, these fluctuations were a prominent driver in species range shifts, often resulting in intraspecific diversification, as has been the case for American pikas (Ochotona princeps). Range shifts and population declines in this thermally-sensitive lagomorph have been linked to historical and contemporary environmental changes across its western North American range, with previous research reconstructing five mitochondrial DNA lineages. Here, we paired genome-wide data (25,244 SNPs) with range-wide sampling to re-examine the number and distribution of intra-specific lineages, and investigate patterns of within- and among-lineage divergence and diversity. Our results provide genomic evidence of O. princeps monophyly, reconstructing six distinct lineages that underwent multiple rounds of divergence (0.809-2.81 mya), including a new Central Rocky Mountain lineage. We further found evidence for population differentiation across multiple spatial scales, and reconstructed levels of standing variation comparable to those found in other small mammals. Overall, our findings demonstrate the influence of past glacial cycles on O. princeps lineage diversification, suggest that current subspecific taxonomy may need to be revisited, and provide an important framework for investigations of American pika adaptive potential in the face of anthropogenic climate change.

Historical dispersal and host-switching formed the evolutionary history of a globally distributed multi-host parasite - The Ligula intestinalis species complex

Studies on parasite biogeography and host spectrum provide insights into the processes driving parasite diversification. Global geographical distribution and a multi-host spectrum make the tapeworm Ligula intestinalis a promising model for studying both the vicariant and ecological modes of speciation in parasites. To understand the relative importance of host association and biogeography in the evolutionary history of this tapeworm, we analysed mtDNA and reduced-represented genomic SNP data for a total of 139 specimens collected from 18 fish-host genera across a distribution range representing 21 countries. Our results strongly supported the existence of at least 10 evolutionary lineages and estimated the deepest divergence at approximately 4.99-5.05 Mya, which is much younger than the diversification of the fish host genera and orders. Historical biogeography analyses revealed that the ancestor of the parasite diversified following multiple vicariance events and was widespread throughout the Palearctic, Afrotropical, and Nearctic between the late Miocene and early Pliocene. Cyprinoids were inferred as the ancestral hosts for the parasite. Later, from the late Pliocene to Pleistocene, new lineages emerged following a series of biogeographic dispersal and host-switching events. Although only a few of the current Ligula lineages show narrow host-specificity (to a single host genus), almost no host genera, even those that live in sympatry, overlapped between different Ligula lineages. Our analyses uncovered the impact of historical distribution shifts on host switching and the evolution of host specificity without parallel host-parasite co-speciation. Historical biogeography reconstructions also found that the parasite colonized several areas (Afrotropical and Australasian) much earlier than was suggested by only recent faunistic data.

Genetic drift drives rapid speciation of an Arctic insular endemic shrew (Sorex pribilofensis)

Episodes of Quaternary environmental change shaped the genomes of extant species, influencing their response to contemporary environments, which are changing rapidly. Island endemics are among the most vulnerable to such change, accounting for a disproportionate number of recent extinctions. To prevent extinctions and conserve island biodiversity it is vital to combine knowledge of species' ecologies with their complex evolutionary histories. The Bering Sea has a history of cyclical island isolation and reconnection, coupled with modern rates of climate change that exceed global averages. The endangered Pribilof Island shrew (Sorex pribilofensis) is endemic to St. Paul Island, Alaska, which was isolated from mainland Beringia ~14,000 years ago by rising sea levels. Using ~11,000 single nucleotide polymorphisms, 17 microsatellites and mitochondrial sequence data, we test predictions about the evolutionary processes driving shrew speciation across Beringia. Our data show considerable differentiation of S. pribilofensis from mainland sibling species, relative to levels of divergence between mainland shrews. We also find a genome-wide loss of diversity and extremely low N_e for S. pribilofensis. We then show that intraspecific genetic diversity is significantly related to interspecific divergence, and that differentiation between S. pribilofensis and other Beringian shrews is highest across loci that are fixed in S. pribilofensis, indicating that strong drift has driven differentiation of this island species. Our findings show that drift as a consequence of Arctic climate cycling can rapidly reshape insular biodiversity. Arctic island species that lack genomic diversity and have evolved in response to past climate may have limited ability to respond to modern environmental changes.

Chromosome-Level Reference Genome Assembly for the American Pika (Ochotona princeps)

The American pika (Ochotona princeps) is an alpine lagomorph found throughout western North America. Primarily inhabiting talus slopes at higher elevations (>2000 m), American pikas are well adapted to cold, montane environments. Warming climates on both historical and contemporary scales have contributed to population declines in American pikas, positioning them as a focal mammalian species for investigating the ecological effects of climate change. To support and expand ongoing research efforts, here, we present a highly contiguous and annotated reference genome assembly for the American pika (OchPri4.0). This assembly was produced using Dovetail de novo proximity ligation methods and annotated through the NCBI Eukaryotic Genome Annotation pipeline. The resulting assembly was chromosome- scale, with a total length of 2.23 Gb across 9350 scaffolds and a scaffold N50 of 75.8 Mb. The vast majority (>97%) of the total assembly length was found within 36 large scaffolds; 33 of these scaffolds correlated to whole autosomes, while the X chromosome was covered by 3 large scaffolds. Additionally, we identified 17 enriched gene ontology terms among American pika-specific genes putatively related to adaptation to high-elevation environments. This high-quality genome assembly will serve as a springboard for exploring the evolutionary underpinnings of behavioral, ecological, and taxonomic diversification in pikas as well as broader-scale eco-evolutionary questions pertaining to cold-adapted species in general.