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

Comparative phylogeography of direct-developing frogs (Anura: Craugastoridae: Pristimantis) in the southern Andes of Colombia

The Andes of South America hosts perhaps the highest amphibian species diversity in the world, and a sizable component of that diversity is comprised of direct-developing frogs of the genus Pristimantis (Anura: Craugastoridae). In order to better understand the initial stages of species formation in these frogs, this study quantified local-scale spatial genetic structuring in three species of Pristimantis. DNA sequences of two mitochondrial gene fragments (16S and COI) were obtained from P. brevifrons, P. palmeri and P. jubatus at different locations in the Cordillera Occidental. We found high levels of genetic diversity in the three species, with highly structured populations (as measured by F(ST)) in P. brevifrons and P. palmeri while P. jubatus showed panmixia. Large effective population sizes, inferred from the high levels of genetic diversity, were found in the three species and two highly divergent lineages were detected within P. jubatus and P. palmeri. Estimated divergence times among populations within P. brevifrons and P. palmeri coincide with the Pleistocene, perhaps due to similar responses to climatic cycling or recent geological history. Such insights have important implications for linking alpha and beta diversity, suggesting regional scale patterns may be associated with local scale processes in promoting differentiation among populations in the Andes.

Hybridization and massive mtDNA unidirectional introgression between the closely related Neotropical toads Rhinella marina and R. schneideri inferred from mtDNA and nuclear markers

Background

The classical perspective that interspecific hybridization in animals is rare has been changing due to a growing list of empirical examples showing the occurrence of gene flow between closely related species. Using sequence data from cyt b mitochondrial gene and three intron nuclear genes (RPL9, c-myc, and RPL3) we investigated patterns of nucleotide polymorphism and divergence between two closely related toad species R. marina and R. schneideri. By comparing levels of differentiation at nuclear and mtDNA levels we were able to describe patterns of introgression and infer the history of hybridization between these species.

Results

All nuclear loci are essentially concordant in revealing two well differentiated groups of haplotypes, corresponding to the morphologically-defined species R. marina and R. schneideri. Mitochondrial DNA analysis also revealed two well-differentiated groups of haplotypes but, in stark contrast with the nuclear genealogies, all R. schneideri sequences are clustered with sequences of R. marina from the right Amazon bank (RAB), while R. marina sequences from the left Amazon bank (LAB) are monophyletic. An Isolation-with-Migration (IM) analysis using nuclear data showed that R. marina and R. schneideri diverged at ≈ 1.69 Myr (early Pleistocene), while R. marina populations from LAB and RAB diverged at ≈ 0.33 Myr (middle Pleistocene). This time of divergence is not consistent with the split between LAB and RAB populations obtained with mtDNA data (≈ 1.59 Myr), which is notably similar to the estimate obtained with nuclear genes between R. marina and R. schneideri. Coalescent simulations of mtDNA phylogeny under the speciation history inferred from nuclear genes rejected the hypothesis of incomplete lineage sorting to explain the conflicting signal between mtDNA and nuclear-based phylogenies.

Conclusions

The cytonuclear discordance seems to reflect the occurrence of interspecific hybridization between these two closely related toad species. Overall, our results suggest a phenomenon of extensive mtDNA unidirectional introgression from the previously occurring R. schneideri into the invading R. marina. We hypothesize that climatic-induced range shifts during the Pleistocene/Holocene may have played an important role in the observed patterns of introgression.

Genetic diversity and biogeographic history inform future conservation management strategies for the rare sunset frog (Spicospina flammocaerulea)

Outlining the distribution of genetic variation, patterns of gene flow and clarifying the biogeographic processes underlying population history are critical components of a comprehensive conservation strategy for endangered or vulnerable species. We provide this information for the vulnerable sunset frog (Spicospina flammocaerulea) using a comprehensive genetic dataset (ND2) with samples from 17 of 22 geographic localities where this species has been found. From genetic, biogeographic and coalescent-based analyses, we document the existing genetic variation, likely movement patterns and explore the biogeographic history of S. flammocaerulea. While catchment-based genetic variation is well documented in other high-rainfall taxa in south-western Australia, a much more complex scenario including dispersal across ridge lines between catchments better explains the distribution of genetic variation and observed patterns of gene flow in S. flammocaerulea. The population history of S. flammocaerulea is strongly indicative of recent population contraction and expansion, which may be related to late Pleistocene climate fluctuations. This suggests that this species can adapt or move in response to fluctuating climates provided suitable habitats or expansion areas are available. However, like many other endemic taxa with limited geographic ranges in south-western Australia, the potential to shift distributions is hampered by being land-locked within an agricultural landscape, limiting management options in the face of climate change.

Late Pleistocene environmental changes lead to unstable demography and population divergence of Anopheles albimanus in the northern Neotropics

We investigated the historical demography of Anopheles albimanus using mosquitoes from five countries and three different DNA regions, the mitochondrial cytochrome oxidase subunit I gene (COI), the single copy nuclear white gene and the ribosomal internal transcribed spacer two (ITS2). All the molecular markers supported the taxonomic status of a single species of An. albimanus. Furthermore, agreement between the COI and the white genes suggested a scenario of Pleistocene geographic fragmentation (i.e., population contraction) and subsequent range expansion across southern Central America.

