A palaeogenomic investigation of overharvest implications in an endemic wild reindeer subspecies.

Overharvest can severely reduce the abundance and distribution of a species and thereby impact its genetic diversity and threaten its future viability. Overharvest remains an ongoing issue for Arctic mammals, which due to climate change now also confront one of the fastest changing environments on Earth. The high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus), endemic to Svalbard, experienced a harvest-induced demographic bottleneck that occurred during the 17-20th centuries. Here, we investigate changes in genetic diversity, population structure, and gene-specific differentiation during and after this overharvesting event. Using whole-genome shotgun sequencing, we generated the first ancient and historical nuclear (n = 11) and mitochondrial (n = 18) genomes from Svalbard reindeer (up to 4000 BP) and integrated these data with a large collection of modern genome sequences (n = 90) to infer temporal changes. We show that hunting resulted in major genetic changes and restructuring in reindeer populations. Near-extirpation followed by pronounced genetic drift has altered the allele frequencies of important genes contributing to diverse biological functions. Median heterozygosity was reduced by 26%, while the mitochondrial genetic diversity was reduced only to a limited extent, likely due to already low pre-harvest diversity and a complex post-harvest recolonization process. Such genomic erosion and genetic isolation of populations due to past anthropogenic disturbance will likely play a major role in metapopulation dynamics (i.e., extirpation, recolonization) under further climate change. Our results from a high-arctic case study therefore emphasize the need to understand the long-term interplay of past, current, and future stressors in wildlife conservation.

The persistent threat of emerging plant disease pandemics to global food security.

Plant disease outbreaks are increasing and threaten food security for the vulnerable in many areas of the world. Now a global human pandemic is threatening the health of millions on our planet. A stable, nutritious food supply will be needed to lift people out of poverty and improve health outcomes. Plant diseases, both endemic and recently emerging, are spreading and exacerbated by climate change, transmission with global food trade networks, pathogen spillover, and evolution of new pathogen lineages. In order to tackle these grand challenges, a new set of tools that include disease surveillance and improved detection technologies including pathogen sensors and predictive modeling and data analytics are needed to prevent future outbreaks. Herein, we describe an integrated research agenda that could help mitigate future plant disease pandemics.

Genomic diversity and differentiation between island and mainland populations of white-tailed eagles (Haliaeetus albicilla).

Divergence in the face of high dispersal capabilities is a documented but poorly understood phenomenon. The white-tailed eagle (Haliaeetus albicilla) has a large geographic dispersal capability and should theoretically be able to maintain genetic homogeneity across its dispersal range. However, following analysis of the genomic variation of white-tailed eagles, from both historical and contemporary samples, clear signatures of ancient biogeographic substructure across Europe and the North-East Atlantic is observed. The greatest genomic differentiation was observed between island (Greenland and Iceland) and mainland (Denmark, Norway and Estonia) populations. The two island populations share a common ancestry from a single mainland population, distinct from the other sampled mainland populations, and despite the potential for high connectivity between Iceland and Greenland they are well separated from each other and are characterized by inbreeding and little variation. Temporal differences also highlight a pattern of regional populations persisting despite the potential for admixture. All sampled populations generally showed a decline in effective population size over time, which may have been shaped by four historical events: (1) Isolation of refugia during the last glacial period 110-115,000 years ago, (2) population divergence following the colonization of the deglaciated areas ~10,000 years ago, (3) human population expansion, which led to the settlement in Iceland ~1100 years ago, and (4) human persecution and exposure to toxic pollutants during the last two centuries.

The evolutionary history of extinct and living lions.

Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions (Panthera leo spelaea), 12 historic lions (Panthera leo leo/Panthera leo melanochaita) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus.

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

Runs of homozygosity in killer whale genomes provide a global record of demographic histories.

