Gene Genealogy and Variance of Interpopulational Nucleotide Differences

A mathematical theory is developed for computing the probability that m genes sampled from one population (species) and n genes sampled from another are derived from l genes that existed at the time of population splitting. The expected time of divergence between the two most closely related genes sampled from two different populations and the time of divergence (coalescence) of all genes sampled are studied by using this theory. It is shown that the time of divergence between the two most closely related genes can be used as an approximate estimate of the time of population splitting (T) only when T identical to t/(2N) is small, where t and N are the number of generations and the effective population size, respectively. The variance of Nei and Li's estimate (d) of the number of net nucleotide differences between two populations is also studied. It is shown that the standard error (Sd) of d is larger than the mean when T is small (T much less than 1). In this case, Sd is reduced considerably by increasing sample size. When T is large (T greater than 1), however, a large proportion of the variance of d is caused by stochastic factors, and increase in the sample size does not help to reduce Sd. To reduce the stochastic variance of d, one must use data from many independent unlinked gene loci.     

Gene Genealogies within Mutant Allelic Classes

A coalescent theory of the gene genealogy within an allelic class that arises by a unique mutational event is developed and analyzed. To interpret this theory it was necessary to expand on existing theory for populations of varying size. Two features of the gene genealogy--the average pairwise distance and the total tree length--within the mutant class and within the nonmutant class are found. An index, I, is proposed that describes the extent to which a genealogy is similar to one from a population of constant size (for which I = 0) or to a star genealogy (for which I = 1). The value of I is positive in growing populations and is generally positive for the gene genealogy for the mutant class. The value of I is negative for a population decreasing in size and for the nonmutant class, if the mutant arose recently. The results are discussed in the context of the infinite sites model of mutation, which is appropriate for nucleotide sequence data, and the generalized stepwise mutation model, which is appropriate for microsatellite loci. The same genealogical methods are used to find the probability of at least one recombination event between a nucleotide that defines an allelic class and a marker at a nearby linked site.     

A Continuous Method for Gene Flow

Most modern population genetics inference methods are based on the coalescence framework. Methods that allow estimating parameters of structured populations commonly insert migration events into the genealogies. For these methods the calculation of the coalescence probability density of a genealogy requires a product over all time periods between events. Data sets that contain populations with high rates of gene flow among them require an enormous number of calculations. A new method, transition probability-structured coalescence (TPSC), replaces the discrete migration events with probability statements. Because the speed of calculation is independent of the amount of gene flow, this method allows calculating the coalescence densities efficiently. The current implementation of TPSC uses an approximation simplifying the interaction among lineages. Simulations and coverage comparisons of TPSC vs. MIGRATE show that TPSC allows estimation of high migration rates more precisely, but because of the approximation the estimation of low migration rates is biased. The implementation of TPSC into programs that calculate quantities on phylogenetic tree structures is straightforward, so the TPSC approach will facilitate more general inferences in many computer programs.     

Mitochondrial DNA phylogeography of Acrossocheilus paradoxus (Cyprinidae) in Taiwan

Nucleotide sequences of 3' end of the cytochrome b gene, tRNA genes, D-loop control region, and the 5' end of the 12S rRNA of mitochondrial DNA (mtDNA) were used to assess the genetic and phylogeographic structure of Acrossocheilus paradoxus populations, a Cyprinidae fish of Taiwan. A hierarchical examination of populations in 12 major streams from three geographical regions using an analysis of molecular variance (AMOVA) indicates high genetic differentiation both among populations (PhiST = 0.511, P < 0.001) and among regions (PhiCT = 0.368, P < 0.001). Limited migration largely contributed to the genetic differentiation. High nucleotide diversity (1.13%) and haplotype diversity (0.80%) were detected among populations. The degree of genetic differentiation was correlated with geographical distance between populations, a result consistent with the one-dimensional stepping stone models. A neighbour-joining tree recovered by (DAMBE) supports the pattern of isolation by distance and reveals a closer relationship between populations of the central and southern regions. A minimum spanning network based on nucleotide substitutions reflected migration routes from populations of the central region to the northern and southern regions, respectively. Postglacial colonization and expansion can explain the phylogeographical pattern. Single and ancient migration events may have allowed the northern region to attain the monophyly of mtDNA alleles. In contrast, most populations within geographical regions are either paraphyletic or polyphyletic due to the relatively shorter time period for coalescence. Both low haplotype number and genetic variability suggest a bottleneck event in the Chingmei population of northern Taiwan. Based on coalescence theory, the monophyly of the Tungkang population of the southern region may be associated with a founder event.     

