The gut microbiota-immune-brain axis in a wild vertebrate: dynamic interactions and health impacts.

The gut microbiota-immune-brain axis is a feedback network which influences diverse physiological processes and plays a pivotal role in overall health and wellbeing. Although research in humans and laboratory mice has shed light into the associations and mechanisms governing this communication network, evidence of such interactions in wild, especially in young animals, is lacking. We therefore investigated these interactions during early development in a population of common buzzards (Buteo buteo) and their effects on individual condition. In a longitudinal study, we used a multi-marker approach to establish potential links between the bacterial and eukaryotic gut microbiota, a panel of immune assays and feather corticosterone measurements as a proxy for long-term stress. Using Bayesian structural equation modeling, we found no support for feedback between gut microbial diversity and immune or stress parameters. However, we did find strong relationships in the feedback network. Immunity was negatively correlated with corticosterone levels, and microbial diversity was positively associated with nestling body condition. Furthermore, corticosterone levels and eukaryotic microbiota diversity decreased with age while immune activity increased. The absence of conclusive support for the microbiota-immune-brain axis in common buzzard nestlings, coupled with the evidence for stress mediated immunosuppression, suggests a dominating role of stress-dominated maturation of the immune system during early development. Confounding factors inherent to wild systems and developing animals might override associations known from adult laboratory model subjects. The positive association between microbial diversity and body condition indicates the potential health benefits of possessing a diverse and stable microbiota.

Genomic and fitness consequences of a near-extinction event in the northern elephant seal.

Understanding the genetic and fitness consequences of anthropogenic bottlenecks is crucial for biodiversity conservation. However, studies of bottlenecked populations combining genomic approaches with fitness data are rare. Theory predicts that severe bottlenecks deplete genetic diversity, exacerbate inbreeding depression and decrease population viability. However, actual outcomes are complex and depend on how a species' unique demography affects its genetic load. We used population genetic and veterinary pathology data, demographic modelling, whole-genome resequencing and forward genetic simulations to investigate the genomic and fitness consequences of a near-extinction event in the northern elephant seal. We found no evidence of inbreeding depression within the contemporary population for key fitness components, including body mass, blubber thickness and susceptibility to parasites and disease. However, we detected a genomic signature of a recent extreme bottleneck (effective population size = 6; 95% confidence interval = 5.0-7.5) that will have purged much of the genetic load, potentially leading to the lack of observed inbreeding depression in our study. Our results further suggest that deleterious genetic variation strongly impacted the post-bottleneck population dynamics of the northern elephant seal. Our study provides comprehensive empirical insights into the intricate dynamics underlying species-specific responses to anthropogenic bottlenecks.

Characterizing offshore polar ocean soundscapes using ecoacoustic intensity and diversity metrics.

Polar offshore environments are considered the last pristine soundscapes, but accelerating climate change and increasing human activity threaten their integrity. In order to assess the acoustic state of polar oceans, there is the need to investigate their soundscape characteristics more holistically. We apply a set of 14 ecoacoustic metrics (EAMs) to identify which metrics are best suited to reflect the characteristics of disturbed and naturally intact polar offshore soundscapes. We used two soundscape datasets: (i) the Arctic eastern Fram Strait (FS), which is already impacted by anthropogenic noise, and (ii) the quasi-pristine Antarctic Weddell Sea (WS). Our results show that EAMs when applied in concert can be used to quantitatively assess soundscape variability, enabling the appraisal of marine soundscapes over broad spatiotemporal scales. The tested set of EAMs was able to show that the eastern FS, which is virtually free from sea ice, lacks seasonal soundscape dynamics and exhibits low acoustic complexity owing to year-round wind-mediated sounds and anthropogenic noise. By contrast, the WS exhibits pronounced seasonal soundscape dynamics with greater soundscape heterogeneity driven in large part by the vocal activity of marine mammal communities, whose composition in turn varies with the prevailing seasonal sea ice conditions.

