Single copy nuclear DNA markers characterized for comparative phylogeography in Australian wet tropics rainforest skinks

Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia

Background

The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia.

Methods

Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data.

Results

Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus’ populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107–0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

Cryptic diversity and ancient diversification in the northern Atlantic Forest Pristimantis (Amphibia, Anura, Craugastoridae)

The Neotropical region is known both for its megadiverse fauna and for the significant deficiency of our knowledge on species limits. The Amazon and Atlantic Forest are the two most diverse and large rainforests in South America, and they harbor many groups of sister species and sister genera. The frog genus Pristimantis is the most speciose genus of terrestrial vertebrates with 546 species, but only three of them occur in the Atlantic forest. Herein, we investigate the diversification history and phylogenetic relationship among the Atlantic Forest Pristimantis lineages in a spatial-temporal framework, using mitochondrial and nuclear genes. Our results reveal a remarkable hidden diversity, with nine highly structure lineages that may correspond to undescribed species, with many cases of sympatry among these divergent lineages. Atlantic forest Pristimantis form a monophyletic group that started to diversify over 40 million years ago. This ancient group shows diversification events that remount the early bursts of the Atlantic Forest diversification history, as well as lineage diversification likely resulting from recent Pleistocene climatic fluctuations. Future work must concentrate in comparing these lineages under an integrative framework including morphology, advertisement calls and other ecology traits to confidently delimit species of Pristimantis in the Atlantic Forest.

Male-male competition and repeated evolution of terrestrial breeding in Atlantic Coastal Forest frogs

Terrestrial breeding is a derived condition in frogs, with multiple transitions from an aquatic ancestor. Shifts in reproductive mode often involve changes in habitat use, and these are typically associated with diversification in body plans, with repeated transitions imposing similar selective pressures. We examine the diversification of reproductive modes, male and female body sizes, and sexual size dimorphism (SSD) in the Neotropical frog genera Cycloramphus and Zachaenus, both endemic to the Atlantic rainforest of Brazil. Species in this clade either breed in rocky streams (saxicolous) or in terrestrial environments, allowing us to investigate reproductive habitat shifts. We constructed a multilocus molecular phylogeny and inferred evolutionary histories of reproductive habitats, body sizes, and SSD. The common ancestor was small, saxicolous, and had low SSD. Terrestrial breeding evolved independently three times and we found a significant association between reproductive habitat and SSD, with shifts to terrestrial breeding evolving in correlation with decreases in male body size, but not female body size. Terrestrial breeding increases the availability of breeding sites and results in concealment of amplexus, egg-laying, and parental care, therefore reducing male-male competition at all stages of reproduction. We conclude that correlated evolution of terrestrial reproduction and small males is due to release from intense male-male competition that is typical of exposed saxicolous breeding.

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae)

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Genetic structure and diversity of the endangered growling grass frog in a rapidly urbanizing region

Two pervasive and fundamental impacts of urbanization are the loss and fragmentation of natural habitats. From a genetic perspective, these impacts manifest as reduced genetic diversity and ultimately reduced genetic viability. The growling grass frog (Litoria raniformis) is listed as vulnerable to extinction in Australia, and endangered in the state of Victoria. Remaining populations of this species in and around the city of Melbourne are threatened by habitat loss, degradation and fragmentation due to urban expansion. We used mitochondrial DNA (mtDNA) and microsatellites to study the genetic structure and diversity of L. raniformis across Melbourne's urban fringe, and also screened four nuclear gene regions (POMC, RAG-1, Rhod and CRYBA1). The mtDNA and nuclear DNA sequences revealed low levels of genetic diversity throughout remnant populations of L. raniformis. However, one of the four regions studied, Cardinia, exhibited relatively high genetic diversity and several unique haplotypes, suggesting this region should be recognized as a separate Management Unit. We discuss the implications of these results for the conservation of L. raniformis in urbanizing landscapes, particularly the potential risks and benefits of translocation, which remains a contentious management approach for this species.

