EVOLUTIONARY HISTORY OF THE NORTHERN LEOPARD FROG: RECONSTRUCTION OF PHYLOGENY, PHYLOGEOGRAPHY, AND HISTORICAL CHANGES IN POPULATION DEMOGRAPHY FROM MITOCHONDRIAL DNA

Balancing selection at a frog antimicrobial peptide locus: fluctuating immune effector alleles?

Balancing selection is common on many defense genes, but it has rarely been reported for immune effector proteins such as antimicrobial peptides (AMPs). We describe genetic diversity at a brevinin-1 AMP locus in three species of leopard frogs (Rana pipiens, Rana blairi, and Rana palustris). Several highly divergent allelic lineages are segregating at this locus. That this unusual pattern results from balancing selection is demonstrated by multiple lines of evidence, including a ratio of nonsynonymous/synonymous polymorphism significantly higher than 1, the ZnS test, incongruence between the number of segregating sites and haplotype diversity, and significant Tajima's D values. Our data are more consistent with a model of fluctuating selection in which alleles change frequencies over time than with a model of stable balancing selection such as overdominance. Evidence for fluctuating selection includes skewed allele frequencies, low levels of synonymous variation, nonneutral values of Tajima's D within allelic lineages, an inverse relationship between the frequency of an allelic lineage and its degree of polymorphism, and divergent allele frequencies among populations. AMP loci could be important sites of adaptive genetic diversity, with consequences for host–pathogen coevolution and the ability of species to resist disease epidemics.

Genetic diversity and structure in Canadian northern leopard frog (Rana pipiens) populations: implications for reintroduction programs

The northern leopard frog ( Rana pipiens Schreber, 1782) underwent a large decline in the western portion of its range and only occurs in 20% of historically occupied sites in Alberta. Its absence may reflect an inability to disperse to these sites because of habitat fragmentation, and human-mediated translocation has been proposed. In this study, we used three criteria to examine the genetic suitability of potential translocation sources: diversity, similarity to area of reintroduction, and evolutionary history. We genotyped 187 samples and sequenced 812 bp of the mitochondrial NADH dehydrogenase 1 gene from 14 Canadian northern leopard frog populations. Nuclear and mitochondrial diversity were highest in Manitoba and western Ontario and declined westward. There was no significant relationship between genetic and geographic distance, suggesting that genetic drift is a driving force affecting the genetic relationships between populations. Regions separated by more than ~50 km were quite differentiated. Therefore, source populations similar to the original inhabitants of an area for reintroduction may be uncommon. Mitochondrial analyses revealed that all populations share a close evolutionary history, belonging to the western haplotype group. While genetic criteria support the use of Manitoba and Ontario as sources, the desirability of environmental similarity to the reintroduction site suggests that ecologically exchangeable Alberta populations should also be considered.

Effective population sizes and temporal stability of genetic structure in Rana pipiens, the northern leopard frog

Although studies of population genetic structure are very common, whether genetic structure is stable over time has been assessed for very few taxa. The question of stability over time is particularly interesting for frogs because it is not clear to what extent frogs exist in dynamic metapopulations with frequent extinction and recolonization, or in stable patches at equilibrium between drift and gene flow. In this study we collected tissue samples from the same five populations of leopard frogs, Rana pipiens, over a 22-30 year time interval (11-15 generations). Genetic structure among the populations was very stable, suggesting that these populations were not undergoing frequent extinction and colonization. We also estimated the effective size of each population from the change in allele frequencies over time. There exist few estimates of effective size for frog populations, but the data available suggest that ranid frogs may have much larger ratios of effective size (Ne) to census size (Nc) than toads (bufonidae). Our results indicate that R. pipiens populations have effective sizes on the order of hundreds to at most a few thousand frogs, and Ne/Nc ratios in the range of 0.1-1.0. These estimates of Ne/Nc are consistent with those estimated for other Rana species. Finally, we compared the results of three temporal methods for estimating Ne. Moment and pseudolikelihood methods that assume a closed population gave the most similar point estimates, although the moment estimates were consistently two to four times larger. Wang and Whitlock's new method that jointly estimates Ne and the rate of immigration into a population (m) gave much smaller estimates of Ne and implausibly large estimates of m. This method requires knowing allele frequencies in the source of immigrants, but was thought to be insensitive to inexact estimates. In our case the method may have failed because we did not know the true source of immigrants for each population. The method may be more sensitive to choice of source frequencies than was previously appreciated, and so should be used with caution if the most likely source of immigrants cannot be identified clearly.