Diving into broad-scale and high-resolution population genomics to decipher drivers of structure and climatic vulnerability in a marine invertebrate.

Species with widespread distributions play a crucial role in our understanding of climate change impacts on population structure. In marine species, population structure is often governed by both high connectivity potential and selection across strong environmental gradients. Despite the complexity of factors influencing marine populations, studying species with broad distribution can provide valuable insights into the relative importance of these factors and the consequences of climate-induced alterations across environmental gradients. We used the northern shrimp Pandalus borealis and its wide latitudinal distribution to identify current drivers of population structure and predict the species' vulnerability to climate change. A total of 1514 individuals sampled across 24° latitude were genotyped at high geographic (54 stations) and genetic (14,331 SNPs) resolutions to assess genetic variation and environmental correlations. Four populations were identified in addition to finer substructure associated with local adaptation. Geographic patterns of neutral population structure reflected predominant oceanographic currents, while a significant proportion of the genetic variation was associated with gradients in salinity and temperature. Adaptive landscapes generated using climate projections suggest a larger genomic offset in the southern extent of the P. borealis range, where shrimp had the largest adaptive standing genetic variation. Our genomic results combined with recent observations point to further deterioration in southern regions and an impending vulnerable status in the regions at higher latitudes for P. borealis. They also provide rare insights into the drivers of population structure and climatic vulnerability of a widespread meroplanktonic species, which is crucial to understanding future challenges associated with invertebrates essential to ecosystem functioning.

Variation in genomic vulnerability to climate change across temperate populations of eelgrass (Zostera marina).

A global decline in seagrass populations has led to renewed calls for their conservation as important providers of biogenic and foraging habitat, shoreline stabilization and carbon storage. Eelgrass (Zostera marina) occupies the largest geographic range among seagrass species spanning a commensurately broad spectrum of environmental conditions. In Canada, eelgrass is managed as a single phylogroup despite occurring across three oceans and a range of ocean temperatures and salinity gradients. Previous research has focused on applying relatively few markers to reveal population structure of eelgrass, whereas a whole-genome approach is warranted to investigate cryptic structure among populations inhabiting different ocean basins and localized environmental conditions. We used a pooled whole-genome re-sequencing approach to characterize population structure, gene flow and environmental associations of 23 eelgrass populations ranging from the Northeast United States to Atlantic, subarctic and Pacific Canada. We identified over 500,000 SNPs, which when mapped to a chromosome-level genome assembly revealed six broad clades of eelgrass across the study area, with pairwise F ST ranging from 0 among neighbouring populations to 0.54 between Pacific and Atlantic coasts. Genetic diversity was highest in the Pacific and lowest in the subarctic, consistent with colonization of the Arctic and Atlantic oceans from the Pacific less than 300 kya. Using redundancy analyses and two climate change projection scenarios, we found that subarctic populations are predicted to be potentially more vulnerable to climate change through genomic offset predictions. Conservation planning in Canada should thus ensure that representative populations from each identified clade are included within a national network so that latent genetic diversity is protected, and gene flow is maintained. Northern populations, in particular, may require additional mitigation measures given their potential susceptibility to a rapidly changing climate.

Application of Omics Tools in Designing and Monitoring Marine Protected Areas For a Sustainable Blue Economy.

A key component of the global blue economy strategy is the sustainable extraction of marine resources and conservation of marine environments through networks of marine protected areas (MPAs). Connectivity and representativity are essential factors that underlie successful implementation of MPA networks, which can safeguard biological diversity and ecosystem function, and ultimately support the blue economy strategy by balancing ocean use with conservation. New "big data" omics approaches, including genomics and transcriptomics, are becoming essential tools for the development and maintenance of MPA networks. Current molecular omics techniques, including population-scale genome sequencing, have direct applications for assessing population connectivity and for evaluating how genetic variation is represented within and among MPAs. Effective baseline characterization and long-term, scalable, and comprehensive monitoring are essential for successful MPA management, and omics approaches hold great promise to characterize the full range of marine life, spanning the microbiome to megafauna across a range of environmental conditions (shallow sea to the deep ocean). Omics tools, such as eDNA metabarcoding can provide a cost-effective basis for biodiversity monitoring in large and remote conservation areas. Here we provide an overview of current omics applications for conservation planning and monitoring, with a focus on metabarcoding, metagenomics, and population genomics. Emerging approaches, including whole-genome sequencing, characterization of genomic architecture, epigenomics, and genomic vulnerability to climate change are also reviewed. We demonstrate that the operationalization of omics tools can enhance the design, monitoring, and management of MPAs and thus will play an important role in a modern and comprehensive blue economy strategy.

