Assessing population structure at high levels of differentiation: microsatellite comparisons of bighorn sheep and large carnivores

Bighorn sheep associations: understanding tradeoffs of sociality and implications for disease transmission

Sociality directly influences mating success, survival rates, and disease, but ultimately likely evolved for its fitness benefits in a challenging environment. The tradeoffs between the costs and benefits of sociality can operate at multiple scales, resulting in different interpretations of animal behavior. We investigated the influence of intrinsic (e.g., relatedness, age) and extrinsic factors (e.g., land cover type, season) on direct contact (simultaneous GPS locations ≤ 25 m) rates of bighorn sheep (Ovis canadensis) at multiple scales near the Waterton-Glacier International Peace Park. During 2002-2012, male and female bighorn were equipped with GPS collars. Indirect contact (GPS locations ≤ 25 m regardless of time) networks identified two major breaks whereas direct contact networks identified an additional barrier in the population, all of which corresponded with prior disease exposure metrics. More direct contacts occurred between same-sex dyads than female-male dyads and between bighorn groups with overlapping summer home ranges. Direct contacts occurred most often during the winter-spring season when bighorn traveled at low speeds and when an adequate number of bighorn were collared in the area. Direct contact probabilities for all dyad types were inversely related to habitat quality, and differences in contact probability were driven by variables related to survival such as terrain ruggedness, distance to escape terrain, and canopy cover. We provide evidence that probabilities of association are higher when there is greater predation risk and that contact analysis provides valuable information for understanding fitness tradeoffs of sociality and disease transmission potential.

Genetic and demographic history define a conservation strategy for earth's most endangered pinniped, the Mediterranean monk seal Monachus monachus

The Mediterranean monk seal (Monachus monachus) is a flagship species for marine conservation, but important aspects of its life history remain unknown. Concerns over imminent extinction motivated a nuclear DNA study of the species in its largest continuous subpopulation in the eastern Mediterranean Sea. Despite recent evidence of partial subpopulation recovery, we demonstrate that there is no reason for complacency, as the species still shares several traits that are characteristic of a critically endangered species: Mediterranean monk seals in the eastern Mediterranean survive in three isolated and genetically depauperate population clusters, with small effective population sizes and high levels of inbreeding. Our results indicated male philopatry over short distances, which is unexpected for a polygynous mammal. Such a pattern may be explained by the species' unique breeding behavior, in which males defend aquatic territories near breeding sites, while females are often forced to search for new pupping areas. Immediate action is necessary to reverse the downward spiral of population decline, inbreeding accumulation and loss of genetic diversity. We propose concrete conservation measures for the Mediterranean monk seal focusing on reducing anthropogenic threats, increasing the population size and genetic diversity, and thus improving the long-term prospects of survival.

Spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis) at the northern limit of their native range

The Canadian Rocky Mountains are one of the few places on Earth where the spatial genetic structure of wide-ranging species has been relatively unaffected by anthropogenic disturbance. We characterized the spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis Shaw, 1804) in the northern portion of their range. Using microsatellites from 1495 individuals and mitochondrial DNA sequences from 188 individuals, we examined both broad- and fine-scale spatial genetic structure, assessed sex-biased gene flow within the northern portion of the species range, and identified geographic patterns of genetic diversity. We found that broad-scale spatial genetic structure was consistent with barriers to movement created by major river valleys. The fine-scale spatial genetic structure was characterized by a strong isolation-by-distance pattern, and analysis of neighborhood size using spatial autocorrelation indicated gene flow frequently occurred over distances of up to 100 km. However, analysis of sex-specific spatial autocorrelation and analysis of mitochondrial haplotype distributions failed to detect any evidence of sex-biased gene flow. Finally, our analyses reveal decreasing genetic diversity with increasing latitude, consistent with patterns of post-glacial recolonization of the Rocky Mountains.

Population genetic structure and Wolbachia infection in an endangered butterfly, Zizina emelina (Lepidoptera, Lycaenidae), in Japan

