Insights into bear evolution from a Pleistocene polar bear genome.

The polar bear (Ursus maritimus) has become a symbol of the threat to biodiversity from climate change. Understanding polar bear evolutionary history may provide insights into apex carnivore responses and prospects during periods of extreme environmental perturbations. In recent years, genomic studies have examined bear speciation and population history, including evidence for ancient admixture between polar bears and brown bears (Ursus arctos). Here, we extend our earlier studies of a 130,000- to 115,000-y-old polar bear from the Svalbard Archipelago using a 10× coverage genome sequence and 10 new genomes of polar and brown bears from contemporary zones of overlap in northern Alaska. We demonstrate a dramatic decline in effective population size for this ancient polar bear's lineage, followed by a modest increase just before its demise. A slightly higher genetic diversity in the ancient polar bear suggests a severe genetic erosion over a prolonged bottleneck in modern polar bears. Statistical fitting of data to alternative admixture graph scenarios favors at least one ancient introgression event from brown bears into the ancestor of polar bears, possibly dating back over 150,000 y. Gene flow was likely bidirectional, but allelic transfer from brown into polar bear is the strongest detected signal, which contrasts with other published work. These findings may have implications for our understanding of climate change impacts: Polar bears, a specialist Arctic lineage, may not only have undergone severe genetic bottlenecks but also been the recipient of generalist, boreal genetic variants from brown bears during critical phases of Northern Hemisphere glacial oscillations.

Ancient bears provide insights into Pleistocene ice age refugia in Southeast Alaska.

During the Late Pleistocene, major parts of North America were periodically covered by ice sheets. However, there are still questions about whether ice-free refugia were present in the Alexander Archipelago along the Southeast (SE) Alaska coast during the last glacial maximum (LGM). Numerous subfossils have been recovered from caves in SE Alaska, including American black (Ursus americanus) and brown (U. arctos) bears, which today are found in the Alexander Archipelago but are genetically distinct from mainland bear populations. Hence, these bear species offer an ideal system to investigate long-term occupation, potential refugial survival and lineage turnover. Here, we present genetic analyses based on 99 new complete mitochondrial genomes from ancient and modern brown and black bears spanning the last ~45,000 years. Black bears form two SE Alaskan subclades, one preglacial and another postglacial, that diverged >100,000 years ago. All postglacial ancient brown bears are closely related to modern brown bears in the archipelago, while a single preglacial brown bear is found in a distantly related clade. A hiatus in the bear subfossil record around the LGM and the deep split of their pre- and postglacial subclades fail to support a hypothesis of continuous occupancy in SE Alaska throughout the LGM for either species. Our results are consistent with an absence of refugia along the SE Alaska coast, but indicate that vegetation quickly expanded after deglaciation, allowing bears to recolonize the area after a short-lived LGM peak.

Genetic drift drives rapid speciation of an Arctic insular endemic shrew (Sorex pribilofensis).

Episodes of Quaternary environmental change shaped the genomes of extant species, influencing their response to contemporary environments, which are changing rapidly. Island endemics are among the most vulnerable to such change, accounting for a disproportionate number of recent extinctions. To prevent extinctions and conserve island biodiversity it is vital to combine knowledge of species' ecologies with their complex evolutionary histories. The Bering Sea has a history of cyclical island isolation and reconnection, coupled with modern rates of climate change that exceed global averages. The endangered Pribilof Island shrew (Sorex pribilofensis) is endemic to St. Paul Island, Alaska, which was isolated from mainland Beringia ~14,000 years ago by rising sea levels. Using ~11,000 single nucleotide polymorphisms, 17 microsatellites and mitochondrial sequence data, we test predictions about the evolutionary processes driving shrew speciation across Beringia. Our data show considerable differentiation of S. pribilofensis from mainland sibling species, relative to levels of divergence between mainland shrews. We also find a genome-wide loss of diversity and extremely low Ne for S. pribilofensis. We then show that intraspecific genetic diversity is significantly related to interspecific divergence, and that differentiation between S. pribilofensis and other Beringian shrews is highest across loci that are fixed in S. pribilofensis, indicating that strong drift has driven differentiation of this island species. Our findings show that drift as a consequence of Arctic climate cycling can rapidly reshape insular biodiversity. Arctic island species that lack genomic diversity and have evolved in response to past climate may have limited ability to respond to modern environmental changes.