Refinement of the Antarctic fur seal (Arctocephalus gazella) reference genome increases continuity and completeness.

The Antarctic fur seal (Arctocephalus gazella) is an important top predator and indicator of the health of the Southern Ocean ecosystem. Although abundant, this species narrowly escaped extinction due to historical sealing and is currently declining as a consequence of climate change. Genomic tools are essential for understanding these anthropogenic impacts and for predicting long-term viability. However, the current reference genome ("arcGaz3") shows considerable room for improvement in terms of both completeness and contiguity. We therefore combined PacBio sequencing, haplotype-aware HiRise assembly and scaffolding based on Hi-C information to generate a refined assembly of the Antarctic fur seal reference genome ("arcGaz4_h1"). The new assembly is 2.53Gb long, has a scaffold N50 of 55.6Mb and includes 18 chromosome-sized scaffolds, which correspond to the 18 chromosomes expected in otariids. Genome completeness is greatly improved, with 23,408 annotated genes and a Benchmarking Universal Single-Copy Orthologs (BUSCO) score raised from 84.7% to 95.2%. We furthermore included the new genome in a reference-free alignment of the genomes of eleven pinniped species to characterize evolutionary conservation across the Pinnipedia using genome-wide Genomic Evolutionary Rate Profiling (GERP). We then implemented Gene Ontology (GO) enrichment analyses to identify biological processes associated with those genes showing the highest levels of either conservation or differentiation between the two major pinniped families, the Otariidae and Phocidae. We show that processes linked to neuronal development, the circulatory system and osmoregulation are overrepresented both in conserved as well as in differentiated regions of the genome.

Runs of homozygosity reveal contrasting histories of inbreeding across global lineages of the edible porcini mushroom, Boletus edulis.

Inbreeding, the mating of individuals that are related through common ancestry, is of central importance in evolutionary and conservation biology due to its impacts on individual fitness and population dynamics. However, while advanced genomic approaches have revolutionised the study of inbreeding in animals, genomic studies of inbreeding are rare in plants and lacking in fungi. We investigated global patterns of inbreeding in the prized edible porcini mushroom Boletus edulis using 225 whole genomes from seven lineages distributed across the northern hemisphere. Genomic inbreeding was quantified using runs of homozygosity (ROHs). We found appreciable variation both among and within lineages, with some individuals having over 20% of their genomes in ROHs. Much of this variation could be explained by a combination of elevation and latitude, and to a lesser extent by predicted habitat suitability during the last glacial maximum. In line with this, the majority of ROHs were short, reflecting ancient common ancestry dating back approximately 200-1700 generations ago, while longer ROHs indicative of recent common ancestry (less than approximately 50 generations ago) were infrequent. Our study reveals the inbreeding legacy of major climatic events in a widely distributed forest mutualist, aligning with prevailing theories and empirical studies of the impacts of historical glaciation events on the dominant forest tree species of the northern hemisphere.

Life-history stage influences immune investment and oxidative stress in response to environmental heterogeneity in Antarctic fur seals.

Immune defenses are crucial for survival but costly to develop and maintain. Increased immune investment is therefore hypothesized to trade-off with other life-history traits. Here, we examined innate and adaptive immune responses to environmental heterogeneity in wild Antarctic fur seals. In a fully crossed, repeated measures design, we sampled 100 pups and their mothers from colonies of contrasting density during seasons of contrasting food availability. Biometric and cortisol data as well as blood for the analysis of 13 immune and oxidative status markers were collected at two key life-history stages. We show that immune responses of pups are more responsive than adults to variation in food availability, but not population density, and are modulated by cortisol and condition. Immune investment is associated with different oxidative status markers in pups and mothers. Our results suggest that early life stages show greater sensitivity to extrinsic and intrinsic effectors, and that immunity may be a strong target for natural selection even in low-pathogen environments such as Antarctica.

Early-life factors shaping the gut microbiota of Common buzzard nestlings.