Persistence and dispersal in a Southern Hemisphere glaciated landscape: the phylogeography of the spotted snow skink (Niveoscincus ocellatus) in Tasmania

Background

The aim of this research was to identify the effects of Pleistocene climate change on the distribution of fauna in Tasmania, and contrast this with biotic responses in other temperate regions in the Northern and Southern Hemisphere that experienced glacial activity during this epoch. This was achieved by examining the phylogeographic patterns in a widely distributed Tasmanian endemic reptile, Niveoscincus ocellatus. 204 individuals from 29 populations across the distributional range of N. ocellatus were surveyed for variation at two mitochondrial genes (ND2, ND4), and two nuclear genes (β-globin, RPS8). Phylogenetic relationships were reconstructed using a range of methods (maximum parsimony, Bayesian inference and haplotype networks), and the demographic histories of populations were assessed (AMOVA, Tajima’s D, Fu’s Fs, mismatch distributions, extended Bayesian skyline plots, and relaxed random walk analyses).

Results

There was a high degree of mitochondrial haplotype diversity (96 unique haplotypes) and phylogeographic structure, where spatially distinct groups were associated with Tasmania’s Northeast and a large area covering Southeast and Central Tasmania. Phylogeographic structure was also present within each major group, but the degree varied regionally, being highest in the Northeast. Only the Southeastern group had a signature of demographic expansion, occurring during the Pleistocene but post-dating the Last Glacial Maximum. In contrast, nuclear DNA had low levels of variation and a lack of phylogeographic structure, and further loci should be surveyed to corroborate the mitochondrial inferences.

Conclusions

The phylogeographic patterns of N. ocellatus indicate Pleistocene range and demographic expansion in N. ocellatus, particularly in the Southeast and Central areas of Tasmania. Expansion in Central and Southeastern areas appears to have been more recent in both demographic and spatial contexts, than in Northeast Tasmania, which is consistent with inferences for other taxa of greater stability and persistence in Northeast Tasmania during the Last Glacial Maximum. These phylogeographic patterns indicate contrasting demographic histories of populations in close proximity to areas directly affected by glaciers in the Southern Hemisphere during the LGM.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0397-y) contains supplementary material, which is available to authorized users.

Testing hypotheses for genealogical discordance in a rainforest lizard

Genealogical discordance, or when different genes tell distinct stories although they evolved under a shared history, often emerges from either coalescent stochasticity or introgression. In this study, we present a strong case of mito-nuclear genealogical discordance in the Australian rainforest lizard species complex of Saproscincus basiliscus and S. lewisi. One of the lineages that comprises this complex, the Southern S. basiliscus lineage, is deeply divergent at the mitochondrial genome but shows markedly less divergence at the nuclear genome. By placing our results in a comparative context and reconstructing the lineages' demography via multilocus and coalescent-based approximate Bayesian computation methods, we test hypotheses for how coalescent variance and introgression contribute to this pattern. These analyses suggest that the observed genealogical discordance likely results from introgression. Further, to generate such strong discordance, introgression probably acted in concert with other factors promoting asymmetric gene flow between the mitochondrial and nuclear genomes, such as selection or sex-biased dispersal. This study offers a framework for testing sources of genealogical discordance and suggests that historical introgression can be an important force shaping the genetic diversity of species and their populations.

Strong selection against hybrids maintains a narrow contact zone between morphologically cryptic lineages in a rainforest lizard