Ecological and genetic divergence between two lineages of middle American túngara frogs Physalaemus (= Engystomops) pustulosus

BACKGROUND

Uncovering how populations of a species differ genetically and ecologically is important for understanding evolutionary processes. Here we combine population genetic methods (microsatellites) with phylogenetic information (mtDNA) to define genetic population clusters of the wide-spread Neotropical túngara frog (Physalaemus pustulosus). We measure gene flow and migration within and between population clusters and compare genetic diversity between population clusters. By applying ecological niche modeling we determine whether the two most divergent genetic groups of the túngara frog (1) inhabit different habitats, and (2) are separated geographically by unsuitable habitat across a gap in the distribution.

RESULTS

Most population structure is captured by dividing all sample localities into two allopatric genetic lineages. The Northern genetic lineage (NW Costa Rica) is genetically homogenous while the Southern lineage (SW Costa Rica and Panama) is sub-divided into three population clusters by both microsatellite and mtDNA analyses. Gene flow is higher within the Northern lineage than within the Southern lineage, perhaps due to increased landscape heterogeneity in the South. Niche modeling reveals differences in suitable habitat between the Northern and Southern lineages: the Northern lineage inhabits dry/pine-oak forests, while the Southern lineage is confined to tropical moist forests. Both lineages seem to have had little movement across the distribution gap, which persisted during the last glacial maximum. The lack of movement was more pronounced for the Southern lineage than for the Northern lineage.

CONCLUSIONS

This study confirms the finding of previous studies that túngara frogs diverged into two allopatric genetic lineages north and south of the gap in the distribution in central Costa Rica several million years ago. The allopatric distribution is attributed to unsuitable habitat and probably other unknown ecological factors present across the distribution gap. Niche conservatism possibly contributes to preventing movements across the gap and gene flow between both groups. Genetic and ecological data indicate that there is the potential for ecological divergence in allopatry between lineages. In this context we discuss whether the Northern and Southern lineages should be recognized as separate species, and we conclude that further studies of pre- and post-zygotic isolation are needed for a final assessment. Identified population clusters should motivate future behavioral and ecological research regarding within-species biodiversity and speciation mechanisms.

Genetic diversification, vicariance, and selection in a polytypic frog

Spatial patterns of heritable phenotypic diversity reflect the relative roles of gene flow and selection in determining geographic variation within a species. We quantified color differentiation and genetic divergence among 20 populations of the red-eyed tree frog (Agalychnis callidryas) in lower Central America. Phylogenetic analyses revealed 5 well-supported mitochondrial DNA clades, and we infer from our phylogeny that geographic barriers have played a large role in structuring populations. Two phenotypic characters varied independently among isolated population groups: Flank coloration distinguished Caribbean from Pacific individuals, whereas leg coloration exhibited a more complex geographic pattern. We detected 3 generalized spatial patterns of genetic and phenotypic diversity: 1) phenotypic differentiation in the presence of historical connectivity, 2) phenotypic uniformity across genetically differentiated regions, and 3) codistribution of genetic and phenotypic characters. These patterns indicate that phenotypic diversification is highly regionalized and can result from spatial variation in localized adaptations, geographic isolation, genetic drift, and/or evolutionary stasis. Although the mode of selection underlying color variation was not the focal objective of this study, we discuss the possible roles of natural and sexual selection in mediating population differentiation. Our study underscores the fact that selection gradients vary across relatively small spatial scales, even in species that occupy relatively homogeneous environments.

Divergence in the face of gene flow: the case of two newts (amphibia: salamandridae)

Understanding the process of divergence requires the quantitative characterization of patterns of gene flow between diverging taxa. New and powerful coalescent-based methods give insight into these processes in unprecedented details by enabling the reconstruction of the temporal distribution of past gene flow. Here, we use sequence variation at eight nuclear markers and mitochondrial DNA (mtDNA) in multiple populations to study diversity, divergence, and gene flow between two subspecies of a salamander, the smooth newt (Lissotriton vulgaris kosswigi and Lissotriton vulgaris vulgaris) in Turkey. The ranges of both subspecies encompass mainly the areas of this important glacial refugial area. Populations in refugia where species have been present for a long time and differentiated in situ should better preserve the record of past gene flow than young populations in postglacial expansion areas. Sequence diversity in both subspecies was substantial (nuclear pi(sil) = 0.69% and 1.31%). We detected long-term demographic stability in these refugial populations with large effective population sizes (N(e)) of the order of 1.5-3 x 10(5) individuals. Gene trees and the isolation with migration (IM) analysis complemented by tests of nested IM models showed that despite deep, pre-Pleistocene divergence of the studied newts, asymmetric introgression from vulgaris to kosswigi has occurred, with signatures of recent gene flow in mtDNA and an anonymous nuclear marker, and evidence for more ancient introgression in nuclear introns. The distribution of migration times raises the intriguing possibility that even the initial divergence may have occurred in the face of gene flow.

Phylogenetic relationships of glassfrogs (Centrolenidae) based on mitochondrial and nuclear genes

