Population differentiation determined from putative neutral and divergent adaptive genetic markers in Eulachon (Thaleichthys pacificus,Osmeridae), an anadromous Pacific smelt

Twelve eulachon (Thaleichthys pacificus, Osmeridae) populations ranging from Cook Inlet, Alaska and along the west coast of North America to the Columbia River were examined by restriction‐site‐associated DNA (RAD) sequencing to elucidate patterns of neutral and adaptive variation in this high geneflow species. A total of 4104 single‐nucleotide polymorphisms (SNPs) were discovered across the genome, with 193 putatively adaptive SNPs as determined by F_ST outlier tests. Estimates of population structure in eulachon with the putatively adaptive SNPs were similar, but provided greater resolution of stocks compared with a putatively neutral panel of 3911 SNPs or previous estimates with 14 microsatellites. A cline of increasing measures of genetic diversity from south to north was found in the adaptive panel, but not in the neutral markers (SNPs or microsatellites). This may indicate divergent selective pressures in differing freshwater and marine environments between regional eulachon populations and that these adaptive diversity patterns not seen with neutral markers could be a consideration when determining genetic boundaries for conservation purposes. Estimates of effective population size (N_e) were similar with the neutral SNP panel and microsatellites and may be utilized to monitor population status for eulachon where census sizes are difficult to obtain. Greater differentiation with the panel of putatively adaptive SNPs provided higher individual assignment accuracy compared to the neutral panel or microsatellites for stock identification purposes. This study presents the first SNPs that have been developed for eulachon, and analyses with these markers highlighted the importance of integrating genome‐wide neutral and adaptive genetic variation for the applications of conservation and management.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Genome‐scale sampling suggests cryptic epigenetic structuring and insular divergence in Canada lynx

Determining the molecular signatures of adaptive differentiation is a fundamental component of evolutionary biology. A key challenge is to identify such signatures in wild organisms, particularly between populations of highly mobile species that undergo substantial gene flow. The Canada lynx (Lynx canadensis) is one species where mainland populations appear largely undifferentiated at traditional genetic markers, despite inhabiting diverse environments and displaying phenotypic variation. Here, we used high‐throughput sequencing to investigate both neutral genetic structure and epigenetic differentiation across the distributional range of Canada lynx. Newfoundland lynx were identified as the most differentiated population at neutral genetic markers, with demographic modelling suggesting that divergence from the mainland occurred at the end of the last glaciation (20–33 KYA). In contrast, epigenetic structure revealed hidden levels of differentiation across the range coincident with environmental determinants including winter conditions, particularly in the peripheral Newfoundland and Alaskan populations. Several biological pathways related to morphology were differentially methylated between populations, suggesting that epigenetic modifications might explain morphological differences seen between geographically peripheral populations. Our results indicate that epigenetic modifications, specifically DNA methylation, are powerful markers to investigate population differentiation in wild and non‐model systems.     

Discovery and characterization of single nucleotide polymorphisms in coho salmon,Oncorhynchus kisutch

Molecular population genetic analyses have become an integral part of ecological investigation and population monitoring for conservation and management. Microsatellites have been the molecular marker of choice for such applications over the last several decades, but single nucleotide polymorphism (SNP) markers are rapidly expanding beyond model organisms. Coho salmon (Oncorhynchus kisutch) is native to the north Pacific Ocean and its tributaries, where it is the focus of intensive fishery and conservation activities. As it is an anadromous species, coho salmon typically migrate across multiple jurisdictional boundaries, complicating management and requiring shared data collection methods. Here, we describe the discovery and validation of a suite of novel SNPs and associated genotyping assays which can be used in the genetic analyses of this species. These assays include 91 that are polymorphic in the species and one that discriminates it from a sister species, Chinook salmon. We demonstrate the utility of these SNPs for population assignment and phylogeographic analyses, and map them against the draft trout genome. The markers constitute a large majority of all SNP markers described for coho salmon and will enable both population‐ and pedigree‐based analyses across the southern part of the species native range.     

Riverscape genomics of a threatened fish across a hydroclimatically heterogeneous river basin