BACKGROUND: Exploring the dynamics of gut microbiome colonisation during early-life stages is important for understanding the potential impact of microbes on host development and fitness. Evidence from model organisms suggests a crucial early-life phase when shifts in gut microbiota can lead to immune dysregulation and reduced host condition. However, our understanding of gut microbiota colonisation in long-lived vertebrates, especially during early development, remains limited. We therefore used a wild population of common buzzard nestlings (Buteo buteo) to investigate connections between the early-life gut microbiota colonisation, environmental and host factors. RESULTS: We targeted both bacterial and eukaryotic microbiota using the 16S and 28S rRNA genes. We sampled the individuals during early developmental stages in a longitudinal design. Our data revealed that age significantly affected microbial diversity and composition. Nest environment was a notable predictor of microbiota composition, with particularly eukaryotic communities differing between habitats occupied by the hosts. Nestling condition and infection with the blood parasite Leucocytozoon predicted microbial community composition. CONCLUSION: Our findings emphasise the importance of studying microbiome dynamics to capture changes occurring during ontogeny. They highlight the role of microbial communities in reflecting host health and the importance of the nest environment for the developing nestling microbiome. Overall, this study contributes to understanding the complex interplay between microbial communities, host factors, and environmental variables, and sheds light on the ecological processes governing gut microbial colonisation during early-life stages.

No evidence for a role of MHC class II genotype in the chemical encoding of heterozygosity and relatedness in Antarctic fur seals.

Despite decades of research, surprisingly little is known about the mechanism(s) by which an individual's genotype is encoded in odour. Many studies have focused on the role of the major histocompatibility complex (MHC) owing to its importance for survival and mate choice. However, the salience of MHC-mediated odours compared to chemicals influenced by the rest of the genome remains unclear, especially in wild populations where it is challenging to quantify and control for the effects of the genomic background. We addressed this issue in Antarctic fur seals by analysing skin swabs together with full-length MHC DQB II exon 2 sequences and data from 41 genome-wide distributed microsatellites. We did not find any effects of MHC relatedness on chemical similarity and there was also no relationship between MHC heterozygosity and chemical diversity. However, multilocus heterozygosity showed a significant positive association with chemical diversity, even after controlling for MHC heterozygosity. Our results appear to rule out a dominant role of the MHC in the chemical encoding of genetic information in a wild vertebrate population and highlight the need for genome-wide approaches to elucidate the mechanism(s) and specific genes underlying genotype-odour associations.

Ninety years of change, from commercial extinction to recovery, range expansion and decline for Antarctic fur seals at South Georgia.

With environmental change, understanding how species recover from overharvesting and maintain viable populations is central to ecosystem restoration. Here, we reconstruct 90 years of recovery trajectory of the Antarctic fur seal at South Georgia (S.W. Atlantic), a key indicator species in the krill-based food webs of the Southern Ocean. After being harvested to commercial extinction by 1907, this population rebounded and now constitutes the most abundant otariid in the World. However, its status remains uncertain due to insufficient and conflicting data, and anthropogenic pressures affecting Antarctic krill, an essential staple for millions of fur seals and other predators. Using integrated population models, we estimated simultaneously the long-term abundance for Bird Island, northwest South Georgia, epicentre of recovery of the species after sealing, and population adjustments for survey counts with spatiotemporal applicability. Applied to the latest comprehensive survey data, we estimated the population at South Georgia in 2007-2009 as 3,510,283 fur seals [95% CI: 3,140,548-3,919,604] (ca. 98% of global population), after 40 years of maximum growth and range expansion owing to an abundant krill supply. At Bird Island, after 50 years of exponential growth followed by 25 years of slow stable growth, the population collapsed in 2009 and has thereafter declined by -7.2% [-5.2, -9.1] per annum, to levels of the 1970s. For the instrumental record, this trajectory correlates with a time-varying relationship between coupled climate and sea surface temperature cycles associated with low regional krill availability, although the effects of increasing krill extraction by commercial fishing and natural competitors remain uncertain. Since 2015, fur seal longevity and recruitment have dropped, sexual maturation has retarded, and population growth is expected to remain mostly negative and highly variable. Our analysis documents the rise and fall of a key Southern Ocean predator over a century of profound environmental and ecosystem change.