Modeling demersal fish and benthic invertebrate assemblages in support of marine conservation planning.

Marine classification schemes based on abiotic surrogates often inform regional marine conservation planning in lieu of detailed biological data. However, these schemes may poorly represent ecologically relevant biological patterns required for effective design and management strategies. We used a community-level modeling approach to characterize and delineate representative mesoscale (tens to thousands of kilometers) assemblages of demersal fish and benthic invertebrates in the Northwest Atlantic. Hierarchical clustering of species occurrence data from four regional annual multispecies trawl surveys revealed three to six groupings (predominant assemblage types) in each survey region, broadly associated with geomorphic and oceanographic features. Indicator analyses identified 3-34 emblematic taxa of each assemblage type. Random forest classifications accurately predicted assemblage distributions from environmental covariates (AUC > 0.95) and identified thermal limits (annual minimum and maximum bottom temperatures) as important predictors of distribution in each region. Using forecasted oceanographic conditions for the year 2075 and a regional classification model, we projected assemblage distributions in the southernmost bioregion (Scotian Shelf-Bay of Fundy) under a high emissions climate scenario (RCP 8.5). Range expansions to the northeast are projected for assemblages associated with warmer and shallower waters of the Western Scotian Shelf over the 21st century as thermal habitat on the relatively cooler Eastern Scotian Shelf becomes more favorable. Community-level modeling provides a biotic-informed approach for identifying broadscale ecological structure required for the design and management of ecologically coherent, representative, well-connected networks of Marine Protected Areas. When combined with oceanographic forecasts, this modeling approach provides a spatial tool for assessing sensitivity and resilience to climate change, which can improve conservation planning, monitoring, and adaptive management.

Eusociality in snapping shrimps is associated with larger genomes and an accumulation of transposable elements.

Despite progress uncovering the genomic underpinnings of sociality, much less is known about how social living affects the genome. In different insect lineages, for example, eusocial species show both positive and negative associations between genome size and structure, highlighting the dynamic nature of the genome. Here, we explore the relationship between sociality and genome architecture in Synalpheus snapping shrimps that exhibit multiple origins of eusociality and extreme interspecific variation in genome size. Our goal is to determine whether eusociality leads to an accumulation of repetitive elements and an increase in genome size, presumably due to reduced effective population sizes resulting from a reproductive division of labor, or whether an initial accumulation of repetitive elements leads to larger genomes and independently promotes the evolution of eusociality through adaptive evolution. Using phylogenetically informed analyses, we find that eusocial species have larger genomes with more transposable elements (TEs) and microsatellite repeats than noneusocial species. Interestingly, different TE subclasses contribute to the accumulation in different species. Phylogenetic path analysis testing alternative causal relationships between sociality and genome architecture is most consistent with the hypothesis that TEs modulate the relationship between sociality and genome architecture. Although eusociality appears to influence TE accumulation, ancestral state reconstruction suggests moderate TE abundances in ancestral species could have fueled the initial transitions to eusociality. Ultimately, we highlight a complex and dynamic relationship between genome and social evolution, demonstrating that sociality can influence the evolution of the genome, likely through changes in demography related to patterns of reproductive skew.

Resolving fine-scale population structure and fishery exploitation using sequenced microsatellites in a northern fish.

The resiliency of populations and species to environmental change is dependent on the maintenance of genetic diversity, and as such, quantifying diversity is central to combating ongoing widespread reductions in biodiversity. With the advent of next-generation sequencing, several methods now exist for resolving fine-scale population structure, but the comparative performance of these methods for genetic assignment has rarely been tested. Here, we evaluate the performance of sequenced microsatellites and a single nucleotide polymorphism (SNP) array to resolve fine-scale population structure in a critically important salmonid in north eastern Canada, Arctic Charr (Salvelinus alpinus). We also assess the utility of sequenced microsatellites for fisheries applications by quantifying the spatial scales of movement and exploitation through genetic assignment of fishery samples to rivers of origin and comparing these results with a 29-year tagging dataset. Self-assignment and simulation-based analyses of 111 genome-wide microsatellite loci and 500 informative SNPs from 28 populations of Arctic Charr in north-eastern Canada identified largely river-specific genetic structure. Despite large differences (~4X) in the number of loci surveyed between panels, mean self-assignment accuracy was similar with the microsatellite loci and the SNP panel (>90%). Subsequent analysis of 996 fishery-collected samples using the microsatellite panel revealed that larger rivers contribute greater numbers of individuals to the fishery and that coastal fisheries largely exploit individuals originating from nearby rivers, corroborating results from traditional tagging experiments. Our results demonstrate the efficacy of sequence-based microsatellite genotyping to advance understanding of fine-scale population structure and harvest composition in northern and understudied species.