Understanding how natural selection generates and maintains adaptive genetic diversity in heterogeneous environments is key to predicting the evolutionary response of populations to rapid environmental change. Detecting selection in complex spatial environments remains challenging, especially for threatened species where the effects of strong genetic drift may overwhelm signatures of selection. We carried out a basinwide riverscape genomic analysis in the threatened southern pygmy perch (Nannoperca australis), an ecological specialist with low dispersal potential. High‐resolution environmental data and 5162 high‐quality filtered SNPs were used to clarify spatial population structure and to assess footprints of selection associated with a steep hydroclimatic gradient and with human disturbance across the naturally and anthropogenically fragmented Murray–Darling Basin (Australia). Our approach included F_ST outlier tests to define neutral loci, and a combination of spatially explicit genotype–environment association analyses to identify candidate adaptive loci while controlling for the effects of landscape structure and shared population history. We found low levels of genetic diversity and strong neutral population structure consistent with expectations based on spatial stream hierarchy and life history. In contrast, variables related to precipitation and temperature appeared as the most important environmental surrogates for putatively adaptive genetic variation at both regional and local scales. Human disturbance also influenced the variation in candidate loci for adaptation, but only at a local scale. Our study contributes to understanding of adaptive evolution along naturally and anthropogenically fragmented ecosystems. It also offers a tangible example of the potential contributions of landscape genomics for informing in situ and ex situ conservation management of biodiversity.     

The population genomic signature of environmental association and gene flow in an ecologically divergent tree species Metrosideros polymorpha (Myrtaceae)

Genomewide markers enable us to study genetic differentiation within a species and the factors underlying it at a much higher resolution than before, which advances our understanding of adaptation in organisms. We investigated genomic divergence in Metrosideros polymorpha, a woody species that occupies a wide range of ecological habitats across the Hawaiian Islands and shows remarkable phenotypic variation. Using 1659 single nucleotide polymorphism (SNP) markers annotated with the genome assembly, we examined the population genetic structure and demographic history of nine populations across five elevations and two ages of substrates on Mauna Loa, the island of Hawaii. The nine populations were differentiated into two genetic clusters distributed on the lower and higher elevations and were largely admixed on the middle elevation. Demographic modelling revealed that the two genetic clusters have been maintained in the face of gene flow, and the effective population size of the high‐altitude cluster was much smaller. A F_ST‐based outlier search among the 1659 SNPs revealed that 34 SNPs (2.05%) were likely to be under divergent selection and the allele frequencies of 21 of them were associated with environmental changes along elevations, such as temperature and precipitation. This study shows a genomic mosaic of M. polymorpha, in which contrasting divergence patterns were found. While most genomic polymorphisms were shared among populations, a small fraction of the genome was significantly differentiated between populations in diverse environments and could be responsible for the dramatic adaptation to a wide range of environments.     

SNP discovery in wild and domesticated populations of blue catfish,Ictalurus furcatus,using genotyping‐by‐sequencing and subsequent SNP validation

Blue catfish, Ictalurus furcatus, are valued in the United States as a trophy fishery for their capacity to reach large sizes, sometimes exceeding 45 kg. Additionally, blue catfish × channel catfish (I. punctatus) hybrid food fish production has recently increased the demand for blue catfish broodstock. However, there has been little study of the genetic impacts and interaction of farmed, introduced and stocked populations of blue catfish. We utilized genotyping‐by‐sequencing (GBS) to capture and genotype SNP markers on 190 individuals from five wild and domesticated populations (Mississippi River, Missouri, D&B, Rio Grande and Texas). Stringent filtering of SNP‐calling parameters resulted in 4275 SNP loci represented across all five populations. Population genetics and structure analyses revealed potential shared ancestry and admixture between populations. We utilized the Sequenom MassARRAY to validate two multiplex panels of SNPs selected from the GBS data. Selection criteria included SNPs shared between populations, SNPs specific to populations, number of reads per individual and number of individuals genotyped by GBS. Putative SNPs were validated in the discovery population and in two additional populations not used in the GBS analysis. A total of 64 SNPs were genotyped successfully in 191 individuals from nine populations. Our results should guide the development of highly informative, flexible genotyping multiplexes for blue catfish from the larger GBS SNP set as well as provide an example of a rapid, low‐cost approach to generate and genotype informative marker loci in aquatic species with minimal previous genetic information.     