Runs of homozygosity (ROH) occur when offspring inherit haplotypes that are identical by descent from each parent. Length distributions of ROH are informative about population history; specifically, the probability of inbreeding mediated by mating system and/or population demography. Here, we investigated whether variation in killer whale (Orcinus orca) demographic history is reflected in genome-wide heterozygosity and ROH length distributions, using a global data set of 26 genomes representative of geographic and ecotypic variation in this species, and two F1 admixed individuals with Pacific-Atlantic parentage. We first reconstructed demographic history for each population as changes in effective population size through time using the pairwise sequential Markovian coalescent (PSMC) method. We found a subset of populations declined in effective population size during the Late Pleistocene, while others had more stable demography. Genomes inferred to have undergone ancestral declines in effective population size, were autozygous at hundreds of short ROH (<1 Mb), reflecting high background relatedness due to coalescence of haplotypes deep within the pedigree. In contrast, longer and therefore younger ROH (>1.5 Mb) were found in low latitude populations, and populations of known conservation concern. These include a Scottish killer whale, for which 37.8% of the autosomes were comprised of ROH >1.5 Mb in length. The fate of this population, in which only two adult males have been sighted in the past five years, and zero fecundity over the last two decades, may be inextricably linked to its demographic history and consequential inbreeding depression.

Will borealization of Arctic tundra herbivore communities be driven by climate warming or vegetation change?

Poleward shifts in species distributions are expected and frequently observed with a warming climate. In Arctic ecosystems, the strong warming trends are associated with increasing greenness and shrubification. Vertebrate herbivores have the potential to limit greening and shrub advance and expansion on the tundra, posing the question of whether changes in herbivore communities could partly mediate the impacts of climate warming on Arctic tundra. Therefore, future changes in the herbivore community in the Arctic tundra will depend on whether the community tracks the changing climates directly (i.e. occurs in response to temperature) or indirectly, in response to vegetation changes (which can be modified by trophic interactions). In this study, we used biogeographic and remotely sensed data to quantify spatial variation in vertebrate herbivore communities across the boreal forest and Arctic tundra biomes. We then tested whether present-day herbivore community structure is determined primarily by temperature or vegetation. We demonstrate that vertebrate herbivore communities are significantly more diverse in the boreal forest than in the Arctic tundra in terms of species richness, phylogenetic diversity and functional diversity. A clear shift in community structure was observed at the biome boundary, with stronger northward declines in diversity in the Arctic tundra. Interestingly, important functional traits characterizing the role of herbivores in limiting tundra vegetation change, such as body mass and woody plant feeding, did not show threshold changes across the biome boundary. Temperature was a more important determinant of herbivore community structure across these biomes than vegetation productivity or woody plant cover. Thus, our study does not support the premise that herbivore-driven limitation of Arctic tundra shrubification or greening would limit herbivore community change in the tundra. Instead, borealization of tundra herbivore communities is likely to result from the direct effect of climate warming.

Killer whale genomes reveal a complex history of recurrent admixture and vicariance.

Reconstruction of the demographic and evolutionary history of populations assuming a consensus tree-like relationship can mask more complex scenarios, which are prevalent in nature. An emerging genomic toolset, which has been most comprehensively harnessed in the reconstruction of human evolutionary history, enables molecular ecologists to elucidate complex population histories. Killer whales have limited extrinsic barriers to dispersal and have radiated globally, and are therefore a good candidate model for the application of such tools. Here, we analyse a global data set of killer whale genomes in a rare attempt to elucidate global population structure in a nonhuman species. We identify a pattern of genetic homogenisation at lower latitudes and the greatest differentiation at high latitudes, even between currently sympatric lineages. The processes underlying the major axis of structure include high drift at the edge of species' range, likely associated with founder effects and allelic surfing during postglacial range expansion. Divergence between Antarctic and non-Antarctic lineages is further driven by ancestry segments with up to four-fold older coalescence time than the genome-wide average; relicts of a previous vicariance during an earlier glacial cycle. Our study further underpins that episodic gene flow is ubiquitous in natural populations, and can occur across great distances and after substantial periods of isolation between populations. Thus, understanding the evolutionary history of a species requires comprehensive geographic sampling and genome-wide data to sample the variation in ancestry within individuals.