Phylogenomics reveals ancient and contemporary gene flow contributing to the evolutionary history of sea ducks (Tribe Mergini).

Insight into complex evolutionary histories continues to build through broad comparative phylogenomic and population genomic studies. In particular, there is a need to understand the extent and scale that gene flow contributes to standing genomic diversity and the role introgression has played in evolutionary processes such as hybrid speciation. Here, we investigate the evolutionary history of the Mergini tribe (sea ducks) by coupling multi-species comparisons with phylogenomic analyses of thousands of nuclear ddRAD-seq loci, including Z-sex chromosome and autosomal linked loci, and the mitogenome assayed across all extant sea duck species in North America. All sea duck species are strongly structured across all sampled marker types (pair-wise species ΦST > 0.2), with clear genetic assignments of individuals to their respective species, and phylogenetic relationships recapitulate known relationships. Despite strong species integrity, we identify at least 18 putative hybrids; with all but one being late generational backcrosses. Most interesting, we provide the first evidence that an ancestral gene flow event between long-tailed ducks (Clangula hyemalis) and true Eiders (Somateria spp.) not only moved genetic material into the former species, but likely generated a novel species - the Steller's eider (Polysticta stelleri) - via hybrid speciation. Despite generally low contemporary levels of gene flow, we conclude that hybridization has and continues to be an important process that shifts novel genetic variation between species within the tribe Mergini. Finally, we outline methods that permit researchers to contrast genomic patterns of contemporary versus ancestral gene flow when attempting to reconstruct potentially complex evolutionary histories.

Isolation and characterization of microsatellite loci in merlins (Falco columbarius) and cross-species amplification in gyrfalcons (F. rusticolus) and peregrine falcons (F. peregrinus).

Merlins, Falco columbarius, breed throughout temperate and high latitude habitats in Asia, Europe, and North America. Like peregrine falcons, F. peregrinus, merlins underwent population declines during the mid-to-late twentieth century, due to organochlorine-based contamination, and have subsequently recovered, at least in North American populations. To better understand levels of genetic diversity and population structuring in contemporary populations and to assess the impact of the twentieth century decline, we used genomic data archived in public databases and constructed genomic libraries to isolate and characterize a suite of 17 microsatellite markers for use in merlins. We also conducted cross-amplification experiments to determine the markers' utility in peregrine falcons and gyrfalcons, F. rusticolus. These markers provide a valuable addition to marker suites that can be used to determine individual identity and conduct genetic analyses on merlins and congeners.

Development and characterization of 12 polymorphic microsatellite loci in the sea sandwort, Honckenya peploides.

Codominant marker systems are better suited to analyze population structure and assess the source of an individual in admixture analyses. Currently, there is no codominant marker system using microsatellites developed for the sea sandwort, Honckenya peploides (L.) Ehrh., an early colonizer in island systems. We developed and characterized novel microsatellite loci from H. peploides, using reads collected from whole genome shotgun sequencing on a 454 platform. The combined output from two shotgun runs yielded a total of 62,669 reads, from which 58 loci were screened. We identified 12 polymorphic loci that amplified reliably and exhibited disomic inheritance. Microsatellite data were collected and characterized for the 12 polymorphic loci in two Alaskan populations of H. peploides: Fossil Beach, Kodiak Island (n = 32) and Egg Bay, Atka Island (n = 29). The Atka population exhibited a slightly higher average number of alleles (3.9) and observed heterozygosity (0.483) than the Kodiak population (3.3 and 0.347, respectively). The overall probability of identity values for both populations was PID = 2.892e-6 and PIDsib = 3.361e-3. We also screened the 12 polymorphic loci in Wilhelmsia physodes (Fisch. ex Ser.) McNeill, the most closely related species to H. peploides, and only one locus was polymorphic. These microsatellite markers will allow future investigations into population genetic and colonization patterns of the beach dune ruderal H. peploides on new and recently disturbed islands.

