Missing the forest for the gene trees: Conservation genetics is more than the identification of distinct population segments.—Zink et al. (2010) reinterpreted Barr et al. (2008) with the apparent agenda of espousing the use of mitochondrial DNA (mtDNA) over microsatellites in conservation genetics. In doing so, Zink et al. (2010) poorly represented both Barr et al. (2008) and the general value of microsatellites in population genetics research. We are compelled to correct some of the misconceptions that may have been created by Zink et al. (2010), and to underscore the value of microsatellite data in the conservation of the endangered Black-capped Vireo (Vireo atricapilla).
We agree with Zink et al. (2010) that mtDNA can be an excellent marker choice for characterizing phylogenetic structure, especially that resulting from disruption of gene flow between groups over many generations. In exhaustively detailing the pitfalls of using microsatellites with very little discussion of the strengths, however, Zink et al. (2010) appear to be making a strawman of our study. Barr et al. (2008) made no attempt to classify reciprocally monophyletic groups, as was the central focus of Zink et al. (2010). The stated purposes of both Zink et al. (2010) and Barr et al. (2008) centered around characterizing population structure in the Black-capped Vireo, whereas the scope of our study was rather more current: to ascertain the level of contemporary gene flow in the songbird given both extensive 20th-century habitat fragmentation and recent significant recovery of census populations. Zink et al. (2010:798) did note that “it is possible to detect fine-scale (i.e., recent) spatial structure with a set of highly polymorphic microsatellite loci,” a statement with which we certainly agree, given the aforementioned purpose of our study.
Zink et al. (2010) asserted that contemporary patterns of isolation and drift would be uniformly better detected by mtDNA because its smaller effective population size (Ne) would impart a shorter coalescence time than that of microsatellites. This is not as clear to us. Though a smaller Ne would theoretically confer faster coalescence, the statistical power gained from using many microsatellite loci versus the inherently single-locus mtDNA likely outweighs the latter marker system's presumed increased sensitivity to genetic drift and gene flow. This is especially important when studying processes that are necessarily very recent and, hence, would not be detected by an mtDNA tree, such as the effects of habitat fragmentation on gene flow. Furthermore, it is questionable whether the differences in coalescence time between mtDNA and microsatellites are as large as suggested by Zink et al. (2010). If, for example, the variance in reproductive success of males is much larger than that of females, the difference between the Ne of mtDNA and that of nuclear DNA (nDNA) will be small (Ballard and Whitlock 2004). Laporte and Charlesworth (2002) also illustrated that when males have higher variance in reproductive success and a lower dispersal rate than females, there is a sharp decrease in the difference in population structure expected between mtDNA and nDNA. Because males are more philopatric than females (Grzybowski 1995) and extrapair copulations are common in songbirds, this scenario is quite likely in Black-capped Vireos. There are, unfortunately, insufficient data to directly compare Ne between microsatellites and mtDNA here, but it would be surprising if the difference is as substantial as argued by Zink et al. (2010).
Zink et al. (2010) assessed how Barr et al. (2008) might be interpreted in the context of distinct population segments (DPSs). Barr et al. (2008) did not attempt to identify DPSs, and the criticisms of microsatellites and their use in this context should not be interpreted as damaging to our study. Barr et al. (2008) reported patterns of differentiation that suggest a recent decrease in gene flow between most of the major extant Black-capped Vireo aggregations in Texas and Oklahoma, as well as an inferred limitation on the Neotropical migrants' dispersal capability. Solely regarding the mtDNA data in the absence of the microsatellite results might lead to the incorrect assumption that gene flow is currently high all across the Black-capped Vireo's range, given that no geographically relevant patterns were detected (see Zink et al. 2010). On the basis of our microsatellite data, we could not agree more with the conservation recommendation by Zink et al. (2010:802) that future effort be made in connecting “high-value existing sites.”
Lost in the lengthy discussion of DPSs by Zink et al. (2010) are some potentially important inferences from the microsatellite data for the conservation of the Black-capped Vireo. For instance, although their populations have dramatically recovered, the level of differentiation we detected suggests that individuals are not freely dispersing between major aggregations (Barr et al. 2008). We do not agree with Zink et al. (2010) that this structure can be explained solely by isolation-by-distance. In fact, we reported that “once distances separating pairs of populations exceed 100 km, there appears to be little evidence that genetic differentiation increases with greater distances” (Barr et al. 2008:3635). This pattern was the result of nearby aggregations having similar allele frequencies, a relationship that breaks down at >100 km, beyond which populations were all significantly differentiated with no evidence of isolation-by-distance (see Barr et al. 2008: figure 2). We suspect this to be a signal of localized dispersal, with genetic drift following founder events or otherwise small Ne obscuring isolation-by-distance at greater distances.
Zink et al. (2010) further criticized our conclusions that biologically significant differentiation exists on the basis of the low values of FST and “subjective” significance tests we reported. Although we appreciate the concerns about the potential problems associated with using P values as indicators of differentiation, the dismissal of small but significant levels of differentiation apparently results from a misunderstanding of the nature of estimating FST from highly polymorphic microsatellites. It is widely recognized that FST decreases with increasing marker polymorphism (Hedrick 1999, 2005; Jost 2008), and correcting for this problem remains a contested subject (see Whitlock 2011). At present, we can say that on the basis of Hedrick (2005) and Jost (2008), the FST values we reported probably underestimated the true differentiation, and that many pairs of Black-capped Vireo aggregations exhibited significant differences in allele frequencies. We also question the complete dismissal by Zink et al. (2010) of significant differentiation in allele frequencies at the level of α < 0.001 in a highly vagile passerine as having no biological relevance. Random sampling error is unlikely, and, as previously discussed, isolation-by-distance breaks down after a short distance. Ignoring these differences risks missing important processes that may be affecting Black-capped Vireo populations, such as habitat fragmentation and restricted dispersal over larger spatial scales.
In order to better explore the assertion that differentiation among sample locations separated by <600 km, as those studied in Barr et al (2008), has no biological significance, we searched the literature for microsatellite-based studies of Neotropical migratory passerines for comparison. We identified eight relevant studies (Gibbs et al. 2000, Clegg et al. 2003, Veit et al. 2005, Ruegg et al. 2006, Davis et al. 2006, Lindsay et al. 2008, Winker et al. 2008, Bull et al. 2010), calculated the average FST between population pairs in these studies separated by ≤600 km, and compared these eight means to that of Black-capped Vireo population pairs using a one-sample t-test. The level of differentiation between major aggregations of other Neotropical migrants studied to date was significantly less than that of the Black-capped Vireo (t = -3.8, df = 7, P = 0.006). Furthermore, FST between Black-capped Vireo aggregations was double, on average, that of other previously studied songbirds (data available on request). Considering these findings, we reject the opinion of Zink et al. (2010) that the level of differentiation we reported is biologically irrelevant.
In conclusion, Barr et al. (2008) analyzed patterns of allele frequencies over space and attempted to ascertain the underlying mechanisms driving or restricting gene flow on the basis of geographic and ecological features—focal areas of study in conservation genetics. While noting that it was neither a goal nor a conclusion of our study, we applaud Zink et al. (2010) for applying mtDNA data to facilitate the management of the Black-capped Vireo by searching for potential DPSs—another, different focal point of conservation genetics. We are pleased that the mtDNA-based analyses, which illustrate little spatial structure over a deeper time-scale, are congruent with our suspicions that the restriction of gene flow between Black-capped Vireo aggregations is likely recent. It is our opinion that future investigators should continue recognizing the merits of both mtDNA and microsatellites and appropriately using both genetic markers in conservation genetics research.
Acknowledgments.
We thank the investigators, especially S. M. Clegg, V. L. Friesen, K. H. Ruegg, and K. Winker, who made their data sets available to us to make comparisons of genetic differentation to their studies. We also thank L. Joseph for his comments on the manuscript of this letter.