Proliferation of East Antarctic Adélie penguins in response to historical deglaciation
BMC Evolutionary Biology volume 15, Article number: 236 (2015)
Major, long-term environmental changes are projected in the Southern Ocean and these are likely to have impacts for marine predators such as the Adélie penguin (Pygoscelis adeliae). Decadal monitoring studies have provided insight into the short-term environmental sensitivities of Adélie penguin populations, particularly to sea ice changes. However, given the long-term nature of projected climate change, it is also prudent to consider the responses of populations to environmental change over longer time scales. We investigated the population trajectory of Adélie penguins during the last glacial-interglacial transition to determine how the species was affected by climate warming over millennia. We focussed our study on East Antarctica, which is home to 30 % of the global population of Adélie penguins.
Using mitochondrial DNA from extant colonies, we reconstructed the population trend of Adélie penguins in East Antarctica over the past 22,000 years using an extended Bayesian skyline plot method. To determine the relationship of East Antarctic Adélie penguins with populations elsewhere in Antarctica we constructed a phylogeny using mitochondrial DNA sequences.
We found that the Adélie penguin population expanded 135-fold from approximately 14,000 years ago. The population growth was coincident with deglaciation in East Antarctica and, therefore, an increase in ice-free ground suitable for Adélie penguin nesting. Our phylogenetic analysis indicated that East Antarctic Adélie penguins share a common ancestor with Adélie penguins from the Antarctic Peninsula and Scotia Arc, with an estimated age of 29,000 years ago, in the midst of the last glacial period. This finding suggests that extant colonies in East Antarctica, the Scotia Arc and the Antarctic Peninsula were founded from a single glacial refuge.
While changes in sea ice conditions are a critical driver of Adélie penguin population success over decadal and yearly timescales, deglaciation appears to have been the key driver of population change over millennia. This suggests that environmental drivers of population trends over thousands of years may differ to drivers over years or decades, highlighting the need to consider millennial-scale trends alongside contemporary data for the forecasting of species’ abundance and distribution changes under future climate change scenarios.
While climate change is a global phenomenon, its environmental effects can vary dramatically in different locations. For example, trends in the extent and duration of sea ice around the Antarctic continent show high spatial heterogeneity [1, 2]. Over a 34 year monitoring period, sea ice extent decreased in the Bellingshausen and Amundsen Seas accompanied by a dramatic shortening of the sea ice season by 100 ± 31 days [1, 3]. Meanwhile, in the Ross Sea, both extent and duration of the sea ice season increased substantially over the same period . Even within Antarctic regions there have been variations in the extent and duration of sea ice. For example, the East Antarctic region, defined here as between 30 and 150°E, has demonstrated considerably more complex trends in sea ice seasonality and extent than the rest of the continent [4, 5]. Since 1980, in some East Antarctic areas (between 95 and 110°E; and isolated pockets between 75 and 150°E), there has been a significant shortening of the sea ice season by up to 93 days . Meanwhile, many other localities (west of 105°E; between 40 and 90°E) have experienced a significant lengthening of the sea ice season . Antarctic sea ice is expected to undergo further declines in the future; in the most extreme climate model scenario (RCP8.5), East Antarctica would be completely free of sea ice in February by 2081–2100, while only small portions of the Weddell and Ross Seas would retain sea ice .
Environmental changes that are currently underway in the Southern Ocean, including changes in the timing of sea ice advance and retreat, receding glaciers, and shifting oceanographic fronts, could lead to major changes in the terrestrial breeding habitats, marine foraging environment and prey availability for higher order predators [7, 8]. The Adélie penguin (Pygoscelis adeliae) is a prime example of a Southern Ocean predator that is likely to be affected by environmental impacts associated with the Anthropocene . The species forms breeding colonies on ice-free land along the Antarctic coastline (Fig. 1) [10, 11] and forages largely in the pack ice zone during the breeding season . Adélie penguin populations are known to be sensitive to changes in sea ice extent, the timing of sea ice retreat [13, 14] and the extent of glaciation . Given the large regional variability in the physical manifestations of climate change, it is likely that geographically distant populations of Adélie penguins will experience very different environmental impacts .
East Antarctica is home to approximately 30 % of the global population of Adélie penguins, with an estimated abundance of 1.14 million breeding pairs . The breeding distribution in this region has expanded over the past several decades, possibly as a result of sustained population growth . Sea ice conditions strongly influence Adélie penguin populations in East Antarctica, although the mechanisms of impact are complex and depend on the nature, extent and timing of the presence of sea ice [13, 18]. The Béchervaise Island population experienced near total reproductive failure in years with extensive near-shore sea ice during the guard stage . The proposed mechanism for this impact was a reduced efficiency of chick provisioning at a crucial time in the breeding season, as more extensive sea ice increases the duration of the adults’ foraging trips and reduces the frequency at which the chicks are fed . The negative impact of unusually extensive sea ice was also felt at the Pétrels Island colony at Pointe Géologie in the 2013/14 breeding season, with emaciated chicks often observed during the summer . In this case, the negative impact of extensive sea ice was compounded by unusually warm air temperatures that caused snow melt and unprecedented amounts of rain . Rain can be fatal to Adélie penguin chicks, as their downy plumage is not waterproof and when wet they may succumb to cold temperatures . In the case of Pétrels Island, the compounded effects of extensive sea ice and warm temperatures were devastating, resulting in 100 % chick mortality in 2013/14 . In the opposite scenario, at Pointe Géologie the size of breeding populations increased six years after a period of low sea ice extent and concentration . As Adélie penguins commence breeding between 5 and 6 years of age, it is likely that sea ice conditions during the fledgling and yearling stages are important, with lower sea ice extent being favourable for young birds . At Béchervaise Island adult penguins were also sensitive to extremes in sea ice concentrations in their winter foraging grounds, with either too much sea ice (>80 % cover) or too little sea ice (<15 % cover) negatively impacting adult survival . Extreme climatic events that alter the sea ice environment may also impact Adélie penguin reproductive success. In 2010 there was a calving of the Mertz Glacier Tongue, which resulted in decreased polynya activity and sea ice production in the area [21, 22] and these local changes in the icescape appear to have negatively impacted Adélie penguin reproductive success at the nearby Pointe Géologie populations in the 2011/2012 and 2012/2013 breeding seasons .
While decadal scale monitoring studies have provided invaluable data on the short-term environmental sensitivities of Adélie penguin populations, given the long-term environmental change projected in the Southern Ocean, it is also prudent to consider the responses of populations to environmental change over longer time scales (e.g. thousands to tens of thousands of years) and during climate regime shifts, for example, during the transition from the last glacial maximum (LGM, 26–19.5 kya; ) to the Holocene warming period (11.7 kya–present) . The responses of Adélie penguins to climate change during the LGM and Holocene have been well-studied in the Ross Sea and Antarctic Peninsula/Scotia Arc regions (Fig. 1) [26, 27]. In both locations, Adélie penguin numbers were much lower during the LGM than they are today, which was thought to be the result of reduced ice-free ground suitable for breeding [28, 29]. Phylogenetic studies have found evidence of two genetic lineages of Adélie penguins that are suggestive of two refuge populations dating to the LGM [29–31]. One of these lineages was comprised solely of individuals from modern Ross Sea colonies, suggesting that a refuge may have been situated somewhere in the vicinity of the Ross Sea during the LGM . The second lineage was comprised of individuals from the Antarctic Peninsula, Scotia Arc, Ross Sea, and East Antarctica [29–31]. However, the samples analysed for the East Antarctic region were limited to two colonies at Gardner and Welch Islands and, given the length of the coastline, may not be representative of the genetic diversity of the broader region. As this second genetic lineage shows no strong geographic affinity, the location of its associated LGM refuge is unknown. Based on a genetic coalescent study, Adélie penguin numbers in the Antarctic Peninsula/Scotia Arc increased during the Holocene warming period from ca. 16 kya , roughly coincident with deglaciation of the region [32, 33]. In the Ross Sea, radiocarbon dated remains suggest that the Adélie penguin distribution expanded from approximately 8 kya, followed by two periods of reduced occupation from 5 to 4 kya and 2 to 1.1 kya .
The impacts of past climate regime shifts on Adélie penguins have been less well studied in East Antarctica. In the Windmill Islands and Vestfold Hills, radiocarbon dating of remains suggests that Adélie penguins were present by 9 kya and 8.5 kya respectively, roughly coincident with local deglaciation and, therefore, an emergence of available ice-free nesting habitat [35, 36]. There is also evidence for a peak in Adélie penguin numbers ca. 4 kya, coinciding with a mid-Holocene warm period [27, 36, 37]. During the LGM, the ancestors of individuals currently breeding at Gardner Island in the Prydz Bay region and Welch Island on the Mawson Coast were restricted to a refuge population along with the ancestors of extant Antarctic Peninsula/Scotia Arc and Ross Sea colonies [29–31]. However, as mentioned previously, the location of this refuge is unknown. Several questions regarding the past dynamics of East Antarctic Adélie penguin populations remain to be answered, including the trends in abundance during and after the LGM, whether an LGM refuge population may have been located in this area, and how populations may have responded to environmental changes across the extensive East Antarctic region. In this study we sought to address these questions using genetic data from multiple extant Adélie penguin colonies across East Antarctica.
We sequenced 56 individuals from six colonies located in East Antarctica (Béchervaise Island, Macey Island, Welch Island, Blakeney Point, Holl Island and Pétrels Island; Fig. 1). A 642 bp fragment of the mitochondrial hypervariable region (HVR) was sequenced for each individual [GenBank: KT932437 - KT932492], and a 902 bp fragment of cytochrome b (CytB) was sequenced for 45 individuals [GenBank: KT932493 - KT932537].
Genetic diversity was high for the HVR, with 85 polymorphic sites in the 642 bp fragment, 55 unique haplotypes recorded from the 56 individuals sequenced, and a mean number of pairwise differences between haplotypes of 8.54 ± 4.01. CytB genetic diversity was lower than for the HVR, with only 11 unique haplotypes among the 45 individuals sequenced. We did not detect any significant genetic structure among the colonies (CytB: FST = 0.023, p = 0.170; HVR: FST = 0.019, p = 0.055). The AMOVA analyses for HVR indicated that 98.07 % of the observed genetic variation occurred within colonies, with only 1.93 % of variation among colonies. The CytB result was similar, with 97.65 % of the observed genetic variation within colonies and 2.35 % of variation among colonies. These findings indicate that the East Antarctic samples can be analysed as a single population with respect to demographic history. The lack of genetic structure among East Antarctic colonies is consistent with a previous study of Adélie penguin genetic structure, which found genetic homogeneity among Adélie penguins around the continent based on seven microsatellite loci .
East Antarctic Adélie penguin abundance over the past 22,000 years
There is evidence of a population expansion in East Antarctic Adélie penguins commencing around 14 kya, towards the end of the glacial-interglacial transition (Fig. 2). Interestingly, during the period of 22 kya to 15 kya, which encompassed the LGM and the majority of the subsequent glacial-interglacial transition period, the N ef of East Antarctic Adélie penguins was less than 1000, and based on the 95 % highest posterior density interval (HPD) may have been zero (Fig. 2). It is therefore possible that Adélie penguins were not present in East Antarctica prior to ca. 15 kya. The population size then began to rapidly increase ca. 14 kya (95 % HPD: 11 kya – 19 kya), with a total increase of approximately 135-fold (Fig. 2). This period of population growth is coincident with increasing temperatures and deglaciation in East Antarctica during the glacial-interglacial transition and early Holocene (Fig. 2).
Antarctic wide mtDNA phylogeny
Our phylogenetic analysis of Adélie penguins from East Antarctica, the Ross Sea  and the Scotia Arc , revealed two strongly supported clades (posterior probability = 1; Fig. 3). One of these clades is comprised solely of penguins sampled at Ross Sea colonies, whereas the other contains individuals from colonies in East Antarctica, the Scotia Arc and the Ross Sea. This confirms the pattern found by Ritchie et al. (2004), who noted the presence of one Ross Sea lineage and one “Antarctic” lineage, which comprised individuals from the Ross Sea, Welch and Gardner Islands in East Antarctica, and Torgersen Island on the Antarctic Peninsula. Clucas et al. (2014) compared Scotia Arc penguins to Ross Sea individuals using a haplotype network, and found that individuals from the Scotia Arc also fell into the “Antarctic” lineage. Our study now shows that Adélie penguins from additional locations in East Antarctica, including Pétrels Island, Holl Island, Macey Island, Blakeney Point and Béchervaise Island, also fall into the “Antarctic” lineage.
Our dated phylogeny indicates that both lineages originated during the last glacial period and probably represent two glacial refugia, with the time to the most recent common ancestor (TMRCA) of the Ross Sea lineage estimated at 33 kya (95 % HPD: 22–48 kya), and of the Antarctic lineage estimated at 29 kya (95 % HPD: 22–39 kya).
Within the two lineages there is a high degree of topological uncertainty, suggesting very little genetic structure among Adélie penguin populations, aside from the division into two lineages. This is consistent with our FST and AMOVA results, as well the published findings of genetic homogeneity among Adélie penguins around much of the continent based on microsatellite DNA markers .
Mitochondrial DNA sequence data has revealed that Adélie penguin numbers in East Antarctica increased 135-fold during the most recent glacial-interglacial transition. This dramatic expansion suggests that the East Antarctic environment is currently much more favourable for Adélie penguins than it was prior to 14 kya. Our genetic analyses shows that the East Antarctic population began to expand between 19 and 11 kya (95 % HPD), with a median estimate of 14 kya. This was coincident with a global change in climate regime, from the LGM (26–19.5 kya) to the warm Holocene (11.7 kya–present) . During this transitional period, glaciers and ice sheets retreated [39–41], Southern Ocean primary productivity increased by between two and five fold [42–44], the winter sea ice field is estimated to have halved in areal extent , and sea ice became seasonal, rather than perennial . These changes in the East Antarctic coastal environment would have created more favourable conditions for Adélie penguins, by increasing the amount of ice-free ground suitable for nesting, by increasing prey abundance as marine productivity rose, and by creating more accessible marine foraging grounds and breeding habitat as sea ice concentrations lessened.
Previous studies based on radiocarbon dated remains, found at the Windmill Islands and Vestfold Hills, placed the occupation of East Antarctica by Adélie penguins at ca. 9 kya [35, 36]; our data from a more extensive set of sites across the region now show that East Antarctica was probably colonised at least 2000 years earlier. Given the confidence interval on our demographic reconstructions (Fig. 2), it is unclear whether Adélie penguins persisted in East Antarctica throughout the LGM in small numbers, or if they colonised the region during the post-glacial period from refugia located elsewhere.
Our phylogenetic analyses (Fig. 3) provide support to the notion that only two Adélie penguin refugia existed in Antarctica during the LGM [30, 31]. One refuge was most likely located near the Ross Sea , whereas the location of the other refuge remains unknown, but was clearly the source population for extant colonies in the Scotia Arc , parts of the Ross Sea [30, 31], and the full breadth of East Antarctica (Fig. 3). Given the high genetic differentiation of this lineage from the Ross Sea lineage, it seems unlikely that the second refuge was located in the vicinity of the Ross Sea, directing attention toward the Scotia Arc or East Antarctica as potential refuge locations.
Our results indicate a median expansion time of the East Antarctic population ca. 14 kya, whereas the Scotia Arc population is estimated to have expanded slightly earlier (ca. 17 kya, median estimate; ), coincident with the earlier deglaciation of that region compared to East Antarctica [32, 33, 41]. Given that the Scotia Arc population expanded first, it may be that the LGM refuge was located in this region and, as the lineage expanded and colonised new areas, it made its way to East Antarctica ca. 14 kya. The South Shetland Islands, in the Scotia Arc, had a similar glacial extent during the LGM as they do today , and could therefore have supported a small refuge population of Adélie penguins prior to widespread deglaciation of the Scotia Arc region from 18 kya. Under this scenario, Adélie penguins from a Scotia Arc refuge could have colonised an ice-free area in East Antarctica, such as the Vestfold Hills , Lützow-Holm Bay  or the Larsemann Hills  ca. 14 kya, subsequently becoming more widespread in East Antarctica as the rest of the region deglaciated from 12 kya . Previous genetic analysis of radiocarbon dated remains indicate that both Adélie penguin lineages were present in the Ross Sea from at least ca. 6 kya, suggesting that individuals from the second refuge had spread to the Ross Sea by this time [30, 31].
It is also possible that an Adélie penguin refuge was located in East Antarctica. While most of the East Antarctic coastline was covered by glacial ice until at least 12 kya, some areas have been ice-free since before the LGM [41, 47] and could, theoretically, have acted as a refuge for small Adélie penguin breeding colonies if local sea ice and foraging conditions were favourable. If this were the case, then these sites were likely to be isolated pockets suitable for breeding Adélie penguin populations rather than broad regions of suitable breeding and marine foraging habitat. Previous studies have indicated that ice-free oases during the LGM were located at several sites within East Antarctica, including at the Vestfold Hills (based on lake sediment records; ), at Lützow-Holm Bay (based on Holocene raised beach deposits; ) and at the Larsemann Hills (based on radiocarbon dated moss deposits; ). The Bunger Hills were also partially ice-free during and since the LGM , however, the region is bounded by ice shelves and is therefore inaccessible to penguins. If Adélie penguins were present in East Antarctica at one of these sites before 14 kya, they may have expanded their range and numbers in the region as ice-free habitat became more plentiful coincident with the widespread retreat of the East Antarctic Ice Sheet from ca. 12 kya . While there is evidence that these sites were ice-free and could have theoretically supported breeding Adélie penguins, the existence of Adélie penguin glacial refugia in East Antarctica or the Scotia Arc could be confirmed by radiocarbon dating of penguin remains from these potential refugia sites. The Vestfold Hills, Lützow-Holm Bay and the South Shetland Islands are all home to penguin breeding colonies today and, while no penguins currently breed at the Larsemann Hills, they may have done so in the past.
While ice-free ground suitable for nesting is a key requirement for the existence of Adélie penguin colonies, the species also requires accessible marine foraging grounds and sufficiently abundant prey to survive. Unlike winter during the LGM, when sea ice extent was double the current winter values, the LGM summer sea ice extent is estimated to have been similar to that seen in summer today . Therefore, adult Adélie penguins provisioning for their chicks during the summer breeding season may have encountered similar sea ice extents as they do today. However, the penguins’ capacity to forage successfully would largely depend on the amount of fast ice present, which can impede their ability to reach foraging grounds [13, 19, 51]. The LGM sea ice records refer to total sea ice extent, with no differentiation between fast and pack ice , therefore it is unknown whether the LGM summer fast ice conditions would have been suitable for Adélie penguin foraging. Coastal polynyas could have facilitated Adélie penguin foraging amidst the sea ice field and also acted as hot spots of primary productivity . Polynyas are known to have existed in several locations in the Weddell and Ross Seas during the LGM [52–56]. There are currently no records of LGM polynyas in East Antarctica, but this is more likely a result of the sparse sediment core record rather than an actual absence of polynyas, which are thought to have been more widespread during the LGM due to increased katabatic winds compared to today [53, 56].
Sea ice extent and seasonality began to shift from LGM conditions ca. 10.4 kya in Prydz Bay, with Holocene sea ice conditions similar to today’s reached between 10 and 9 kya . This decline in sea ice occurred after the initial increase in abundance of East Antarctic Adélie penguins (ca. 14 kya; Fig. 2), suggesting that sea ice conditions were not the primary driver of population expansion. Recent studies in East Antarctica have shown that sea ice variation is a key driver of Adélie penguin population dynamics and key demographic parameters over yearly and decadal time scales [13, 18, 20, 57, 58], and sea ice declines are predicted to result in decreasing numbers of Adélie penguins in the most northerly latitudes of their breeding range over the coming decades . It appears, however, that changes in sea ice extent and seasonality during the glacial-interglacial transition were not the key driver of East Antarctic Adélie penguin population expansion. This suggests that environmental drivers of population trends over thousands of years may differ to drivers over years or decades.
Increases in primary productivity did not commence until between 12 and 10 kya off MacRobertson Land [60, 61], and from ca. 10 kya off Adélie Land , indicating that changes in primary productivity were not an initial driver of the Adélie penguin abundance increase either. As Adélie penguins began increasing from 14 kya, this suggests that prey abundance was already sufficiently high for penguin survival prior to the increase in East Antarctic primary productivity between 12 and 10 kya. This is supported by evidence from another meso-predator, the Weddell seal, which persisted in East Antarctica throughout the LGM and post-glacial period in similar numbers to today . For the Weddell seal’s population size in the region to be unchanged, the species must have had sufficiently abundant prey over this period; therefore, it is plausible that prey resources for the Adélie penguin, which breeds at the same time of year as the Weddell seal, were also sufficient prior to 14 kya. The overall indication, based on the timing of the East Antarctic Adélie penguin population increase, is that deglaciation, leading to increased ice-free area for nesting, was the initial environmental variable that changed sufficiently to allow the post-glacial abundance increase. Subsequent changes to sea ice and primary productivity after the initial population expansion may have sustained this trend in unison with increasing breeding habitat availability.
The observed post-glacial expansion of Adélie penguins is common to many penguin species, including emperors , kings , gentoos and chinstraps , which all had populations smaller in size and restricted in range during the LGM, and which expanded post-glacially. Interestingly, the increase in abundance of East Antarctic Adélie penguins began earlier and was far greater than that of the sympatric, closely related emperor penguins , which were also restricted to refugia during the LGM . The East Antarctic emperor penguin population increased in abundance during the Holocene; however, the expansion was only 5.7 fold and did not commence until ca. 10 kya , approximately 4000 years after the Adélie penguin expansion. This suggests that the two species were influenced by different environmental drivers during the post-glacial period. While both species responded positively to declining ice, for the Adélie penguins it is likely that terrestrial ice sheet retreat was the key factor, whereas the expansion of emperor penguins was more closely coupled with reductions in sea ice extent .
Our study has shown that East Antarctic Adélie penguins responded similarly to Scotia Arc Adélie penguins during climate warming following the LGM, with both populations undergoing an increase in population size coincident with an expansion of ice-free breeding habitat. Increases in Adélie penguin numbers in these two regions occurred asynchronously, in line with local deglaciation, indicating that Adélie penguins are sensitive to local glacier and ice sheet retreat. As climate change progresses, glaciers and ice sheets in Antarctica are expected to retreat further. This study highlights the possibility that, in regions where sea ice and prey conditions remain favourable, Adélie penguin numbers may expand in line with this deglaciation as additional breeding sites become exposed. In the Ross Sea, such an expansion has already occurred in response to receding glacial ice at the Beaufort Island colony, with an increase of 84 % in the Adélie penguin population between 1983 and 2010, concurrent with a 543 m retreat of the glacier field . While the future trends in Adélie penguin abundance remain uncertain, our study suggests that ice sheet retreat and the availability of ice-free breeding areas may be overriding factors in determining the millennial scale abundance trends of this species. However, for Adélie penguin populations to expand in line with increasing breeding habitat, prey must be abundant and accessible enough to meet the requirements of the expanding population. Whether this will be the case in the future remains to be seen, as the impacts of climate change on Adélie penguin prey species, such as Antarctic krill (Euphausia superba), are currently ambiguous [8, 66, 67].
Our study also demonstrates that the environmental drivers underlying changes in population size over yearly, decadal and millennial timescales are not necessarily the same. In the case of East Antarctic Adélie penguins, sea ice conditions are a critical factor in population success over decadal and yearly timescales, however, deglaciation appears to have been the key driver of population change over millennia. This finding has important implications for the forecasting of species’ abundance and distribution changes under future climate change scenarios, and highlights the need to consider millennial scale trends alongside contemporary data for understanding species’ responses to climate change.
Specimens of muscle tissue were collected between 2012 and 2014 from the carcasses of recently deceased Adélie penguins at Blakeney Point (n = 13), Holl Island (n = 7), Macey Island (n = 2), Welch Island (n = 13), Béchervaise Island (n = 11) and Pétrels Island (n = 10) (Fig. 1). We refer to these collections as the East Antarctic region samples. Tissue was transported and stored at −20 °C.
Genomic DNA was extracted using the QIAGEN DNeasy Blood and Tissue Kit following the manufacturer’s protocols. The mitochondrial hypervariable region (HVR) and cytochrome B (CytB) were sequenced, as these markers have been successfully used to reconstruct the demographic history of the closely related emperor penguin , and the use of HVR allowed for comparisons with published Adélie penguin datasets from the Ross Sea [30, 31] and Antarctic Peninsula/Scotia Arc .
HVR was amplified and sequenced for all individuals using primers AP1STR (5′-CCACCCTATACATACAATTCCCCTCCC-3′)  and H-A650 (5′-CTGACATAGGAACCAGAGGCGC-3′) [29, 31, 68]. CytB was amplified and sequenced for individuals from Béchervaise Island (n = 11), Macey Island (n = 1), Welch Island (n = 11), Blakeney Point (n = 12), and Pétrels Island (n = 10) using primers CytB-F1 (5′-ACTGCAGACACAACCCTAGC-3′)  and CytB-R1 (5′-GGGAAGAGGATCAGGAGGGT-3′) . For both HVR and CytB, reaction mixes consisted of 7.5 μL of GoTaq Green Master Mix (Promega), 0.2 μM of each primer, and 5–10 ng of gDNA, made up to 15 μL with ddH2O. Annealing temperatures for HVR and CytB PCRs were 52.5 and 60 °C, respectively. Bi-directional Sanger sequencing using the PCR primer pairs was carried out at the Australian Genome Research Facility (AGRF). Occasional heteroplasmic sites were present in the HVR dataset, as expected for Adélie penguin HVR [29, 69], and these were re-scored according to IUPAC ambiguity codes when the secondary peak was >40 % of the height of the primary peak in both forward and reverse sequences.
jModeltest  was used to estimate the optimal nucleotide substitution model for each dataset and Arlequin v3.5  was used to calculate summary statistics for the CytB and HVR datasets and to perform analyses of molecular variance (AMOVA).
East Antarctic demographic reconstructions
The demographic history of East Antarctic Adélie penguins over the past 22,000 years was reconstructed using the coalescent extended Bayesian skyline plot method  within BEAST v2.1.3 . The nucleotide substitution model for both HVR and CytB were specified as HKY  with four gamma categories, which was selected by jModelTest  as the optimal model in both cases. A strict molecular clock was used with lognormal substitution rate priors specified for HVR (mean = 0.55 substitutions/site/Myr, SD = 0.15) to reflect the published substitution rate ; and for CytB (mean = 0.039 substitutions/site/Myr, SD = 1.5) to reflect the substitution rate of CytB in the closely related emperor penguin .
The posterior distribution of effective population size through time was generated using the Markov chain Monte Carlo (MCMC) sampling procedure, which was run for 120 million generations with samples drawn every 5000 steps. Tracer v1.5 was used to visualise the sampling trace and to check the effective sample size values (ESS) to confirm convergence, with most ESS values exceeding 1000, and all values exceeding 200. Three independent BEAST analyses with different random number seeds were performed to ensure reproducibility of the posterior distribution.
The effective population size (N e )  is the number of individuals in an ideal population (i.e. a population with an equal sex ratio, random mating and no variation in reproductive success) that undergoes random genetic change at the same rate as the real population . Because real populations do not usually adhere to these idealised constraints, N e is almost always considerably smaller than the census size of a population (N); a recent review across a wide range of taxonomic groups found a median N e/N ratio of 0.14 . Because mitochondrial DNA was used for this study, our estimate is for female effective population size (N ef ). To convert the population size parameter of the demographic model (N ef *tau) to N ef , we divided the parameter by the generation length of Adélie penguins, which we have estimated at 13.3 years using the formula tau = A + [S/(1-S)], where A = age of first breeding (estimated at 4.85 for females, ) and S = yearly survival rate after reaching breeding age (estimated at 0.894, ). Generation length can be difficult to assess accurately for penguins, as both age at first breeding and annual survival rates may differ by location and be influenced by anomalous environmental conditions. It should therefore be noted that any variance in our estimate of tau would affect the absolute value of N ef in our results, but has no bearing on either the timing or magnitude of the abundance increase reported. It should also be noted that the measure of N ef applies to the total breeding population and is not necessarily constrained to our study area, which may be inside the boundaries of a larger panmictic breeding population.
Antarctic wide phylogeny
In order to determine the phylogenetic placement of East Antarctic Adélie penguins within the global population we constructed a phylogeny for a dataset of HVR sequences of Adélie penguins from East Antarctica (this study), from the Scotia Arc and Antarctic Peninsula (hereinafter referred to as Scotia Arc) , and from the Ross and Somov Seas (hereinafter referred to as Ross Sea) . The Scotia Arc dataset consisted of 36 individuals from the South Shetland Islands, the South Orkney Islands, the South Sandwich Islands, and Lagotellerie Island (Fig. 1) [GenBank PopSet: 634224762; ]. The Ross Sea dataset consisted of 49 individuals, randomly chosen from within GenBank PopSet: 45443792, which contains Adélie penguins sampled from Cape Adare, Port Martin, Adélie Cove, Edmonson Point and the Balleny Islands (Fig. 1) .
Bayesian phylogenetic analyses were performed on this dataset using BEAST v2.1.3 . The nucleotide substitution model and substitution rate prior for HVR were specified as described above for the demographic reconstruction, and the coalescent Bayesian skyline tree prior was used . The posterior distribution of phylogenetic trees was generated using the MCMC sampling procedure, which was run for 200 million iterations with trees logged every 40,000 steps. Tracer v.1.5 was used to check effective sample size (ESS) values to confirm convergence with all values >200. Four independent BEAST analyses, from different random number seeds, were performed to ensure reproducibility of the posterior distribution. The maximum clade credibility tree (after a burn in of 5 %) was selected using TreeAnnotator v2.1.3 and then visualised in FigTree v1.4. The heights specified for each node in the tree are median values.
Availability of supporting data
The nucleic acid sequences supporting the results of this article are available in the GenBank repository, accession numbers: KT932437 - KT932537.
Vaughan DG, Comiso JC, Allison I, Carrasco J, Kaser G, Kwok R, et al. Observations: Cryosphere. In: Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, editors. Climate change 2013: The physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change. Cambridge, UK and New York, USA: Cambridge University Press; 2013. p. 317–82.
Zwally HJ, Comiso JC, Parkinson CL, Cavalieri DJ, Gloersen P. Variability of Antarctic sea ice 1979–1998. J Geophys Res. 2002;107:9–19.
Stammerjohn S, Massom R, Rind D, Martinson D. Regions of rapid sea ice change: An inter-hemispheric seasonal comparison. Geophys Res Lett. 2012;39:L06501.
Massom R, Reid P, Stammerjohn S, Raymond B, Fraser A, Ushio S. Change and variability in East Antarctic sea ice seasonality, 1979/80–2009/10. PLoS One. 2013;8(5):e64756. doi:10.1371/journal.pone.0064756.
Fraser AD, Massom RA, Michael KJ, Galton-Fenzi BK, Lieser JL. East Antarctic landfast sea ice distribution and variability, 2000–08. J Climate. 2012;25:1137–56.
Collins M, Knutti R, Arblaster J, Dufresne J-L, Fichefet T, Friedlingstein P, et al. Long-term climate change: Projections, commitments and irreversibility. In: Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, editors. Climate change 2013: The physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change. Cambridge, UK and New York, USA: Cambridge University Press; 2013. p. 1029–136.
Smetacek V, Nicol S. Polar ocean ecosystems in a changing world. Nature. 2005;437:362–8.
Constable AJ, Melbourne-Thomas J, Corney SP, Arrigo KR, Barbraud C, Barnes DK, et al. Climate change and Southern Ocean ecosystems I: How changes in physical habitats directly affect marine biota. Glob Chang Biol. 2014;20:3004–25.
Lewis SL, Maslin MA. Defining the anthropocene. Nature. 2015;519:171–80.
Trathan P, Ballard G. Adélie penguin. In: Borboroglu PG, Boersma PD, editors. Penguins: Natural history and conservation. Seattle: University of Washington Press; 2013.
Lynch H, LaRue M. First global census of the Adélie penguin. Auk. 2014;131:457–66.
Clarke J, Emmerson LM, Otahal P. Environmental conditions and life history constraints determine foraging range in breeding Adélie penguins. Mar Ecol Prog Ser. 2006;310:247–61.
Emmerson L, Southwell C. Sea ice cover and its influence on Adélie penguin reproductive performance. Ecology. 2008;89:2096–102.
Emmerson L, Pike R, Southwell C. Reproductive consequences of environment-driven variation in Adélie penguin breeding phenology. Mar Ecol Prog Ser. 2011;440:203–16.
Baroni C, Orombelli G. Abandoned penguin rookeries as Holocene paleoclimatic indicators in Antarctica. Geology. 1994;22:23–6.
Forcada J, Trathan P, Reid K, Murphy E, Croxall J. Contrasting population changes in sympatric penguin species in association with climate warming. Glob Chang Biol. 2006;12:411–23.
Southwell C, Emmerson L. Large-scale occupancy surveys in East Antarctica discover new Adélie penguin breeding sites and reveal an expanding breeding distribution. Antarct Sci. 2013;25:531–5.
Kato A, Ropert-Coudert Y, Naito Y. Changes in Adélie penguin breeding populations in Lützow-Holm Bay, Antarctica, in relation to sea-ice conditions. Polar Biol. 2002;25:934–8.
Ropert-Coudert Y, Kato A, Meyer X, Pellé M, MacIntosh AJ, Angelier F, et al. A complete breeding failure in an Adélie penguin colony correlates with unusual and extreme environmental events. Ecography. 2015;38:111–3.
Jenouvrier S, Barbraud C, Weimerskirch H. Sea ice affects the population dynamics of Adélie penguins in Terre Adélie. Polar Biol. 2006;29:413–23.
Tamura T, Williams GD, Fraser AD, Ohshima KI. Potential regime shift in decreased sea ice production after the Mertz Glacier calving. Nat Commun. 2012;3:826.
Dragon A-C, Houssais M-N, Herbaut C, Charrassin J-B. A note on the intraseasonal variability in an Antarctic polynia: Prior to and after the Mertz Glacier calving. J Mar Syst. 2014;130:46–55.
Raymond B, Lea MA, Patterson T, Andrews-Goff V, Sharples R, Charrassin JB, et al. Important marine habitat off East Antarctica revealed by two decades of multi-species predator tracking. Ecography. 2014;38:121–9.
Clark PU, Dyke AS, Shakun JD, Carlson AE, Clark J, Wohlfarth B, et al. The last glacial maximum. Science. 2009;325:710–4.
Jouzel J, Masson-Delmotte V, Cattani O, Dreyfus G, Falourd S, Hoffmann G, et al. Orbital and millennial Antarctic climate variability over the past 800,000 years. Science. 2007;317:793–6.
Millar CD, Subramanian S, Heupink TH, Swaminathan S, Baroni C, Lambert DM. Adélie penguins and temperature changes in Antarctica: a long-term view. Integrative zoology. 2012;7:113–20.
Emslie S, Polito M, Brasso R, Patterson W, Sun L. Ornithogenic soils and the paleoecology of pygoscelid penguins in Antarctica. Quat Int. 2014;352:4–15.
Li C, Zhang Y, Li J, Kong L, Hu H, Pan H, et al. Two Antarctic penguin genomes reveal insights into their evolutionary history and molecular changes related to the Antarctic environment. GigaScience. 2014;3:27.
Clucas GV, Dunn MJ, Dyke G, Emslie SD, Levy H, Naveen R, et al. A reversal of fortunes: Climate change ‘winners’ and ‘losers’ in Antarctic Peninsula penguins. Sci Rep. 2014;4:5024.
Lambert DM, Ritchie PA, Millar CD, Holland B, Drummond AJ, Baroni C. Rates of evolution in ancient DNA from Adélie penguins. Science. 2002;295:2270–3.
Ritchie PA, Millar CD, Gibb GC, Baroni C, Lambert DM. Ancient DNA enables timing of the Pleistocene origin and Holocene expansion of two Adélie penguin lineages in Antarctica. Mol Biol Evol. 2004;21:240–8.
Ó Cofaigh C, Davies BJ, Livingstone SJ, Smith JA, Johnson JS, Hocking EP, et al. Reconstruction of ice-sheet changes in the Antarctic Peninsula since the Last Glacial Maximum. Quat Sci Rev. 2014;100:87–110.
Hodgson DA, Graham AG, Roberts SJ, Bentley MJ, Ó Cofaigh C, Verleyen E, et al. Terrestrial and submarine evidence for the extent and timing of the Last Glacial Maximum and the onset of deglaciation on the maritime-Antarctic and sub-Antarctic islands. Quat Sci Rev. 2014;100:137–58.
Emslie SD, Coats L, Licht K. A 45,000 yr record of Adélie penguins and climate change in the Ross Sea. Antarctica Geology. 2007;35:61–4.
Huang T, Sun L, Wang Y, Zhu R. Penguin occupation in the Vestfold Hills. Antarct Sci. 2009;21:131–4.
Emslie SD, Woehler EJ. A 9000-year record of Adélie penguin occupation and diet in the Windmill Islands, East Antarctica. Antarct Sci. 2005;17:57–66.
Huang T, Sun L, Wang Y, Liu X, Zhu R. Penguin population dynamics for the past 8500 years at Gardner Island, Vestfold Hills. Antarct Sci. 2009;21:571–8.
Roeder AD, Marshall RK, Mitchelson AJ, Visagathilagar T, Ritchie PA, Love DR, et al. Gene flow on the ice: genetic differentiation among Adélie penguin colonies around Antarctica. Mol Ecol. 2001;10:1645–56.
Bentley MJ, Ó Cofaigh C, Anderson JB, Conway H, Davies B, Graham AG, et al. A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum. Quat Sci Rev. 2014;100:1–9.
Barbara L, Crosta X, Massé G, Ther O. Deglacial environments in eastern Prydz Bay, East Antarctica. Quat Sci Rev. 2010;29:2731–40.
Mackintosh AN, Verleyen E, O’Brien PE, White DA, Jones RS, McKay R, et al. Retreat history of the East Antarctic Ice sheet since the last glacial maximum. Quat Sci Rev. 2014;100:10–30.
Denis D, Crosta X, Schmidt S, Carson DS, Ganeshram RS, Renssen H et al. Holocene productivity changes off Adélie Land (East Antarctica). Paleoceanography. 2009; doi:10.1029/2008pa001689.
Kohfeld KE, Le Quéré C, Harrison SP, Anderson RF. Role of marine biology in glacial-interglacial CO2 cycles. Science. 2005;308:74–8.
Hillenbrand C-D, Cortese G. Polar stratification: A critical view from the Southern Ocean. Palaeogeogr Palaeoclimatol Palaeoecol. 2006;242:240–52.
Gersonde R, Crosta X, Abelmann A, Armand L. Sea-surface temperature and sea ice distribution of the southern ocean at the EPILOG last glacial maximum—a circum-antarctic view based on siliceous microfossil records. Quat Sci Rev. 2005;24:869–96.
Crosta X, Denis D, Ther O. Sea ice seasonality during the Holocene, Adélie Land, East Antarctica. Mar Micropaleontol. 2008;66:222–32.
Gibson JA, Paterson KS, White CA, Swadling KM. Evidence for the continued existence of Abraxas Lake, Vestfold Hills, East Antarctica during the Last Glacial Maximum. Antarct Sci. 2009;21:269–78.
Hayashi M, Yoshida Y. Holocene raised beaches in the Lützow-Holm Bay region, East Antarctica. Memoirs of National Institute of Polar Research Special issue. 1994;50:49–84.
Burgess J, Spate A, Shevlin J. The onset of deglaciation in the Larsemann Hills, Eastern Antarctica. Antarct Sci. 1994;6:491–5.
Gore D, Rhodes E, Augustinus P, Leishman M, Colhoun E, Rees-Jones J. Bunger Hills, East Antarctica: Ice free at the last glacial maximum. Geology. 2001;29:1103–6.
Ainley DG. The Adélie penguin: Bellwether of climate change. New York: Columbia University Press; 2002.
Thatje S, Hillenbrand C-D, Mackensen A, Larter R. Life hung by a thread: Endurance of Antarctic fauna in glacial periods. Ecology. 2008;89:682–92.
Smith JA, Hillenbrand C-D, Pudsey CJ, Allen CS, Graham AG. The presence of polynyas in the Weddell Sea during the Last Glacial Period with implications for the reconstruction of sea-ice limits and ice sheet history. Earth Planet Sci Lett. 2010;296:287–98.
Brambati A, Melis R, Quaia T, Salvi G. Late Quaternary climatic changes in the Ross Sea area, Antarctica. In: Gamble JA, Skinner DNB, Henry SA, editors. Antarctica at the close of a millennium. Wellington, New Zealand: Royal Society New Zealand; 2002. p. 359–64.
Mackensen A, Grobe H, Hubberten H-W, Kuhn G. Benthic foraminiferal assemblages and the d13C-signal in the Atlantic sector of the Southern Ocean: Glacial-to-interglacial contrasts. In: Zahn R, Pedersen TF, Kaminski MA, Labeyrie L, editors. Carbon cycling in the glacial ocean: Constraints on the ocean’s role in global change. Berlin, Germany: Springer; 1994. p. 105–44.
Sprenk D, Weber M, Kuhn G, Wennrich V, Hartmann T, Seelos K. Seasonal changes in glacial polynya activity inferred from Weddell Sea varves. Clim Past. 2014;10:1239–51.
Emmerson L, Southwell C. Adélie penguin survival: Age structure, temporal variability and environmental influences. Oecologia. 2011;167:951–65.
Southwell C, Emmerson L, McKinlay J, Newbery K, Takahashi A, Kato A, et al. Spatially extensive standardized surveys reveal widespread, multi-decadal increase in East Antarctic Adélie penguin populations. PLoS ONE. 2015; doi:10.1371/journal.pone.0139877.
Ainley D, Russel J, Jenouvrier S, Woehler E, Lyver PB, Fraser W, et al. Antarctic penguin response to habitat change as Earth’s troposphere reaches 2 °C above preindustrial levels. Ecol Monogr. 2010;80:49–66.
Anderson R, Ali S, Bradtmiller L, Nielsen S, Fleisher M, Anderson B, et al. Wind-driven upwelling in the Southern Ocean and the deglacial rise in atmospheric CO2. Science. 2009;323:1443–8.
Mackintosh A, Golledge N, Domack E, Dunbar R, Leventer A, White D, et al. Retreat of the East Antarctic ice sheet during the last glacial termination. Nat Geosci. 2011;4:195–202.
Younger J. Contrasting responses to a climate regime change by sympatric meso-predators. In: Palaeoecological changes in populations of Antarctic ice-dependent predators and their environmental drivers. Hobart: University of Tasmania PhD; 2015. p. 59–85.
Younger JL, Clucas GV, Kooyman G, Wienecke B, Rogers AD, Trathan PN, et al. Too much of a good thing: Sea ice extent may have forced emperor penguins into refugia during the last glacial maximum. Glob Chang Biol. 2015;21:2215–26.
Trucchi E, Gratton P, Whittington JD, Cristofari R, Le Maho Y, Stenseth NC et al. King penguin demography since the last glaciation inferred from genome-wide data. Proceedings of the Royal Society B: Biological Sciences. 2014; doi:10.1098/rspb.2014.0528.
LaRue MA, Ainley DG, Swanson M, Dugger KM, Lyver POB, Barton K et al. Climate change winners: Receding ice fields facilitate colony expansion and altered dynamics in an Adélie penguin metapopulation. PLoS ONE. 2013; doi:10.1371/journal.pone.0060568.
Nicol S, Brierley AS. Through a glass less darkly—new approaches for studying the distribution, abundance and biology of euphausiids. Deep-Sea Res II Top Stud Oceanogr. 2010;57:496–507.
Atkinson A, Nicol S, Kawaguchi S, Pakhomov E, Quetin L, Ross R, et al. Fitting Euphausia superba into Southern Ocean food-web models: A review of data sources and their limitations. CCAMLR Science. 2012;19:219–45.
Roeder AD, Ritchie PA, Lambert DM. New DNA markers for penguins. Conserv Genet. 2002;3:341–4.
Millar CD, Dodd A, Anderson J, Gibb GC, Ritchie PA, Baroni C, et al. Mutation and evolutionary rates in Adélie penguins from the Antarctic. PLoS Genet. 2008;4:1000209.
Posada D. jModelTest: Phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.
Excoffier L, Lischer HEL. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.
Heled J, Drummond AJ. Bayesian inference of population size history from multiple loci. BMC Evol Biol. 2008;8:289.
Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu C-H, Xie D, et al. BEAST 2: A software platform for Bayesian evolutionary analysis. PLoS Comput Biol. 2014;10(4):e1003537. doi:10.1371/journal.pcbi.1003537.
Hasegawa M, Kishino H, Yano T-A. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985;22:160–74.
Wright S. Evolution in Mendelian populations. Genetics. 1931;16:97.
Crow J, Kimura M. An introduction to population genetics. Minneapolis: Burgess Publishing Co; 1970.
Palstra FP, Ruzzante DE. Genetic estimates of contemporary effective population size: What can they tell us about the importance of genetic stochasticity for wild population persistence? Mol Ecol. 2008;17:3428–47.
Ainley DG, LeReseche RE, Sladen WJ. Breeding biology of the Adelie penguin. Berkeley: University of California Press; 1983.
Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005;22:1185–92.
We are very thankful to the logistics team at the Australian Antarctic Division for enabling our fieldwork in a particularly difficult season. We also gratefully acknowledge Laura Morrissey, Christian Gallagher, Dan Laban, Shane Kern and Steve Hankins for their assistance in the field at Casey; Helen Achurch, Julie McInnes and Cecilia Carrea for collecting specimens at Mawson; and Thierry Raclot for collecting specimens at Pétrels Island. This research was funded by the Australian Antarctic Science Program Project 4184 (KM, LE, JY) and the Sea World Research and Rescue Foundation (KM, JY). Sample collections were made through AAS project #4088 and IPEV project #1091. All genetic sampling was from already deceased penguins, so no live birds were disturbed. Fieldwork in the vicinity of penguin colonies in the Australian Antarctic Territory was approved by the Australian Antarctic Division’s animal ethics committee.
The authors declare that they have no competing interests.
JY carried out the fieldwork, participated in the molecular laboratory work, analysed the data, interpreted the data, drafted the manuscript and participated in conceiving and designing the study. LE coordinated the fieldwork, carried out fieldwork and participated in conceiving and designing the study, and interpreting the data. CS carried out fieldwork and participated in designing the study and interpreting the data. PL carried out the molecular laboratory work and participated in drafting the manuscript and analysing the data. KM conceived of the study, coordinated the study and participated in interpreting the data. All authors read and approved the final manuscript.
About this article
Cite this article
Younger, J., Emmerson, L., Southwell, C. et al. Proliferation of East Antarctic Adélie penguins in response to historical deglaciation. BMC Evol Biol 15, 236 (2015). https://doi.org/10.1186/s12862-015-0502-2