Organic matter quantity and source affects microbial community structure and function following volcanic eruption on Kasatochi Island, Alaska.

In August 2008, Kasatochi volcano erupted and buried a small island in pyroclastic deposits and fine ash; since then, microbes, plants and birds have begun to re-colonize the initially sterile surface. Five years post-eruption, bacterial 16S rRNA gene and fungal internal transcribed spacer (ITS) copy numbers and extracellular enzyme activity (EEA) potentials were one to two orders of magnitude greater in pyroclastic materials with organic matter (OM) inputs relative to those without, despite minimal accumulation of OM (< 0.2%C). When normalized by OM levels, post-eruptive surfaces with OM inputs had the highest ß-glucosidase, phosphatase, NAGase and cellobiohydrolase activities, and had microbial population sizes approaching those in reference soils. In contrast, the strongest factor determining bacterial community composition was the dominance of plants versus birds as OM input vectors. Although soil pH ranged from 3.9 to 7.0, and %C ranged 100×, differentiation between plant- and bird-associated microbial communities suggested that cell dispersal or nutrient availability are more likely drivers of assembly than pH or OM content. This study exemplifies the complex relationship between microbial cell dispersal, soil geochemistry, and microbial structure and function; and illustrates the potential for soil microbiota to be resilient to disturbance.

Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change.

Polar bears (PBs) are superbly adapted to the extreme Arctic environment and have become emblematic of the threat to biodiversity from global climate change. Their divergence from the lower-latitude brown bear provides a textbook example of rapid evolution of distinct phenotypes. However, limited mitochondrial and nuclear DNA evidence conflicts in the timing of PB origin as well as placement of the species within versus sister to the brown bear lineage. We gathered extensive genomic sequence data from contemporary polar, brown, and American black bear samples, in addition to a 130,000- to 110,000-y old PB, to examine this problem from a genome-wide perspective. Nuclear DNA markers reflect a species tree consistent with expectation, showing polar and brown bears to be sister species. However, for the enigmatic brown bears native to Alaska's Alexander Archipelago, we estimate that not only their mitochondrial genome, but also 5-10% of their nuclear genome, is most closely related to PBs, indicating ancient admixture between the two species. Explicit admixture analyses are consistent with ancient splits among PBs, brown bears and black bears that were later followed by occasional admixture. We also provide paleodemographic estimates that suggest bear evolution has tracked key climate events, and that PB in particular experienced a prolonged and dramatic decline in its effective population size during the last ca. 500,000 years. We demonstrate that brown bears and PBs have had sufficiently independent evolutionary histories over the last 4-5 million years to leave imprints in the PB nuclear genome that likely are associated with ecological adaptation to the Arctic environment.

Where east meets west: Phylogeography of the high Arctic North American brant goose.