Zizina emelina (de l'Orza) is listed on Japan's Red Data List as an endangered species because of loss of its principal food plant and habitat. We compared parts of the mitochondrial and nuclear genes of this species to investigate the level of genetic differentiation among the 14 extant populations. We also examined infection of the butterfly with the bacterium Wolbachia to clarify the bacterium's effects on the host population's genetic structure. Mitochondrial and nuclear DNA analyses revealed that haplotype composition differed significantly among most of the populations, and the fixation index F ST was positively correlated with geographic distance. In addition, we found three strains of Wolbachia, one of which was a male killer; these strains were prevalent in several populations. There was linkage between some host mitochondrial haplotypes and the three Wolbachia strains, although no significant differences were found in a comparison of host mitochondrial genetic diversity with nuclear genetic diversity in Wolbachia-infected or -uninfected populations. These genetic analyses and Wolbachia infection findings show that Z. emelina has little migratory activity and that little gene flow occurs among the current populations.

Physiological vagility and its relationship to dispersal and neutral genetic heterogeneity in vertebrates

Vagility is the inherent power of movement by individuals. Vagility and the available duration of movement determine the dispersal distance individuals can move to interbreed, which affects the fine-scale genetic structure of vertebrate populations. Vagility and variation in population genetic structure are normally explained by geographic variation and not by the inherent power of movement by individuals. We present a new, quantitative definition for physiological vagility that incorporates aerobic capacity, body size, body temperature and the metabolic cost of transport, variables that are independent of the physical environment. Physiological vagility is the speed at which an animal can move sustainably based on these parameters. This meta-analysis tests whether this definition of physiological vagility correlates with empirical data for maximal dispersal distances and measured microsatellite genetic differentiation with distance {[F(ST)/[1-F(ST))]/ln distance} for amphibians, reptiles, birds and mammals utilizing three locomotor modes (running, flying, swimming). Maximal dispersal distance and physiological vagility increased with body mass for amphibians, reptiles and mammals utilizing terrestrial movement. The relative slopes of these relationships indicate that larger individuals require longer movement durations to achieve maximal dispersal distances. Both physiological vagility and maximal dispersal distance were independent of body mass for flying vertebrates. Genetic differentiation with distance was greatest for terrestrial locomotion, with amphibians showing the greatest mean and variance in differentiation. Flying birds, flying mammals and swimming marine mammals showed the least differentiation. Mean physiological vagility of different groups (class and locomotor mode) accounted for 98% of the mean variation in genetic differentiation with distance in each group. Genetic differentiation with distance was not related to body mass. The physiological capacity for movement (physiological vagility) quantitatively predicts genetic isolation by distance in the vertebrates examined.

Genomic consequences of genetic rescue in an insular population of bighorn sheep (Ovis canadensis)

Genetic rescue is a management intervention whereby a small population is supplemented with individuals from other populations in an attempt to reverse the effects of inbreeding and increased genetic load. One such rescue was recently documented in the population of bighorn sheep (Ovis canadensis) within the National Bison Range wildlife refuge (Montana, USA). Here, we examine the locus-specific effects of rescue in this population using a newly developed genome-wide set of 195 microsatellite loci and first-generation linkage map. We found that the rate of introgression varied among loci and that 111 loci, 57% of those examined, deviated from patterns of neutral inheritance. The most common deviation was an excess of homozygous genotypes relative to neutral expectations, indicative of directional selection. As in previous study of this rescue, individuals with more introduced alleles had higher reproductive success and longevity. In addition, we found 30 loci, distributed throughout the genome, which seem to have individual effects on these life history traits. Although the potential for outbreeding depression is a major concern when translocating individuals between populations, we found no evidence of such effects in this population.

Translating effects of inbreeding depression on component vital rates to overall population growth in endangered bighorn sheep

Evidence of inbreeding depression is commonly detected from the fitness traits of animals, yet its effects on population growth rates of endangered species are rarely assessed. We examined whether inbreeding depression was affecting Sierra Nevada bighorn sheep (Ovis canadensis sierrae), a subspecies listed as endangered under the U.S. Endangered Species Act. Our objectives were to characterize genetic variation in this subspecies; test whether inbreeding depression affects bighorn sheep vital rates (adult survival and female fecundity); evaluate whether inbreeding depression may limit subspecies recovery; and examine the potential for genetic management to increase population growth rates. Genetic variation in 4 populations of Sierra Nevada bighorn sheep was among the lowest reported for any wild bighorn sheep population, and our results suggest that inbreeding depression has reduced adult female fecundity. Despite this population sizes and growth rates predicted from matrix-based projection models demonstrated that inbreeding depression would not substantially inhibit the recovery of Sierra Nevada bighorn sheep populations in the next approximately 8 bighorn sheep generations (48 years). Furthermore, simulations of genetic rescue within the subspecies did not suggest that such activities would appreciably increase population sizes or growth rates during the period we modeled (10 bighorn sheep generations, 60 years). Only simulations that augmented the Mono Basin population with genetic variation from other subspecies, which is not currently a management option, predicted significant increases in population size. Although we recommend that recovery activities should minimize future losses of genetic variation, genetic effects within these endangered populations-either negative (inbreeding depression) or positive (within subspecies genetic rescue)-appear unlikely to dramatically compromise or stimulate short-term conservation efforts. The distinction between detecting the effects of inbreeding depression on a component vital rate (e.g., fecundity) and the effects of inbreeding depression on population growth underscores the importance of quantifying inbreeding costs relative to population dynamics to effectively manage endangered populations.

The genetic impact of translocations and habitat fragmentation in chamois (Rupicapra) spp

The chamois is a useful species with which to investigate the combined genetic impact of habitat fragmentation, over hunting, and translocations. Genetic variation within and between chamois (genus Rupicapra) populations was analyzed in 259 individuals from 16 sampling sites located in Italy, Spain, Slovakia, and the Czech Republic. Two mitochondrial DNA markers (control region and cytochrome b) and 11 nuclear microsatellites were typed. The principal results of this study can be summarized as follows: 1) high and significant differentiation between almost all chamois populations is observed even on a microgeographical scale, probably caused by the patchy distribution of this species, sharp geographical barriers to gene flow, and drift effects related to recent bottlenecks; 2) historical translocation events have left a clear genetic signature, including interspecific hybridization in some Alpine localities; 3) the Apennine subspecies of chamois, Rupicapra pyrenaica ornata, shows a high and similar level of divergence (about 1.5 My) from the Pyrenean (Rupicapra pyrenaica pyrenaica) and the Alpine (Rupicapra rupicapra) chamois; therefore, the specific status of these taxa should be revised. These results confirm the potential of population genetic analyses to dissect and interpret complex patterns of diversity in order to define factors important to conservation and management.