Herbarium collection-based phylogenetics of the ragweeds (Ambrosia, Asteraceae).

Ambrosia (Asteraceae) is a taxonomically difficult genus of weedy, wind-pollinated plants with an apparent center of diversity in the Sonoran Desert of North America. Determining Ambrosia's evolutionary relationships has been the subject of much interest, with numerous studies using morphological characters, cytology, comparative phytochemistry, and chloroplast restriction site variation to produce conflicting accounts the relationships between Ambrosia species, as well as the classification of their close relatives in Franseria and Hymenoclea. To resolve undetermined intra-generic relationships within Ambrosia, we used DNA extracted from tissues obtained from seed banks and herbarium collections to generate multi-locus genetic data representing nearly all putative species, including four from South America. We performed Bayesian and Maximum-Likelihood phylogenetic analyses of six chloroplast-genome and two nuclear-genome markers, enabling us to infer monophyly for the genus, resolve major infra-generic species clusters, as well as to resolve open questions about the evolutionary relationships of several Ambrosia species and former members of Franseria. We also provide molecular data supporting the hypothesis that A. sandersonii formed through the hybridization of A. eriocentra and A. salsola. The topology of our chloroplast DNA phylogeny is almost entirely congruent with the most recent molecular work based on chloroplast restriction site variation of a much more limited sampling of 14 North American species of Ambrosia, although our improved sampling of global Ambrosia diversity enables us to draw additional conclusions. As our study is the first direct DNA sequence-based phylogenetic analyses of Ambrosia, we analyze the data in relation to previous taxonomic studies and discuss several instances of chloroplast/nuclear incongruence that leave the precise geographic center of origin of Ambrosia in question.

Genetic polymorphisms of GRIN2A and GRIN2B modify the neurobehavioral effects of low-level lead exposure in children.

Lead (Pb) is neurotoxic and children are highly susceptible to this effect, particularly within the context of continuous low-level Pb exposure. A current major challenge is identification of children who may be uniquely susceptible to Pb toxicity because of genetic predisposition. Learning and memory are among the neurobehavioral processes that are most notably affected by Pb exposure, and modification of N-methyl-D-aspartate receptors (NMDAR) that regulate these processes during development are postulated to underlie these adverse effects of Pb. We examined the hypothesis that polymorphic variants of genes encoding glutamate receptor, ionotropic, NMDAR subunits 2A and 2B, GRIN2A and GRIN2B, exacerbate the adverse effects of Pb exposure on these processes in children. Participants were subjects who participated as children in the Casa Pia Dental Amalgam Clinical Trial and for whom baseline blood Pb concentrations and annual neurobehavioral test results over the 7 year course of the clinical trial were available. Genotyping assays were performed for variants of GRIN2A (rs727605 and rs1070503) and GRIN2B (rs7301328 and rs1806201) on biological samples acquired from 330 of the original 507 trial participants. Regression modeling strategies were employed to evaluate the association between genotype status, Pb exposure, and neurobehavioral test outcomes. Numerous significant adverse interaction effects between variants of both GRIN2A and GRIN2B, individually and in combination, and Pb exposure were observed particularly among boys, preferentially within the domains of Learning & Memory and Executive Function. In contrast, very few interaction effects were observed among similarly genotyped girls with comparable Pb exposure. These findings support observations of an essential role of GRIN2A and GRIN2B on developmental processes underlying learning and memory as well as other neurological functions in children and demonstrate, further, modification of Pb effects on these processes by specific variants of both GRIN2A and GRIN2B genes. These observations highlight the importance of genetic factors in defining susceptibility to Pb neurotoxicity and may have important public health implications for future strategies aimed at protecting children and adolescents from potential health risks associated with low-level Pb exposure.

Significant upregulation of U1 and U4 spliceosomal snRNAs by ATP nanoliposomes explains acceleration of wound healing, due to increased pre-mRNA processing to functional mRNA.