Genetic variation in Arctic species is often influenced by vicariance during the Pleistocene, as ice sheets fragmented the landscape and displaced populations to low- and high-latitude refugia. The formation of secondary contact or suture zones during periods of ice sheet retraction has important consequences on genetic diversity by facilitating genetic connectivity between formerly isolated populations. Brant geese (Branta bernicla) are a maritime migratory waterfowl (Anseriformes) species that almost exclusively uses coastal habitats. Within North America, brant geese are characterized by two phenotypically distinct subspecies that utilize disjunct breeding and wintering areas in the northern Pacific and Atlantic. In the Western High Arctic of Canada, brant geese consist of individuals with an intermediate phenotype that are rarely observed nesting outside this region. We examined the genetic structure of brant geese populations from each subspecies and areas consisting of intermediate phenotypes using mitochondrial DNA (mtDNA) control region sequence data and microsatellite loci. We found a strong east-west partition in both marker types consistent with refugial populations. Within subspecies, structure was also observed at mtDNA while microsatellite data suggested the presence of only two distinct genetic clusters. The Western High Arctic (WHA) appears to be a secondary contact zone for both Atlantic and Pacific lineages as mtDNA and nuclear genotypes were assigned to both subspecies, and admixed individuals were observed in this region. The mtDNA sequence data outside WHA suggests no or very restricted intermixing between Atlantic and Pacific wintering populations which is consistent with published banding and telemetry data. Our study indicates that, although brant geese in the WHA are not a genetically distinct lineage, this region may act as a reservoir of genetic diversity and may be an area of high conservation value given the potential of low reproductive output in this species.

Powassan virus in mammals, Alaska and New Mexico, U.S.A., and Russia, 2004-2007.

Powassan virus is endemic to the United States, Canada, and the Russian Far East. We report serologic evidence of circulation of this virus in Alaska, New Mexico, and Siberia. These data support further studies of viral ecology in rapidly changing Arctic environments.

Complete mitochondrial genome of a Pleistocene jawbone unveils the origin of polar bear.

The polar bear has become the flagship species in the climate-change discussion. However, little is known about how past climate impacted its evolution and persistence, given an extremely poor fossil record. Although it is undisputed from analyses of mitochondrial (mt) DNA that polar bears constitute a lineage within the genetic diversity of brown bears, timing estimates of their divergence have differed considerably. Using next-generation sequencing technology, we have generated a complete, high-quality mt genome from a stratigraphically validated 130,000- to 110,000-year-old polar bear jawbone. In addition, six mt genomes were generated of extant polar bears from Alaska and brown bears from the Admiralty and Baranof islands of the Alexander Archipelago of southeastern Alaska and Kodiak Island. We show that the phylogenetic position of the ancient polar bear lies almost directly at the branching point between polar bears and brown bears, elucidating a unique morphologically and molecularly documented fossil link between living mammal species. Molecular dating and stable isotope analyses also show that by very early in their evolutionary history, polar bears were already inhabitants of the Artic sea ice and had adapted very rapidly to their current and unique ecology at the top of the Arctic marine food chain. As such, polar bears provide an excellent example of evolutionary opportunism within a widespread mammalian lineage.

A climate for speciation: rapid spatial diversification within the Sorex cinereus complex of shrews.

The cyclic climate regime of the late Quaternary caused dramatic environmental change at high latitudes. Although these events may have been brief in periodicity from an evolutionary standpoint, multiple episodes of allopatry and divergence have been implicated in rapid radiations of a number of organisms. Shrews of the Sorex cinereus complex have long challenged taxonomists due to similar morphology and parapatric geographic ranges. Here, multi-locus phylogenetic and demographic assessments using a coalescent framework were combined to investigate spatiotemporal evolution of 13 nominal species with a widespread distribution throughout North America and across Beringia into Siberia. For these species, we first test a hypothesis of recent differentiation in response to Pleistocene climate versus more ancient divergence that would coincide with pre-Pleistocene perturbations. We then investigate the processes driving diversification over multiple continents. Our genetic analyses highlight novel diversity within these morphologically conserved mammals and clarify relationships between geographic distribution and evolutionary history. Demography within and among species indicates both regional stability and rapid expansion. Ancestral ecological differentiation coincident with early cladogenesis within the complex enabled alternating and repeated episodes of allopatry and expansion where successive glacial and interglacial phases each promoted divergence. The Sorex cinereus complex constitutes a valuable model for future comparative assessments of evolution in response to cyclic environmental change.

Persistence and diversification of the Holarctic shrew, Sorex tundrensis (Family Soricidae), in response to climate change.