Multiple merger gene genealogies in two species: Monophyly, paraphyly, and polyphyly for two examples of Lambda coalescents

Probabilities of monophyly, paraphyly, and polyphyly of two-species gene genealogies are computed for modest sample sizes and compared for two different Λ coalescent processes. Coalescent processes belonging to the Λ coalescent family admit asynchronous multiple mergers of active ancestral lineages. Assigning a timescale to the time of divergence becomes a central issue when different populations have different coalescent processes running on different timescales. Clade probabilities in single populations are also computed, which can be useful for testing for taxonomic distinctiveness of an observed set of monophyletic lineages. The coalescence rates of multiple merger coalescent processes are functions of coalescent parameters. The effect of coalescent parameters on the probabilities studied depends on the coalescent process, and if the population is ancestral or derived. The probability of reciprocal monophyly tends to be somewhat lower, when associated with a Λ coalescent, under the null hypothesis that two groups come from the same population. However, even for fairly recent divergence times, the probability of monophyly tends to be higher as a function of the number of generations for coalescent processes that admit multiple mergers, and is sensitive to the parameter of one of the example processes.     

Phylogenetic relationships in the genus mus,based on paternally,maternally, and biparentally inherited characters

Several species in the rodent genus Mus are used as model research organisms, but comparative studies of these mice have been hampered by the lack of a well-supported phylogeny. We used DNA sequences from six genes representing paternally, maternally, and biparentally inherited regions of the genome to infer phylogenetic relationships among 10 species of Mus commonly used in laboratory research. Our sample included seven species from the subgenus Mus; one species each from the subgenera Pyromys, Coelomys, and Nannomys; and representatives from three additional murine genera, which served as outgroups in the phylogenetic analyses. Although each of the six genes yielded a unique phylogeny, several clades were supported by four or more gene trees. Nodes that conflicted between trees were generally characterized by weak support for one or both of the alternative topologies, thus providing no compelling evidence that any individual gene, or part of the genome, was misleading with respect to the evolutionary history of these mice. Analysis of the combined data resulted in a fully resolved tree that strongly supports monophyly of the genus Mus, monophyly of the subgenus Mus, division of the subgenus Mus into Palearctic (M. musculus, M. macedonicus, M. spicilegus, and M. spretus) and Asian (M. cervicolor, M. cookii, and M. caroli) clades, monophyly of the house mice (M. m. musculus, "M. m. molossinus," M. m. castaneus, and M. m. domesticus), and a sister-group relationship between M. macedonicus and M. spicilegus. Other clades that were strongly supported by one or more gene partitions were not strongly supported by the combined data. This appears to reflect a localized homoplasy in one partition obscuring the phylogenetic signal from another, rather than differences in gene or genome histories.     

The probability of preservation of a newly arisen gene duplicate.

Newly emerging data from genome sequencing projects suggest that gene duplication, often accompanied by genetic map changes, is a common and ongoing feature of all genomes. This raises the possibility that differential expansion/contraction of various genomic sequences may be just as important a mechanism of phenotypic evolution as changes at the nucleotide level. However, the population-genetic mechanisms responsible for the success vs. failure of newly arisen gene duplicates are poorly understood. We examine the influence of various aspects of gene structure, mutation rates, degree of linkage, and population size (N) on the joint fate of a newly arisen duplicate gene and its ancestral locus. Unless there is active selection against duplicate genes, the probability of permanent establishment of such genes is usually no less than 1/(4N) (half of the neutral expectation), and it can be orders of magnitude greater if neofunctionalizing mutations are common. The probability of a map change (reassignment of a key function of an ancestral locus to a new chromosomal location) induced by a newly arisen duplicate is also generally >1/(4N) for unlinked duplicates, suggesting that recurrent gene duplication and alternative silencing may be a common mechanism for generating microchromosomal rearrangements responsible for postreproductive isolating barriers among species. Relative to subfunctionalization, neofunctionalization is expected to become a progressively more important mechanism of duplicate-gene preservation in populations with increasing size. However, even in large populations, the probability of neofunctionalization scales only with the square of the selective advantage. Tight linkage also influences the probability of duplicate-gene preservation, increasing the probability of subfunctionalization but decreasing the probability of neofunctionalization.     

