Evolutionary factors affecting the cross‐species utility of newly developed microsatellite markers in seabirds

Microsatellite loci are ideal for testing hypotheses relating to genetic segregation at fine spatio‐temporal scales. They are also conserved among closely related species, making them potentially useful for clarifying interspecific relationships between recently diverged taxa. However, mutations at primer binding sites may lead to increased nonamplification, or disruptions that may result in decreased polymorphism in nontarget species. Furthermore, high mutation rates and constraints on allele size may also with evolutionary time, promote an increase in convergently evolved allele size classes, biasing measures of interspecific genetic differentiation. Here, we used next‐generation sequencing to develop microsatellite markers from a shotgun genome sequence of the sub‐Antarctic seabird, the thin‐billed prion (Pachyptila belcheri), that we tested for cross‐species amplification in other Pachyptila and related sub‐Antarctic species. We found that heterozygosity decreased and the proportion of nonamplifying loci increased with phylogenetic distance from the target species. Surprisingly, we found that species trees estimated from interspecific F_ST provided better approximations of mtDNA relationships among the studied species than those estimated using D_C, even though F_ST was more affected by null alleles. We observed a significantly nonlinear second order polynomial relationship between microsatellite and mtDNA distances. We propose that the loss of linearity with increasing mtDNA distance stems from an increasing proportion of homoplastic allele size classes that are identical in state, but not identical by descent. Therefore, despite high cross‐species amplification success and high polymorphism among the closely related Pachyptila species, we caution against the use of microsatellites in phylogenetic inference among distantly related taxa.