Glassfrogs (family Centrolenidae) represent an exceptionally diverse group among Neotropical anurans, but their evolutionary relationships never have been assessed from a molecular perspective. Mitochondrial and nuclear markers were used to develop a novel hypothesis of centrolenid phylogeny. Ingroup sampling included 100 terminals, with 78 (53%) of the named species in the family, representing most of the phenotypic diversity described for the group. Thirty-five species representing taxa traditionally associated with glassfrogs were used as outgroups. Gene sampling consisted of complete or partial sequences of three mitochondrial (12S, 16S, ND1) and three nuclear markers (c-myc exon 2, RAG1, POMC) for a total of approximately 4362bp. Phylogenies were estimated using maximum parsimony, maximum likelihood, and Bayesian analyses for individual genes and combined datasets. The separate analysis of mitochondrial and nuclear datasets allowed us to clarify the relationships within glassfrogs; also, we corroborate the sister-group relationship between Allophryne ruthveni and glassfrogs. The new phylogeny differs significantly from all previous morphology-based hypotheses of relationships, and shows that hypotheses based on few traits are likely to misrepresent evolutionary history. Traits previously hypothesized as unambiguous synapomorphies are shown to be homoplastic, and all genera in the current taxonomy (Centrolene, Cochranella, Hyalinobatrachium, Nymphargus) are found to be poly- or paraphyletic. The new topology implies a South American origin of glassfrogs and reveals allopatric speciation as the most important speciation mechanism. The phylogeny profoundly affects the traditional interpretations of glassfrog taxonomy, character evolution, and biogeography-topics that now require more extensive evaluation in future studies.

A dated phylogeny of the palm tribe Chamaedoreeae supports Eocene dispersal between Africa, North and South America

The palm tribe Chamaedoreeae reaches its higher diversity in Central America, however, its distribution ranges from the north eastern part of Mexico to Bolivia with a disjunction to the Mascarene Islands in the Indian Ocean. The disjunct distribution of Chamaedoreeae is generally considered a result of Gondwana vicariance and extinction from Africa and/or Madagascar. However, latitudinal migrations and their role in shaping the distribution of this tribe in the Americas have been largely overlooked. In this study we used seven plastid and two nuclear DNA regions to investigate the phylogenetic relationships and biogeography of the Chamaedoreeae. The resulting phylogeny fully resolved the generic relationships within the tribe. The exact area of origin of the tribe remains uncertain, but dating analyses indicated an initial diversification of the Chamaedoreeae during the Early Eocene, followed by long distance dispersion to the Mascarene Islands in the late Miocene. The radiation of Hyophorbe could have taking place on islands in the Indian Ocean now submerged, but its former presence in Africa or Madagascar cannot be ruled out. At least two independent migrations between North and South America predating the rise of the Panama isthmus need to be postulated to explain the distribution of Chamaedoreeae, one during the Middle Eocene and a second during the Miocene. Whereas the traditional interpretation of distribution of Chamaedoreeae species assumes a west Gondwana origin of the group, the findings presented in this paper make it equally possible to interpret the group as a primarily boreotropical element.

Phylogeography of the Pygmy Rain Frog (Pristimantis ridens) across the lowland wet forests of isthmian Central America

We used a phylogeographic approach to elucidate the evolutionary history of a lineage of frogs, known as Pristimantis (formerly Eleutherodactylus) ridens (Anura: Brachycephalidae), restricted to the wet forests occurring along the Caribbean versant of isthmian Central America as well as the disjunct wet forest on the Pacific slope of Costa Rica. We placed our phylogeographic study of P. ridens within a larger molecular phylogenetic analysis of Central American Pristimantis. All phylogenetic inferences were based on a 1455 base pair fragment of mitochondrial DNA, containing the complete ND2 gene and five flanking tRNA genes. Our reconstruction of the intraspecific phylogeny of P. ridens yielded a basal trichotomy dating to an estimated 12+ million years ago (Ma), consisting of central Panama, western Panama, and Costa Rica plus Honduras. Thus, the presence of P. ridens appears to predate the completion of the Isthmus 3.1Ma. Using a parametric bootstrap (SOWH) test, we evaluated four a priori zoogeographic hypotheses for the origin and spread of P. ridens. This analysis suggested that the P. ridens populations on the Caribbean versant of Costa Rica were established by Pacific versant ancestors only recently, in contrast to the very old lineages found in Panama. Our results support a model of Miocene colonization, long-term geographic stasis, followed by rapid dispersal across the Caribbean lowlands during the Pliocene or Pleistocene.

Cryptic diversity and deep divergence in an upper Amazonian leaflitter frog, Eleutherodactylus ockendeni

Background

The forests of the upper Amazon basin harbour some of the world's highest anuran species richness, but to date we have only the sparsest understanding of the distribution of genetic diversity within and among species in this region. To quantify region-wide genealogical patterns and to test for the presence of deep intraspecific divergences that have been documented in some other neotropical anurans, we developed a molecular phylogeny of the wide-spread terrestrial leaflitter frog Eleutherodactylus ockendeni (Leptodactylidae) from 13 localities throughout its range in Ecuador using data from two mitochondrial genes (16S and cyt b; 1246 base pairs). We examined the relation between divergence of mtDNA and the nuclear genome, as sampled by five species-specific microsatellite loci, to evaluate indirectly whether lineages are reproductively isolated where they co-occur. Our extensive phylogeographic survey thus assesses the spatial distribution of E. ockendeni genetic diversity across eastern Ecuador.