Phenotypically cryptic lineages comprise an important yet understudied part of biodiversity; in particular, we have much to learn about how these lineages are formed and maintained. To better understand the evolutionary significance of such lineages, we studied a hybrid zone between two morphologically cryptic phylogeographic lineages in the rainforest lizard, Lampropholis coggeri. Analyzing a multilocus genetic dataset through cline inference, individual-based methods and population measures of disequilibrium and using simulations to explore our genetic results in context of theoretical expectations, we inferred the processes maintaining this hybrid zone. We find that these lineages meet in a hybrid zone that is narrow (≈400 m) relative to inferred dispersal rate. Further, the hybrid zone exhibits substantial genetic disequilibrium and sharply coincident and largely concordant clines. Based on our knowledge about the region's biogeography, the species' natural history, and our simulation results, we suggest that strong selection against hybrids structures this system. As all clines show a relatively narrow range of introgression, we posit that this hybrid zone might not yet be in equilibrium. Nonetheless, our results clearly show that phylogeographic lineages can evolve substantial reproductive isolation without concomitant morphological diversification, suggesting that such lineages can constitute a significant component of evolutionary diversity.

What are the consequences of combining nuclear and mitochondrial data for phylogenetic analysis? Lessons from Plethodon salamanders and 13 other vertebrate clades

Background

The use of mitochondrial DNA data in phylogenetics is controversial, yet studies that combine mitochondrial and nuclear DNA data (mtDNA and nucDNA) to estimate phylogeny are common, especially in vertebrates. Surprisingly, the consequences of combining these data types are largely unexplored, and many fundamental questions remain unaddressed in the literature. For example, how much do trees from mtDNA and nucDNA differ? How are topological conflicts between these data types typically resolved in the combined-data tree? What determines whether a node will be resolved in favor of mtDNA or nucDNA, and are there any generalities that can be made regarding resolution of mtDNA-nucDNA conflicts in combined-data trees? Here, we address these and related questions using new and published nucDNA and mtDNA data for Plethodon salamanders and published data from 13 other vertebrate clades (including fish, frogs, lizards, birds, turtles, and mammals).

Results

We find widespread discordance between trees from mtDNA and nucDNA (30-70% of nodes disagree per clade), but this discordance is typically not strongly supported. Despite often having larger numbers of variable characters, mtDNA data do not typically dominate combined-data analyses, and combined-data trees often share more nodes with trees from nucDNA alone. There is no relationship between the proportion of nodes shared between combined-data and mtDNA trees and relative numbers of variable characters or levels of homoplasy in the mtDNA and nucDNA data sets. Congruence between trees from mtDNA and nucDNA is higher on branches that are longer and deeper in the combined-data tree, but whether a conflicting node will be resolved in favor mtDNA or nucDNA is unrelated to branch length. Conflicts that are resolved in favor of nucDNA tend to occur at deeper nodes in the combined-data tree. In contrast to these overall trends, we find that Plethodon have an unusually large number of strongly supported conflicts between data types, which are generally resolved in favor of mtDNA in the combined-data tree (despite the large number of nuclear loci sampled).

Conclusions

Overall, our results from 14 vertebrate clades show that combined-data analyses are not necessarily dominated by the more variable mtDNA data sets. However, given cases like Plethodon, there is also the need for routine checking of incongruence between mtDNA and nucDNA data and its impacts on combined-data analyses.

Comparative multi-locus phylogeography confirms multiple vicariance events in co-distributed rainforest frogs

Though Pleistocene refugia are frequently cited as drivers of species diversification, comparisons of molecular divergence among sister species typically indicate a continuum of divergence times from the Late Miocene, rather than a clear pulse of speciation events at the Last Glacial Maximum. Community-scale inference methods that explicitly test for multiple vicariance events, and account for differences in ancestral effective population size and gene flow, are well suited for detecting heterogeneity of species' responses to past climate fluctuations. We apply this approach to multi-locus sequence data from five co-distributed frog species endemic to the Wet Tropics rainforests of northeast Australia. Our results demonstrate at least two episodes of vicariance owing to climate-driven forest contractions: one in the Early Pleistocene and the other considerably older. Understanding how repeated cycles of rainforest contraction and expansion differentially affected lineage divergence among co-distributed species provides a framework for identifying evolutionary processes that underlie population divergence and speciation.

