Evolutionary history of the northern leopard frog: reconstruction of phylogeny, phylogeography, and historical changes in population demography from mitochondrial DNA

This study uses a combined methodological approach including phylogenetic, phylogeographic, and demographic analyses to understand the evolutionary history of the northern leopard frog, Rana pipiens. We tested hypotheses concerning how (or if) known geological events and key features of the species biology influenced the contemporary geographic and genetic distribution of R. pipiens. We assayed mitochondrial DNA variation from 389 individuals within 35 populations located throughout the species range. Our a priori expectations for patterns and processes influencing the current genetic structure of R. pipiens were supported by the data. However, our analyses revealed specific aspects of R. pipiens evolutionary history that were unexpected. The phylogenetic analysis indicated that R. pipiens is split into populations containing discrete eastern or western haplotypes, with the Mississippi River and Great Lakes region dividing the geographic ranges. Nested clade analysis indicated that the biological process most often invoked to explain the pattern of haplotype position is restricted gene flow with isolation by distance. Demographic analyses showed evidence of both historical bottlenecks and population expansions. Surprisingly, the genetic evidence indicated that the western haplotypes had significantly reduced levels of genetic diversity relative to the eastern haplotypes and that major range expansions occurred in both regions well before the most recent glacial retreat. This study provides a detailed history of how a widespread terrestrial vertebrate responded to episodic Pleistocene glacial events in North America. Moreover, this study illustrates how complementary methods of data analysis can be used to disentangle recent and ancient effects on the genetic structure of a species.     

Living on the edge: the role of geography and environment in structuring genetic variation in the southernmost populations of a tropical oak

Understanding the factors determining genetic diversity and structure in peripheral populations is a long‐standing goal of evolutionary biogeography, yet little empirical information is available for tropical species. In this study, we combine information from nuclear microsatellite markers and niche modelling to analyse the factors structuring genetic variation across the southernmost populations of the tropical oak Quercus segoviensis. First, we tested the hypothesis that genetic variability decreases with population isolation and increases with local habitat suitability and stability since the Last Glacial Maximum (LGM). Second, we employed a recently developed multiple matrix regression with randomisation (MMRR) approach to study the factors associated with genetic divergence among the studied populations and test the relative contribution of environmental and geographic isolation to contemporary patterns of genetic differentiation. We found that genetic diversity was negatively correlated with average genetic differentiation with other populations, indicating that isolation and limited gene flow have contributed to erode genetic variability in some populations. Considering the relatively small size of the study area (<120 km), analyses of genetic structure indicate a remarkable inter‐population genetic differentiation. Environmental dissimilarity and differences in current and past climate niche suitability and their additive effects were not associated with genetic differentiation after controlling for geographic distance, indicating that local climate does not contribute to explain spatial patterns of genetic structure. Overall, our data indicate that geographic isolation, but not current or past climate, is the main factor determining contemporary patterns of genetic diversity and structure within the southernmost peripheral populations of this tropical oak.     

Vicariance and marine migration in continental island populations of a frog endemic to the Atlantic Coastal forest

The theory of island biogeography is most often studied in the context of oceanic islands where all island inhabitants are descendants from founding events involving migration from mainland source populations. Far fewer studies have considered predictions of island biogeography in the case of continental islands, where island formation typically splits continuous populations and thus vicariance also contributes to the diversity of island populations. We examined one such case on continental islands in southeastern Brazil, to determine how classic island biogeography predictions and past vicariance explain the population genetic diversity of Thoropa taophora, a frog endemic to the Atlantic Coastal Forest. We used nuclear microsatellite markers to examine the genetic diversity of coastal and island populations of this species. We found that island isolation has a role in shaping the genetic diversity of continental island species, with island populations being significantly less diverse than coastal populations. However, area of the island and distance from coast had no significant effect on genetic diversity. We also found no significant differences between migration among coastal populations and migration to and from islands. We discuss how vicariance and the effects of continued migration between coastal and island populations interact to shape evolutionary patterns on continental islands.     

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub‐legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west‐to‐east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.     