Results

We identified three distinct and well-supported clades within the Ecuadorean range of E. ockendeni: an uplands clade spanning north to south, a northeastern and central lowlands clade, and a central and southeastern clade, which is basal. Clades are separated by 12% to 15% net corrected p-distance for cytochrome b, with comparatively low sequence divergence within clades. Clades marginally overlap in some geographic areas (e.g., Napo River basin) but are reproductively isolated, evidenced by diagnostic differences in microsatellite PCR amplification profiles or DNA repeat number and coalescent analyses (in MDIV) best modelled without migration. Using Bayesian (BEAST) and net phylogenetic estimates, the Southeastern Clade diverged from the Upland/Lowland clades in the mid-Miocene or late Oligocene. Lowland and Upland clades speciated more recently, in the early or late Miocene.

Conclusion

Our findings uncover previously unsuspected cryptic species diversity within the common leaflitter frog E. ockendeni, with at least three different species in Ecuador. While these clades are clearly geographically circumscribed, they do not coincide with any existing landscape barriers. Divergences are ancient, from the Miocene, before the most dramatic mountain building in the Ecuadorean Andes. Therefore, this diversity is not a product of Pleistocene refuges. Our research coupled with other studies suggests that species richness in the upper Amazon is drastically underestimated by current inventories based on morphospecies.

The role of tropical dry forest as a long-term barrier to dispersal: a comparative phylogeographical analysis of dry forest tolerant and intolerant frogs

We used a comparative phylogeographical approach to investigate the origins of the disjunct wet forest biota of the Golfo Dulce region along the Pacific slope of Costa Rica. This region is isolated by Pacific dry forests north and south and isolated from Caribbean wet forests by mountains. We studied three sympatric lowland frog species in the Craugastor fitzingeri species group that prefer wet forest but differ in their response to dry habitats. In dry forest, C. fitzingeri can survive along streams while C. crassidigitus and C. talamancae are entirely absent. We collected samples from across the ranges of all three species, and obtained mitochondrial DNA sequence data from the COI and cytochrome b genes. We observed significant phylogeographical structure in C. crassidigitus and C. talamancae, but much less in C. fitzingeri, demonstrating that mountain barriers and dry forest habitat have reduced mitochondrial gene flow in the strictly wet-forest species. Additionally, we discovered that the Golfo Dulce and Central Panama populations of C. crassidigitus appear to have diverged in the Pliocene or earlier, suggesting that the dry forest separating these populations is old. Our phylogenetic analysis of 12 of approximately 16 species of the C. fitzingeri species group suggests that the three lowland species are each other's closest relatives. Because of this shared phylogenetic history, we attribute the striking differences in phylogeographical structure to the different ecologies of the frogs. In summary, we find that what appear to be minor differences in the natural history of these three closely related species may profoundly impact the potential for dispersal, range size, and cladogenesis.

Impact of Plio-Pleistocene arid cycling on the population history of a southwestern Australian frog

Southwestern Australia is regarded as a global biodiversity hotspot. The region contains a high number of endemic species, ranging from Gondwanan relicts to much more recently evolved plant and animal species. Myobatrachid frogs are diverse in southwestern Australia, and while we know they have speciated in situ in the southwest, we know little about the temporal and geographical patterning of speciation events. Crinia georgiana is an ideal subject to test hypotheses concerning the effect of climatic history on southwestern Australian anurans, as it is an old lineage with a broad distribution covering the entire region. We compiled an extensive phylogeographical data set based on 1085 bp of the mitochondrial gene ND2 for 68 individuals from 18 sites across the species' range. Two major genetic clades were identified which were largely confined to the high rainfall and southeast coastal biogeographical zones, respectively. The clades appear to have diverged around the Plio-Pleistocene border (1.26-1.72 million years ago), concordant with increasing intensity and frequency of arid climate cycles. Subsequent phylogeographical structure appears to have developed primarily during the Pleistocene climatic fluctuations that also have been integral in generating species diversity in the endemic southwestern Australian flora. Phylogeographical analyses identified several dispersal routes, possible refugial areas within the range of the species and also regions of secondary contact. Dispersal routes identified may now be closed to the species because of habitat destruction and salinity problems in inland regions, posing concerns about the evolutionary potential of the species in light of predicted climate change.

Applying new inter-individual approaches to assess fine-scale population genetic diversity in a neotropical frog, Eleutherodactylus ockendeni

