Limited geographic genetic structure detected in a widespread Palearctic corvid, Nucifraga caryocatactes

Introduction

In Eurasia, phylogeographic studies of many widespread vertebrate species have revealed a variety of geographical patterns of population structure influenced by current and historical landscapes, with little overall consensus among species. Using mitochondrial DNA, east–west splits have been documented for a variety of vertebrates including bats (Flanders et al., 2009), and several avian species (e.g., Eurasian magpie (Pica pica; Kryukov et al., 2004), rook (Corvus frugilegus; Haring, Gamauf & Kryukov, 2007), and red-breasted flycatcher (Ficedula parva; Zink et al., 2008)). For other species, multiple splits have occurred (e.g., root vole (Microtus oeconomus; Brunhoff et al., 2003) and reed bunting (Emberiza schoeniclus; Zink et al., 2008)), or peninsula populations are isolated (e.g., great bustard (Otis tarda; Pitra, Lieckfeldt & Alonso, 2000)). In contrast, little population structure has been detected in some widespread species, such as otters (Lutra lutra; Ferrando et al., 2004) and several avifauna species (e.g., great spotted woodpecker (Dendrocopos major; Zink, Drovetski & Rohwer, 2002), common sandpiper (Actitis hypoleucos; Zink et al., 2008), and Eurasian magpie (Pica pica; Zhang et al., 2012)). Some of the observed phylogeographic patterns have been explained by post-glacial colonization from single or multiple refugia, but may also be influenced by barriers to dispersal, such as mountain ranges (e.g., Ural Mountains), large areas of inhospitable habitat (e.g., Tibetan Plateau), or large bodies of water.

The Eurasian nutcracker (Nucifraga caryocatactes, Linnaeus, 1758) is a corvid with a widespread Palearctic distribution. Although generally classified as a resident species of continental coniferous forests, nutcrackers are known to irruptively disperse to take advantage of mast conifer seed crops (Haring, Gamauf & Kryukov, 2007), similar to its North American sister species, Clark’s nutcracker (N. columbiana; Tomback, 1998). Strong geographic genetic structure has not been found in Clark’s nutcracker, despite numerous potential physical barriers to dispersal and thus gene flow (Dohms & Burg, 2013). A previous study by Haring and colleagues (2007) found no population structure in N. caryocatactes throughout Eurasia. However, Haring, Gamauf & Kryukov (2007) only used 11 specimens, thus additional data may shed further light on nutcracker population genetic structure.

In this study, we use a highly variable and rapidly evolving mitochondrial DNA marker, the control region, to further investigate population structure of N. caryocatactes in Eurasia. Based on the ecology of this species, we predict little range-wide population genetic structure.

Materials & Methods

Samples

Tissue samples (n = 62) collected throughout the Eurasian nutcracker’s range (Fig. 1) were acquired from the Burke Museum of Natural History and Culture at the University of Washington (Table S1). Previously published control region (CR) sequences (n = 11) were obtained from GenBank (EU070770 and EU070886–EU070895; Haring, Gamauf & Kryukov, 2007).

Results

Genetic analyses

We sequenced and aligned the 924 bp control region (CR) sequences for 62 individuals from 12 populations (Table 1; GenBank accession nos. KJ999615–KJ999676). We aligned 11 additional GenBank sequences (Haring, Gamauf & Kryukov, 2007) with sequences from our samples and obtained a 305 bp area of overlap. Statistical parsimony networks did not suggest strong geographic structure for the 924 bp sequence (Fig. 2), nor for the larger dataset using the overlapping 305 bp fragment for all 73 individuals (Fig. 3). Ncarcar5 could not be connected to the network in the larger dataset, which was also found by Haring, Gamauf & Kryukov (2007).

Table 1:

Diversity of a 924 bp mitochondrial DNA sequence.

Haplotype diversity for a 924 bp fragment of mtDNA control region from 12 populations of Nucifraga caryocatactes throughout Eurasia.

Population	Location	n	Hn	Hd	π
BK	Badzhal’skiy Krebet, Russia	3	3	1.000	0.00361
GA	Gorno-Altaysk, Russia	4	4	1.000	0.00328
KA	Kamchatka, Russia	6	4	0.800	0.00369
KD	Irkutsk Oblast, Russia	5	5	1.000	0.00523
MA	Magadanskaya Oblast, Russia	8	8	1.000	0.00397
OK	Kyzyl, Russia	1	1	-	-
OR	Ola River headwaters, Russia	2	1	0.000	0.00000
PK	Primorsky Kray, Russia	11	11	1.000	0.00993
SO	Sakhalinksya Oblast, Russia	5	4	0.900	0.00195
TO	Tyumenskaya Oblast, Russia	7	6	0.952	0.00622
TU	Irkutsk Oblast, Russia	5	4	0.900	0.00326
UL	Ulaanbaatar, Mongolia	5	4	0.900	0.00611
Overall		62	45	0.967	0.00537

DOI: 10.7717/peerj.371/table-1

Notes:

n

