Global diversity and oceanic divergence of humpback whales (Megaptera novaeangliae).

Humpback whales (Megaptera novaeangliae) annually undertake the longest migrations between seasonal feeding and breeding grounds of any mammal. Despite this dispersal potential, discontinuous seasonal distributions and migratory patterns suggest that humpbacks form discrete regional populations within each ocean. To better understand the worldwide population history of humpbacks, and the interplay of this species with the oceanic environment through geological time, we assembled mitochondrial DNA control region sequences representing approximately 2700 individuals (465 bp, 219 haplotypes) and eight nuclear intronic sequences representing approximately 70 individuals (3700 bp, 140 alleles) from the North Pacific, North Atlantic and Southern Hemisphere. Bayesian divergence time reconstructions date the origin of humpback mtDNA lineages to the Pleistocene (880 ka, 95% posterior intervals 550-1320 ka) and estimate radiation of current Northern Hemisphere lineages between 50 and 200 ka, indicating colonization of the northern oceans prior to the Last Glacial Maximum. Coalescent analyses reveal restricted gene flow between ocean basins, with long-term migration rates (individual migrants per generation) of less than 3.3 for mtDNA and less than 2 for nuclear genomic DNA. Genetic evidence suggests that humpbacks in the North Pacific, North Atlantic and Southern Hemisphere are on independent evolutionary trajectories, supporting taxonomic revision of M. novaeangliae to three subspecies.

Multigene analysis suggests ecological speciation in the fungal pathogen Claviceps purpurea.

Claviceps purpurea is an important pathogen of grasses and source of novel chemical compounds. Three groups within this species (G1, G2 and G3) have been recognized based on habitat association, sclerotia and conidia morphology, as well as alkaloid production. These groups have further been supported by Random Amplification of Polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP) markers, suggesting this species may be more accurately described as a species complex. However, all divergent ecotypes can coexist in sympatric populations with no obvious physical barriers to prevent gene flow. In this study, we used both phylogenetic and population genetic analyses to test for speciation within C. purpurea using DNA sequences from ITS, a RAS-like locus, and a portion of beta-tubulin. The G1 types are significantly divergent from the G2/G3 types based on each of the three loci and the combined dataset, whereas the G2/G3 types are more integrated with one another. Although the G2 and G3 lineages have not diverged as much as the G1 lineage based on DNA sequence data, the use of three DNA loci does reliably separate the G2 and G3 lineages. However, the population genetic analyses strongly suggest little to no gene flow occurring between the different ecotypes, and we argue that this process is driven by adaptations to ecological habitats; G1 isolates are associated with terrestrial grasses, G2 isolates are found in wet and shady environments, and G3 isolates are found in salt marsh habitats.

Widespread unidirectional transfer of mitochondrial DNA: a case in western Palaearctic water frogs.

Interspecies transfer of mitochondrial (mt) DNA is a common phenomenon in plants, invertebrates and vertebrates, normally linked with hybridization of closely related species in zones of sympatry or parapatry. In central Europe, in an area north of 48 degrees N latitude and between 8 degrees and 22 degrees E longitude, western Palaearctic water frogs show massive unidirectional introgression of mtDNA: 33.7% of 407 Rana ridibunda possessed mtDNA specific for Rana lessonae. By contrast, no R. lessonae with R. ridibunda mtDNA was observed. That R. ridibunda with introgressed mitochondrial genomes were found exclusively within the range of the hybrid Rana esculenta and that most hybrids had lessonae mtDNA (90.4% of 335 individuals investigated) is evidence that R. esculenta serves as a vehicle for transfer of lessonae mtDNA into R. ridibunda. Such introgression has occurred several times independently. The abundance and wide distribution of individuals with introgressed mitochondrial genomes show that R. lessonae mt genomes work successfully in a R. ridibunda chromosomal background despite their high sequence divergence from R. ridibunda mtDNAs (14.2-15.2% in the ND2/ND3 genes). Greater effectiveness of enzymes encoded by R. lessonae mtDNA may be advantageous to individuals of R. ridibunda and probably R. esculenta in the northern parts of their ranges.

Strength and tempo of directional selection in the wild.

Directional selection is a major force driving adaptation and evolutionary change. However, the distribution, strength, and tempo of phenotypic selection acting on quantitative traits in natural populations remain unclear across different study systems. We reviewed the literature (1984-1997) that reported the strength of directional selection as indexed by standardized linear selection gradients (beta). We asked how strong are viability and sexual selection, and whether strength of selection is correlated with the time scale over which it was measured. Estimates of the magnitude of directional selection (absolute value of beta) were exponentially distributed, with few estimates greater than 0.50 and most estimates less than 0.15. Sexual selection (measured by mating success) appeared stronger than viability selection (measured by survival). Viability selection that was measured over short periods (days) was typically stronger than selection measured over longer periods (months and years), but the strength of sexual selection did not vary with duration of selection episodes; as a result, sexual selection was stronger than viability selection over longer time scales (months and years), but not over short time scales (days).

Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach.