We assess patterns of genetic diversity of a neotropical leaflitter frog, Eleutherodactylus ockendeni, in the upper Amazon of Ecuador without a priori delineation of biological populations and with sufficiently intensive sampling to assess inter-individual patterns. We mapped the location of each collected frog across a 5.4 x 1 km landscape at the Jatun Sacha Biological Station, genotyped 185 individuals using five species-specific DNA microsatellite loci, and sequenced a fragment of mitochondrial cytochrome b for a subset of 51 individuals. The microsatellites were characterized by high allelic diversity and homozygote excess across all loci, suggesting that when pooled the sample is not a panmictic population. We conclude that the lack of panmixia is not attributable to the influence of null alleles or biased sampling of consanguineous family groups. Multiple methods of population cluster analysis, using both Bayesian and maximum likelihood approaches, failed to identify discrete genetic clusters across the sampled area. Using multivariate spatial autocorrelation, kinship coefficients and relatedness coefficients, we identify a continuous isolation by distance population structure, with a first patch size of ca. 260 m and apparently large population sizes. Analysis of mtDNA corroborates the observation of high genetic diversity at fine scales: there are multiple haplotypes, they are non-randomly distributed and a binary haplotype correlogram shows significant spatial genetic autocorrelation. We demonstrate the utility of inter-individual genetic methods and caution against making a priori assumptions about population genetic structure based simply on arbitrary or convenient patterns of sampling.

Biogeography and speciation of a direct developing frog from the coastal arid zone of Western Australia

Within the southwestern Australian biodiversity hotspot, the Shark Bay region displays high levels of plant and animal endemism, particularly in the herpetofauna. The region has been subjected to dramatic climatic fluctuations and has been geologically active from the Late Miocene to the present. The myobatrachid frog Arenophryne rotunda, a Shark Bay endemic, provides an ideal opportunity to examine the relative effects of fluctuating climates and geological activity on the biota of Shark Bay. A comprehensive phylogeographic analysis of A. rotunda, based on data comprising 1154 bp of the mitochondrial gene ND2, is presented. My results demonstrate a major genetic break that divides this species at the northern edge of the Victorian Plateau into northern and southern species lineages, dating to the Late Miocene, with a further division of the southern species lineage across the Murchison Gorge dating to the Plio-Pleistocene border. Both of these periods are related to prominent geological activity and climatic shifts in the Shark Bay region. Interpretation of phylogeographic results point to the prominent role of fluctuating Pleistocene climates and associated coastal landscape evolution in the generation of phylogeographic structure within the distinct A. rotunda species lineages. Similar processes have been invoked to explain the diversity of other Shark Bay biota.

Demographic and genetic estimates of effective population and breeding size in the amphibian Rana temporaria

Genetic methods for estimating effective population size ( Ne) or the effective number of breeders ( Nb) have become popular, but comparisons of these estimates with demographic estimates of Ne and Nb are rare, especially in anurans. We used three genetic (linkage disequilibrium, temporal moments, Bayesian coalescent-based method) and three demographic models, the latter considering number of breeding individuals, sex ratio, reproductive skew, and other demographic data, to estimate Ne and Nb in two subarctic populations (T and P) of the common frog Rana temporaria, subject to long-term capture-recapture studies. Demographic estimates of Ne based on total population size ( Ne ([T])= 44.5-56.9; Ne ([P])= 68.8-93.7) deviated markedly from the genetic estimates obtained using the linkage disequilibrium method ( Ne ([T])= 97.1; Ne ([P])= 13.2). The demographic estimates of Nb, taking into consideration sex ratio and variance in reproductive success ( Nb ([T])= 10.1-39.7; Nb ([P])= 3.9-21.3), were higher than the genetic estimates ( Nb ([T])= 3.7-5.4; Nb ([P])= 3.5-3.9). The main factors affecting the effective size estimates were sex ratio and reproductive skew. The discrepancies between corresponding Ne and Nb estimates highlight the sensitivity of both demographic and genetic estimates on their underlying assumptions. Yet the ratios of effective or breeding effective size to the census population size were similar to those reported earlier for anurans, reinforcing the view that the discrepancy between actual and effective breeding sizes in anuran populations is typically very large.

Gene sampling strategies for multi-locus population estimates of genetic diversity (theta)

Background

Theoretical work suggests that data from multiple nuclear loci provide better estimates of population genetic parameters than do single loci, but just how many loci are needed and how much sequence is required from each has been little explored.

Methodology/Principle Findings

To investigate how much data is required to estimate the population genetic parameter θ (4N_eμ) accurately under ideal circumstances, we simulated datasets of DNA sequences under three values of θ per site (0.1, 0.01, 0.001), varying in both the total number of base pairs sequenced per individual and the number of equal-length loci. From these datasets we estimated θ using the maximum likelihood coalescent framework implemented in the computer program Migrate. Our results corroborated the theoretical expectation that increasing the number of loci impacted the accuracy of the estimate more than increasing the sequence length at single loci. However, when the value of θ was low (0.001), the per-locus sequence length was also important for estimating θ accurately, something that has not been emphasized in previous work.

Conclusions/Significance

Accurate estimation of θ required data from at least 25 independently evolving loci. Beyond this, there was little added benefit in terms of decreasing the squared coefficient of variation of the coalescent estimates relative to the extra effort required to sample more loci.

Ecological correlates of population genetic structure: a comparative approach using a vertebrate metacommunity

