Genomic evidence of geographically widespread effect of gene flow from polar bears into brown bears

Polar bears are an arctic, marine adapted species that is closely related to brown bears. Genome analyses have shown that polar bears are distinct and genetically homogeneous in comparison to brown bears. However, these analyses have also revealed a remarkable episode of polar bear gene flow into the population of brown bears that colonized the Admiralty, Baranof and Chichagof islands (ABC islands) of Alaska. Here, we present an analysis of data from a large panel of polar bear and brown bear genomes that includes brown bears from the ABC islands, the Alaskan mainland and Europe. Our results provide clear evidence that gene flow between the two species had a geographically wide impact, with polar bear DNA found within the genomes of brown bears living both on the ABC islands and in the Alaskan mainland. Intriguingly, while brown bear genomes contain up to 8.8% polar bear ancestry, polar bear genomes appear to be devoid of brown bear ancestry, suggesting the presence of a barrier to gene flow in that direction.     

Remote biopsy darting and marking of polar bears

Remote biopsy darting of polar bears (Ursus maritimus) is less invasive and time intensive than physical capture and is therefore useful when capture is challenging or unsafe. We worked with two manufacturers to develop a combination biopsy and marking dart for use on polar bears. We had an 80% success rate of collecting a tissue sample with a single biopsy dart and collected tissue samples from 143 polar bears on land, in water, and on sea ice. Dye marks ensured that 96% of the bears were not resampled during the same sampling period, and we recovered 96% of the darts fired. Biopsy heads with 5 mm diameters collected an average of 0.12 g of fur, tissue, and subcutaneous adipose tissue, while biopsy heads with 7 mm diameters collected an average of 0.32 g. Tissue samples were 99.3% successful (142 of 143 samples) in providing a genetic and sex identification of individuals. We had a 64% success rate collecting adipose tissue and we successfully examined fatty acid signatures in all adipose samples. Adipose lipid content values were lower compared to values from immobilized or harvested polar bears, indicating that our method was not suitable for quantifying adipose lipid content.     

Genetic turnovers and northern survival during the last glacial maximum in European brown bears

The current phylogeographic pattern of European brown bears (Ursus arctos) has commonly been explained by postglacial recolonization out of geographically distinct refugia in southern Europe, a pattern well in accordance with the expansion/contraction model. Studies of ancient DNA from brown bear remains have questioned this pattern, but have failed to explain the glacial distribution of mitochondrial brown bear clades and their subsequent expansion across the European continent. We here present 136 new mitochondrial sequences generated from 346 remains from Europe, ranging in age between the Late Pleistocene and historical times. The genetic data show a high Late Pleistocene diversity across the continent and challenge the strict confinement of bears to traditional southern refugia during the last glacial maximum (LGM). The mitochondrial data further suggest a genetic turnover just before this time, as well as a steep demographic decline starting in the mid‐Holocene. Levels of stable nitrogen isotopes from the remains confirm a previously proposed shift toward increasing herbivory around the LGM in Europe. Overall, these results suggest that in addition to climate, anthropogenic impact and inter‐specific competition may have had more important effects on the brown bear's ecology, demography, and genetic structure than previously thought.     

Canadian polar bear population structure using genome‐wide markers

Predicting the consequences of environmental changes, including human‐mediated climate change on species, requires that we quantify range‐wide patterns of genetic diversity and identify the ecological, environmental, and historical factors that have contributed to it. Here, we generate baseline data on polar bear population structure across most Canadian subpopulations (n = 358) using 13,488 genome‐wide single nucleotide polymorphisms (SNPs) identified with double‐digest restriction site‐associated DNA sequencing (ddRAD). Our ddRAD dataset showed three genetic clusters in the sampled Canadian range, congruent with previous studies based on microsatellites across the same regions; however, due to a lack of sampling in Norwegian Bay, we were unable to confirm the existence of a unique cluster in that subpopulation. These data on the genetic structure of polar bears using SNPs provide a detailed baseline against which future shifts in population structure can be assessed, and opportunities to develop new noninvasive tools for monitoring polar bears across their range.     