number of individuals in population

H_n

number of haplotypes

H_d

haplotype diversity

π

nucleotide diversity within the population

For the 62 individuals we sequenced, we found 45 unique haplotypes and high levels of genetic diversity in most populations (Table 1). We found 57 polymorphic sites within the 924 bp sequence and 22 within the 305 bp sequence. Haplotype diversity for the 924 bp sequence varied from 0.000 (Ola River headwaters (OR)) to 1.000 (five populations), and all but the OR population had a haplotype diversity equal to or greater than 0.800 (Table 1). Nucleotide diversity ranged from 0.00000 (OR) to 0.00993 (PK; Table 1). Overall haplotype diversity (H_d) 0.967 and nucleotide diversity (π) was 0.00537.

Discussion

As predicted, analyses of Eurasian nutcracker mitochondrial DNA control region sequences did not detect significant population genetic structure. All populations except nutcrackers from the Ola River headwaters (OR; n = 2) exhibited high haplotype diversity and relatively high nucleotide diversity. No geographic clustering was observed in statistical parsimony networks, even when integrating samples from the western part of the range. Despite potential barriers to dispersal for this species, such as isolation on an island (e.g., Sakhalin Oblast (SO)) or peninsula (e.g., Kamchatka (KA)), most populations of N. caryocatactes do not appear to be geographically differentiated from each other, likely due to gene flow during irruptive dispersal in search of mast pine seed crops. Overall, our work supports that done by Haring, Gamauf & Kryukov (2007) where no significant split was seen between the east and west for nutcrackers and is similar to the pattern found in N. caryocatactes sister species, N. columbiana (Dohms & Burg, 2013).

Compared to Haring, Gamauf & Kryukov (2007), our study found a higher level of haplotype diversity (H_d = 0.967 vs 0.844 and π = 0.00537 vs 0.00279). This may be due to the portion of control region sequenced and the larger sample sizes used in this study. The sequences obtained from our samples are predominantly composed of domains I and II of the mtDNA control region (Saunders & Edwards, 2000), whereas Haring, Gamauf & Kryukov (2007) sequenced primarily domain II, which is considered less variable (Ruokonen & Kvist, 2002).

Haring, Gamauf & Kryukov (2007) state that low genetic diversity may suggest a bottleneck in this species and a single glacial refugium. With high levels of nucleotide and haplotype diversity, our findings do not support a historical bottleneck but rather point toward two possible scenarios: multiple refugia with gene flow or a single refugium with large population size during expansion. Expansion from multiple refugia with gene flow after colonization can produce similar genetic patterns to those in species that expanded slowly from a single refugium with a large population size, retaining high genetic diversity and limited geographic structuring of populations (Hewitt, 2004). Given the dispersal potential of nutcrackers, it is possible that a large population expanded out of a single refugium, but it is equally plausible that expansion occurred out of multiple refugia with subsequent gene flow between geographically distinct populations due to irruptive dispersal events. The multiple refugia scenario could have produced the large number of haplotypes, often with high levels of sequence divergence pattern seen here. For example, individuals from Primorsky Kray (PK) are found scattered throughout the parsimony network (Fig. 2), in some cases with a large number of mutations between PK individuals and other haplotypes, yet found clustered with geographically distant individuals from Irkutsk Oblast (KD) and KA. This level of divergence is often associated with isolation in and subsequent colonization from multiple refugia (Hewitt, 2004). With the high haplotype diversity across all populations, it is not possible to determine which population(s), if any sampled here, may be in the location of the original refugium or refugia. Without additional present day samples from the Alps and Himalaya Mountain ranges, it is difficult to tell using genetic signatures if these areas served as refugia for nutcrackers during the LGM.

Our findings do not support a single refugium in the Altai Mountains of southern Mongolia, as postulated by Haring, Gamauf & Kryukov (2007). Rather, our data show highly divergent haplotypes which could be the result of prolonged isolation in multiple refugia. Scots pine (Pinus sylvestris), an important source of food for nutcrackers, is thought to have survived in refugia near the Alps (Naydenov et al., 2007) and in the east, unglaciated portions of the Himalayas could have served as a refugium for high latitude species (Zhuo, Baoyin & Petit-Maire, 1998; Owen, Finkel & Caffee, 2002). Alternatively, high levels of plant endemism have been found in the mountains of southern and eastern China, suggestive of long-term suitable habitats (Zhuo, Baoyin & Petit-Maire, 1998). Nutcrackers may have survived in these bands of suitable habitat in the southwest and southeast areas of Eurasia and expanded northward from multiple refugia as glaciers retreated.

Conclusions

Overall, Eurasian nutcrackers exhibit limited geographic genetic structure throughout their range, as might be expected from a resident bird with irruptive dispersal patterns. Our study found high genetic diversity, which suggests that a population bottleneck has not occurred in this species as previously hypothesized. A more detailed phylogeographical study could include additional genetic sampling from northern and southern parts of N. caryocatactes’ range to further investigate structure across the range of this species.

Supplemental Information