Findings {#Sec1}
========

Common raven (*Corvus corax*) is a widespread, Holarctic bird that is listed as least concern on the IUCN Red List of Threatened Species (<http://www.iucnredlist.org/search>). Although common ravens in Alaska are not morphologically diverse (2 subspecies) \[[@CR1]\], genetic assessments using mitochondrial DNA markers have shown differences among lineages \[[@CR2]\] and have suggested that isolated, endemic populations are found in the Aleutian Islands of Alaska \[[@CR3]\]. Island populations were possibly isolated in glacial refugia as has been suggested in other studies of Aleutian landbirds \[[@CR3]--[@CR6]\]. Alternatively, populations might be post-glacial in origin having colonized after the last glacial maximum (\~15,000 ybp) \[[@CR7]\]. Given the recent timing of these events, a quickly evolving molecular marker is needed to assess the genetic diversity of populations, connectivity among populations, and to understand the history of ravens throughout the Holarctic. We developed eight microsatellite loci that will be used to assess population structure, rates of gene flow, colonization history, and to identify isolated populations of common raven.

Total genomic DNA was extracted from the muscle tissue of two common ravens from Alaska using a QiAamp DNA Mini Kit (Qiagen, Valencia, CA). The protocols of Glenn and Schable \[[@CR8]\] and Katzinel et al. \[[@CR9]\] were used to isolate microsatellite loci. Briefly, restriction enzymes were used to digest DNA, adapters were ligated onto fragments, microsatellite-containing fragments were separated using biotinylated probes and streptavidin beads, and fragments yields were increased using PCR. Fragments were run on a Roche 454 FLX with Titanium chemistry (454 Life Sciences, a Roche Company, Branford, CT) with similarly treated libraries, each with combinations of unique adapter and Roche MID-tags. We used demuxipy (<https://github.com/faircloth-lab/demuxipy/>) to demultiplex fragments, and used a version of MSATCOMMANDER \[[@CR10]\] to find 591 sequences containing microsatellites from 1467 sequence reads. Primer 3 \[[@CR11]\] was used to design primers with the CAG primer sequence (5′-CAGTCGGGCGTCATCA-3′) prepended to the 5′ end of one primer of each pair. This enables the use of a fluorescently labeled primer with the CAG sequence during PCR to tag amplified fragments \[[@CR12], [@CR13]\]. A short sequence (GTTT) was prepended to the primer without the CAG sequence to encourage adenylation of the amplicons for accurate genotyping \[[@CR14], [@CR15]\].

Fifty primer pairs were tested for amplification and polymorphism using the DNA from 10 individuals. Amplifications were in 20 µl volumes \[250 µg/mL BSA, 10× Buffer B (Fisher Scientific), 25 mM MgCl~2~, 5 µM unlabeled primer, 0.5 µM tag labeled primer, 5 µM universal dye-labeled primer, 2.5 mM dNTPs, 0.5 units Taq DNA polymerase (Fisher Scientific), and 20 ng DNA template\] using a BioRad MyCycler thermal cycler. Touchdown cycling conditions were used to amplify DNA and for attachment of the universal dye-labeled primer with the CAG tag. PCR parameters included an initial denaturation step of 2 min 30 s at 95 °C then 20 cycles of 95 °C for 20 s, 65 °C to 50 °C annealing temperature for 20 s (decreasing 0.5 °C per cycle), and extension step of 72 °C for 30 s followed by 15 cycles of 95 °C for 20 s, 55 °C to 45 °C for 20 s, and 72 °C for 30 s. Cycles were followed with a final extension step of 72 °C for 10 min. An ABI3730XL sequencer (Applied Biosystems) was used to determine genotypes. We found that fifteen primers amplified a product and eight of these primers exhibited polymorphism.

Initially, we tested for polymorphism using ten individuals from seven locations in Alaska including locations from the Aleutian Islands, the mainland of Alaska, and islands in southeast Alaska. Herein, we report the statistics associated with polymorphic loci using 18--20 individuals from Fairbanks, Alaska to avoid any Wahlund effect caused by population structure \[[@CR16]\]. We also tested the polymorphism of one to three individuals each for seven other species in the family Corvidae (Tables [1](#Tab1){ref-type="table"}, [2](#Tab2){ref-type="table"}) to assess cross-species amplification. We used Gene Mapper software (Applied Biosystems) to score alleles and to determine the number of alleles per locus (K), observed heterozygosity (H~o~), and expected heterozygosity (H~e~), using Arlequin ver. 3.5 \[[@CR17]\]. We also tested for Hardy--Weinberg equilibrium and linkage disequilibrium using Arlequin and used Microchecker \[[@CR18]\] to test for genotyping errors. We found that the number of alleles per locus ranged from 3 to 8, H~o~ ranged from 0.25 to 0.80, and H~e~ ranged from 0.30 to 0.80 for common ravens (Table [1](#Tab1){ref-type="table"}). All loci were in Hardy--Weinberg equilibrium and linkage equilibrium after correction for multiple tests and we found no evidence of null alleles, stuttering, or large-allele dropout.Table 1Characteristics of eight microsatellite loci in the common raven *Corvus corax*LocusPrimer sequence (5′--3′)Primer labelRepeat motifT~a~ (°C)nHWEKH~o~H~e~GenBank accession no.Coco6F: \*AACCAAGCTGTCAAAGGAGAVIC(ATCT)^13^60/50200.39680.550.54KT288258R: TTTCCTGGTGCGGAGAAGTACoco12F: \*GAGGCTTTCAGCGATATGGG6FAM(AGAT)^7^60/50180.17350.720.72KT288259R: GCTGCACTTGCCTGTTGTTTACoco30F: \*GTGCAAACCACTCCAGAAT6FAM(GGAT)^16^60/50200.91060.700.80KT288260R: GCTGTCACCAGTGGCTTTAATCoco31F: \*TGGCTACTTACCCGGAATCTCNED(ATGT)^8^60/50200.47030.250.30KT288261R: CTTCTCTACCGCCCATCACACoco32F: ATGCAAATCTCCGTGTATCANED(ATCT)^9^60/50200.42280.800.80KT288262R: \*CCTCCCTTTGTCCTTTCTGGCoco36F: CCTTCTGCGCTGTGTTCAAANED(CTT)^8^60/50200.00850.400.71KT288263R: \*CTGAGGTGCAAGCTTTCCTGCoco45F: TGACCAACACAGAGCAGATATT6FAM(GT)^11^60/50190.35330.530.54KT288264R: \*GGAACCTCTGAGTCCAAACoco50F: \*ATCCCATCAAAGGTCCACGVIC(AC)^13^60/50200.15750.500.64KT288265R: CCACGACAGCTGATGAGAGAAll loci were in Hardy--Weinberg and linkage equilibrium after Bonferroni correction for multiple tests*Asterisks* show location of CAG tag for each primer pair, *T* ~*a*~ annealing temperatures for touchdown protocol for each locus, *n* number of individuals genotyped, *K* number of alleles, *H* ~*o*~ observed heterozygosity, *H* ~*e*~ expected heterozygosityTable 2Number of alleles in eight microsatellite loci in seven Corvidae speciesSpeciesnNumber of alleles per locusCoco6Coco12Coco30Coco31Coco32Coco36Coco45Coco50Fish crow (*Corvus ossifragus*)241313221Northwestern crow (*Corvus caurinus*)2----311331American crow (*Corvus brachyrhynchos*)1----221211Steller's jay (*Cyanocitta stelleri*)3----21--213Blue jay (*Cyanocitta cristata*)2----32--212Gray jay (*Perisoreus canadensis*)2----314112Florida scrub-jay (*Aphelocoma coerulescens*)3----------------*n* number of individuals attempted to amplify and genotype, -- individuals failed to amplify at locus

Preliminary data from nine Alaska populations of common raven showed that the loci developed in this study are useful for population-level assessments with pairwise FST values ranging from 0.011 to 0.51. In addition, these loci are polymorphic in these populations (H~e~ 0.40--0.74; alleles per locus, 2.88--6.75). Eight loci that were developed for other Corvidae have shown polymorphism in common raven \[[@CR19]\]. By combining our newly developed markers with those from other taxa, it is possible that researchers could assess very fine-scale population structure.

In other Corvidae taxa, we found that the majority of loci amplified and were polymorphic in six of the seven species (Table [2](#Tab2){ref-type="table"}). All loci amplified in the fish crow (*Corvus ossifragus*) and five of eight loci were polymorphic. Loci that were monomorphic in other taxa might be polymorphic if a larger sample size was assessed. We suggest that these loci could be successfully used for population-level assessments in these species (Table [2](#Tab2){ref-type="table"}) and possibly in other members of the genera *Corvus* and *Perisoreus* but are unlikely to be useful in *Aphelocoma* (Table [2](#Tab2){ref-type="table"}).

Availability of the supporting data {#Sec2}
===================================

Microsatellite sequences were deposited in the National Center for Biotechnology Information (<http://www.ncbi.nlm.nih.gov>). They are accessible via GenBank accession numbers listed in Table [1](#Tab1){ref-type="table"}.

SL, TG, and BF developed and identified the microsatellite loci. LW, TL, and CS identified polymorphic loci and genotyped samples used in this manuscript. KW performed field work and provided DNA samples of ravens. CLP analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements {#FPar1}
================

This work was supported by the Florida Institute of Technology, University of Alaska Museum, and by the U. S. Department of Energy under award number DE-FC09-07SR22506 to the University of Georgia Research Foundation.

Competing interests {#FPar2}
===================

The authors declare that they have no competing interests.

Ethics statement {#FPar3}
================

Methods for obtaining samples for this research used the procedures outlined in "Guidelines to the Use of Wild Birds in Research" published by The Ornithological Council (<http://www.nmnh.si.edu/BIRDNET/guide/index.html>). These guidelines are the appropriate standard for Animal Welfare Act compliance. Bird sampling was approved by the University of Alaska Institutional Animal Care and Use Committee and in compliance with permits from state and federal wildlife agencies.