Evolutionary and ecological traps for brown bears Ursus arctos in human‐modified landscapes

相似文献   

Landscape models for nuclear genetic diversity and genetic structure in white-footed mice (Peromyscus leucopus)

Dramatic changes in the North American landscape over the last 12 000 years haveshaped the genomes of the small mammals, such as the white-footed mouse (Peromyscusleucopus), which currently inhabit the region. However, very recent interactions ofpopulations with each other and the environment are expected to leave the most pronouncedsignature on rapidly evolving nuclear microsatellite loci. We analyzed landscapecharacteristics and microsatellite markers of P. leucopus populations along atransect from southern Ohio to northern Michigan, in order to evaluate hypotheses aboutthe spatial distribution of genetic heterogeneity. Genetic diversity increased to thenorth and was best approximated by a single-variable model based on habitat availabilitywithin a 0.5-km radius of trapping sites. Interpopulation differentiation measured byclustering analysis was highly variable and not significantly related to latitude orhabitat availability. Interpopulation differentiation measured as F_STvalues and chord distance was correlated with the proportion of habitat intervening, butwas best explained by agricultural distance and by latitude. The observed gradients indiversity and interpopulation differentiation were consistent with recent habitatavailability being the major constraint on effective population size in this system, andcontradicted the predictions of both the postglacial expansion and core-peripheryhypotheses.     

Thirty‐five years of spatial autocorrelation analysis in population genetics: an essay in honour of Robert Sokal (1926–2012)

We could not start this review, literally from the beginning, without expressing our sadness over the passing of Professor Robert R. Sokal. We are sure, nevertheless, that the importance of his scientific achievements will ensure he is long remembered. In this modest tribute to Professor Sokal, we highlight his contributions to the field of population genetics and spatial statistical methods. Specifically, we discuss how two papers, co‐authored with Professor N. L. Oden and published in the pages of the Biological Journal of the Linnean Society in 1978, revolutionized the field of analytical population genetics. In these papers, Sokal and Oden created an elegant framework for inferring evolutionary processes (e.g. isolation‐by‐distance, demic diffusion, selection gradients, genetic drift) from the spatial autocorrelation analysis of genetic variation patterns. We also highlight the pivotal importance of Sokal's work to the development of emerging fields (e.g. landscape and conservation genetics). We hope this virtual issue containing the papers that Professor Sokal published in BJLS, and later, related papers by other researchers, will help to remember his work and maintain his legacy of spatial analysis in genetics, ecology, and evolutionary biology. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ??, ??–??.     

Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Speciation occurs when populations diverge and become reproductively isolated from each other. Natural selection is commonly accepted to play a large role in this process, and it has been widely assumed that reproductive isolation often results as a by‐product of divergence driven by adaptation in allopatry. When such populations come into secondary contact, reinforcement can act to strengthen reproductive isolation, but the frequency and importance of this process are still unknown. Here, we explored genomic signatures of selection in allopatry and sympatry for loci associated with reproductive isolation using a natural primate hybrid zone. By analysing reduced‐representation sequencing data, we quantified admixture and population structure across a howler monkey hybrid zone and examined the relationship between locus‐specific differentiation and introgression. We detected extensive admixture that was mostly limited to the narrow contact zone. Loci with reduced introgression into the heterospecific genomic background (the pattern expected for loci associated with reproductive isolation due to selection against hybrids) were significantly more differentiated between allopatric parental populations than loci with neutral and increased introgression, supporting the hypothesis that reproductive isolation is a by‐product of divergence in allopatry. Further, loci with reduced introgression showed greater differentiation in sympatry than in allopatry, suggesting a role for reinforcement. Thus, our results reflect multiple forms of selection that have shaped reproductive isolation in this system. We conclude that reproductive isolation may have initially been driven by divergence in allopatry, but later reinforced by divergent selection in sympatry.     

Mitogenetic structure of brown bears (Ursus arctos L.) in northeastern Europe and a new time frame for the formation of European brown bear lineages