PREDICTING NUCLEAR GENE COALESCENCE FROM MITOCHONDRIAL DATA: THE THREE-TIMES RULE

Abstract.— Coalescence theory predicts when genetic drift at nuclear loci will result in fixation of sequence differences to produce monophyletic gene trees. However, the theory is difficult to apply to particular taxa because it hinges on genetically effective population size, which is generally unknown. Neutral theory also predicts that evolution of monophyly will be four times slower in nuclear than in mitochondrial genes primarily because genetic drift is slower at nuclear loci. Variation in mitochondrial DNA (mtDNA) within and between species has been studied extensively, but can these mtDNA data be used to predict coalescence in nuclear loci? Comparison of neutral theories of coalescence of mitochondrial and nuclear loci suggests a simple rule of thumb. The “three‐times rule” states that, on average, most nuclear loci will be monophyletic when the branch length leading to the mtDNA sequences of a species is three times longer than the average mtDNA sequence diversity observed within that species. A test using mitochondrial and nuclear intron data from seven species of whales and dolphins suggests general agreement with predictions of the three‐times rule. We define the coalescence ratio as the mitochondrial branch length for a species divided by intraspecific mtDNA diversity. We show that species with high coalescence ratios show nuclear monophyly, whereas species with low ratios have polyphyletic nuclear gene trees. As expected, species with intermediate coalescence ratios show a variety of patterns. Especially at very high or low coalescence ratios, the three‐times rule predicts nuclear gene patterns that can help detect the action of selection. The three‐times rule may be useful as an empirical benchmark for evaluating evolutionary processes occurring at multiple loci.     

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.     

The distribution of the coalescence time and the number of pairwise nucleotide differences in a model of population divergence or speciation with an initial period of gene flow

This paper is concerned with a model of “isolation with an initial period of migration”, where a panmictic ancestral population split into n descendant populations which exchanged migrants symmetrically at a constant rate for a period of time and subsequently became completely isolated. In the limit as the population split occurred an infinitely long time ago, the model becomes an “isolation after migration” model, describing completely isolated descendant populations which arose from a subdivided ancestral population. The probability density function of the coalescence time of a pair of genes and the probability distribution of the number of pairwise nucleotide differences are derived for both models. Whilst these are theoretical results of interest in their own right, they also give an exact analytical expression for the likelihood, for data consisting of the numbers of nucleotide differences between pairs of DNA sequences where each pair is at a different, independent locus. The behaviour of the distribution of the number of pairwise nucleotide differences under these models is illustrated and compared to the corresponding distributions under the “isolation with migration” and “complete isolation” models. It is shown that the distribution of the number of nucleotide differences between a pair of DNA sequences from different descendant populations in the model of “isolation with an initial period of migration” can be quite different from that under the “isolation with migration model”, even if the average migration rate over time (and hence the total number of migrants) is the same in both scenarios. It is also illustrated how the results can be extended to other demographic scenarios that can be described by a combination of isolated panmictic populations and “symmetric island” models.     

Gene Genealogy in Three Related Populations: Consistency Probability between Gene and Population Trees

A genealogical relationship among genes at a locus (gene tree) sampled from three related populations was examined with special reference to population relatedness (population tree). A phylogenetically informative event in a gene tree constructed from nucleotide differences consists of interspecific coalescences of genes in each of which two genes sampled from different populations are descended from a common ancestor. The consistency probability between gene and population trees in which they are topologically identical was formulated in terms of interspecific coalescences. It was found that the consistency probability thus derived substantially increases as the sample size of genes increases, unless the divergence time of populations is very long compared to population sizes. Hence, there are cases where large samples at a locus are very useful in inferring a population tree.     