Identifying ecological factors associated with population genetic differentiation is important for understanding microevolutionary processes and guiding the management of threatened populations. We identified ecological correlates of several population genetic parameters for three interacting species (two garter snakes and an anuran) that occupy a common landscape. Using multiple regression analysis, we found that species interactions were more important in explaining variation in population genetic parameters than habitat and nearest-neighbour characteristics. Effective population size was best explained by census size, while migration was associated with differences in species abundance. In contrast, genetic distance was poorly explained by the ecological correlates that we tested, but geographical distance was prominent in models for all species. We found substantially different population dynamics for the prey species relative to the two predators, characterized by larger effective sizes, lower gene flow and a state of migration-drift equilibrium. We also identified an escarpment formed by a series of block faults that serves as a barrier to dispersal for the predators. Our results suggest that successful landscape-level management should incorporate genetic and ecological data for all relevant species, because even closely associated species can exhibit very different population genetic dynamics on the same landscape.

Refugial isolation and secondary contact in the dyeing poison frog Dendrobates tinctorius

Recent palaeoclimactic research suggests that fluctuating environmental conditions throughout the Pleistocene of Amazonia occurred with previously unrecognized frequency. This has resulted in a theoretical shift from glacially influenced fluctuations to those driven by precessional rhythms. This theoretical revolution has a profound impact on expectations of biotic diversity within biogeographical regions that have long been based on the idea of large-scale landscape fragmentation associated with increased aridity and glacial cycles. Generally speaking, this shifts phylogeographical expectations from that of (i) large areas of sympatry of closely related (but not sister) species whose origins lie in separate refugia, and current distributions are the results of cyclic connectivity of those two refugia (refuge hypothesis), to that of (ii) fine scale genetic structure, often associated with elevation, and divergence well below expected speciation levels [disturbance-vicariance (DV) hypothesis]. To date there have been few tests of the expectations of the DV hypothesis based on empirical studies of Neotropical floral and faunal communities. Herein we examine phylogeographical structure of Dendrobates tinctorius, an amphibian species endemic to the uplands of the eastern Guiana Shield, based on sampling of 114 individuals from 24 localities. Phylogenetic, nested clade, and dispersal-vicariance (DIVA) analyses of cytochrome b sequence data reveal the presence of two mitochondrial lineages that are associated with previously identified western and eastern uplands of this area. The geographical distribution of mitochondrial haplotypes and the results of DIVA and coalescent analyses suggest that there has been extensive secondary contact between these lineages indicating a complex history of connectivity between these western and eastern highlands, supporting the predictions of the DV hypothesis.

Fine-scale population differentiation and gene flow in a terrestrial salamander (Plethodon cinereus) living in continuous habitat

Several recent studies have shown that amphibian populations may exhibit high genetic subdivision in areas with recent fragmentation and urban development. Less is known about the potential for genetic differentiation in continuous habitats. We studied genetic differentiation of red-backed salamanders (Plethodon cinereus) across a 2-km transect through continuous forest in Virginia, USA. Mark-recapture studies suggest very little dispersal for this species, whereas homing experiments and post-Pleistocene range expansion both suggest greater dispersal abilities. We used six microsatellite loci to examine genetic population structure and differentiation between eight subpopulations of red-backed salamanders at distances from 200 m to 2 km. We also used several methods to extrapolate dispersal frequencies and test for sex-biased dispersal. We found small, but detectable differentiation among populations, even at distances as small as 200 m. Differentiation was closely correlated with distance and both Mantel tests and assignment tests were consistent with an isolation-by-distance model for the population. Extrapolations of intergenerational variance in spatial position (sigma(2)<15 m(2)) and pair-wise dispersal frequencies (4 Nm < 25 for plots separated by 300 m) both suggest limited gene flow. Additionally, tests for sex-biased dispersal imply that dispersal frequency is similarly low for both sexes. We suggest that these low levels of gene flow and the infrequent dispersal observed in mark-recapture studies may be reconciled with homing ability and range expansion if dispersing animals rarely succeed in breeding in saturated habitats, if dispersal is flexible depending on the availability of habitat, or if dispersal frequency varies across the geographic range of red-backed salamanders.

Hylid frog phylogeny and sampling strategies for speciose clades

How should characters and taxa be sampled to resolve efficiently the phylogeny of ancient and highly speciose groups? We addressed this question empirically in the treefrog family Hylidae, which contains > 800 species and may be nonmonophyletic with respect to other anuran families. We sampled 81 species (54 hylids and 27 outgroups) for two mitochondrial genes (12S, ND1), two nuclear genes (POMC, c-myc), and morphology (144 characters) in an attempt to resolve higher-level relationships. We then added 117 taxa to the combined data set, many of which were sampled for only one gene (12S). Despite the relative incompleteness of the majority of taxa, the resulting trees placed all taxa in the expected higher-level clades with strong support, despite some taxa being > 90% incomplete. Furthermore, we found no relationship between the completeness of a taxon and the support (parsimony bootstrap or Bayesian posterior probabilities) for its localized placement on the tree. Separate analysis of the data set with the most taxa (12S) gives a somewhat problematic estimate of higher-level relationships, suggesting that data sets scored only for some taxa (ND1, nuclear genes, morphology) are important in determining the outcome of the combined analysis. The results show that hemiphractine hylids are not closely related to other hylids and should be recognized as a distinct family. They also show that the speciose genus Hyla is polyphyletic, but that its species can be arranged into three monophyletic genera. A new classification of hylid frogs is proposed. Several potentially misleading signals in the morphological data are discussed.