We estimated the phylogenetic relationships of brown bear maternal haplotypes from countries of northeastern Europe (Estonia, Finland and European Russia), using sequences of mitochondrial DNA (mtDNA) control region of 231 bears. Twenty-five mtDNA haplotypes were identified. The brown bear population in northeastern Europe can be divided into three haplogroups: one with bears from all three countries, one with bears from Finland and Russia, and the third composed almost exclusively of bears from European Russia. Four haplotypes from Finland and European Russia matched exactly with haplotypes from Slovakia, suggesting the significance of the current territory of Slovakia in ancient demographic processes of brown bears. Based on the results of this study and those from the recent literature, we hypothesize that the West Carpathian Mountains have served either as one of the northernmost refuge areas or as an important movement corridor for brown bears of the Eastern lineage towards northern Europe during or after the last ice age. Bayesian analyses were performed to investigate the temporal framework of brown bear lineages in Europe. The molecular clock was calibrated using Beringian brown bear sequences derived from radiocarbon-dated ancient samples, and the estimated mutation rate was 29.8% (13.3%-47.6%) per million years. The whole European population and Western and Eastern lineages formed about 175,000, 70,000 and 25,000 years before present, respectively. Our approach to estimating the time frame of brown bear evolution demonstrates the importance of using an appropriate mutation rate, and this has implications for other studies of Pleistocene populations.     

Diversification in a biodiversity hot spot: landscape correlates of phylogeographic patterns in the African spotted reed frog

The Eastern Afromontane Biodiversity Hotspot is known for microendemism and exceptional population genetic structure. The region's landscape heterogeneity is thought to limit gene flow between fragmented populations and create opportunities for regional adaptation, but the processes involved are poorly understood. Using a combination of phylogeographic analyses and circuit theory, I investigate how characteristics of landscape heterogeneity including regional distributions of slope, rivers and streams, habitat and hydrological basins (drainages) impact genetic distance among populations of the endemic spotted reed frog (Hyperolius substriatus), identifying corridors of connectivity as well as barriers to dispersal. Results show that genetic distance among populations is most strongly correlated to regional and local hydrologic structure and the distribution of suitable habitat corridors, not isolation by distance. Contrary to expectations, phylogeographic structure is not coincident with the two montane systems, but instead corresponds to the split between the region's two major hydrological basins (Zambezi and East Central Coastal). This results in a paraphyletic relationship for the Malawian Highlands populations with respect to the Eastern Arc Mountains and implies that the northern Malawian Highlands are the diversity centre for H. substriatus. Although the Malawian Highlands collectively hold the greatest genetic diversity, individual populations have lower diversity than their Eastern Arc counterparts, with an overall pattern of decreasing population diversity from north to south. Through the study of intraspecific differentiation across a mosaic of ecosystem and geographic heterogeneity, we gain insight into the processes of diversification and a broader understanding of the role of landscape in evolution.     

Comparative analyses of plastid and AFLP data suggest different colonization history and asymmetric hybridization between Betula pubescens and B. nana

Birches (Betula spp.) hybridize readily, confounding genetic signatures of refugial isolation and postglacial migration. We aimed to distinguish hybridization from range‐shift processes in the two widespread and cold‐adapted species B. nana and B. pubescens, previously shown to share a similarly east–west‐structured variation in plastid DNA (pDNA). We sampled the two species throughout their ranges and included reference samples of five other Betula species and putative hybrids. We analysed 901 individual plants using mainly nuclear high‐resolution markers (amplified fragment length polymorphisms; AFLPs); a subset of 64 plants was also sequenced for two pDNA regions. Whereas the pDNA variation as expected was largely shared between B. nana and B. pubescens, the two species were distinctly differentiated at AFLP loci. In B. nana, both the AFLP and pDNA results corroborated the former pDNA‐based hypothesis that it expanded from at least two major refugia in Eurasia, one south of and one east of the North European ice sheets. In contrast, B. pubescens showed a striking lack of geographic structuring of its AFLP variation. We identified a weak but significant increase in nuclear (AFLP) gene flow from B. nana into B. pubescens with increasing latitude, suggesting hybridization has been most frequent at the postglacial expansion front of B. pubescens and that hybrids mainly backcrossed to B. pubescens. Incongruence between pDNA and AFLP variation in B. pubescens can be explained by efficient expansion from a single large refugium combined with leading‐edge hybridization and plastid capture from B. nana during colonization of new territory already occupied by this more cold‐tolerant species.     