Mathematical consequences of the genealogical species concept

A genealogical species is defined as a basal group of organisms whose members are all more closely related to each other than they are to any organisms outside the group ("exclusivity"), and which contains no exclusive group within it. In practice, a pair of species is so defined when phylogenies of alleles from a sample of loci shows them to be reciprocally monophyletic at all or some specified fraction of the loci. We investigate the length of time it takes to attain this status when an ancestral population divides into two descendant populations of equal size with no gene exchange, and when genetic drift and mutation are the only evolutionary forces operating. The number of loci used has a substantial effect on the probability of observing reciprocal monophyly at different times after population separation, with very long times needed to observe complete reciprocal monophyly for a large number of loci. In contrast, the number of alleles sampled per locus has a relatively small effect on the probability of reciprocal monophyly. Because a single mitochondrial or chloroplast locus becomes reciprocally monophyletic much faster than does a single nuclear locus, it is not advisable to use mitochondrial and chloroplast DNA to recognize genealogical species for long periods after population divergence. Using a weaker criterion of assigning genealogical species status when more than 50% of sampled nuclear loci show reciprocal monophyly, genealogical species status depends much less on the number of sampled loci, and is attained at roughly 4-7 N generations after populations are isolated, where N is the historically effective population size of each descendant. If genealogical species status is defined as more than 95% of sampled nuclear loci showing reciprocal monophyly, this status is attained after roughly 9-12 N generations.     

The complete mitochondrial genomes of largemouth bass of the northern subspecies (Micropterus salmoides salmoides) and Florida subspecies (Micropterus salmoides floridanus) and their applications in the identification of largemouth bass species

The largemouth bass belongs to the family Centrarchidae, which includes two subspecies: the northern subspecies, Micropterus salmoides salmoides, and the Florida subspecies, Micropterus salmoides floridanus. In this study, the complete mitochondrial genomes of the two subspecies were sequenced, and their genetic differences were identified. The mitogenomes of M. s. salmoides and M. s. floridanus are 16,486 and 16,479?bp in length, respectively. The two subspecies consisted of 37 genes (13 protein-coding genes, 2 ribosomal RNA, and 22 transfer RNA), which are typical for vertebrate mtDNA. Phylogenetic analysis provided statistical support for the monophyly of the family Centrarchidae. Comparison of the two subspecies' mitogenomes revealed a relatively high number (450) of single nucleotide polymorphisms (SNPs) in protein-coding genes. We characterized SNPs in the partial cytochrome c oxidase subunit 1 gene of different individuals from three cultured populations, one wild northern subspecies population, and one wild Florida subspecies population. Twenty-eight SNPs were fixed with alternative nucleotides in the two subspecies, which could be used for differentiating them. Based on this gene, phylogenetic tree and genetic distance analyses supported that cultured largemouth bass in China belongs to the northern subspecies.     

Inbreeding coefficients and coalescence times.

This paper describes the relationship between probabilities of identity by descent and the distribution of coalescence times. By using the relationship between coalescence times and identity probabilities, it is possible to extend existing results for inbreeding coefficients in regular systems of mating to find the distribution of coalescence times and the mean coalescence times. It is also possible to express Sewall Wright's FST as the ratio of average coalescence times of different pairs of genes. That simplifies the analysis of models of subdivided populations because the average coalescence time can be found by computing separately the time it takes for two genes to enter a single subpopulation and time it takes for two genes in the same subpopulation to coalesce. The first time depends only on the migration matrix and the second time depends only on the total number of individuals in the population. This approach is used to find FST in the finite island model and in one- and two-dimensional stepping-stone models. It is also used to find the rate of approach of FST to its equilibrium value. These results are discussed in terms of different measures of genetic distance. It is proposed that, for the purposes of describing the amount of gene flow among local populations, the effective migration rate between pairs of local populations, M, which is the migration rate that would be estimated for those two populations if they were actually in an island model, provides a simple and useful measure of genetic similarity that can be defined for either allozyme or DNA sequence data.     