Phylogeographic structure of the fire ant Solenopsis invicta in its native South American range: roles of natural barriers and habitat connectivity

We generated mitochondrial DNA (mtDNA) sequence data from 402 individuals of the fire ant Solenopsis invicta collected from 11 native populations and analyzed these data using a combination of demographic, phylogenetic, and phylogeographic methods to infer features of the evolutionary history of this species. Prior expectations regarding high levels of genetic structure and isolation by distance among populations were supported by the data, but we also discovered several unanticipated patterns. Our analyses revealed a major genetic break between S. invicta mtDNA haplotypes that coincides with the Mesopotamia wetlands region of South America, resulting in two higher level nested clade groupings. In addition, we identified contrasting patterns of genetic differentiation within these two major groups, which may reflect differences in connectivity of suitable habitat in different parts of the native range of S. invicta. Our study represents the first attempt to understand the phylogeographic history of S. invicta across its native range.     

The genetic structure of Drosophila ananassae populations from Asia, Australia and Samoa

Information about genetic structure and historical demography of natural populations is central to understanding how natural selection changes genomes. Drosophila ananassae is a widespread species occurring in geographically isolated or partially isolated populations and provides a unique opportunity to investigate population structure and molecular variation. We assayed microsatellite repeat-length variation among 13 populations of D. ananassae to assess the level of structure among the populations and to make inferences about their ancestry and historic biogeography. High levels of genetic structure are apparent among all populations, particularly in Australasia and the South Pacific, and patterns are consistent with the hypothesis that the ancestral populations are from Southeast Asia. Analysis of population structure and use of F-statistics and Bayesian analysis suggest that the range expansion of the species into the Pacific is complex, with multiple colonization events evident in some populations represented by lineages that show no evidence of recent admixture. The demographic patterns show isolation by distance among populations and population expansion within all populations. A morphologically distinct sister species, D. pallidosa, collected in Malololelei, Samoa, appears to be more closely related to some of the D. ananassae populations than many of the D. ananassae populations are to one another. The patterns of genotypic diversity suggest that many of the individuals that we sampled may be morphologically indistinguishable nascent species.     

In situ genetic differentiation in a Hispaniolan lizard (Ameiva chrysolaema): a multilocus perspective

A previous phylogeographic study of mitochondrial haplotypes for the Hispaniolan lizard Ameiva chrysolaema revealed deep genetic structure associated with seawater inundation during the late Pliocene/early Pleistocene and evidence of subsequent population expansion into formerly inundated areas. We revisit hypotheses generated by our previous study using increased geographic sampling of populations and analysis of three nuclear markers (alpha-enolase intron 8, alpha-cardiac-actin intron 4, and beta-actin intron 3) in addition to mitochondrial haplotypes (ND2). Large genetic discontinuities correspond spatially and temporally with historical barriers to gene flow (sea inundations). NCPA cross-validation analysis and Bayesian multilocus analyses of divergence times (IMa and MCMCcoal) reveal two separate episodes of fragmentation associated with Pliocene and Pleistocene sea inundations, separating the species into historically separate Northern, East-Central, West-Central, and Southern population lineages. Multilocus Bayesian analysis using IMa indicates asymmetrical migration from the East-Central to the West-Central populations following secondary contact, consistent with expectations from the more pervasive sea inundation in the western region. The West-Central lineage has a genetic signature of population growth consistent with the expectation of geographic expansion into formerly inundated areas. Within each lineage, significant spatial genetic structure indicates isolation by distance at comparable temporal scales. This study adds to the growing body of evidence that vicariant speciation may be the prevailing source of lineage accumulation on oceanic islands. Thus, prior theories of island biogeography generally underestimate the role and temporal scale of intra-island vicariant processes.     

Phylogeography of a facultatively migratory dragonfly, Libellula quadrimaculata (Odonata: Anisoptera)