Integrative approaches to guide conservation decisions: Using genomics to define conservation units and functional corridors

Climate change and increasing habitat loss greatly impact species survival, requiring range shifts, phenotypic plasticity and/or evolutionary change for long‐term persistence, which may not readily occur unaided in threatened species. Therefore, defining conservation actions requires a detailed assessment of evolutionary factors. Existing genetic diversity needs to be thoroughly evaluated and spatially mapped to define conservation units (CUs) in an evolutionary context, and we address that here. We also propose a multidisciplinary approach to determine corridors and functional connectivity between CUs by including genetic diversity in the modelling while controlling for isolation by distance and phylogeographic history. We evaluate our approach on a Near Threatened Iberian endemic rodent by analysing genotyping‐by‐sequencing (GBS) genomic data from 107 Cabrera voles (Microtus cabrerae), screening the entire species distribution to define categories of CUs and their connectivity: We defined six management units (MUs) which can be grouped into four evolutionarily significant units (ESUs) and three (putatively) adaptive units (AUs). We demonstrate that the three different categories of CU can be objectively defined using genomic data, and their characteristics and connectivity can inform conservation decision‐making. In particular, we show that connectivity of the Cabrera vole is very limited in eastern Iberia and that the pre‐Pyrenean and part of the Betic geographic nuclei contribute the most to the species genetic diversity. We argue that a multidisciplinary framework for CU definition is essential and that this framework needs a strong evolutionary basis.     

Chloroplast DNA diversity associated with protected slopes and valleys for hybridizing Eucalyptus species on isolated ranges in south‐eastern Australia

Aim To relate genetic diversity to topographic features and to investigate genetic interactions between Eucalyptus species in a local centre of endemism and diversity in south‐eastern Australia. Location Grampian Ranges, Victoria, Australia. Methods We documented chloroplast DNA (cpDNA) variation for a group of endemic Eucalyptus species (E. serraensis, E. verrucata and E. victoriana) that dominate rocky, high‐elevation ridgelines of the Grampian Ranges and for one closely‐related, widespread species (E. baxteri) occupying flanking slopes and valleys. We documented genetic patterns across the landscape using cpDNA microsatellites, and related them to topographic features (exposed west‐facing versus protected east‐facing slopes and valleys). We also determined the extent of local haplotype sharing between populations of endemic species and neighbouring E. baxteri downslope with cpDNA microsatellites, and haplotype sharing between the endemic group and more distantly related species (E. obliqua, E. pauciflora and E. willisii) with sequences of the J_LA+ chloroplast region. Results We detected 26 cpDNA microsatellite haplotypes in a relatively small area of c. 20 km × 50 km. Populations of E. baxteri on east‐facing slopes and valleys had greater cpDNA microsatellite diversity than E. baxteri and endemic species on exposed west‐facing slopes. Endemic species frequently shared chloroplast haplotypes with E. baxteri downslope. Sharing of J_LA+ haplotypes with species outside the endemic group was mostly restricted to E. victoriana, which had cpDNA more similar to the species from other sections of Eucalyptus (E. obliqua, E. willisii and E. pauciflora). Main conclusions Intensive sampling of related species on small isolated mountain ranges allowed us to relate genetic diversity to fine‐scale habitats and to document extensive local haplotype sharing between species. This study contributes to a general understanding of the environmental conditions that enable plant population persistence by linking concentrations of genetic diversity to particular habitats.     

Assessing the utility of whole genome amplified DNA for next‐generation molecular ecology