Genetic diversity and demographic history of the leopard seal: A Southern Ocean top predator.

Leopard seals (Hydrurga leptonyx) are top predators that can exert substantial top-down control of their Antarctic prey species. However, population trends and genetic diversity of leopard seals remain understudied, limiting our understanding of their ecological role. We investigated the genetic diversity, effective population size and demographic history of leopard seals to provide fundamental data that contextualizes their predatory influence on Antarctic ecosystems. Ninety leopard seals were sampled from the northern Antarctic Peninsula during the austral summers of 2008-2019 and a 405bp segment of the mitochondrial control region was sequenced for each individual. We uncovered moderate levels of nucleotide (π = 0.013) and haplotype (Hd = 0.96) diversity, and the effective population size was estimated at around 24,000 individuals (NE = 24,376; 95% CI: 16,876-33,126). Consistent with findings from other ice-breeding pinnipeds, Bayesian skyline analysis also revealed evidence for population expansion during the last glacial maximum, suggesting that historical population growth may have been boosted by an increase in the abundance of sea ice. Although leopard seals can be found in warmer, sub-Antarctic locations, the species' core habitat is centered on the Antarctic, making it inherently vulnerable to the loss of sea ice habitat due to climate change. Therefore, detailed assessments of past and present leopard seal population trends are needed to inform policies for Antarctic ecosystems.

Effects of hunting on genetic diversity, inbreeding and dispersal in Finnish black grouse (Lyrurus tetrix).

Intensive hunting activities such as commercial fishing and trophy hunting can have profound influences on natural populations. However, less intensive recreational hunting can also have subtle effects on animal behaviour, habitat use and movement, with implications for population persistence. Lekking species such as the black grouse (Lyrurus tetrix) may be especially prone to hunting as leks are temporally and spatially predictable, making them easy targets. Furthermore, inbreeding in black grouse is mainly avoided through female-biased dispersal, so any disruptions to dispersal caused by hunting could lead to changes in gene flow, increasing the risk of inbreeding. We therefore investigated the impact of hunting on genetic diversity, inbreeding and dispersal on a metapopulation of black grouse in Central Finland. We genotyped 1065 adult males and 813 adult females from twelve lekking sites (six hunted, six unhunted) and 200 unrelated chicks from seven sites (two hunted, five unhunted) at up to thirteen microsatellite loci. Our initial confirmatory analysis of sex-specific fine-scale population structure revealed little genetic structure in the metapopulation. Levels of inbreeding did not differ significantly between hunted and unhunted sites in neither adults nor chicks. However, immigration rates into hunted sites were significantly higher among adults compared to immigration into unhunted sites. We conclude that the influx of migrants into hunted sites may compensate for the loss of harvested individuals, thereby increasing gene flow and mitigating inbreeding. Given the absence of any obvious barriers to gene flow in Central Finland, a spatially heterogeneous matrix of hunted and unhunted regions may be crucial to ensure sustainable harvests into the future.

Heavy metal contamination in pristine environments: lessons from the Juan Fernandez fur seal (Arctocephalus philippii philippii).

Heavy metals, including mercury (Hg) and cadmium (Cd), occur naturally or anthropogenically and are considered toxic to the environment and human health. However, studies on heavy metal contamination focus on locations close to industrialized settlements, while isolated environments with little human activity are often ignored due to perceived low risk. This study reports heavy metal exposure in Juan Fernandez fur seals (JFFS), a marine mammal endemic to an isolated and relatively pristine archipelago off the coast of Chile. We found exceptionally high concentrations of Cd and Hg in JFFS faeces. Indeed, they are among the highest reported for any mammalian species. Following analysis of their prey, we concluded that diet is the most likely source of Cd contamination in JFFS. Furthermore, Cd appears to be absorbed and incorporated into JFFS bones. However, it was not associated with mineral changes observed in other species, suggesting Cd tolerance/adaptations in JFFS bones. The high levels of silicon found in JFFS bones may counteract the effects of Cd. These findings are relevant to biomedical research, food security and the treatment of heavy metal contamination. It also contributes to understanding the ecological role of JFFS and highlights the need for surveillance of apparently pristine environments.

