Carnivore Contact: A Species Fracture Zone Delineated Amongst Genetically Structured North American Marten Populations (Martes americana and Martes caurina)

• Published in: Frontiers in Genetics
• Main Authors: Lucid, Michael, Cushman, Sam, Robinson, Lacy, Kortello, Andrea, Hausleitner, Doris, Mowat, Garth, Ehlers, Shannon, Gillespie, Sara, Svancara, Leona K., Sullivan, Jack, Rankin, Andrew, Paetkau, David
• Format: Text
• Language: English
• Published: Frontiers Media S.A. 2020
• Subjects: Genetics; Canada; Pacific; American marten; Martes americana
• Online Access: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC7370953/; http://www.ncbi.nlm.nih.gov/pubmed/32754203; https://doi.org/10.3389/fgene.2020.00735

North American martens are forest dependent, influenced by human activity, and climate vulnerable. They have long been managed and harvested throughout their range as the American marten (Martes americana). Recent work has expanded evidence for the original description of two species in North America — M. americana and the Pacific Coast marten, M. caurina — but the geographic boundary between these groups has not been described in detail. From 2010 to 2016 we deployed 734 multi-taxa winter bait stations across a 53,474 km(2) study area spanning seven mountain ranges within the anticipated contact zone along the border of Canada and the United States. We collected marten hair samples and developed genotypes for 15 polymorphic microsatellite loci for 235 individuals, and 493 base-pair sequences of the mtDNA gene COI for 175 of those individuals. Both nuclear and mitochondrial genetic structure identified a sharp break across the Clark Fork Valley, United States with M. americana and M. caurina occurring north and south of the break, respectively. We estimated global effective population size (N(e)) for each mountain range, clinal genetic neighborhood sizes (NS), calculated observed (H(o)) and expected (H(e)) heterozygosity, fixation index (F(ST)()), and clinal measures of allelic richness (Ar), H(o), and inbreeding coefficient (F(IS)). Despite substantial genetic structure, we detected hybridization along the fracture zone with both contemporary (nuclear DNA) and historic (mtDNA) gene flow. Marten populations in our study area are highly structured and the break across the fracture zone being the largest documented in North America (F(ST) range 0.21–0.34, mean = 0.27). With the exception of the Coeur d’Alene Mountains, marten were well distributed across higher elevation portions of our sampling area. Clinal NS values were variable suggesting substantial heterogeneity in marten density and movement. For both M. americana and M. caurina, elevationaly dependent gene flow and high genetic population structure suggest ...