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.

Whole-genome survey reveals extensive variation in genetic diversity and inbreeding levels among peregrine falcon subspecies.

In efforts to prevent extinction, resource managers are often tasked with increasing genetic diversity in a population of concern to prevent inbreeding depression or improve adaptive potential in a changing environment. The assumption that all small populations require measures to increase their genetic diversity may be unwarranted, and limited resources for conservation may be better utilized elsewhere. We test this assumption in a case study focused on the peregrine falcon (Falco peregrinus), a cosmopolitan circumpolar species with 19 named subspecies. We used whole-genome resequencing to generate over two million single nucleotide polymorphisms (SNPs) from multiple individuals of all peregrine falcon subspecies. Our analyses revealed extensive variation among subspecies, with many island-restricted and nonmigratory populations possessing lower overall genomic diversity, elevated inbreeding coefficients (F ROH)-among the highest reported, and extensive runs of homozygosity (ROH) compared to mainland and migratory populations. Similarly, the majority of subspecies that are either nonmigratory or restricted to islands show a much longer history of low effective population size (N e). While mutational load analyses indicated an increased proportion of homozygous-derived deleterious variants (i.e., drift load) among nonmigrant and island populations compared to those that are migrant or reside on the mainland, no significant differences in the proportion of heterozygous deleterious variants (i.e., inbreeding load) was observed. Our results provide evidence that high levels of inbreeding may not be an existential threat for some populations or taxa. Additional factors such as the timing and severity of population declines are important to consider in management decisions about extinction potential.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina).

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

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.