New view of population genetics of zooplankton: RAD‐seq analysis reveals population structure of the North Atlantic planktonic copepod Centropages typicus

Detection of population genetic structure of zooplankton at medium‐to‐small spatial scales in the absence of physical barriers has remained challenging and controversial. The large population sizes and high rates of gene flow characteristic of zooplankton have made resolution of geographical differentiation very difficult, especially when using few genetic markers and assuming equilibrium conditions. Next‐generation sequencing now allows simultaneous sampling of hundreds to thousands of genetic markers; new analytical approaches allow studies under nonequilibrium conditions and directional migration. Samples of the North Atlantic Ocean planktonic copepod, Centropages typicus, were analysed using restriction site‐associated DNA (RAD) sequencing on a PROTON platform. Although prior studies revealed no genetic differentiation of populations across the geographical range of the species, analysis of RAD tags showed significant structure across the North Atlantic Ocean. We also compared the likelihood for models of connectivity among NW Atlantic populations under various directional flow scenarios that replicate oceanographic conditions of the sampled domain. High‐density marker sampling with RAD sequencing markedly outperformed other technical and analytical approaches in detection of population genetic structure and characterization of connectivity of this high geneflow zooplankton species.     

Characterizing genic and nongenic molecular markers: comparison of microsatellites and SNPs

The implications of transitioning to single nucleotide polymorphism (SNPs) from microsatellite markers (MSs) have been investigated in a number of population genetics studies, but the effect of genomic location on the amount of information each type of marker reveals has not been explored in detail. We developed novel SNP markers flanking 1 kb regions of 13 genic (within gene or <1 kb away from gene) and 13 nongenic (>10 kb from annotated gene) MSs in the threespine stickleback genome to obtain comparable data for both types of markers. We analysed patterns of genetic diversity and divergence on various geographic scales after converting the SNP loci within each genomic region into haplotypes. Marker type (SNP haplotype or MS) and location (genic or nongenic) significantly affected most estimates of population diversity and divergence. Between‐lineage divergence was significantly higher in SNP haplotypes (genic and nongenic), however, within‐lineage divergence was similar between marker types. Most divergence and diversity measures were uncorrelated between markers, except for population differentiation which was correlated between MSs and SNP haplotypes (both genic and nongenic). Broad‐scale population structure and assignment were similarly resolved by both marker types, however, only the MSs were able to delimit fine‐scale population structuring, particularly when genic and nongenic markers were combined. These results demonstrate that estimates of genetic variability and differentiation among populations can be strongly influenced by marker type, their genomic location in relation to genes and by the interaction of these two factors. This highlights the importance of having an awareness of the inherent strengths and limitations associated with different molecular tools to select the most appropriate methods for accurately addressing various ecological and evolutionary questions.     

Standing genetic diversity and selection at functional gene loci are associated with differential invasion success in two non‐native fish species

Invasive species are expected to experience a unique combination of high genetic drift due to demographic factors while also experiencing strong selective pressures. The paradigm that reduced genetic diversity should limit the evolutionary potential of invasive species, and thus, their potential for range expansion has received little empirical support, possibly due to the choice of genetic markers. Our goal was to test for effects of genetic drift and selection at functional genetic markers as they relate to the invasion success of two paired invasive goby species, one widespread (successful) and one with limited range expansion (less successful). We genotyped fish using two marker types: single nucleotide polymorphisms (SNPs) in known‐function, protein‐coding genes and microsatellites to contrast the effects of neutral genetic processes. We identified reduced allelic variation in the invaded range for the less successful tubenose goby. SNPs putatively under selection were responsible for the observed differences in population structure between marker types for round goby (successful) but not tubenose goby (less successful). A higher proportion of functional loci experienced divergent selection for round goby, suggesting increased evolutionary potential in invaded ranges may be associated with round goby's greater invasion success. Genes involved in thermal tolerance were divergent for round goby populations but not tubenose goby, consistent with the hypothesis that invasion success for fish in temperate regions is influenced by capacity for thermal tolerance. Our results highlight the need to incorporate functional genetic markers in studies to better assess evolutionary potential for the improved conservation and management of species.     