Phylogeography and demographic history of Lacerta lepida in the Iberian Peninsula: multiple refugia, range expansions and secondary contact zones

BACKGROUND

The Iberian Peninsula is recognized as an important refugial area for species survival and diversification during the climatic cycles of the Quaternary. Recent phylogeographic studies have revealed Iberia as a complex of multiple refugia. However, most of these studies have focused either on species with narrow distributions within the region or species groups that, although widely distributed, generally have a genetic structure that relates to pre-Quaternary cladogenetic events. In this study we undertake a detailed phylogeographic analysis of the lizard species, Lacerta lepida, whose distribution encompasses the entire Iberian Peninsula. We attempt to identify refugial areas, recolonization routes, zones of secondary contact and date demographic events within this species.

RESULTS

Results support the existence of 6 evolutionary lineages (phylogroups) with a strong association between genetic variation and geography, suggesting a history of allopatric divergence in different refugia. Diversification within phylogroups is concordant with the onset of the Pleistocene climatic oscillations. The southern regions of several phylogroups show a high incidence of ancestral alleles in contrast with high incidence of recently derived alleles in northern regions. All phylogroups show signs of recent demographic and spatial expansions. We have further identified several zones of secondary contact, with divergent mitochondrial haplotypes occurring in narrow zones of sympatry.

CONCLUSIONS

The concordant patterns of spatial and demographic expansions detected within phylogroups, together with the high incidence of ancestral haplotypes in southern regions of several phylogroups, suggests a pattern of contraction of populations into southern refugia during adverse climatic conditions from which subsequent northern expansions occurred. This study supports the emergent pattern of multiple refugia within Iberia but adds to it by identifying a pattern of refugia coincident with the southern distribution limits of individual evolutionary lineages. These areas are important in terms of long-term species persistence and therefore important areas for conservation.

Microsatellite analysis supports mitochondrial phylogeography of the hellbender (Cryptobranchus alleganiensis)

We investigated genetic diversity of the hellbender (Cryptobranchus alleganiensis) throughout its range in the eastern US using nuclear markers and compared our results to a previously published mitochondrial analysis. A variety of nuclear markers, including protein-coding gene introns and microsatellites were tested but only microsatellites were variable enough for population level analysis. Microsatellite loci showed moderate among population sharing of alleles, in contrast to the reciprocal monophyly exhibited by mitochondrial DNA. However, analyses using F-statistics and Bayesian clustering algorithms showed considerable population subdivision and clustered hellbender populations into the same major groups as the mtDNA. The microsatellites combined with the mtDNA data suggest that gene flow is severely restricted or non-existent among eight major groups, and potentially among populations (rivers) within groups. The combined mtDNA and microsatellite data suggest that the currently recognized hellbender subspecies are paraphyletic. We suggest that the eight independent groups identified in our study should be managed as such, rather than basing conservation decisions on the two named subspecies of hellbender.

Patterns of persistence and isolation indicate resilience to climate change in montane rainforest lizards

Globally, montane tropical diversity is characterized by extraordinary local endemism that is not readily explained by current environmental variables indicating a strong imprint of history. Montane species often exist as isolated populations under current climatic conditions and may have remained isolated throughout recent climatic cycles, leading to substantial genetic and phenotypic divergence. Alternatively, populations may have become contiguous during colder climates resulting in less divergence. Here we compare responses to historical climate fluctuation in a montane specialist skink, Lampropholis robertsi, and its more broadly distributed congener, L. coggeri, both endemic to rainforests of northeast Australia. To do so, we combine spatial modelling of potential distributions under representative palaeoclimates, multi-locus phylogeography and analyses of phenotypic variation. Spatial modelling of L. robertsi predicts strong isolation among disjunct montane refugia during warm climates, but with potential for localized exchange during the most recent glacial period. In contrast, predicted stable areas are more widespread and connected in L. coggeri. Both species exhibit pronounced phylogeographic structuring for mitochondrial and nuclear genes, attesting to low dispersal and high persistence across multiple isolated regions. This is most prominent in L. robertsi, for which coalescent analyses indicate that most populations persisted in isolation throughout the climate cycles of the Pleistocene. Morphological divergence, principally in body size, is more evident among isolated populations of L. robertsi than L. coggeri. These results highlight the biodiversity value of isolated montane populations and support the general hypothesis that tropical montane regions harbour high levels of narrow-range taxa because of their resilience to past climate change.