Molecular ecological approaches to studying the evolutionary impact of selective harvesting in wildlife

Harvesting of wildlife populations by humans is usually targeted by sex, age or phenotypic criteria, and is therefore selective. Selective harvesting has the potential to elicit a genetic response from the target populations in several ways. First, selective harvesting may affect population demographic structure (age structure, sex ratio), which in turn may have consequences for effective population size and hence genetic diversity. Second, wildlife-harvesting regimes that use selective criteria based on phenotypic characteristics (e.g. minimum body size, horn length or antler size) have the potential to impose artificial selection on harvested populations. If there is heritable genetic variation for the target characteristic and harvesting occurs before the age of maturity, then an evolutionary response over time may ensue. Molecular ecological techniques offer ways to predict and detect genetic change in harvested populations, and therefore have great utility for effective wildlife management. Molecular markers can be used to assess the genetic structure of wildlife populations, and thereby assist in the prediction of genetic impacts by delineating evolutionarily meaningful management units. Genetic markers can be used for monitoring genetic diversity and changes in effective population size and breeding systems. Tracking evolutionary change at the phenotypic level in the wild through quantitative genetic analysis can be made possible by genetically determined pedigrees. Finally, advances in genome sequencing and bioinformatics offer the opportunity to study the molecular basis of phenotypic variation through trait mapping and candidate gene approaches. With this understanding, it could be possible to monitor the selective impacts of harvesting at a molecular level in the future. Effective wildlife management practice needs to consider more than the direct impact of harvesting on population dynamics. Programs that utilize molecular genetic tools will be better positioned to assess the long-term evolutionary impact of artificial selection on the evolutionary trajectory and viability of harvested populations.

Gene flow in complex landscapes: testing multiple hypotheses with causal modeling

Predicting population-level effects of landscape change depends on identifying factors that influence population connectivity in complex landscapes. However, most putative movement corridors and barriers have not been based on empirical data. In this study, we identify factors that influence connectivity by comparing patterns of genetic similarity among 146 black bears (Ursus americanus), sampled across a 3,000-km(2) study area in northern Idaho, with 110 landscape-resistance hypotheses. Genetic similarities were based on the pairwise percentage dissimilarity among all individuals based on nine microsatellite loci (average expected heterozygosity=0.79). Landscape-resistance hypotheses describe a range of potential relationships between movement cost and land cover, slope, elevation, roads, Euclidean distance, and a putative movement barrier. These hypotheses were divided into seven organizational models in which the influences of barriers, distance, and landscape features were statistically separated using partial Mantel tests. Only one of the competing organizational models was fully supported: patterns of genetic structure are primarily related to landscape gradients of land cover and elevation. The alternative landscape models, isolation by barriers and isolation by distance, are not supported. In this black bear population, gene flow is facilitated by contiguous forest cover at middle elevations.

Genetic rescue of an insular population of large mammals