Breaking RAD: an evaluation of the utility of restriction site‐associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site‐associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome‐scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq‐based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD‐tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best‐case modelling scenario, we found that RADseq using six‐ or eight‐base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD‐tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Is MHC diversity a better marker for conservation than neutral genetic diversity? A case study of two contrasting dolphin populations

Genetic diversity is essential for populations to adapt to changing environments. Measures of genetic diversity are often based on selectively neutral markers, such as microsatellites. Genetic diversity to guide conservation management, however, is better reflected by adaptive markers, including genes of the major histocompatibility complex (MHC). Our aim was to assess MHC and neutral genetic diversity in two contrasting bottlenose dolphin (Tursiops aduncus) populations in Western Australia—one apparently viable population with high reproductive output (Shark Bay) and one with lower reproductive output that was forecast to decline (Bunbury). We assessed genetic variation in the two populations by sequencing the MHC class II DQB, which encompasses the functionally important peptide binding regions (PBR). Neutral genetic diversity was assessed by genotyping twenty‐three microsatellite loci. We confirmed that MHC is an adaptive marker in both populations. Overall, the Shark Bay population exhibited greater MHC diversity than the Bunbury population—for example, it displayed greater MHC nucleotide diversity. In contrast, the difference in microsatellite diversity between the two populations was comparatively low. Our findings are consistent with the hypothesis that viable populations typically display greater genetic diversity than less viable populations. The results also suggest that MHC variation is more closely associated with population viability than neutral genetic variation. Although the inferences from our findings are limited, because we only compared two populations, our results add to a growing number of studies that highlight the usefulness of MHC as a potentially suitable genetic marker for animal conservation. The Shark Bay population, which carries greater adaptive genetic diversity than the Bunbury population, is thus likely more robust to natural or human‐induced changes to the coastal ecosystem it inhabits.     

Levels and patterns of population genetic diversity in the red seaweed Chondrus crispus (Florideophyceae): a direct comparison of single nucleotide polymorphisms and microsatellites

Single nucleotide polymorphisms (SNPs) are predicted to supersede microsatellites as the marker of choice for population genetic studies in the near future. To date, however, very few studies have directly compared both marker systems in natural populations, particularly in non‐model organisms. In the present study, we compared the utility of SNPs and microsatellites for population genetic analysis of the red seaweed Chondrus crispus (Florideophyceae). Six SNP loci yielded very different patterns of intrapopulation genetic diversity compared to those obtained using seven moderately (mean 5.2 alleles) polymorphic microsatellite loci, although Bayesian clustering analysis gave largely congruent results between the two marker classes. A weak but significant pattern of isolation‐by‐distance was observed across scales from a few hundred metres to approximately 200 km using the combined SNP and microsatellite data set of 13 loci. Over larger scales, however, there was little correlation between genetic divergence and geographical distance. Our findings suggest that even a moderate number of SNPs is sufficient to determine patterns of genetic diversity across natural populations, and also highlight the fact that patterns of genetic variation in seaweeds arise through a complex interplay of short‐ and long‐term natural processes, as well as anthropogenic influence. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 108 , 251–262.     

Population genomic analysis suggests strong influence of river network on spatial distribution of genetic variation in invasive saltcedar across the southwestern United States

Understanding the complex influences of landscape and anthropogenic elements that shape the population genetic structure of invasive species provides insight into patterns of colonization and spread. The application of landscape genomics techniques to these questions may offer detailed, previously undocumented insights into factors influencing species invasions. We investigated the spatial pattern of genetic variation and the influences of landscape factors on population similarity in an invasive riparian shrub, saltcedar (Tamarix L.) by analysing 1,997 genomewide SNP markers for 259 individuals from 25 populations collected throughout the southwestern United States. Our results revealed a broad‐scale spatial genetic differentiation of saltcedar populations between the Colorado and Rio Grande river basins and identified potential barriers to population similarity along both river systems. River pathways most strongly contributed to population similarity. In contrast, low temperature and dams likely served as barriers to population similarity. We hypothesize that large‐scale geographic patterns in genetic diversity resulted from a combination of early introductions from distinct populations, the subsequent influence of natural selection, dispersal barriers and founder effects during range expansion.     