The biogeography of a widely distributed dragonfly (Libellula quadrimaculata) was examined via a phylogenetic analysis of 416 bp of the mitochondrial cytochrome oxidase I subunit (COI). Phylogenetic analyses under parsimony and minimum evolution produced trees with similar topologies, and revealed strong support for three clades corresponding to populations in Asia, Europe and North America. However, resolution was poor within clades, and genetic distances between populations within continents was quite low (1–2%). Several populations of this species are known to engage in periodic mass migrations, and it is possible that introgression from gene flow due to the mobility of this species has obscured phylogenetic patterns within continents. I was unable to test for phylogenetic patterns coincident with historical glacial refugia given the lack of phylogenetic patterns within continents. However, given that some sequence divergence was observed between populations within continents, it is possible that phylogenetic patterns exist, and subsequent studies should make use of larger data sets, and molecular data from faster evolving genes. Despite the propensity for periodic, short distance migrations in L. quadrimaculata, gene flow appears to be limited and does not influence the phylogenetic relationships of populations between continents.     

Species richness and evolutionary niche dynamics: a spatial pattern-oriented simulation experiment

Evolutionary processes underlying spatial patterns in species richness remain largely unexplored, and correlative studies lack the theoretical basis to explain these patterns in evolutionary terms. In this study, we develop a spatially explicit simulation model to evaluate, under a pattern-oriented modeling approach, whether evolutionary niche dynamics (the balance between niche conservatism and niche evolution processes) can provide a parsimonious explanation for patterns in species richness. We model the size, shape, and location of species' geographical ranges in a multivariate heterogeneous environmental landscape by simulating an evolutionary process in which environmental fluctuations create geographic range fragmentation, which, in turn, regulates speciation and extinction. We applied the model to the South American domain, adjusting parameters to maximize the correspondence between observed and predicted patterns in richness of about 3,000 bird species. Predicted spatial patterns, which closely resemble observed ones (r2=0.795), proved sensitive to niche dynamics processes. Our simulations allow evaluation of the roles of both evolutionary and ecological processes in explaining spatial patterns in species richness, revealing the enormous potential of the link between ecology and historical biogeography under integrated theoretical and methodological frameworks.     

Linking genetic and ecological differentiation in an ungulate with a circumpolar distribution

Genetic differentiation among populations may arise from the disruption of gene flow due to local adaptation to distinct environments and/or neutral accumulation of mutations and genetic drift resulted from geographical isolation. Quantifying the role of these processes in determining the genetic structure of natural populations remains challenging. Here, we analyze the relative contribution of isolation‐by‐resistance (IBR), isolation‐by‐environment (IBE), genetic drift and historical isolation in allopatry during Pleistocene glacial cycles on shaping patterns of genetic differentiation in caribou/reindeer populations Rangifer tarandus across the entire distribution range of the species. Our study integrates analyses at range‐wide and regional scales to partial out the effects of historical and contemporary isolation mechanisms. At the circumpolar scale, our results indicate that genetic differentiation is predominantly explained by IBR and historical isolation. At a regional scale, we found that IBR, IBE and population size significantly explained the spatial distribution of genetic variation among populations belonging to the Euro‐Beringian lineage within North America. In contrast, genetic differentiation among populations within the North American lineage was predominantly explained by IBR and population size, but not IBE. We also found discrepancies between genetic and ecotype designation across the Holarctic species distribution range. Overall, these results indicate that multiple isolating mechanisms have played roles in shaping the spatial distribution of genetic variation across the distribution range of a large mammal with high potential for gene flow. Considering multiple spatial scales and simultaneously testing a comprehensive suite of potential isolating mechanisms, our study contributes to understand the ecological and evolutionary processes underlying organism–landscape interactions.     

Genetic differentiation and gene flow among populations of the alpine butterfly, Parnassius smintheus, vary with landscape connectivity