A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source-sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.

The strength of phenotypic selection in natural populations.

How strong is phenotypic selection on quantitative traits in the wild? We reviewed the literature from 1984 through 1997 for studies that estimated the strength of linear and quadratic selection in terms of standardized selection gradients or differentials on natural variation in quantitative traits for field populations. We tabulated 63 published studies of 62 species that reported over 2,500 estimates of linear or quadratic selection. More than 80% of the estimates were for morphological traits; there is very little data for behavioral or physiological traits. Most published selection studies were unreplicated and had sample sizes below 135 individuals, resulting in low statistical power to detect selection of the magnitude typically reported for natural populations. The absolute values of linear selection gradients |beta| were exponentially distributed with an overall median of 0.16, suggesting that strong directional selection was uncommon. The values of |beta| for selection on morphological and on life-history/phenological traits were significantly different: on average, selection on morphology was stronger than selection on phenology/life history. Similarly, the values of |beta| for selection via aspects of survival, fecundity, and mating success were significantly different: on average, selection on mating success was stronger than on survival. Comparisons of estimated linear selection gradients and differentials suggest that indirect components of phenotypic selection were usually modest relative to direct components. The absolute values of quadratic selection gradients |gamma| were exponentially distributed with an overall median of only 0.10, suggesting that quadratic selection is typically quite weak. The distribution of gamma values was symmetric about 0, providing no evidence that stabilizing selection is stronger or more common than disruptive selection in nature.

Usefulness of single nucleotide polymorphism data for estimating population parameters.

Single nucleotide polymorphism (SNP) data can be used for parameter estimation via maximum likelihood methods as long as the way in which the SNPs were determined is known, so that an appropriate likelihood formula can be constructed. We present such likelihoods for several sampling methods. As a test of these approaches, we consider use of SNPs to estimate the parameter Theta = 4N(e)micro (the scaled product of effective population size and per-site mutation rate), which is related to the branch lengths of the reconstructed genealogy. With infinite amounts of data, ML models using SNP data are expected to produce consistent estimates of Theta. With finite amounts of data the estimates are accurate when Theta is high, but tend to be biased upward when Theta is low. If recombination is present and not allowed for in the analysis, the results are additionally biased upward, but this effect can be removed by incorporating recombination into the analysis. SNPs defined as sites that are polymorphic in the actual sample under consideration (sample SNPs) are somewhat more accurate for estimation of Theta than SNPs defined by their polymorphism in a panel chosen from the same population (panel SNPs). Misrepresenting panel SNPs as sample SNPs leads to large errors in the maximum likelihood estimate of Theta. Researchers collecting SNPs should collect and preserve information about the method of ascertainment so that the data can be accurately analyzed.

Perspective: gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies.

Molecular methods as applied to the biogeography of single species (phylogeography) or multiple codistributed species (comparative phylogeography) have been productively and extensively used to elucidate common historical features in the diversification of the Earth's biota. However, only recently have methods for estimating population divergence times or their confidence limits while taking into account the critical effects of genetic polymorphism in ancestral species become available, and earlier methods for doing so are underutilized. We review models that address the crucial distinction between the gene divergence, the parameter that is typically recovered in molecular phylogeographic studies, and the population divergence, which is in most cases the parameter of interest and will almost always postdate the gene divergence. Assuming that population sizes of ancestral species are distributed similarly to those of extant species, we show that phylogeographic studies in vertebrates suggest that divergence of alleles in ancestral species can comprise from less than 10% to over 50% of the total divergence between sister species, suggesting that the problem of ancestral polymorphism in dating population divergence can be substantial. The variance in the number of substitutions (among loci for a given species or among species for a given gene) resulting from the stochastic nature of DNA change is generally smaller than the variance due to substitutions along allelic lines whose coalescence times vary due to genetic drift in the ancestral population. Whereas the former variance can be reduced by further DNA sequencing at a single locus, the latter cannot. Contrary to phylogeographic intuition, dating population divergence times when allelic lines have achieved reciprocal monophyly is in some ways more challenging than when allelic lines have not achieved monophyly, because in the former case critical data on ancestral population size provided by residual ancestral polymorphism is lost. In the former case differences in coalescence time between species pairs can in principle be explained entirely by differences in ancestral population size without resorting to explanations involving differences in divergence time. Furthermore, the confidence limits on population divergence times are severely underestimated when those for number of substitutions per site in the DNA sequences examined are used as a proxy. This uncertainty highlights the importance of multilocus data in estimating population divergence times; multilocus data can in principle distinguish differences in coalescence time (T) resulting from differences in population divergence time and differences in T due to differences in ancestral population sizes and will reduce the confidence limits on the estimates. We analyze the contribution of ancestral population size (theta) to T and the effect of uncertainty in theta on estimates of population divergence (tau) for single loci under reciprocal monophyly using a simple Bayesian extension of Takahata and Satta's and Yang's recent coalescent methods. The confidence limits on tau decrease when the range over which ancestral population size theta is assumed to be distributed decreases and when tau increases; they generally exclude zero when tau/(4Ne) > 1. We also apply a maximum-likelihood method to several single and multilocus data sets. With multilocus data, the criterion for excluding tau = 0 is roughly that l tau/(4Ne) > 1, where l is the number of loci. Our analyses corroborate recent suggestions that increasing the number of loci is critical to decreasing the uncertainty in estimates of population divergence time.

Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach.