Evolution of the germline mutation rate across vertebrates.

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself1. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies2. Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis3. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals.

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

Signatures of Selection on Mitonuclear Integrated Genes Uncover Hidden Mitogenomic Variation in Fur Seals.

Nuclear copies of mitochondrial genes (numts) are commonplace in vertebrate genomes and have been characterized in many species. However, relatively little attention has been paid to understanding their evolutionary origins and to disentangling alternative sources of insertions. Numts containing genes with intact mitochondrial reading frames represent good candidates for this purpose. The sequences of the genes they contain can be compared with their mitochondrial homologs to characterize synonymous to nonsynonymous substitution rates, which can shed light on the selection pressures these genes have been subjected to. Here, we characterize 25 numts in the Antarctic fur seal (Arctocephalus gazella) genome. Among those containing genes with intact mitochondrial reading frames, three carry multiple substitutions in comparison to their mitochondrial homologs. Our analyses reveal that one represents a historic insertion subjected to strong purifying selection since it colonized the Otarioidea in a genomic region enriched in retrotransposons. By contrast, the other two numts appear to be more recent and their large number of substitutions can be attributed to noncanonical insertions, either the integration of heteroplasmic mtDNA or hybridization. Our study sheds new light on the evolutionary history of pinniped numts and uncovers the presence of hidden sources of mitonuclear variation.

Low heritability and high phenotypic plasticity of salivary cortisol in response to environmental heterogeneity in a wild pinniped.

Individuals are unique in how they interact with and respond to their environment. Correspondingly, unpredictable challenges or environmental stressors often produce an individualized response of the hypothalamic-pituitary-adrenal (HPA) axis and its downstream effector cortisol. We used a fully crossed, repeated measures design to investigate the factors shaping individual variation in baseline cortisol in Antarctic fur seal pups and their mothers. Saliva samples were collected from focal individuals at two breeding colonies, one with low and the other with high density, during two consecutive years of contrasting food availability. Mothers and pups were sampled concurrently at birth and shortly before weaning, while pups were additionally sampled every 20 days. We found that heritability was low for baseline cortisol, while within-individual repeatability and among-individual variability were high. A substantial proportion of the variation in baseline cortisol could be explained in pups and mothers by a combination of intrinsic and extrinsic factors including sex, weight, day, season, and colony of birth. Our findings provide detailed insights into the individualization of endocrine phenotypes and their genetic and environmental drivers in a wild pinniped. Furthermore, the strong associations between cortisol and life history traits that we report in fur seals could have important implications for understanding the population dynamics of species impacted by environmental change.

Demographic Reconstruction of Antarctic Fur Seals Supports the Krill Surplus Hypothesis.

Much debate surrounds the importance of top-down and bottom-up effects in the Southern Ocean, where the harvesting of over two million whales in the mid twentieth century is thought to have produced a massive surplus of Antarctic krill. This excess of krill may have allowed populations of other predators, such as seals and penguins, to increase, a top-down hypothesis known as the 'krill surplus hypothesis'. However, a lack of pre-whaling population baselines has made it challenging to investigate historical changes in the abundance of the major krill predators in relation to whaling. Therefore, we used reduced representation sequencing and a coalescent-based maximum composite likelihood approach to reconstruct the recent demographic history of the Antarctic fur seal, a pinniped that was hunted to the brink of extinction by 18th and 19th century sealers. In line with the known history of this species, we found support for a demographic model that included a substantial reduction in population size around the time period of sealing. Furthermore, maximum likelihood estimates from this model suggest that the recovered, post-sealing population at South Georgia may have been around two times larger than the pre-sealing population. Our findings lend support to the krill surplus hypothesis and illustrate the potential of genomic approaches to shed light on long-standing questions in population biology.