Environmental processes govern demography, species movements, community turnover and diversification and yet in many respects these dynamics are still poorly understood at high latitudes. We investigate the combined effects of climate change and geography through time for a widespread Holarctic shrew, Sorex tundrensis. We include a comprehensive suite of closely related outgroup taxa and three independent loci to explore phylogeographic structure and historical demography. We then explore the implications of these findings for other members of boreal communities. The tundra shrew and its sister species, the Tien Shan shrew (Sorex asper), exhibit strong geographic population structure across Siberia and into Beringia illustrating local centres of endemism that correspond to Late Pleistocene refugia. Ecological niche predictions for both current and historical distributions indicate a model of persistence through time despite dramatic climate change. Species tree estimation under a coalescent process suggests that isolation between populations has been maintained across timeframes deeper than the periodicity of Pleistocene glacial cycling. That some species such as the tundra shrew have a history of persistence largely independent of changing climate, whereas other boreal species shifted their ranges in response to climate change, highlights the dynamic processes of community assembly at high latitudes.

Adverse foraging conditions may impact body mass and survival of a high Arctic seabird.

Tradeoffs between current reproduction and future survival are widely recognized, but may only occur when food is limited: when foraging conditions are favorable, parents may be able to reproduce without compromising their own survival. We investigated these tradeoffs in the little auk (Alle alle), a small seabird with a single-egg clutch. During 2005-2007, we examined the relationship between body mass and survival of birds breeding under contrasting foraging conditions at two Arctic colonies. We used corticosterone levels of breeding adults as a physiological indicator of the foraging conditions they encountered during each reproductive season. We found that when foraging conditions were relatively poor (as reflected in elevated levels of corticosterone), parents ended the reproductive season with low body mass and suffered increased post-breeding mortality. A positive relationship between body mass and post-breeding survival was found in one study year; light birds incurred higher survival costs than heavy birds. The results of this study suggest that reproducing under poor foraging conditions may affect the post-breeding survival of long-lived little auks. They also have important demographic implications because even a small change in adult survival may have a large effect on populations of long-lived species.

Using next generation sequencing of alpine plants to improve fecal metabarcoding diet analysis for Dall's sheep.

OBJECTIVES: Dall's sheep (Ovis dalli dalli) are important herbivores in the mountainous ecosystems of northwestern North America, and recent declines in some populations have sparked concern. Our aim was to improve capabilities for fecal metabarcoding diet analysis of Dall's sheep and other herbivores by contributing new sequence data for arctic and alpine plants. This expanded reference library will provide critical reference sequence data that will facilitate metabarcoding diet analysis of Dall's sheep and thus improve understanding of plant-animal interactions in a region undergoing rapid climate change. DATA DESCRIPTION: We provide sequences for the chloroplast rbcL gene of 16 arctic-alpine vascular plant species that are known to comprise the diet of Dall's sheep. These sequences contribute to a growing reference library that can be used in diet studies of arctic herbivores.

Population structure and plumage polymorphism: The intraspecific evolutionary relationships of a polymorphic raptor, Buteo jamaicensis harlani.

BACKGROUND: Phenotypic and molecular genetic data often provide conflicting patterns of intraspecific relationships confounding phylogenetic inference, particularly among birds where a variety of environmental factors may influence plumage characters. Among diurnal raptors, the taxonomic relationship of Buteo jamaicensis harlani to other B. jamaicensis subspecies has been long debated because of the polytypic nature of the plumage characteristics used in subspecies or species designations. RESULTS: To address the evolutionary relationships within this group, we used data from 17 nuclear microsatellite loci, 430 base pairs of the mitochondrial control region, and 829 base pairs of the melanocortin 1 receptor (Mc1r) to investigate molecular genetic differentiation among three B. jamaicensis subspecies (B. j. borealis, B. j. calurus, B. j. harlani). Bayesian clustering analyses of nuclear microsatellite loci showed no significant differences between B. j. harlani and B. j. borealis. Differences observed between B. j. harlani and B. j. borealis in mitochondrial and microsatellite data were equivalent to those found between morphologically similar subspecies, B. j. borealis and B. j. calurus, and estimates of migration rates among all three subspecies were high. No consistent differences were observed in Mc1r data between B. j. harlani and other B. jamaicensis subspecies or between light and dark color morphs within B. j. calurus, suggesting that Mc1r does not play a significant role in B. jamaicensis melanism. CONCLUSIONS: These data suggest recent interbreeding and gene flow between B. j. harlani and the other B. jamaicensis subspecies examined, providing no support for the historical designation of B. j. harlani as a distinct species.