Natural populations worldwide are increasingly fragmented by habitat loss. Isolation at small population size is thought to reduce individual and population fitness via inbreeding depression. However, little is known about the time-scale over which adverse genetic effects may develop in natural populations or the number and types of traits likely to be affected. The benefits of restoring gene flow to isolates are therefore also largely unknown. In contrast, the potential costs of migration (e.g. disease spread) are readily apparent. Management for ecological connectivity has therefore been controversial and sometimes avoided. Using pedigree and life-history data collected during 25 years of study, we evaluated genetic decline and rescue in a population of bighorn sheep founded by 12 individuals in 1922 and isolated at an average size of 42 animals for 10-12 generations. Immigration was restored experimentally, beginning in 1985. We detected marked improvements in reproduction, survival and five fitness-related traits among descendants of the 15 recent migrants. Trait values were increased by 23-257% in maximally outbred individuals. This is the first demonstration, to our knowledge, of increased male and female fitness attributable to outbreeding realized in a fully competitive natural setting. Our findings suggest that genetic principles deserve broader recognition as practical management tools with near-term consequences for large-mammal conservation.

Population structure of California coyotes corresponds to habitat-specific breaks and illuminates species history

Little is known about the relationship between animal movements and the emergent structure of populations, especially for species occupying large continuous distributions. Some such mammals disperse disproportionately into habitat similar to their natal habitat, a behavioural bias that might be expected to lead to habitat-conforming genetic structure. We hypothesized that coyotes (Canis latrans) would exhibit such natal-biased dispersal, and used 13 microsatellite loci to test, correspondingly, whether genetic structure conformed to major habitat breaks. First, we used a model-based approach to assign coyote genotypes to distinct genetic clusters irrespective of geographical location. Visualization on a geographical information system revealed a strong concordance between the locations of cluster assignments and habitat bioregions, not explainable in terms of physical dispersal barriers or intervening low-quality habitat. Next, we used a multiple Mantel test, which controlled for effects of geographical distance (itself, marginally significant; P = 0.06), to statistically determine that genetic distance was indeed higher between than within bioregions (P < 0.001). Whereas previously published examples of landscape effects on gene flow have typically been explainable in terms of species-wide habitat affinities or dispersal barriers, our finding that genetic subdivisions were associated with unobstructed boundaries between contiguous habitats suggests a role for intraspecific variability in habitat affinities as a factor underlying genetic structure. In addition, our data combined with previously published data suggest a pattern of genetic isolation-by-distance throughout western North America, consistent with independent evidence that the western half of the coyote range predates European settlement.

Phylogenetic analysis of the tribe Bovini using microsatellites

The objective of the present study was to determine if the generally accepted phylogenetic relationships in the tribe Bovini correspond to a phylogenetic scheme derived from polymorphisms at 20 bovine microsatellite loci. This study comprises 17 representative populations: eight Bos taurus, two Bos indicus, one Poëphagus, one Bibos, one Bison, three Bubalus and one Syncerus. Phylogenetic analyses using (delta mu)2 and chord (DC) distances revealed substantial divergence among species. Neighbor-joining trees with both distance measures showed only minor differences. Bos taurus and Bos indicus grouped first, followed by Bos frontalis (Mithan) and Bos grunniens (Yak), Bison bison branched off next and Bubalus bubalis and Syncerus caffer emerged as the two most divergent species from the Bos clade. These findings would suggest that Bos, Poëphagus, and Bibos should be integrated into the Bos genus with each group classified as a subgenus. On the other hand, Bison, Bubalus and Syncerus should each be considered a separate genus. Direct estimates of the divergence times were calculated using the (delta mu)2 genetic distance. Bos taurus and Bos indicus were estimated to have diverged 0.31-0.82 MYA, Bos and Poëphagus: 0.57-1.53 MYA, Bos and Bibos: 0.57-1.52 MYA, Bos and Bison: 0.46-1.23 MYA, Bos and Bubalus: 1.85-4.93 MYA and Bos and Syncerus: 0.98-2.61 MYA.

New methods employing multilocus genotypes to select or exclude populations as origins of individuals

A new method for assigning individuals of unknown origin to populations, based on the genetic distance between individuals and populations, was compared to two existing methods based on the likelihood of multilocus genotypes. The distribution of the assignment criterion (genetic distance or genotype likelihood) for individuals of a given population was used to define the probability that an individual belongs to the population. Using this definition, it becomes possible to exclude a population as the origin of an individual, a useful extension of the currently available assignment methods. Using simulated data based on the coalescent process, the different methods were evaluated, varying the time of divergence of populations, the mutation model, the sample size, and the number of loci. Likelihood-based methods (especially the Bayesian method) always performed better than distance methods. Other things being equal, genetic markers were always more efficient when evolving under the infinite allele model than under the stepwise mutation model, even for equal values of the differentiation parameter F(st). Using the Bayesian method, a 100% correct assignment rate can be achieved by scoring ca. 10 microsatellite loci (H approximately 0.6) on 30-50 individuals from each of 10 populations when the F(st) is near 0.1.