Levels of gene flow among populations vary both inter- and intraspecifically, and understanding the ecological bases of variation in levels of gene flow represents an important link between the ecological and evolutionary dynamics of populations. The effects of habitat spatial structure on gene flow have received considerable attention; however, most studies have been conducted at a single spatial scale and without background data on how individual movement is affected by landscape features. We examined the influence of habitat connectivity on inferred levels of gene flow in a high-altitude, meadow-dwelling butterfly, Parnassius smintheus. For this species, we had background data on the effects of landscape structure on both individual movement and on small-scale population genetic differentiation. We compared genetic differentiation and patterns of isolation by distance, based on variation at seven microsatellite loci, among three regions representing two levels of connectivity of high-altitude, nonforested habitats. We found that reduced connectivity of habitats, resulting from more forest cover at high altitudes, was associated with greater genetic differentiation among populations (higher estimated FST), a breakdown of isolation by distance, and overall lower levels of inferred gene flow. These observed differences were consistent with expectations based on our knowledge of the movement behaviour of this species and on previous population genetic analyses conducted at the smaller spatial scale. Our results indicate that the role of gene flow may vary among groups of populations depending on the interplay between individual movement and the structure of the surrounding landscape.     

Revisiting Amazonian phylogeography: insights into diversification hypotheses and novel perspectives

The Amazon Basin harbors one of the richest biotas on Earth, such that a number of diversification hypotheses have been formulated to explain patterns of Amazonian biodiversity and biogeography. For nearly two decades, phylogeographic approaches have been applied to better understand the underlying causes of genetic differentiation and geographic structure among Amazonian organisms. Although this research program has made progress in elucidating several aspects of species diversification in the region, recent methodological and theoretical developments in the discipline of phylogeography will provide new perspectives through more robust hypothesis testing. Herein, we outline central aspects of Amazonian geology and landscape evolution as well as climate and vegetation dynamics through the Neogene and Quaternary to contextualize the historical settings considered by major hypotheses of diversification. We address each of these hypotheses by reviewing key phylogeographic papers and by expanding their respective predictions. We also propose future directions for devising and testing hypotheses. Specifically, combining the exploratory power of phylogeography with the statistical rigor of coalescent methods will greatly expand analytical inferences on the evolutionary history of Amazonian biota. Incorporation of non-genetic data from Earth science disciplines into the phylogeographic approach is key to a better understanding of the influence of climatic and geophysical events on patterns of Amazonian biodiversity and biogeography. In addition, achieving such an integrative enterprise must involve overcoming issues such as limited geographic and taxonomic sampling. These future challenges likely will be accomplished by a combination of extensive collaborative research and incentives for conducting basic inventories.     

Exploring demographic, physical, and historical explanations for the genetic structure of two lineages of Greater Antillean bats

Observed patterns of genetic structure result from the interactions of demographic, physical, and historical influences on gene flow. The particular strength of various factors in governing gene flow, however, may differ between species in biologically relevant ways. We investigated the role of demographic factors (population size and sex-biased dispersal) and physical features (geographic distance, island size and climatological winds) on patterns of genetic structure and gene flow for two lineages of Greater Antillean bats. We used microsatellite genetic data to estimate demographic characteristics, infer population genetic structure, and estimate gene flow among island populations of Erophylla sezekorni/E. bombifrons and Macrotus waterhousii (Chiroptera: Phyllostomidae). Using a landscape genetics approach, we asked if geographic distance, island size, or climatological winds mediate historical gene flow in this system. Samples from 13 islands spanning Erophylla's range clustered into five genetically distinct populations. Samples of M. waterhousii from eight islands represented eight genetically distinct populations. While we found evidence that a majority of historical gene flow between genetic populations was asymmetric for both lineages, we were not able to entirely rule out incomplete lineage sorting in generating this pattern. We found no evidence of contemporary gene flow except between two genetic populations of Erophylla. Both lineages exhibited significant isolation by geographic distance. Patterns of genetic structure and gene flow, however, were not explained by differences in relative effective population sizes, island area, sex-biased dispersal (tested only for Erophylla), or surface-level climatological winds. Gene flow among islands appears to be highly restricted, particularly for M. waterhousii, and we suggest that this species deserves increased taxonomic attention and conservation concern.     

Population Structure and Dispersal Patterns within and between Atlantic and Mediterranean Populations of a Large-Range Pelagic Seabird