A stable foraging polymorphism buffers Galápagos sea lions against environmental change.

Understanding the ability of animals to cope with a changing environment is critical in a world affected by anthropogenic disturbance.1 Individual foraging strategies may influence the coping ability of entire populations, as these strategies can be adapted to contrasting conditions, allowing populations with foraging polymorphisms to be more resilient toward environmental change.2,3 However, environmentally dependent fitness consequences of individual foraging strategies and their effects on population dynamics have not been conclusively documented.4,5 Here, we use biologging data from endangered Galápagos sea lion females (Zalophus wollebaeki) to show that benthically foraging individuals dig after sand-dwelling prey species while pelagic foragers hunt in more open waters. These specialized foraging behaviors result in distinct and temporally stable patterns of vibrissae abrasion. Using vibrissae length as a visual marker for the benthic versus pelagic foraging strategies, we furthermore uncovered an environment-dependent fitness trade-off between benthic and pelagic foragers, suggesting that the foraging polymorphism could help to buffer the population against the negative effects of climate change. However, demographic projections suggest that this buffering effect is unlikely to be sufficient to reverse the ongoing population decline of the past four decades.6 Our study shows how crucial a deeper understanding of behavioral polymorphisms can be for predicting how populations cope within a rapidly changing world. VIDEO ABSTRACT.

Origin and expansion of the world's most widespread pinniped: Range-wide population genomics of the harbour seal (Phoca vitulina).

The harbour seal (Phoca vitulina) is the most widely distributed pinniped, occupying a wide variety of habitats and climatic zones across the Northern Hemisphere. Intriguingly, the harbour seal is also one of the most philopatric seals, raising questions as to how it colonized its current range. To shed light on the origin, remarkable range expansion, population structure and genetic diversity of this species, we used genotyping-by-sequencing to analyse ~13,500 biallelic single nucleotide polymorphisms from 286 individuals sampled from 22 localities across the species' range. Our results point to a Northeast Pacific origin of the harbour seal, colonization of the North Atlantic via the Canadian Arctic, and subsequent stepping-stone range expansions across the North Atlantic from North America to Europe, accompanied by a successive loss of genetic diversity. Our analyses further revealed a deep divergence between modern North Pacific and North Atlantic harbour seals, with finer-scale genetic structure at regional and local scales consistent with strong philopatry. The study provides new insights into the harbour seal's remarkable ability to colonize and adapt to a wide range of habitats. Furthermore, it has implications for current harbour seal subspecies delineations and highlights the need for international and national red lists and management plans to ensure the protection of genetically and demographically isolated populations.

Heterogeneous Genomic Divergence Landscape in Two Commercially Important European Scallop Species.

Two commercially important scallop species of the genus Pecten are found in Europe: the north Atlantic Pecten maximus and the Mediterranean Pecten jacobaeus whose distributions abut at the Almeria-Orán front. Whilst previous studies have quantified genetic divergence between these species, the pattern of differentiation along the Pecten genome is unknown. Here, we mapped RADseq data from 235 P. maximus and 27 P. jacobaeus to a chromosome-level reference genome, finding a heterogeneous landscape of genomic differentiation. Highly divergent genomic regions were identified across 14 chromosomes, while the remaining five showed little differentiation. Demographic and comparative genomics analyses suggest that this pattern resulted from an initial extended period of isolation, which promoted divergence, followed by differential gene flow across the genome during secondary contact. Single nucleotide polymorphisms present within highly divergent genomic regions were located in areas of low recombination and contrasting patterns of LD decay were found between the two species, hinting at the presence of chromosomal inversions in P. jacobaeus. Functional annotations revealed that highly differentiated regions were enriched for immune-related processes and mRNA modification. While future work is necessary to characterize structural differences, this study provides new insights into the speciation genomics of P. maximus and P. jacobaeus.