Elucidation of cryptic diversity in a widespread nearctic treefrog reveals episodes of mitochondrial gene capture as frogs diversified across a dynamic landscape

We investigate the evolutionary history of the wide-ranging Nearctic treefrog Hyla arenicolor through the integration of extensive range-wide sampling, phylogenetic analyses of multilocus genetic data, and divergence dating. Previous phylogeographic studies of this frog documented a potential signature of introgressive hybridization from an ecologically and morphologically divergent sister species. Based on our Bayesian phylogenetic analyses of mitochondrial DNA, we inferred strong phylogeographic structure in H. arenicolor as indicated by seven well-supported clades, five of which correspond to well-defined biogeographic regions. Clades from the Balsas Basin and southwestern Central Mexican Plateau in Mexico, and the Grand Canyon of Arizona, group with the morphologically, behaviorally, and ecologically divergent mountain treefrogs in the H. eximia group, rendering H. arenicolor as paraphyletic. The phylogenetic position of at least two of these three H. arenicolor clades within the H. eximia group, however, is most likely the result of several episodes of introgressive hybridization and subsequent mitochondrial gene capture separated in time and space, as supported by evidence from the nuclear genes. Hyla arenicolor from the Balsas Basin appear to be deeply divergent from other H. arenicolor and represent a distinctly different species. Results suggests that introgressive hybridization events, both ancient and contemporary, coupled with late Neogene vicariance and Pleistocene climate-driven range shifts, have all played a role in the historical diversification of H. arenicolor.

The evolution of two partner LINE/SINE families and a full-length chromodomain-containing Ty3/Gypsy LTR element in the first reptilian genome of Anolis carolinensis

Transposable elements have been characterized in a number of vertebrates, including whole genomes of mammals, birds, and fishes. The Anolis carolinensis draft assembly provides the first opportunity to study retroposons in a reptilian genome. Here, we identified and reconstructed a number of retroposons based on database searches: Five Sauria short interspersed element (SINE) subfamilies, 5S-Sauria SINE chimeras, Anolis Bov-B long interspersed element (LINE), Anolis SINE 2, Anolis LINE 2, Anolis LINE 1, Anolis CR 1, and a chromodomain-containing Ty3/Gypsy LTR element. We focused on two SINE families (Anolis Sauria SINE and Anolis SINE 2) and their partner LINE families (Anolis Bov-B LINE and Anolis LINE 2). We demonstrate that each SINE/LINE pair is distributed similarly and predict that the retrotransposition of evolutionarily younger Sauria SINE members is via younger Bov-B LINE members while a correlation also exists between their respective evolutionarily older SINE/LINE members. The evolutionarily youngest Sauria SINE sequences evolved as part of novel rolling-circle transposons. The evolutionary time frame when Bov-B LINEs and Sauria SINEs were less active in their retrotransposition is characterized by a high retrotransposition burst of Anolis SINE 2 and Anolis LINE 2 elements. We also characterized the first full-length chromoviral LTR element in amniotes (Amn-ichi). This newly identified chromovirus is widespread in the Anolis genome and has been very well preserved, indicating that it is still active. Transposable elements in the Anolis genome account for approximately 20% of the total DNA sequence, whereas the proportion is more than double that in many mammalian genomes in which such elements have important biological functions. Nevertheless, 20% transposable element coverage is sufficient to predict that Anolis retroposons and other mobile elements also may have biologically and evolutionarily relevant functions. The new SINEs and LINEs and other ubiquitous genomic elements characterized in the Anolis genome will prove very useful for studies in comparative genomics, phylogenetics, and functional genetics.