Dispersal is critically linked to the demographic and evolutionary trajectories of populations, but in most seabird species it may be difficult to estimate. Using molecular tools, we explored population structure and the spatial dispersal pattern of a highly pelagic but philopatric seabird, the Cory''s shearwater Calonectris diomedea. Microsatellite fragments were analysed from samples collected across almost the entire breeding range of the species. To help disentangle the taxonomic status of the two subspecies described, the Atlantic form C. d. borealis and the Mediterranean form C. d. diomedea, we analysed genetic divergence between subspecies and quantified both historical and recent migration rates between the Mediterranean and Atlantic basins. We also searched for evidence of isolation by distance (IBD) and addressed spatial patterns of gene flow. We found a low genetic structure in the Mediterranean basin. Conversely, strong genetic differentiation appeared in the Atlantic basin. Even if the species was mostly philopatric (97%), results suggest recent dispersal between basins, especially from the Atlantic to the Mediterranean (aprox. 10% of migrants/generation across the last two generations). Long-term gene flow analyses also suggested an historical exchange between basins (about 70 breeders/generation). Spatial analysis of genetic variation indicates that distance is not the main factor in shaping genetic structure in this species. Given our results we recommend gathering more data before concluded whether these taxa should be treated as two species or subspecies.     

Bridging the gap between community ecology and historical biogeography: niche conservatism and community structure in emydid turtles

Historical (phylogenetic) biogeography and community ecology were once integrated as part of the broader study of organismal diversity, but in recent decades have become largely separate disciplines. This is unfortunate because many patterns studied by community ecologists may originate through processes studied by historical biogeographers and vice versa. In this study, we explore the causes of a geographic pattern of community structure (habitat use) in the emydid turtle assemblages of eastern North America, with more semi-terrestrial species of the subfamily Emydinae in the north and more aquatic species of Deirochelyinae in the south. Specifically, we address the factors that prevent northern emydines from invading southern communities. We test for competitive exclusion by examining patterns of range overlap, and test for the role of niche conservatism using analyses of climatic and physiological data based on a multilocus molecular phylogeny. We find no support for competitive exclusion, whereas several lines of evidence support the idea that niche conservatism has prevented northern emydines from dispersing into southern communities. Our results show how understanding the causes of patterns of historical biogeography may help explain patterns of community structure.     