Effective population size,genetic diversity,and coalescence time in subdivided populations

A formula for the effective population size for the finite island model of subdivided populations is derived. The formula indicates that the effective size can be substantially greater than the actual number of individuals in the entire population when the migration rate among subpopulations is small. It is shown that the mean nucleotide diversity, coalescence time, and heterozygosity for genes sampled from the entire population can be predicted fairly well from the theory for randomly mating populations if the effective population size for the finite island model is used.     

Mitochondrial 16S DNA variation in populations of Triatoma infestans from Argentina

Abstract.  Variation in the mtDNA 16S ribosomal RNA gene in populations of Triatoma infestans (Klug) was surveyed. DNA sequence comparisons yielded 18 haplotypes among 130 individuals from 16 localities that represent a large proportion of the range of T. infestans in Argentina. The most common genotype in all populations was found in 76.9% of individuals and two other haplotypes were shared among different populations. The remaining 15 haplotypes were present exclusively in one of the populations, suggesting currently low levels of genetic exchange. Analysis of mtDNA 16S sequences uncovered substantial genetic variation among T. infestans populations. Haplotype and nucleotide diversities varied among populations, from 0% to 0.84% and 0% to 0.29%, respectively. It appears that this locus has a low mutation rate. Uncorrected pairwise differences of T. infestans haplotypes ranged from 0% to 1.2%. The molecular phylogeny supported the monophyly of T. infestans haplotypes and clustered two different pairs of haplotypes with a moderate degree of bootstrap support (∼ 60%). Mitochondrial DNA phylogeographic differentiation was not evident, suggesting a recent rapid spread of the species. Analysis of molecular variance showed hierarchical structure in the data. Considerably less variation was found among T. infestans populations from the northwest and northeast regions than among those belonging to the central area. Such a lack of variation may be indicative of one or more past population bottlenecks.     

Relationship between DNA Polymorphism and Fixation Time

When there is no recombination among nucleotide sites in DNA sequences, DNA polymorphism and fixation of mutants at nucleotide sites are mutually related. Using the method of gene genealogy, the relationship between the DNA polymorphism and the fixation of mutant nucleotide was quantitatively investigated under the assumption that mutants are selectively neutral, that there is no recombination among nucleotide sites, and that the population is a random mating population with N diploid individuals. The results obtained indicate that the expected number of nucleotide differences between two DNA sequences randomly sampled from the population is 42% less when a mutant at a particular nucleotide site reaches fixation than at a random time, and that heterozygosity is also expected to be less when fixation takes place than at a random time, but the amount of reduction depends on the value of 4Nv in this case, where v is the mutation rate per DNA sequence per generation. The formula for obtaining the expected number of nucleotide differences between the two DNA sequences for a given fixation time is also derived, and indicates that, even when it takes a large number of generations for a mutant to reach fixation, this number is 33% less than at a random time. The computer simulation conducted suggests that the expected number of nucleotide differences between the two DNA sequences at the time when an advantageous mutant becomes fixed is essentially the same as that of neutral mutant if the fixation time is the same. The effect of recombination on the amount of DNA polymorphism was also investigated by using computer simulation.     

Phylogenetic relationships in European erpobdellid leeches (Hirudinea: Erpobdellidae) inferred from restriction-site data of the 18S ribosomal gene and ITS2 region

The PCR-amplified 185 rRNA gene and ITS2 region were used for a restriction-site analysis to infer phylogenetic relationships among European representatives of the leech family Erpobdellidae and to test the phylogenetic reliability of the morphological characters used in this family till now. Fragment patterns produced by 11 restriction endonucleases in the 18S rRNA gene and by nine in the ITS2 region showed a higher interspecific variation in the latter than in the former genome component. No intraspecific variation could be detected, even among geographically very distant populations, except for the two supposed subspecies Trocheta b. bykowskii and T. b. krasense . The estimated overall nucleotide divergence (d) was lower in the genus Erpobdella (d = 0.015–0.061) than in the genera Dina and Trocheta (d = 0.049–0.090). E. testacea and E. monostriata do not seem to be as closely related to each other as previously thought; their pairwise nucleotide divergence (0.061) is highest within the genus Erpobdella . The nucleotide divergence between T. b. bykowskii and T. b. krasense (0.090), which greatly resemble each other morphologically, was surprisingly high. Phylogenetic trees were inferred via neighbour-joining and maximum parsimony methods. The trees support the monophyly of Erpobdellidae and the subdivision into two groups according to the pattern of annulation (the Erpobdella clade and the Dina-Trocheta clade). Further subdivision based on annulation patterns and other morphological characters is not supported. There was no molecular evidence for the monophyly of the genus Dina and the monophyly of the genus Trocheta .     