Conservation Genomics in a Changing Arctic.

Although logistically challenging to study, the Arctic is a bellwether for global change and is becoming a model for questions pertinent to the persistence of biodiversity. Disruption of Arctic ecosystems is accelerating, with impacts ranging from mixing of biotic communities to individual behavioral responses. Understanding these changes is crucial for conservation and sustainable economic development. Genomic approaches are providing transformative insights into biotic responses to environmental change, but have seen limited application in the Arctic due to a series of limitations. To meet the promise of genome analyses, we urge rigorous development of biorepositories from high latitudes to provide essential libraries to improve the conservation, monitoring, and management of Arctic ecosystems through genomic approaches.

Identifying reliable indicators of fitness in polar bears.

Animal structural body size and condition are often measured to evaluate individual health, identify responses to environmental change and food availability, and relate food availability to effects on reproduction and survival. A variety of condition metrics have been developed but relationships between these metrics and vital rates are rarely validated. Identifying an optimal approach to estimate the body condition of polar bears is needed to improve monitoring of their response to decline in sea ice habitat. Therefore, we examined relationships between several commonly used condition indices (CI), body mass, and size with female reproductive success and cub survival among polar bears (Ursus maritimus) measured in two subpopulations over three decades. To improve measurement and application of morphometrics and CIs, we also examined whether CIs are independent of age and structural size-an important assumption for monitoring temporal trends-and factors affecting measurement precision and accuracy. Maternal CIs and mass measured the fall prior to denning were related to cub production. Similarly, maternal CIs, mass, and length were related to the mass of cubs or yearlings that accompanied her. However, maternal body mass, but not CIs, measured in the spring was related to cub production and only maternal mass and length were related to the probability of cub survival. These results suggest that CIs may not be better indicators of fitness than body mass in part because CIs remove variation associated with body size that is important in affecting fitness. Further, CIs exhibited variable relationships with age for growing bears and were lower for longer bears despite body length being related to cub survival and female reproductive success. These results are consistent with findings from other species indicating that body mass is a useful metric to link environmental conditions and population dynamics.

Evidence that pairing with genetically similar mates is maladaptive in a monogamous bird.

BACKGROUND: Evidence of multiple genetic criteria of mate choice is accumulating in numerous taxa. In many species, females have been shown to pair with genetically dissimilar mates or with extra-pair partners that are more genetically compatible than their social mates, thereby increasing their offsprings' heterozygosity which often correlates with offspring fitness. While most studies have focused on genetically promiscuous species, few studies have addressed genetically monogamous species, in which mate choice tends to be mutual. RESULTS: Here, we used microsatellite markers to assess individual global heterozygosity and genetic similarity of pairs in a socially and genetically monogamous seabird, the black-legged kittiwake Rissa tridactyla. We found that pairs were more genetically dissimilar than expected by chance. We also identified fitness costs of breeding with genetically similar partners: (i) genetic similarity of pairs was negatively correlated with the number of chicks hatched, and (ii) offspring heterozygosity was positively correlated with growth rate and survival. CONCLUSION: These findings provide evidence that breeders in a genetically monogamous species may avoid the fitness costs of reproducing with a genetically similar mate. In such species that lack the opportunity to obtain extra-pair fertilizations, mate choice may therefore be under high selective pressure.