Genetic isolation of fragmented populations is exacerbated by drift and selection

Reduced genetic variation at marker loci in small populations has been well documented, whereas the relationship between quantitative genetic variation and population size has attracted little empirical investigation. Here we demonstrate that both neutral and quantitative genetic variation are reduced in small populations of a fragmented plant metapopulation, and that both drift and selective change are enhanced in small populations. Measures of neutral genetic differentiation (F(ST)) and quantitative genetic differentiation (Q(ST)) in two traits were higher among small demes, and Q(ST) between small populations exceeded that expected from drift alone. This suggests that fragmented populations experience both enhanced genetic drift and divergent selection on phenotypic traits, and that drift affects variation in both neutral markers and quantitative traits. These results highlight the need to integrate natural selection into conservation genetic theory, and suggests that small populations may represent reservoirs of genetic variation adaptive within a wide range of environments.     

Unraveling hierarchical genetic structure in a marine metapopulation: A comparison of three high‐throughput genotyping approaches

Marine metapopulations often exhibit subtle population structure that can be difficult to detect. Given recent advances in high‐throughput sequencing, an emerging question is whether various genotyping approaches, in concert with improved sampling designs, will substantially improve our understanding of genetic structure in the sea. To address this question, we explored hierarchical patterns of structure in the coral reef fish Elacatinus lori using a high‐resolution approach with respect to both genetic and geographic sampling. Previously, we identified three putative E. lori populations within Belize using traditional genetic markers and sparse geographic sampling: barrier reef and Turneffe Atoll; Glover's Atoll; and Lighthouse Atoll. Here, we systematically sampled individuals at ~10 km intervals throughout these reefs (1,129 individuals from 35 sites) and sequenced all individuals at three sets of markers: 2,418 SNPs; 89 microsatellites; and 57 nonrepetitive nuclear loci. At broad spatial scales, the markers were consistent with each other and with previous findings. At finer spatial scales, there was new evidence of genetic substructure, but our three marker sets differed slightly in their ability to detect these patterns. Specifically, we found subtle structure between the barrier reef and Turneffe Atoll, with SNPs resolving this pattern most effectively. We also documented isolation by distance within the barrier reef. Sensitivity analyses revealed that the number of loci (and alleles) had a strong effect on the detection of structure for all three marker sets, particularly at small spatial scales. Taken together, these results illustrate empirically that high‐throughput genotyping data can elucidate subtle genetic structure at previously‐undetected scales in a dispersive marine fish.     

Genetic variation across species' geographical ranges: the central-marginal hypothesis and beyond

There is growing interest in quantifying genetic population structure across the geographical ranges of species to understand why species might exhibit stable range limits and to assess the conservation value of peripheral populations. However, many assertions regarding peripheral populations rest on the long-standing but poorly tested supposition that peripheral populations exhibit low genetic diversity and greater genetic differentiation as a consequence of smaller effective population size and greater geographical isolation relative to geographically central populations. We reviewed 134 studies representing 115 species that tested for declines in within-population genetic diversity and/or increases in among-population differentiation towards range margins using nuclear molecular genetic markers. On average, 64.2% of studies detected the expected decline in diversity, 70.2% of those that tested for it showed increased differentiation and there was a positive association between these trends. In most cases, however, the difference in genetic diversity between central and peripheral population was not large. Although these results were consistent across plants and animals, strong taxonomic and biogeographical biases in the available studies call for a cautious generalization of these results. Despite the large number of studies testing these simple predictions, very few attempted to test possible mechanisms causing reduced peripheral diversity or increased differentiation. Almost no study incorporated a phylogeographical framework to evaluate historical influences on contemporary genetic patterns. Finally, there has been little effort to test whether these geographical trends in putatively neutral variation at marker loci are reflected by quantitative genetic trait variation, which is likely to influence the adaptive potential of populations across the geographical range.     

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.)

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well‐known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.