A new method for the estimation of migration rates and effective population sizes is described. It uses a maximum-likelihood framework based on coalescence theory. The parameters are estimated by Metropolis-Hastings importance sampling. In a two-population model this method estimates four parameters: the effective population size and the immigration rate for each population relative to the mutation rate. Summarizing over loci can be done by assuming either that the mutation rate is the same for all loci or that the mutation rates are gamma distributed among loci but the same for all sites of a locus. The estimates are as good as or better than those from an optimized FST-based measure. The program is available on the World Wide Web at http://evolution.genetics. washington.edu/lamarc.html/.

Spontaneous heterosis in larval life-history traits of hemiclonal frog hybrids.

European water frog hybrids Rana esculenta (Rana ridibunda x Rana lessonae) reproduce hemiclonally, transmitting only their ridibunda genome to gametes. We compared fitness-related larval life-history traits of natural R. esculenta from Poland with those of the two sympatric parental species and of newly generated F1 hybrids. Compared with either parental species, F1 hybrid offspring had higher survival, higher early growth rates, a more advanced developmental stage by day 49, and earlier metamorphosis, but similar mass at metamorphosis. R. esculenta from natural lineages had trait values intermediate between those of F1 offspring and of the two parental species. The data support earlier observations on natural R. esculenta that had faster larval growth, earlier metamorphosis, and higher resistance to hypoxic conditions compared with either parental species. Observing larval heterosis in F1 hybrids in survival, growth rate, and time to metamorphosis, however, at an even higher degree than in hybrids from natural lineages, demonstrates that heterosis is spontaneous and results from hybridity per se rather than from subsequent interclonal selection; in natural lineages the effects of hybridity and of clonal history are confounded. This is compelling evidence for spontaneous heterosis in hybrid clonals. Results on hemiclonal fish hybrids (Poeciliopsis) showed no spontaneous heterosis; thus, our frog data are not applicable to all hybrid clonals. Our data do show, however, that heterosis is an important potential source for the extensively observed ecological success of hybrid clonals. We suggest that heterosis and interclonal selection together shape fitness of natural R. esculenta lineages.

Mitochondrial DNA reveals formation of nonhybrid frogs by natural matings between hemiclonal hybrids.

The European water frog Rana esculenta (RL), a natural hybrid between R. ridibunda (RR) and R. lessonae (LL), reproduces by hybridogenesis: haploid gametes usually contain an intact chromosome set of R. ridibunda (R); the lessonae nuclear genome (L) is lost from the germ line. Hybridity is restored in the next generation, via fertilization by syntopic R. lessonae. Matings between two hybrids (RL x RL) usually give inviable R. ridibunda (RR) progeny. The adult R. ridibunda subpopulation of Trubeschloo, a gravel pit in northern Switzerland, consists only of females. Fragment patterns for mitochondrial DNA (mtDNA) of these R. ridibunda were identical with those of syntopic R. esculenta and of local populations of R. lessonae; they differed from the patterns in eastern European populations of R. lessonae and of R. ridibunda mtDNAs (3.7% and 9.3% estimated sequence divergence, respectively). In contrast, mtDNAs of two R. ridibunda from an introduced Swiss population with both sexes, although different (2.7% divergence) from each other, were typical R. ridibunda rather than R. lessonae mtDNAs. These data, together with unisexuality, demonstrate conclusively that the all-female R. ridibunda population at Trubeschloo originated from matings between two R. esculenta. The formation of independently reproducing R. ridibunda populations via such hybrid x hybrid matings is precluded because progeny of these matings are unisexual. Recombination in the regenerated fertile R. ridibunda females, followed by matings with R. lessonae, nevertheless provides a mechanism for meiotic reshuffling of genetic material in ridibunda haplotypes that is not typically available in hemiclonal lineages.