Integrating climate adaptation and biodiversity conservation in the global ocean.

The impacts of climate change and the socioecological challenges they present are ubiquitous and increasingly severe. Practical efforts to operationalize climate-responsive design and management in the global network of marine protected areas (MPAs) are required to ensure long-term effectiveness for safeguarding marine biodiversity and ecosystem services. Here, we review progress in integrating climate change adaptation into MPA design and management and provide eight recommendations to expedite this process. Climate-smart management objectives should become the default for all protected areas, and made into an explicit international policy target. Furthermore, incentives to use more dynamic management tools would increase the climate change responsiveness of the MPA network as a whole. Given ongoing negotiations on international conservation targets, now is the ideal time to proactively reform management of the global seascape for the dynamic climate-biodiversity reality.

Global marine biodiversity in the context of achieving the Aichi Targets: ways forward and addressing data gaps.

In 2010, the Conference of the Parties of the Convention on Biological Diversity agreed on the Strategic Plan for Biodiversity 2011-2020 in Aichi Prefecture, Japan. As this plan approaches its end, we discussed whether marine biodiversity and prediction studies were nearing the Aichi Targets during the 4th World Conference on Marine Biodiversity held in Montreal, Canada in June 2018. This article summarises the outcome of a five-day group discussion on how global marine biodiversity studies should be focused further to better understand the patterns of biodiversity. We discussed and reviewed seven fundamental biodiversity priorities related to nine Aichi Targets focusing on global biodiversity discovery and predictions to improve and enhance biodiversity data standards (quantity and quality), tools and techniques, spatial and temporal scale framing, and stewardship and dissemination. We discuss how identifying biodiversity knowledge gaps and promoting efforts have and will reduce such gaps, including via the use of new databases, tools and technology, and how these resources could be improved in the future. The group recognised significant progress toward Target 19 in relation to scientific knowledge, but negligible progress with regard to Targets 6 to 13 which aimed to safeguard and reduce human impacts on biodiversity.

Fine-scale temperature-associated genetic structure between inshore and offshore populations of sea scallop (Placopecten magellanicus).

In the northwest Atlantic Ocean, sea scallop (Placopecten magellanicus) has been characterized by a latitudinal genetic cline with a breakpoint between northern and southern genetic clusters occurring at ~45°N along eastern Nova Scotia, Canada. Using 96 diagnostic single-nucleotide polymorphisms (SNPs) capable of discriminating between northern and southern clusters, we examined fine-scale genetic structure of scallops among 27 sample locations, spanning the largest geographic range evaluated in this species to date (~37-51°N). Here, we confirmed previous observations of northern and southern groups, but we show that the boundary between northern and southern clusters is not a discrete latitudinal break. Instead, at latitudes near the previously described boundary, we found unexpected patterns of fine-scale genetic structure occurring between inshore and offshore sites. Scallops from offshore sites, including St. Pierre Bank and the eastern Scotian Shelf, clustered with southern stocks, whereas inshore sites at similar latitudes clustered with northern stocks. Our analyses revealed significant genetic divergence across small spatial scales (i.e., 129-221 km distances), and that spatial structure over large and fine scales was strongly associated with temperature during seasonal periods of thermal minima. Clear temperature differences between inshore and offshore locations may explain the fine-scale structuring observed, such as why southern lineages of scallop occur at higher latitudes in deeper, warmer offshore waters. Our study supports growing evidence that fine-scale population structure in marine species is common, often environmentally associated, and that consideration of environmental and genomic data can significantly enhance the identification of marine diversity and management units.

Temporal dynamics of genetic clines of invasive European green crab (Carcinus maenas) in eastern North America.

