The Effect of Extirpation and Reintroduction on Genetic Variability of the Gray Wolf B. vonholdt, D. Stahler, D. Earl, D. Smith, J. Pollinger and R. K. Wayne Ecology and Evolutionary Biology, University of California, Los Angeles
TALK OUTLINE I) Introduction A) Relationships & genetic diversity 1. High gene flow but ecology matters II) Extinct wolves: historic wolves of the West A) Loss of genetic variation and relationships B) Population size III) Natural and human assisted reintroduction A) Yellowstone Natural Park 1. Genealogical reconstruction 2. Levels of variation and inbreeding
NJ coyote s 0.005 substitutions/site Crufus388 la4 la13 wolves lu33 la14 la11 la2 la10 la3 la5 la7 la1 la6 la12 la9 CSI CAU la8 lu1 lu17 lu26 lu23 lu2 lu3 lu14 lu21 lu30 lu13 lu24 lu15 lu16 lu25 lu31 lu7 lu10 lu11 lu12 lu29 lu5 lu34 lu6 lu18 lu19 lu20 lu22 lu4 lu8 lu9 lu28 lu32 lu33 lu27 Mexican wolf RELATIONSHIPS Vila et al., Mol. Ecol., 1999
TYPICAL SPECIES LARGE CARNIVORES 10 Nm 1 0.1 0.01 175 500 1000 5000 10000 20000 Distance (Km) Several subspecies Few subspecies
Ecology Matters!! Geffen et al., 2004
Distance Based Redundancy analysis (dbrda) Marginal tests Variable set F P %var Fst (microsatellite) Distance 1.822 0.0878 31.29 25.39 11.20 7.05 15.65 78.81 70.53 Fst (mtdna RFLP) Distance 4.929 0.0094 45.10 17.77 2.43 22.54 9.68 73.16 75.51 Anderson and Legendre, 1999
Ecological Imprinting Dispersal Cervus elaphus Alces gigas X X Rangifer tarandus
Genetically Defined Ecotypes Carmichael et al. 2006 Carmichael, 2006
The Great Extermination: Historic wolves of The American West Leonard et al., 2005
occidentalis arctos lycaon nubilis baileyi
Genetic Variation Historic wolves have about three times the variation of modern wolves (θ of 0.027 vs 0.010) and at least 43% of the variation has been lost. White -extant Black - historic only Red -Mexican wolf Samp le Haplotype Subsp ecies Localit y JAL 471 lu33 C.l. baileyi Sierra Madre, Chihua hua, Mexico JAL 473 lu33 C.l. baileyi Colonia Garcia, Chihuahua, Mexico JAL 474 lu60 C.l. baileyi Colonia Juarez, Chihuahua, Mexico JAL 475 lu32 C.l. baileyi Tanks, Ari zona JAL 476 lu33 C.l. baileyi Ft. Bowie, Arizona JAL 477 lu47 C.l. baileyi Fort Massachuse tts, New Mexico JAL 478 lu32 C.l. baileyi Fort Massachuse tts, New Mexico JAL 545 lu33 C.l. baileyi Cloverdale, New Mexico JAL 515 lu32 C.l.nub ilus Platt River, Colorado JAL 516 lu48 C.l.nub ilus Bent County, Colorado JAL 517 lu38 C.l.nub ilus Kansas JAL 518 lu28 C.l.nub ilus Fort Hark er, Kansas JAL 519 lu49 C.l.nub ilus Grove City, Kansas JAL 520 lu38 C.l.nub ilus Grove City, Kansas JAL 521 lu50 C.l.nub ilus Fort Kearny, Nebraska JAL 522 lu28 C.l.nub ilus Fort Kearny, Nebraska JAL 523 lu38 C.l.nub ilus Platt River, Nebraska JAL 524 lu50 C.l.nub ilus Guyo Canyon, New Mexico JAL 525 lu28 C.l.nub ilus Santa Rosa, New Mexico JAL 526 lu50 C.l.nub ilus Santa Rosa, New Mexico JAL 527 lu51 C.l.nub ilus Carthage, New Mexico JAL 528 lu28 C.l.nub ilus Medora, North Dakota JAL 529 lu50 C.l.nub ilus Witchita Forest Reserve, Oklahoma JAL 530 lu52 C.l.nub ilus Afton, Oklahoma JAL 539 lu50 C.l.nub ilus Haynes, New Mexico JAL 540 lu28 C.l.nub ilus La Sol, Utah JAL 541 lu50 C.l.nub ilus Box Elder County, Utah JAL 542 lu53 C.l.nub ilus Heart Draw, Utah JAL 543 lu28 C.l.nub ilus Laramie, Wyom ing JAL 544 lu28 C.l.nub ilus Wagon Creek, Wyoming JAL 480 lu32 C.l.nub ilus Porcup ine, Labrado r JAL 557 lu54 C.l.nub ilus Labrad or Peninsula JAL 558 lu32 C.l.nub ilus Labrad or Peninsula JAL 559 lu54 C.l.nub ilus Labrad or Peninsula
100 --Eurasian wolves --NA wolves * Historic wolves *Extant NA wolves 1 change 59 99 lu53 lu32 lu54 61 IWC * IWA 59 IWB lu61 59 lu1 lu17 lu26 lu23 lu30 lu31 lu22 lu8 lu9 lu37 lu12 * lu29 72* lu7 60 lu10 lu11 91 lu28 * lu52 * * lu38 lu48 * * * lu49 * * * * * * IWD 63 100 * * lu2 lu3 lu14 lu21 54 lu15 lu13 lu24 lu5 53 lu16 lu18 84 lu19 lu20 59 lu33 52 lu4 lu27 52 76 * * HWD HWA HWC HWE lu6 lu34 lu47 la11 lu60 lu25 * HWB * Mexican wolf lu50 lu51 * Southern clade
Genetic Variation White -extant Black - historic only Red -Mexican wolf Clade (lu 33, 47,50,51) Sample Haplotype Subspecies Localit y JAL 471 lu33 C.l. baileyi Sierra Madre, Chihuahua, Mexico JAL 473 lu33 C.l. baileyi Colonia Garcia, Chihuahua, Mex ico JAL 474 lu60 C.l. baileyi Colonia Juarez, Chihuahua, Mexico JAL 475 lu32 C.l. baileyi Tanks, Arizona JAL 476 lu33 C.l. baileyi Ft. B owie, Arizona JAL 477 lu47 C.l. baileyi Fort Massachusetts, New Mexico JAL 478 lu32 C.l. baileyi Fort Massachusetts, New Mexico JAL 545 lu33 C.l. baileyi Cloverdale, Ne w Mexico JAL 515 lu32 C.l.nubilus Platt River, Colorado JAL 516 lu48 C.l.nubilus Bent County, Colorado JAL 517 lu38 C.l.nubilus Kansas JAL 518 lu28 C.l.nubilus Fort Harker, Kansas JAL 519 lu49 C.l.nubilus Grove City, Kansas JAL 520 lu38 C.l.nubilus Grove City, Kansas JAL 521 lu50 C.l.nubilus Fort Kearny, Nebraska JAL 522 lu28 C.l.nubilus Fort Kearny, Nebraska JAL 523 lu38 C.l.nubilus Platt River, Nebraska JAL 524 lu50 C.l.nubilus Guyo Canyo n, New Mexico JAL 525 lu28 C.l.nubilus Santa Rosa, N ew Mexico JAL 526 lu50 C.l.nubilus Santa Rosa, N ew Mexico JAL 527 lu51 C.l.nubilus Carthage, Ne w Mexico JAL 528 lu28 C.l.nubilus Medora, North Dakota JAL 529 lu50 C.l.nubilus Witchita Forest Reserve, Oklahoma JAL 530 lu52 C.l.nubilus Afton, Oklahoma JAL 539 lu50 C.l.nubilus Haynes, New Mexico JAL 540 lu28 C.l.nubilus La Sol, Utah JAL 541 lu50 C.l.nubilus Box Elder County, Utah JAL 542 lu53 C.l.nubilus Heart Draw, Utah JAL 543 lu28 C.l.nubilus Laramie, Wyoming JAL 544 lu28 C.l.nubilus Wagon Creek, Wyoming JAL 480 lu32 C.l.nubilus Porcupine, Labrado r JAL 557 lu54 C.l.nubilus Labrador Peninsula JAL 558 lu32 C.l.nubilus Labrador Peninsula JAL 559 lu54 C.l.nubilus Labrador Peninsula Nebraska Oklahoma Utah
NJ coyotes 0.005 substitutions/site Crufus388 la4 la13 wolves lu33 la14 la11 la2 la10 la3 la5 la7 la1 la6 la12 la9 CSI CAU la8 lu1 lu17 lu26 lu23 lu2 lu3 lu14 lu21 lu30 lu13 lu24 lu15 lu16 lu25 lu31 lu7 lu10 lu11 lu12 lu29 lu5 lu34 lu6 lu18 lu19 lu20 lu22 lu4 lu8 lu9 lu28 lu32 lu27 Genetic Variation & population size value of θ of 0.027 implies a population size of 840,000 wolves (assuming 60% of adults reproduce, 50% of the population is adult and a 3 year generation time) Vila et al., Mol. Ecol., 1999, Leonard et al., 2005
Consider recovery involving both ecotypes over a large area N>300 N>700
Yellowstone Wolf Recovery 1995-1996 Abundant prey populations and protection led to rapid recovery
The Trophic Cascade
Yellowstone National Park Reintroduction event 1995: 14 wolves from Hinton, Alberta 1996: 17 wolves from Fort St. John, B.C. Genetic Analysis 30 canid microsatellite markers Sample collection (N = 344) mostly blood/tissue samples; serum, hair, teeth and scat Pedigree construction (N>200) Parentage analysis, relatedness simulation, observational pedigree likelihood estimation Bridgett vonholdt Dan Stahler
Pedigree Example: the Druid Peak pack since its establishment in 1996 Assembly rules for reintroduction Pack size = 6; Number of breeding individuals = 0
Druid Peak pack 1997 Pack size = 10; Number of breeding individuals = 3
Druid Peak pack 1998 Pack size = 10; Number of breeding individuals = 2
Druid Peak pack 1999 Pack size = 9; Number of breeding individuals = 2
Druid Peak pack 2000 Pack size = 18; Number of breeding individuals = 4
Druid Peak pack 2001 Pack size = 18; Number of breeding individuals = 4
Pack size = 18; Number of breeding individuals = 4 Druid Peak pack 2002 Pack splitting
Druid Peak pack 2003 Pack splitting events Pack size = 16 Number of breeding individuals = 4
Druid Peak pack 2004 Pack size = 14 Number of breeding individuals = 4
Pedigree of eleven wolf packs Assembly rules for reintroduction
Reintroduction Genetics Bensch et al. 2006
No apparent difference in the distribution of heterozygosity of F IS Proportion observed 0.50 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05 0.00 Breeding Population Nonbreeding Population 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 Heterozygosity Heterozygosity Breeding Nonbreeding Mean 0.7417 0.7405 Std Error 0.0093 0.0076 Median 0.7333 0.7586 Mode 0.8276 0.7667 Std Deviation 0.0761 0.1100 Sample Variance 0.0058 0.0121 Kurtosis -0.3738 2.6960 Skewness -0.1057-1.0330 N 67 208 F IS 0.008 0.011
Heterozygosity compared annually among breeding status 0.8400 0.8200 0.8000 0.7800 Heterozygosity 0.7600 0.7400 0.7200 0.7000 0.6800 0.6600 0.6400 0.6200 0.6000 Non-breeding population Breeding population 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 Year
Litter size and inbreeding No relationship between F and litter size 12 10 y = -8.3164x + 5.4079 R 2 adj = - 0.002 n =45 observed high count litter size during den emergence period 8 6 4 2 0-0.10-0.05 0.00 0.05 0.10 inbreeding coefficient (F) of the pups
Pup survival and inbreeding No relationship between F and pup survival 100 90 proportion of pups in litter surviving until winter 80 70 60 50 40 30 20 y = -139.44x + 78.125 R 2 adj = 0.016 n = 45 10 0-0.10-0.05 0.00 0.05 0.10 inbreeding coefficient (F) of the pups
Observed and predicted heterozygosity 0.8000 0.7500 0.7000 0.6500 0.6000 0.5500 0.5000 0.4500 0.4000 Predicted heterozygosity Observed heterozygosity 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 N e estimates from SPARKS and PM2000
Annual observed F IS and predicted F 0.2 0.15 0.1 Inbreeding 0.05 0-0.05-0.1 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 Year Projected FIS FIS observed
Heterozygosity of the adult breeding pools 0.8500 Observed Heterozygosity 0.8000 0.7500 0.7000 0.6500 0.6000 Observed Random Managed 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 Year Mate Suitability Index (MSI) = increase genetic diversity, decreases mean kinship, minimize inbreeding and exclude unknown lineages, MateRx
Black and White matings more common that pure colors (negative assortative mating) MALES KK Kk kk FEMALES KK 0 (0) 0 (0) 0 (0) 0 Kk 0 (0) 7 (13) 31 (25) 38 kk 0 (0) 48 (42) 70 (76) 118 0 55 101 156 pvalue = 0.01250471
CONCLUSIONS 1) The American West once contained a genetically diverse populations of wolves, numbering > 700,000 individuals. 2) Reintroduction should involve a large area and more individuals than planned of both wolf ecotypes to their appropriate environments. Mixing of ecotypes is not bad. 3) Large scale reintroductions such as occurred in Yellowstone even if highly isolated, can preserve variation and avoid inbreeding over the short term. 4) Long term preservation should involve gene flow.
Acknowledgements Yellowstone National Park, NPS Dr. Steve Fain, USFWS Forensics Lab, Ashland, OR Ed Bangs; USFWS MT Carter Niemeyer, USFWS ID Joe Fontaine, USFWS MT Roger Stradley (pilot), Gallatin Flying Services, MT Yellowstone Park Foundation Nez Perce Tribe, ID Deb Guernsey, Rick McIntyre, and numerous Yellowstone Wolf Project field technicians Photography provided by: Dan Stahler/NPS, Bridgett vonholdt and Carl Swaboda
Annual population F IS 200 0.03 180 0.02 160 0.01 Population N 140 120 100 80 60 0-0.01-0.02-0.03 FIS 40 20 Predicted ~0.20-0.04-0.05 0 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 Year -0.06 Population N FIS
Annual population heterozygosity Population N 200 180 160 140 120 100 80 60 0.7700 0.7600 0.7500 0.7400 0.7300 Heterozygosity 40 20 0 Predicted ~0.62 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 Year 0.7200 0.7100 Population N Heterozygosity
Differential Survivorship