Delayed wound healing is one of the hallmarks of diabetic complications and certain autoimmune inflammatory diseases. Extensive wound healing studies in rabbits have indicated that the delivery of ATP encapsulated in unilamellar nanoliposomes causes rapid cell proliferation and fast tracks the wound healing process. In the current study, we explored the possible molecular mechanism underlying this response by comparing gene expression in cultured rabbit kidney cells treated with either ATP nanoliposomes (containing 1 mg Mg-ATP/ml formulation) or control nanoliposomes (containing 1 mg/ml unmetabolisable gamma-thio-ATP/ml formulation). High-quality total RNA was isolated 24 h from the cells and subjected to RNA seq technology, which revealed significant overexpression of specific noncoding RNAs. The U1 spliceosomal RNA, U1 snRNA, was upregulated more than 250-fold following treatment with ATP nanoliposomes. This multifunctional U1 spliceosomal RNA may function in transcription by speeding up the critical splicing step, thereby facilitating faster processing of pre-mRNA to translatable mRNA.

Genomic Characterization of a South American Phytophthora Hybrid Mandates Reassessment of the Geographic Origins of Phytophthora infestans.

As the oomycete pathogen causing potato late blight disease, Phytophthora infestans triggered the famous 19th-century Irish potato famine and remains the leading cause of global commercial potato crop destruction. But the geographic origin of the genotype that caused this devastating initial outbreak remains disputed, as does the New World center of origin of the species itself. Both Mexico and South America have been proposed, generating considerable controversy. Here, we readdress the pathogen's origins using a genomic data set encompassing 71 globally sourced modern and historical samples of P. infestans and the hybrid species P. andina, a close relative known only from the Andean highlands. Previous studies have suggested that the nuclear DNA lineage behind the initial outbreaks in Europe in 1845 is now extinct. Analysis of P. andina's phased haplotypes recovered eight haploid genome sequences, four of which represent a previously unknown basal lineage of P. infestans closely related to the famine-era lineage. Our analyses further reveal that clonal lineages of both P. andina and historical P. infestans diverged earlier than modern Mexican lineages, casting doubt on recent claims of a Mexican center of origin. Finally, we use haplotype phasing to demonstrate that basal branches of the clade comprising Mexican samples are occupied by clonal isolates collected from wild Solanum hosts, suggesting that modern Mexican P. infestans diversified on Solanum tuberosum after a host jump from a wild species and that the origins of P. infestans are more complex than was previously thought.

Determination of genetic relatedness from low-coverage human genome sequences using pedigree simulations.

We develop and evaluate methods for inferring relatedness among individuals from low-coverage DNA sequences of their genomes, with particular emphasis on sequences obtained from fossil remains. We suggest the major factors complicating the determination of relatedness among ancient individuals are sequencing depth, the number of overlapping sites, the sequencing error rate and the presence of contamination from present-day genetic sources. We develop a theoretical model that facilitates the exploration of these factors and their relative effects, via measurement of pairwise genetic distances, without calling genotypes, and determine the power to infer relatedness under various scenarios of varying sequencing depth, present-day contamination and sequencing error. The model is validated by a simulation study as well as the analysis of aligned sequences from present-day human genomes. We then apply the method to the recently published genome sequences of ancient Europeans, developing a statistical treatment to determine confidence in assigned relatedness that is, in some cases, more precise than previously reported. As the majority of ancient specimens are from animals, this method would be applicable to investigate kinship in nonhuman remains. The developed software grups (Genetic Relatedness Using Pedigree Simulations) is implemented in Python and freely available.

Epigenetic divergence as a potential first step in darter speciation.