Two genetically distinct lineages of European green crabs (Carcinus maenas) were independently introduced to eastern North America, the first in the early 19th century and the second in the late 20th century. These lineages first came into secondary contact in southeastern Nova Scotia, Canada (NS), where they hybridized, producing latitudinal genetic clines. Previous studies have documented a persistent southward shift in the clines of different marker types, consistent with existing dispersal and recruitment pathways. We evaluated current clinal structure by quantifying the distribution of lineages and fine-scale hybridization patterns across the eastern North American range (25 locations, ~39 to 49°N) using informative single nucleotide polymorphisms (SNPs; n = 96). In addition, temporal changes in the genetic clines were evaluated using mitochondrial DNA and microsatellite loci (n = 9-11) over a 15-year period (2000-2015). Clinal structure was consistent with prior work demonstrating the existence of both northern and southern lineages with a hybrid zone occurring between southern New Brunswick (NB) and southern NS. Extensive later generation hybrids were detected in this region and in southeastern Newfoundland. Temporal genetic analysis confirmed the southward progression of clines over time; however, the rate of this progression was slower than predicted by forecasting models, and current clines for all marker types deviated significantly from these predictions. Our results suggest that neutral and selective processes contribute to cline dynamics, and ultimately, highlight how selection, hybridization, and dispersal can collectively influence invasion success.

Extensive hybridization following a large escape of domesticated Atlantic salmon in the Northwest Atlantic.

Domestication is rife with episodes of interbreeding between cultured and wild populations, potentially challenging adaptive variation in the wild. In Atlantic salmon, Salmo salar, the number of domesticated individuals far exceeds wild individuals, and escape events occur regularly, yet evidence of the magnitude and geographic scale of interbreeding resulting from individual escape events is lacking. We screened juvenile Atlantic salmon using 95 single nucleotide polymorphisms following a single, large aquaculture escape in the Northwest Atlantic and report the landscape-scale detection of hybrid and feral salmon (27.1%, 17/18 rivers). Hybrids were reproductively viable, and observed at higher frequency in smaller wild populations. Repeated annual sampling of this cohort revealed decreases in the presence of hybrid and feral offspring over time. These results link previous observations of escaped salmon in rivers with reports of population genetic change, and demonstrate the potential negative consequences of escapes from net-pen aquaculture on wild populations.

Genomewide evidence of environmentally mediated secondary contact of European green crab (Carcinus maenas) lineages in eastern North America.

Genetic-environment associations are increasingly revealed through population genomic data and can occur through a number of processes, including secondary contact, divergent natural selection, or isolation by distance. Here, we investigate the influence of the environment, including seasonal temperature and salinity, on the population structure of the invasive European green crab (Carcinus maenas) in eastern North America. Green crab populations in eastern North America are associated with two independent invasions, previously shown to consist of distinct northern and southern ecotypes, with a contact zone in southern Nova Scotia, Canada. Using a RAD-seq panel of 9,137 genomewide SNPs, we detected 41 SNPs (0.49%) whose allele frequencies were highly correlated with environmental data. A principal components analysis of 25 environmental variables differentiated populations into northern, southern, and admixed sites in concordance with the observed genomic spatial structure. Furthermore, a spatial principal components analysis conducted on genomic and geographic data revealed a high degree of global structure (p < .0001) partitioning a northern and southern ecotype. Redundancy and partial redundancy analyses revealed that among the environmental variables tested, winter sea surface temperature had the strongest association with spatial structuring, suggesting that it is an important factor defining range and expansion limits of each ecotype. Understanding environmental thresholds associated with intraspecific diversity will facilitate the ability to manage current and predict future distributions of this aquatic invasive species.

Genotyping-by-sequencing of genome-wide microsatellite loci reveals fine-scale harvest composition in a coastal Atlantic salmon fishery.

Individual assignment and genetic mixture analysis are commonly utilized in contemporary wildlife and fisheries management. Although microsatellite loci provide unparalleled numbers of alleles per locus, their use in assignment applications is increasingly limited. However, next-generation sequencing, in conjunction with novel bioinformatic tools, allows large numbers of microsatellite loci to be simultaneously genotyped, presenting new opportunities for individual assignment and genetic mixture analysis. Here, we scanned the published Atlantic salmon genome to identify 706 microsatellite loci, from which we developed a final panel of 101 microsatellites distributed across the genome (average 3.4 loci per chromosome). Using samples from 35 Atlantic salmon populations (n = 1,485 individuals) from coastal Labrador, Canada, a region characterized by low levels of differentiation in this species, this panel identified 844 alleles (average of 8.4 alleles per locus). Simulation-based evaluations of assignment and mixture identification accuracy revealed unprecedented resolution, clearly identifying 26 rivers or groups of rivers spanning 500 km of coastline. This baseline was used to examine the stock composition of 696 individuals harvested in the Labrador Atlantic salmon fishery and revealed that coastal fisheries largely targeted regional groups (<300 km). This work suggests that the development and application of large sequenced microsatellite panels presents great potential for stock resolution in Atlantic salmon and more broadly in other exploited anadromous and marine species.