Local origin and diversification among Lampornis hummingbirds: A Mesoamerican taxon

The huge biodiversity found in Mesoamerica is often explained by its geographic situation as a natural bridge between two large biogeographic regions. Often overlooked, however, are the high levels of speciation and diversification in the area. Here we assess the phylogenetic relationships within a Mesoamerican group of hummingbirds (Lampornis). We sequenced both mtDNA (1,143 bp of cyt b and 727 bp of ND5) and nuclear genes (505 bp of AK-5 intron and 567 bp of c-mos) for each of the seven recognised species and outgroups. We find two or three clades of similar age within this genus: L. clemenciae and L. amethystinus (singly or as each other's sister taxa) and a Central American clade. This Central-American clade presents a clear bipartition between northern (L. viridipallens and L. sybillae) and southern Mesoamerica, which is shared with many other Mesoamerican organisms. Our analyses suggest that L. hemileucus does not belong in the genus Lampornis. While we refrain to apply a time-scale to our data because of the lack of an appropriate calibration, our results indicate that the genus Lampornis predates the uprising of the Panama land-bridge, and that diversification among the isthmian species (L. castaneoventris and L. calolaema) is a very recent event. Our results strongly suggest a local Mesoamerican origin for this genus.

Asynchronous Colonization of Madagascar by the Four Endemic Clades of Primates, Tenrecs, Carnivores, and Rodents as Inferred from Nuclear Genes

Madagascar harbors four large adaptive radiations of endemic terrestrial mammals: lemurs, tenrecs, carnivorans, and rodents. These rank among the most spectacular examples of evolutionary diversification, but their monophyly and origins are debated. The lack of Tertiary fossils from Madagascar leaves molecular studies as most promising to solve these controversies. We provide a simultaneous reconstruction of phylogeny and age of the four radiations based on a 3.5-kb data set from three nuclear genes (ADRA2B, vWF, and AR). The analysis supports each as a monophyletic clade, sister to African taxa, and thereby identifies four events of colonization out of Africa. To infer the time windows for colonization, we take into account both the divergence from the closest non-insular sister group and the initial intra-insular radiation, which is a novel but conservative approach in studies of the colonization history of Madagascar. We estimate that lemurs colonized Madagascar between 60 million years ago (Mya) (split from lorises) and 50 Mya (lemur radiation) (70-41 Mya taking 95% credibility intervals into account), tenrecs between 42 and 25 Mya (50-20 Mya), carnivorans between 26 and 19 Mya (33-14 Mya), and rodents between 24 and 20 Mya (30-15 Mya). These datings suggest at least two asynchronous colonization events: by lemurs in the Late Cretaceous-Middle Eocene, and by carnivorans and rodents in the Early Oligocene-Early Miocene. The colonization by tenrecs may have taken place simultaneously with either of these two events, or in a third event in the Late Eocene-Oligocene. Colonization by at least lemurs, rodents, and carnivorans appears to have occurred by overseas rafting rather than via a land bridge hypothesized to have existed between 45 and 26 Mya, but the second scenario cannot be ruled out if credibility intervals are taken into account.

Reconciling nuclear microsatellite and mitochondrial marker estimates of population structure: breeding population structure of Chesapeake Bay striped bass (Morone saxatilis)

Comparative analyses of nuclear and organelle genetic markers may help delineate evolutionarily significant units or management units, although population differentiation estimates from multiple genomes can also conflict. Striped bass (Morone saxatilis) are long-lived, highly migratory anadromous fish recently recovered from a severe decline in population size. Previous studies with protein, nuclear DNA and mitochondrial DNA (mtDNA) markers produced discordant results, and it remains uncertain if the multiple tributaries within Chesapeake Bay constitute distinct management units. Here, 196 young-of-the-year (YOY) striped bass were sampled from Maryland's Choptank, Potomac and Nanticoke Rivers and the north end of Chesapeake Bay in 1999 and from Virginia's Mataponi and Rappahannock Rivers in 2001. A total of 10 microsatellite loci exhibited between two and 27 alleles per locus with observed heterozygosities between 0.255 and 0.893. The 10-locus estimate of R(ST) among the six tributaries was -0.0065 (95% confidence interval -0.0198 to 0.0018). All R(ST) and all but one theta estimates for pairs of populations were not significantly different from zero. Reanalysis of Chesapeake Bay striped bass mtDNA data from two previous studies estimated population differentiation between theta=-0.002 and 0.160, values generally similar to mtDNA population differentiation predicted from microsatellite R(ST) after adjusting for reduced effective population size and uniparental inheritance in organelle genomes. Based on mtDNA differentiation, breeding sex ratios or gene flow may have been slightly male biased in some years. The results reconcile conflicting past studies based on different types of genetic markers, supporting a single Chesapeake Bay management unit encompassing a panmictic striped bass breeding population.

Cenozoic biogeography and evolution in direct-developing frogs of Central America (Leptodactylidae: Eleutherodactylus) as inferred from a phylogenetic analysis of nuclear and mitochondrial genes

We report the first phylogenetic analysis of DNA sequence data for the Central American component of the genus Eleutherodactylus (Anura: Leptodactylidae: Eleutherodactylinae), one of the most ubiquitous, diverse, and abundant components of the Neotropical amphibian fauna. We obtained DNA sequence data from 55 specimens representing 45 species. Sampling was focused on Central America, but also included Bolivia, Brazil, Jamaica, and the USA. We sequenced 1460 contiguous base pairs (bp) of the mitochondrial genome containing ND2 and five neighboring tRNA genes, plus 1300 bp of the c-myc nuclear gene. The resulting phylogenetic inferences were broadly concordant between data sets and among analytical methods. The subgenus Craugastor is monophyletic and its initial radiation was potentially rapid and adaptive. Within Craugastor, the earliest splits separate three northern Central American species groups, milesi, augusti, and alfredi, from a clade comprising the rest of Craugastor. Within the latter clade, the rhodopis group as formerly recognized comprises three deeply divergent clades that do not form a monophyletic group; we therefore restrict the content of the rhodopis group to one of two northern clades, and use new names for the other northern (mexicanus group) and one southern clade (bransfordii group). The new rhodopis and bransfordii groups together form the sister taxon to a clade comprising the biporcatus, fitzingeri, mexicanus, and rugulosus groups. We used a Bayesian MCMC approach together with geological and biogeographic assumptions to estimate divergence times from the combined DNA sequence data. Our results corroborated three independent dispersal events for the origins of Central American Eleutherodactylus: (1) an ancestor of Craugastor entered northern Central America from South American in the early Paleocene, (2) an ancestor of the subgenus Syrrhophus entered northern Central America from the Caribbean at the end of the Eocene, and (3) a wave of independent dispersal events from South America coincided with formation of the Isthmus of Panama during the Pliocene. We elevate the subgenus Craugastor to the genus rank.

Long-distance gene flow and cross-Andean dispersal of lowland rainforest bees (Apidae: Euglossini) revealed by comparative mitochondrial DNA phylogeography

Euglossine bees (Apidae; Euglossini) exclusively pollinate hundreds of orchid species and comprise up to 25% of bee species richness in neotropical rainforests. As one of the first studies of comparative phylogeography in a neotropical insect group, we performed a mitochondrial DNA (mtDNA)-based analysis of 14 euglossine species represented by populations sampled across the Andes and/or across the Amazon basin. The mtDNA divergences within species were consistently low; across the 12 monophyletic species the mean intraspecific divergence among haplotypes was 0.9% (range of means, 0-1.9%). The cytochrome oxidase 1 (CO1) divergence among populations separated by the Andes (N = 11 species) averaged 1.1% (range 0.0-2.0%). The mtDNA CO1 data set displayed homogeneous rates of nucleotide substitution, permitting us to infer dispersal across the cordillera long after the final Andean uplift based on arthropod molecular clocks of 1.2-1.5% divergence per million years. Gene flow across the 3000-km breadth of the Amazon basin was inferred from identical cross-Amazon haplotypes found in five species. Although mtDNA haplotypes for 12 of the 14 euglossine species were monophyletic, a reticulate CO1 phylogeny was recovered in Euglossa cognata and E. mixta, suggesting large ancestral populations and recent speciation. Reference to closely related outgroups suggested recent speciation for the majority of species. Phylogeographical structure across a broad spatial scale is weaker in euglossine bees than in any neotropical group previously examined, and may derive from a combination of Quaternary speciation, population expansion and/or long-distance gene flow.

Genetic variation coincides with geographic structure in the common bush-tanager (Chlorospingus ophthalmicus) complex from Mexico

Cloud forests are distributed in the Neotropics, from northern Mexico to Argentina, under very specific ecological conditions, namely slopes with high humidity input from clouds and mist. Its distribution in Mesoamerica is highly fragmented, similar to an archipelago, and taxa are thus frequently represented as sets of isolated populations, each restricted to particular mountain ranges and often showing a high degree of divergence, both morphologically and genetically. The common bush-tanager (Chlorospingus ophthalmicus, Aves: Thraupidae) inhabits cloud forests from eastern and southern Mexico south to northwestern Argentina. Here we use 676bp of mtDNA (around the ATPase 8 gene) to explore the genetic variation and phylogeographic structure of the Mexican populations of C. ophthalmicus. Phylogenetic analyses of mtDNA sequences indicate deep genetic structure. Five major clades, which segregate according to geographic breaks, are identified (starting from the deepest one in the phylogeny): (1) Southern Chiapas and Northern Central America, (2) Tuxtlas massif, (3) Sierra Madre del Sur, (4) Eastern Oaxaca and Northern Chiapas, and (5) Sierra Madre Oriental. The long history of isolation undergone by each clade, as suggested by the phylogeny, implies that the species status of each of them should be revised.