Recent studies show that epigenetic variation in the form of DNA methylation may serve as a substrate for selection. Theory suggests that heritable epigenetic marks that increase fitness should increase in frequency in a population, and these changes may result in novel morphology, behaviour, or physiology, and ultimately reproductive isolation. Therefore, epigenetic variation might provide the first substrate for selection during the course of evolutionary divergence. This hypothesis predicts that populations in the earliest stages of divergence will differentiate in their methylome prior to any genetic differentiation. While several studies have investigated natural epigenetic variation, empirical studies that test predictions about its role in speciation are surprisingly scarce. Here, we investigate DNA methylation variation using an isoschizomeric digest method, Methyl-Sensitive Amplified Polymorphism, across multiple stages of evolutionary divergence in natural populations of North American stream fishes. We show that epigenetic differentiation between methylomes is greater than genetic divergence among closely related populations across two river drainages. Additionally, we demonstrate that epigenetic divergence is a stronger predictor of the strength of behavioural reproductive isolation and suggest that changes in the methylome could influence the evolution of reproductive isolation between species. Our findings suggest a role for epigenetics not only in the initiation of divergence, but also in the maintenance of species boundaries over greater evolutionary timescales.

Persistence of the mitochondrial lineage responsible for the Irish potato famine in extant new world phytophthora infestans.

The plant pathogen Phytophthora infestans emerged in Europe in 1845, triggering the Irish potato famine and massive European potato crop losses that continued until effective fungicides were widely employed in the 20th century. Today the pathogen is ubiquitous, with more aggressive and virulent strains surfacing in recent decades. Recently, complete P. infestans mitogenome sequences from 19th-century herbarium specimens were shown to belong to a unique lineage (HERB-1) predicted to be rare or extinct in modern times. We report 44 additional P. infestans mitogenomes: four from 19th-century Europe, three from 1950s UK, and 37 from modern populations across the New World. We use phylogenetic analyses to identify the HERB-1 lineage in modern populations from both Mexico and South America, and to demonstrate distinct mitochondrial haplotypes were present in 19th-century Europe, with this lineage initially diversifying 75 years before the first reports of potato late blight.

Mutation-order divergence by sexual selection: diversification of sexual signals in similar environments as a first step in speciation.

The origin of species remains a central question, and recent research focuses on the role of ecological differences in promoting speciation. Ecological differences create opportunities for divergent selection (i.e. 'ecological' speciation), a Darwinian hypothesis that hardly requires justification. In contrast, 'mutation-order' speciation proposes that, instead of adapting to different environments, populations find different ways to adapt to similar environments, implying that speciation does not require ecological differences. This distinction is critical as it provides an alternative hypothesis to the prevailing view that ecological differences drive speciation. Speciation by sexual selection lies at the centre of debates about the importance of ecological differences in promoting speciation; here, we present verbal and mathematical models of mutation-order divergence by sexual selection. We develop three general cases and provide a two-locus population genetic model for each. Results indicate that alternative secondary sexual traits can fix in populations that initially experience similar natural and sexual selection and that divergent traits and preferences can remain stable in the face of low gene flow. This stable divergence can facilitate subsequent divergence that completes or reinforces speciation. We argue that a mutation-order process could explain widespread diversity in secondary sexual traits among closely related, allopatric species.

Herbarium specimens reveal a historical shift in phylogeographic structure of common ragweed during native range disturbance.

Invasive plants provide ample opportunity to study evolutionary shifts that occur after introduction to novel environments. However, although genetic characters pre-dating introduction can be important determinants of later success, large-scale investigations of historical genetic structure have not been feasible. Common ragweed (Ambrosia artemisiifolia L.) is an invasive weed native to North America that is known for its allergenic pollen. Palynological records from sediment cores indicate that this species was uncommon before European colonization of North America, and ragweed populations expanded rapidly as settlers deforested the landscape on a massive scale, later becoming an aggressive invasive with populations established globally. Towards a direct comparison of genetic structure now and during intense anthropogenic disturbance of the late 19th century, we sampled 45 natural populations of common ragweed across its native range as well as historical herbarium specimens collected up to 140 years ago. Bayesian clustering analyses of 453 modern and 473 historical samples genotyped at three chloroplast spacer regions and six nuclear microsatellite loci reveal that historical ragweed's spatial genetic structure mirrors both the palaeo-record of Ambrosia pollen deposition and the historical pattern of agricultural density across the landscape. Furthermore, for unknown reasons, this spatial genetic pattern has changed substantially in the intervening years. Following on previous work relating morphology and genetic expression between plants collected from eastern North America and Western Europe, we speculate that the cluster associated with humans' rapid transformation of the landscape is a likely source of these aggressive invasive populations.