DNA quantity can be a hindrance in ecological and evolutionary research programmes due to a range of factors including endangered status of target organisms, available tissue type, and the impact of field conditions on preservation methods. A potential solution to low‐quantity DNA lies in whole genome amplification (WGA) techniques that can substantially increase DNA yield. To date, few studies have rigorously examined sequence bias that might result from WGA and next‐generation sequencing of nonmodel taxa. To address this knowledge deficit, we use multiple displacement amplification (MDA) and double‐digest RAD sequencing on the grey mouse lemur (Microcebus murinus) to quantify bias in genome coverage and SNP calls when compared to raw genomic DNA (gDNA). We focus our efforts in providing baseline estimates of potential bias by following manufacturer's recommendations for starting DNA quantities (>100 ng). Our results are strongly suggestive that MDA enrichment does not introduce systematic bias to genome characterization. SNP calling between samples when genotyping both de‐novo and with a reference genome are highly congruent (>98%) when specifying a minimum threshold of 20X stack depth to call genotypes. Relative genome coverage is also similar between MDA and gDNA, and allelic dropout is not observed. SNP concordance varies based on coverage threshold, with 95% concordance reached at ~12X coverage genotyping de‐novo and ~7X coverage genotyping with the reference genome. These results suggest that MDA may be a suitable solution for next‐generation molecular ecological studies when DNA quantity would otherwise be a limiting factor.     

Distinct patterns of hybridization across a suture zone in a coral reef fish (Dascyllus trimaculatus)

Hybrid zones are natural laboratories for investigating the dynamics of gene flow, reproductive isolation, and speciation. A predominant marine hybrid (or suture) zone encompasses Christmas Island (CHR) and Cocos (Keeling) Islands (CKE), where 15 different instances of interbreeding between closely related species from Indian and Pacific Oceans have been documented. Here, we report a case of hybridization between genetically differentiated Pacific and Indian Ocean lineages of the three‐spot dascyllus, Dascyllus trimaculatus (Rüppell, 1829). Field observations indicate there are subtle color differences between Pacific and Indian Ocean lineages. Most importantly, population densities of color morphs and genetic analyses (mitochondrial DNA and SNPs obtained via RADSeq) suggest that the pattern of hybridization within the suture zone is not homogeneous. At CHR, both color morphs were present, mitochondrial haplotypes of both lineages were observed, and SNP analyses revealed both pure and hybrid genotypes. Meanwhile, in CKE, the Indian Ocean color morphs were prevalent, only Indian Ocean mitochondrial haplotypes were observed, and SNP analysis showed hybrid individuals with a large proportion (~80%) of their genotypes assigning to the Indian Ocean lineage. We conclude that CHR populations are currently receiving an influx of individuals from both ocean basins, with a greater influence from the Pacific Ocean. In contrast, geographically isolated CKE populations appear to be self‐recruiting and with more influx of individuals from the Indian Ocean. Our research highlights how patterns of hybridization can be different at scales of hundreds of kilometers, due to geographic isolation and the history of interbreeding between lineages.     

Isolation by instability: Historical climate change shapes population structure and genomic divergence of treefrogs in the Neotropical Cerrado savanna

Although the impact of Pleistocene glacial cycles on the diversification of the tropical biota was once dismissed, increasing evidence suggests that Pleistocene climatic fluctuations greatly affected the distribution and population divergence of tropical organisms. Landscape genomic analyses coupled with paleoclimatic distribution models provide a powerful way to understand the consequences of past climate changes on the present‐day tropical biota. Using genome‐wide SNP data and mitochondrial DNA, combined with projections of the species distribution across the late Quaternary until the present, we evaluate the effect of paleoclimatic shifts on the genetic structure and population differentiation of Hypsiboas lundii, a treefrog endemic to the South American Cerrado savanna. Our results show a recent and strong genetic divergence in H. lundii across the Cerrado landscape, yielding four genetic clusters that do not seem congruent with any current physical barrier to gene flow. Isolation by distance (IBD) explains some of the population differentiation, but we also find strong support for past climate changes promoting range shifts and structuring populations even in the presence of IBD. Post‐Pleistocene population persistence in four main areas of historical stable climate in the Cerrado seems to have played a major role establishing the present genetic structure of this treefrog. This pattern is consistent with a model of reduced gene flow in areas with high climatic instability promoting isolation of populations, defined here as “isolation by instability,” highlighting the effects of Pleistocene climatic fluctuations structuring populations in tropical savannas.