Evidence for differential assortative female preference in association with refugial isolation of rainbow skinks in Australia's tropical rainforests

BACKGROUND

Divergence driven by female preference can give rise to pre-mating isolation more rapidly than post-mating isolation can evolve through the accumulation of allelic incompatibilities. Moreover pre-mating isolation may be more effective at maintaining morphological differentiation between divergent populations. In the context of Australian rainforest endemic skinks that were historically subjected to refugial isolation, this study examined the following predictions: 1) that assortative female preference is associated with more recent divergence of southern C. rubrigularis (S-RED) and C. rhomboidalis (BLUE), but not with deeply divergent S-RED and northern C. rubrigularis (N-RED); and 2) that upon secondary contact, morphological differentiation is maintained between S-RED and BLUE, whereas N-RED and S-RED remain morphogically indistinguishable.

PRINCIPAL FINDINGS

Female preference trials found no evidence for assortative female preference between N-RED and S-RED, supporting a previous genetic hybrid zone study which inferred post-mating but no pre-mating isolation. In contrast there is evidence for assortative female preference between S-RED and BLUE, with BLUE females preferring to associate with BLUE males, but S-RED females showing no preference. Multi-locus coalescent analyses, used to estimate post-divergence gene-flow between proximally located S-RED and BLUE populations, rejected zero gene-flow from BLUE to S-RED and thus RED and BLUE have maintained morphological differentiation despite secondary contact. Morphometric analyses confirmed a lack of morphological divergence between N-RED and S-RED and established that BLUE is morphologically divergent from RED in traits other than throat colour.

CONCLUSIONS/SIGNIFICANCE

Long-term isolation has been sufficient to generate post-mating isolation but no morphological divergence between N-RED and S-RED. In contrast, greater morphological differentiation is associated with evidence for assortative female preference between more recently diverged S-RED and BLUE. Combined with previous estimates of lineage-wide gene flow, these results are consistent with the suggestion that assortative female preference is more effective than post-mating isolation in maintaining morphological differentiation between divergent populations.

Molecular phylogenetics of the lizard genus Microlophus (squamata:tropiduridae): aligning and retrieving indel signal from nuclear introns

We use a multigene data set (the mitochondrial locus and nine nuclear gene regions) to test phylogenetic relationships in the South American "lava lizards" (genus Microlophus) and describe a strategy for aligning noncoding sequences that accounts for differences in tempo and class of mutational events. We focus on seven nuclear introns that vary in size and frequency of multibase length mutations (i.e., indels) and present a manual alignment strategy that incorporates insertions and deletions (indels) for each intron. Our method is based on mechanistic explanations of intron evolution that does not require a guide tree. We also use a progressive alignment algorithm (Probabilistic Alignment Kit; PRANK) and distinguishes insertions from deletions and avoids the "gapcost" conundrum. We describe an approach to selecting a guide tree purged of ambiguously aligned regions and use this to refine PRANK performance. We show that although manual alignment is successful in finding repeat motifs and the most obvious indels, some regions can only be subjectively aligned, and there are limits to the size and complexity of a data matrix for which this approach can be taken. PRANK alignments identified more parsimony-informative indels while simultaneously increasing nucleotide identity in conserved sequence blocks flanking the indel regions. When comparing manual and PRANK with two widely used methods (CLUSTAL, MUSCLE) for the alignment of the most length-variable intron, only PRANK recovered a tree congruent at deeper nodes with the combined data tree inferred from all nuclear gene regions. We take this concordance as an objective function of alignment quality and present a strongly supported phylogenetic hypothesis for Microlophus relationships. From this hypothesis we show that (1) a coded indel data partition derived from the PRANK alignment contributed significantly to nodal support and (2) the indel data set permitted detection of significant conflict between mitochondrial and nuclear data partitions, which we hypothesize arose from secondary contact of distantly related taxa, followed by hybridization and mtDNA introgression.