Genetic polymorphisms of catechol-O-methyltransferase modify the neurobehavioral effects of mercury in children.

Mercury (Hg) is neurotoxic and children may be particularly susceptible to this effect. A current major challenge is identification of children who may be uniquely susceptible to Hg toxicity because of genetic disposition. This study examined the hypothesis that genetic variants of catechol-O-methyltransferase (COMT) that are reported to alter neurobehavioral functions that are also affected by Hg in adults might modify the adverse neurobehavioral effects of Hg exposure in children. Five hundred and seven children, 8-12 yr of age at baseline, participated in a clinical trial to evaluate the neurobehavioral effects of Hg from dental amalgam tooth fillings. Subjects were evaluated at baseline and at seven subsequent annual intervals for neurobehavioral performance and urinary Hg levels. Following the clinical trial, genotyping assays were performed for single-nucleotide polymorphisms (SNPs) of COMT rs4680, rs4633, rs4818, and rs6269 on biological samples provided by 330 of the trial participants. Regression-modeling strategies were employed to evaluate associations between allelic status, Hg exposure, and neurobehavioral test outcomes. Similar analysis was performed using haplotypes of COMT SNPs. Among girls, few interactions for Hg exposure and COMT variants were found. In contrast, among boys, numerous gene-Hg interactions were observed between individual COMT SNPs, as well as with a common COMT haplotype affecting multiple domains of neurobehavioral function. These findings suggest increased susceptibility to the adverse neurobehavioral effects of Hg among children with common genetic variants of COMT, and may have important implications for strategies aimed at protecting children from the potential health risks associated with Hg exposure.

Ancient reindeer mitogenomes reveal island-hopping colonisation of the Arctic archipelagos.

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Host-gut microbiota interactions shape parasite infections in farmed Atlantic salmon.

Animals and their associated microbiota share long evolutionary histories. However, it is not always clear how host genotype and microbiota interact to affect phenotype. We applied a hologenomic approach to explore how host-microbiota interactions shape lifetime growth and parasite infection in farmed Atlantic salmon (Salmo salar). Multi-omics data sets were generated from the guts of 460 salmon, 82% of which were naturally infected with an intestinal cestode. A single Mycoplasma bacterial strain, MAG01, dominated the gut metagenome of large, non-parasitized fish, consistent with previous studies showing high levels of Mycoplasma in the gut microbiota of healthy salmon. While small and/or parasitized salmon also had high abundance of MAG01, we observed increased alpha diversity in these individuals, driven by increased frequency of low-abundance Vibrionaceae and other Mycoplasma species that carried known virulence genes. Colonization by one of these cestode-associated Mycoplasma strains was associated with host individual genomic variation in long non-coding RNAs. Integrating the multi-omic data sets revealed coordinated changes in the salmon gut mRNA transcriptome and metabolome that correlated with shifts in the microbiota of smaller, parasitized fish. Our results suggest that the gut microbiota of small and/or parasitized fish is in a state of dysbiosis that partly depends on the host genotype, highlighting the value of using a hologenomic approach to incorporate the microbiota into the study of host-parasite dynamics.IMPORTANCEStudying host-microbiota interactions through the perspective of the hologenome is gaining interest across all life sciences. Intestinal parasite infections are a huge burden on human and animal health; however, there are few studies investigating the role of the hologenome during parasite infections. We address this gap in the largest multi-omics fish microbiota study to date using natural cestode infection of farmed Atlantic salmon. We find a clear association between cestode infection, salmon lifetime growth, and perturbation of the salmon gut microbiota. Furthermore, we provide the first evidence that the genetic background of the host may partly determine how the gut microbiota changes during parasite-associated dysbiosis. Our study therefore highlights the value of a hologenomic approach for gaining a more in-depth understanding of parasitism.