Phylogeography of Hemibarbus labeo (Cyprinidae): secondary contact of ancient lineages of mtDNA

Lin, C.-J., Lin, H.-D., Wang, J.-P., Chao, S.-C. & Chiang, T.-Y. (2007). Phylogeography of Hemibarbus labeo (Cyprinidae): secondary contact of ancient lineages of mtDNA. —Zoologica Scripta, 39, 23–35. Nucleotide sequences of the D-loop control region of mtDNA were used to assess the genetic structure and phylogeography of Hemibarbus labeo in Mainland China and Taiwan. A hierarchical analysis of molecular variance of populations in 11 major streams from three geographical regions revealed significant structuring among populations and geographical regions. High levels of nucleotide diversity (π = 1.88%) and haplotype diversity (h = 96 ± 0.009%) suggest a large effective population size. A maximum likelihood tree based on mtDNA variation identified two ancient mtDNA lineages, which split approximately 3.39 million years ago. Overlapping distribution of the major lineages displayed low correspondence with geographical regions and reflects a scenario of secondary mergence after long isolation. Gene genealogy further revealed a unidirectional migration. Nevertheless, there existed a phylogeographical structure that mostly agreed with a biogeographical hypothesis. That is, within each lineage, a close phylogeny between populations of the River-Campagna and East-Pacific regions was supported by the mtDNA gene genealogy, although monophyly of each geographical region was not supported. The degree of genetic differentiation was correlated with geographical distances between populations, displaying a pattern of ‘isolation by distance’. Gene genealogy of mtDNA revealed that Yangtzejiang population may act as a divergence centre of H. labeo. In addition, Taiwan population was colonized via a recent a founder event, likely from population Yangtzejiang River about 131 000 years before present. Low haplotype number and genetic variability also suggested possible bottleneck events in the Rongjiang and Dazhangjiang populations.     

Huge populations and old species of Costa Rican and Panamanian dirt frogs inferred from mitochondrial and nuclear gene sequences

Molecular genetic data were used to investigate population sizes and ages of Eleutherodactylus (Anura: Leptodactylidae), a species-rich group of small leaf-litter frogs endemic to Central America. Population genetic structure and divergence was investigated for four closely related species surveyed across nine localities in Costa Rica and Panama. DNA sequence data were collected from a mitochondrial gene (ND2) and a nuclear gene (c-myc). Phylogenetic analyses yielded concordant results between loci, with reciprocal monophyly of mitochondrial DNA haplotypes for all species and of c-myc haplotypes for three of the four species. Estimates of genetic differentiation among populations (FST) based upon mitochondrial data were always higher than nuclear-based FST estimates, even after correcting for the expected fourfold lower effective population size (Ne) of the mitochondrial genome. Comparing within-population variation and the relative mutation rates of the two genes revealed that the Ne of the mitochondrial genome was 15-fold lower than the estimate of the nuclear genome based on c-myc. Nuclear FST estimates were approximately 0 for the most proximal pairs of populations, but ranged from 0.5 to 1.0 for all other pairs, even within the same nominal species. The nuclear locus yielded estimates of Ne within localities on the order of 105. This value is two to three orders of magnitude larger than any previous Ne estimate from frogs, but is nonetheless consistent with published demographic data. Applying a molecular clock model suggested that morphologically indistinguishable populations within one species may be 107 years old. These results demonstrate that even a geologically young and dynamic region of the tropics can support very old lineages that harbour great levels of genetic diversity within populations. The association of high nucleotide diversity within populations, large divergence between populations, and high species diversity is also discussed in light of neutral community models.