Historical demography, selection, and coalescence of mitochondrial and nuclear (genes in Prochilodus species of northern South America

Fishes of the genus Prochilodus are ecologically and commercially important, ubiquitous constituents of large river biota in South America. Recent ecologic and demographic studies indicate that these fishes exist in large, stable populations with adult census numbers exceeding one million individuals. Abundance data present a stark contrast to very low levels of genetic diversity (theta) and small effective population sizes (Ne) observed in a mitochondrial (mt) DNA dataset obtained for two species, Prochilodus mariae, and its putative sister taxon, Prochilodus rubrotaeniatus. Both species occupy major river drainages (Orinoco, Essequibo, and Negro) of northeastern South America. Disparity between expectations based on current abundance and life history information and observed genetic data in these lineages could result from historical demographic bottlenecks, or alternatively, natural selection (i.e., a mtDNA selective sweep). To ascertain underlying processes that affect mtDNA diversity in these species we compared theta and Ne estimates obtained from two, unlinked nuclear loci (calmodulin intron-4 and elongation factor-1alpha intron-6) using an approach based on coalescent theory. Genetic diversity and Ne estimated from mtDNA and nuclear sequences were uniformly low in P. rubrotaeniatus from the Rio Negro, suggesting that this population has encountered a historical bottleneck. For all P. mariae populations, theta and Ne based on nuclear sequences were comparable to expectations based on current adult census numbers and were significantly greater than mtDNA estimates, suggesting that a selective mtDNA sweep has occurred in this species. Comparative genetic analysis indicates that a suite of evolutionary processes involving historical demography and natural selection have influenced patterns of genetic variation and speciation in this important Neotropical fish group.     

Retention of gene diversity during the spread of a non‐native plant species

Spatial expansion, which is a crucial stage in the process to successful biological invasion, is anticipated to profoundly affect the magnitude and spatial distribution of genetic diversity in novel colonized areas. Here, we show that, contrasting common expectations, Pyrenean rocket (Sisymbrium austriacum), retained SNP diversity as this introduced plant species descended in the Meuse River Basin. Allele frequencies did not mirror between‐population distances along the predominant expansion axis. Reconstruction of invasion history based on the genotypes of historical herbarium specimens indicated no influence of additional introductions or multiple points of entry on this nongradual pattern. Assignment analysis suggested the admixture of distant upstream sources in recently founded downstream populations. River dynamics seem to have facilitated occasional long‐distance dispersal which brought diversity to the expansion front and so maintained evolutionary potential. Our findings highlight the merit of a historical framework in interpreting extant patterns of genetic diversity in introduced species and underscore the need to integrate long‐distance dispersal events in theoretical work on the genetic consequences of range expansion.     

Landscape genetics of a specialized grasshopper inhabiting highly fragmented habitats: a role for spatial scale

Aim The study of geographical discontinuities in the distribution of genetic variability in natural populations is a central topic in both evolutionary and conservation research. In this study, we aimed to analyse (1) the factors associated with genetic diversity at the landscape spatial scale in the highly specialized grasshopper Mioscirtus wagneri and (2) to identify the relative contribution of alternative factors to the observed patterns of genetic structure in this species. Location La Mancha region, Central Spain. Methods We sampled 28 populations of the grasshopper M. wagneri and genotyped 648 individuals at seven microsatellite loci. We employed a causal modelling approach to identify the most influential variables associated with genetic differentiation within a multiple hypothesis‐testing framework. Results We found that genetic diversity differs among populations located in different river basins and decreases with population isolation. Causal modelling analyses showed variability in the relative influence of the studied landscape features across different spatial scales. When a highly isolated population is considered, the analyses suggested that geographical distance is the only factor explaining the genetic differentiation between populations. When that population is excluded, the causal modelling analysis revealed that elevation and river basins are also relevant factors contributing to explaining genetic differentiation between the studied populations. Main conclusions These results indicate that the spatial scale considered and the inclusion of outlier populations may have important consequences on the inferred contribution of alternative landscape factors on the patterns of genetic differentiation even when all populations are expected to similarly respond to landscape structure. Thus, a multiscale perspective should also be incorporated into the landscape genetics framework to avoid biased conclusions derived from the spatial scale analysed and/or the geographical distribution of the studied populations.     

Widespread geographical distribution of mitochondrial haplotypes in rock-dwelling cichlid fishes from Lake Tanganyika

The spectacularly diverse cichlid fish species flocks of the East African Rift Lakes have elicited much debate on the potential evolutionary mechanisms responsible for the origin of these adaptive radiations. An historical perspective on population structure may offer insights into the processes driving population differentiation and possibly speciation. Here, we examine mitochondrial DNA (mtDNA) sequence variation in two endemic species of rock-dwelling cichlids, Simochromis babaulti and S. diagramma , from Lake Tanganyika. Phylogeographic analyses were used to infer what factors might have been important in the genetic structuring of Simochromis populations. Patterns of mtDNA differentiation in Simochromis were compared to those of other rock-dwelling cichlids to distinguish between competing hypotheses concerning the processes underlying their evolution. In striking contrast to previous findings, populations of Simochromis , even those separated by up to 300 km, were found to share mitochondrial DNA haplotypes. There is no correspondence between mtDNA genealogies and the geographical distribution of populations. Only S. babaulti , but not S. diagramma was found to have a significant association between genetic and geographical distance. These phylogeographic patterns suggest that the evolutionary effects of abiotic and biotic factors shaping population genetic structure may differ substantially even among closely related species of rock-dwelling cichlids. Physical events and barriers to gene flow that are believed to have had a major impact on the geographical distribution and intralacustrine speciation of Tropheus do not seem to have equally strongly affected its close relative Simochromis . These findings emphasize that no single mechanism can be responsible for the formation of population structure, speciation, and the adaptive radiation of all cichlid fishes.