Comparing alignment methods for inferring the history of the new world lizard genus Mabuya (Squamata: Scincidae)

The rapid increase in the ability to generate molecular data, and the focus on model-based methods for tree reconstruction have greatly advanced the use of phylogenetics in many fields. The recent flurry of new analytical techniques has focused almost solely on tree reconstruction, whereas alignment issues have received far less attention. In this paper, we use a diverse sampling of gene regions from lizards of the genus Mabuya to compare the impact, on phylogeny estimation, of new maximum likelihood alignment algorithms with more widely used methods. Sequences aligned under different optimality criteria are analyzed using partitioned Bayesian analysis with independent models and parameter settings for each gene region, and the most strongly supported phylogenetic hypothesis is then used to test the hypothesis of two colonizations of the New World by African scincid lizards. Our results show that the consistent use of model-based methods in both alignment and tree reconstruction leads to trees with more optimal likelihood scores than the use of independent criteria in alignment and tree reconstruction. We corroborate and extend earlier evidence for two independent colonizations of South America by scincid lizards. Relationships within South American Mabuya are found to be in need of taxonomic revision, specifically complexes under the names M. heathi, M. agilis, and M. bistriata (sensu, M.T. Rodrigues, Papeis Avulsos de Zoologia 41 (2000) 313).

Reinforcement drives rapid allopatric speciation

Allopatric speciation results from geographic isolation between populations. In the absence of gene flow, reproductive isolation arises gradually and incidentally as a result of mutation, genetic drift and the indirect effects of natural selection driving local adaptation. In contrast, speciation by reinforcement is driven directly by natural selection against maladaptive hybridization. This gives individuals that choose the traits of their own lineage greater fitness, potentially leading to rapid speciation between the lineages. Reinforcing natural selection on a population of one of the lineages in a mosaic contact zone could also result in divergence of the population from the allopatric range of its own lineage outside the zone. Here we test this with molecular data, experimental crosses, field measurements and mate choice experiments in a mosaic contact zone between two lineages of a rainforest frog. We show that reinforcing natural selection has resulted in significant premating isolation of a population in the contact zone not only from the other lineage but also, incidentally, from the closely related main range of its own lineage. Thus we show the potential for reinforcement to drive rapid allopatric speciation.

Comparative phylogeographic summary statistics for testing simultaneous vicariance

Testing for simultaneous vicariance across comparative phylogeographic data sets is a notoriously difficult problem hindered by mutational variance, the coalescent variance, and variability across pairs of sister taxa in parameters that affect genetic divergence. We simulate vicariance to characterize the behaviour of several commonly used summary statistics across a range of divergence times, and to characterize this behaviour in comparative phylogeographic datasets having multiple taxon-pairs. We found Tajima's D to be relatively uncorrelated with other summary statistics across divergence times, and using simple hypothesis testing of simultaneous vicariance given variable population sizes, we counter-intuitively found that the variance across taxon pairs in Nei and Li's net nucleotide divergence (pi(net)), a common measure of population divergence, is often inferior to using the variance in Tajima's D across taxon pairs as a test statistic to distinguish ancient simultaneous vicariance from variable vicariance histories. The opposite and more intuitive pattern is found for testing more recent simultaneous vicariance, and overall we found that depending on the timing of vicariance, one of these two test statistics can achieve high statistical power for rejecting simultaneous vicariance, given a reasonable number of intron loci (> 5 loci, 400 bp) and a range of conditions. These results suggest that components of these two composite summary statistics should be used in future simulation-based methods which can simultaneously use a pool of summary statistics to test comparative the phylogeographic hypotheses we consider here.