hybriddetective: A workflow and package to facilitate the detection of hybridization using genomic data in r.

The ability to detect and characterize hybridization in nature has long been of interest to many fields of biology and often has direct implications for wildlife management and conservation. The capacity to identify the presence of hybridization, and quantify the numbers of individuals belonging to different hybrid classes, permits inference on the magnitude of, and timescale over which, hybridization has been or is occurring. Here, we present an r package and associated workflow developed for the detection, with estimates of efficiency and accuracy, of multigenerational hybrid individuals using genetic or genomic data in conjunction with the program newhybrids. This package includes functions for the identification and testing of diagnostic panels of markers, the simulation of multigenerational hybrids, and the quantification and visualization of the efficiency and accuracy with which hybrids can be detected. Overall, this package delivers a streamlined hybrid analysis platform, providing improvements in speed, ease of use and repeatability over current ad hoc approaches. The latest version of the package and associated documentation are available on GitHub (https://github.com/bwringe/hybriddetective).

The Genome Sizes of Ostracod Crustaceans Correlate with Body Size and Evolutionary History, but not Environment.

Within animals, a positive correlation between genome size and body size has been detected in several taxa but not in others, such that it remains unknown how pervasive this pattern may be. Here, we provide another example of a positive relationship in a group of crustaceans whose genome sizes have not previously been investigated. We analyze genome size estimates for 46 species across the 2 most diverse orders of Class Ostracoda, commonly known as seed shrimps, including 29 new estimates made using Feulgen image analysis densitometry and flow cytometry. Genome sizes in this group range ~80-fold, a level of variability that is otherwise not seen in crustaceans with the exception of some malacostracan orders. We find a strong positive correlation between genome size and body size across all species, including after phylogenetic correction. We additionally detect evidence of XX/XO sex determination in 3 species of marine ostracods where male and female genome sizes were estimated. On average, genome sizes are larger but less variable in Order Myodocopida than in Order Podocopida, and marine ostracods have larger genomes than freshwater species, but this appears to be explained by phylogenetic inertia. The relationship between phylogeny, genome size, body size, and habitat is complex in this system and provides a baseline for future studies examining the interactions of these biological traits.

RAD sequencing reveals genomewide divergence between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic.

Genomic studies of invasive species can reveal both invasive pathways and functional differences underpinning patterns of colonization success. The European green crab (Carcinus maenas) was initially introduced to eastern North America nearly 200 years ago where it expanded northwards to eastern Nova Scotia. A subsequent invasion to Nova Scotia from a northern European source allowed further range expansion, providing a unique opportunity to study the invasion genomics of a species with multiple invasions. Here, we use restriction-site-associated DNA sequencing-derived SNPs to explore fine-scale genomewide differentiation between these two invasions. We identified 9137 loci from green crab sampled from 11 locations along eastern North America and compared spatial variation to mitochondrial COI sequence variation used previously to characterize these invasions. Overall spatial divergence among invasions was high (pairwise FST ~0.001 to 0.15) and spread across many loci, with a mean FST ~0.052 and 52% of loci examined characterized by FST values >0.05. The majority of the most divergent loci (i.e., outliers, ~1.2%) displayed latitudinal clines in allele frequency highlighting extensive genomic divergence among the invasions. Discriminant analysis of principal components (both neutral and outlier loci) clearly resolved the two invasions spatially and was highly correlated with mitochondrial divergence. Our results reveal extensive cryptic intraspecific genomic diversity associated with differing patterns of colonization success and demonstrates clear utility for genomic approaches to delineating the distribution and colonization success of aquatic invasive species.

Nuclear DNA content correlates with depth, body size, and diversification rate in amphipod crustaceans from ancient Lake Baikal, Russia.

Lake Baikal in Russia is a large, ancient lake that has been the site of a major radiation of amphipod crustaceans. Nearly 400 named species are known in this single lake, and it is thought that many more await description. The size and depth of Lake Baikal, in particular, may have contributed to the radiation of endemic amphipods by providing a large number of microhabitats for species to invade and subsequently experience reproductive isolation. Here we investigate the possibility that large-scale genomic changes have also accompanied diversification in these crustaceans. Specifically, we report genome size estimates for 36 species of Baikal amphipods, and examine the relationship between genome size, body size, and the maximum depths at which the amphipods are found in the lake. Genome sizes ranged nearly 8-fold in this sample of amphipod species, from 2.15 to 16.63 pg, and there were significant, positive, phylogenetically corrected relationships between genome size, body size, maximum depth, and diversification rate among these species. Our results suggest that major genomic changes, including transposable element proliferation, have accompanied speciation that was driven by selection for differences in body size and habitat preference in Lake Baikal amphipods.

Range-wide parallel climate-associated genomic clines in Atlantic salmon.

Clinal variation across replicated environmental gradients can reveal evidence of local adaptation, providing insight into the demographic and evolutionary processes that shape intraspecific diversity. Using 1773 genome-wide single nucleotide polymorphisms we evaluated latitudinal variation in allele frequency for 134 populations of North American and European Atlantic salmon (Salmo salar). We detected 84 (4.74%) and 195 (11%) loci showing clinal patterns in North America and Europe, respectively, with 12 clinal loci in common between continents. Clinal single nucleotide polymorphisms were evenly distributed across the salmon genome and logistic regression revealed significant associations with latitude and seasonal temperatures, particularly average spring temperature in both continents. Loci displaying parallel clines were associated with several metabolic and immune functions, suggesting a potential basis for climate-associated adaptive differentiation. These climate-based clines collectively suggest evidence of large-scale environmental associated differences on either side of the North Atlantic. Our results support patterns of parallel evolution on both sides of the North Atlantic, with evidence of both similar and divergent underlying genetic architecture. The identification of climate-associated genomic clines illuminates the role of selection and demographic processes on intraspecific diversity in this species and provides a context in which to evaluate the impacts of climate change.

Identifying patterns of dispersal, connectivity and selection in the sea scallop, Placopecten magellanicus, using RADseq-derived SNPs.

Understanding patterns of dispersal and connectivity among marine populations can directly inform fisheries conservation and management. Advances in high-throughput sequencing offer new opportunities for estimating marine connectivity. We used restriction-site-associated DNA sequencing to examine dispersal and realized connectivity in the sea scallop Placopecten magellanicus, an economically important marine bivalve. Based on 245 individuals sampled rangewide at 12 locations from Newfoundland to the Mid-Atlantic Bight, we identified and genotyped 7163 single nucleotide polymorphisms; 112 (1.6%) were identified as outliers potentially under directional selection. Bayesian clustering revealed a discontinuity between northern and southern samples, and latitudinal clines in allele frequencies were observed in 42.9% of the outlier loci and in 24.6% of neutral loci. Dispersal estimates derived using these clines and estimates of linkage disequilibrium imply limited dispersal; 373.1 ± 407.0 km (mean ± SD) for outlier loci and 641.0 ± 544.6 km (mean ± SD) for neutral loci. Our analysis suggests restricted dispersal compared to the species range (>2000 km) and that dispersal and effective connectivity differ. These observations support the hypothesis that limited effective dispersal structures scallop populations along eastern North America. These findings can help refine the appropriate scale of management and conservation in this commercially valuable species.

parallelnewhybrid: an R package for the parallelization of hybrid detection using newhybrids.

Hybridization among populations and species is a central theme in many areas of biology, and the study of hybridization has direct applicability to testing hypotheses about evolution, speciation and genetic recombination, as well as having conservation, legal and regulatory implications. Yet, despite being a topic of considerable interest, the identification of hybrid individuals, and quantification of the (un)certainty surrounding the identifications, remains difficult. Unlike other programs that exist to identify hybrids based on genotypic information, newhybrids is able to assign individuals to specific hybrid classes (e.g. F1 , F2 ) because it makes use of patterns of gene inheritance within each locus, rather than just the proportions of gene inheritance within each individual. For each comparison and set of markers, multiple independent runs of each data set should be used to develop an estimate of the hybrid class assignment accuracy. The necessity of analysing multiple simulated data sets, constructed from large genomewide data sets, presents significant computational challenges. To address these challenges, we present parallelnewhybrid, an r package designed to decrease user burden when undertaking multiple newhybrids analyses. parallelnewhybrid does so by taking advantage of the parallel computational capabilities inherent in modern computers to efficiently and automatically execute separate newhybrids runs in parallel. We show that parallelization of analyses using this package affords users several-fold reductions in time over a traditional serial analysis. parallelnewhybrid consists of an example data set, a readme and three operating system-specific functions to execute parallel newhybrids analyses on each of a computer's c cores. parallelnewhybrid is freely available on the long-term software hosting site github (www.github.com/bwringe/parallelnewhybrid).