Phylogenetic analyses of mtDNA sequences reveal three cryptic lineages in the widespread neotropical frog Leptodactylus fuscus (Schneider, 1799) (Anura, Leptodactylidae)

Leptodactylus fuscus is a neotropical frog ranging from Panamá to Argentina, to the east of the Andes mountains, and also inhabiting Margarita, Trinidad, and the Tobago islands. We performed phylogenetic analyses of 12S rRNA, 16S rRNA, tRNA-Leu, and ND1 mitochondrial (mt) DNA sequences from specimens collected across the geographic distribution of L. fuscus to examine two alternative hypotheses: (i) L. fuscus is a single, widely distributed species, or (ii) L. fuscus is a species complex. We tested statistically for geographic association and partitioning of genetic variation among mtDNA clades. The mtDNA data supported the hypothesis of several cryptic species within L. fuscus. Unlinked mtDNA and nuclear markers supported independently the distinctness of a 'northern' phylogenetic unit. In addition, the mtDNA data divided the southern populations into two clades that showed no sister relationship to each other, consistent with high differentiation and lack of gene flow among southern populations as suggested by allozyme data. Concordance between mtDNA and allozyme patterns suggests that cryptic speciation has occurred in L. fuscus without morphological or call differentiation. This study illustrates a case in which lineage splitting during the speciation process took place without divergence in reproductive isolation mechanisms (e.g. advertisement call in frogs), contrary to expectations predicted using a biological species framework.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society, 2006, 87 , 325–341. No claim to original US government works.     

Isolation and characterization of 20 polymorphic microsatellite loci in the migratory freshwater fish Leporinus obtusidens (Characiformes: Anostomidae) using 454 shotgun pyrosequencing

Twenty polymorphic microsatellite markers were developed for the Neotropical fish Leporinus obtusidens using a next generation sequencing approach and tested in two other characifomes species, Schizodon platae and Prochilodus lineatus. Microsatellite loci alleles in L. obtusidens ranged between 2 and 20 alleles per locus (mean = 5·7), with expected heterozygosity values ranging from 0·097 to 0·956 (mean = 0·578) and observed heterozygosity values ranging from 0·000 to 0·800 (mean = 0·400) in a sample of 20 specimens from the lower Paraná River (Argentina). Most of these markers will be a valuable tool for captive breeding and stocking programmes, as well as for analyses of population connectivity and genetic structure in this broadly distributed Neotropical migratory fish.     

Higher genetic diversity in recolonized areas than in refugia of Alnus glutinosa triggered by continent‐wide lineage admixture

Genetic admixture is supposed to be an important trigger of species expansions because it can create the potential for selection of genotypes suitable for new climatic conditions. Up until now, however, no continent‐wide population genetic study has performed a detailed reconstruction of admixture events during natural species expansions. To fill this gap, we analysed the postglacial history of Alnus glutinosa, a keystone species of European swamp habitats, across its entire distribution range using two molecular markers, cpDNA and nuclear microsatellites. CpDNA revealed multiple southern refugia located in the Iberian, Apennine, Balkan and Anatolian Peninsulas, Corsica and North Africa. Analysis of microsatellites variation revealed three main directions of postglacial expansion: (i) from the northern part of the Iberian Peninsula to Western and Central Europe and subsequently to the British Isles, (ii) from the Apennine Peninsula to the Alps and (iii) from the eastern part of the Balkan Peninsula to the Carpathians followed by expansion towards the Northern European plains. This challenges the classical paradigm that most European populations originated from refugial areas in the Carpathians. It has been shown that colonizing lineages have met several times and formed secondary contact zones with unexpectedly high population genetic diversity in Central Europe and Scandinavia. On the contrary, limited genetic admixture in southern refugial areas of A. glutinosa renders rear‐edge populations in the Mediterranean region more vulnerable to extinction due to climate change.     

The aggregate site frequency spectrum for comparative population genomic inference

Understanding how assemblages of species responded to past climate change is a central goal of comparative phylogeography and comparative population genomics, an endeavour that has increasing potential to integrate with community ecology. New sequencing technology now provides the potential to perform complex demographic inference at unprecedented resolution across assemblages of nonmodel species. To this end, we introduce the aggregate site frequency spectrum (aSFS), an expansion of the site frequency spectrum to use single nucleotide polymorphism (SNP) data sets collected from multiple, co‐distributed species for assemblage‐level demographic inference. We describe how the aSFS is constructed over an arbitrary number of independent population samples and then demonstrate how the aSFS can differentiate various multispecies demographic histories under a wide range of sampling configurations while allowing effective population sizes and expansion magnitudes to vary independently. We subsequently couple the aSFS with a hierarchical approximate Bayesian computation (hABC) framework to estimate degree of temporal synchronicity in expansion times across taxa, including an empirical demonstration with a data set consisting of five populations of the threespine stickleback (Gasterosteus aculeatus). Corroborating what is generally understood about the recent postglacial origins of these populations, the joint aSFS/hABC analysis strongly suggests that the stickleback data are most consistent with synchronous expansion after the Last Glacial Maximum (posterior probability = 0.99). The aSFS will have general application for multilevel statistical frameworks to test models involving assemblages and/or communities, and as large‐scale SNP data from nonmodel species become routine, the aSFS expands the potential for powerful next‐generation comparative population genomic inference.     

Phylogeographic analysis of Pimoidae (Arachnida: Araneae) inferred from mitochondrial cytochrome c oxidase subunit I and nuclear 28S rRNA gene regions

Using mitochondrial DNA cytochrome c oxidase subunit I and nuclear DNA 28S rRNA data, we explored the phylogenetic relationships of the family Pimoidae (Arachnida: Araneae) and tested the North America to Asia dispersal hypothesis. Sequence data were analysed using maximum parsimony and Bayesian inference. A phylogenetic analysis suggested that vicariance, instead of dispersal, better explained the present distribution pattern of Pimoidae. Times of divergence events were estimated using penalized likelihood method. The dating analysis suggested that the emergence time of Pimoidae was approximately 140 million years ago (Ma). The divergence time of the North American and Asian species of Pimoa was approximately 110 Ma. Our phylogenetic hypothesis supports the current morphology‐based taxonomy and suggests that the cave dwelling might have played an important role in the speciation of pimoids in arid areas.     

Molecular evidence for Pleistocene glacial cycles driving diversification of a North American desert spider, Agelenopsis aperta

The influence of historical climatic vs. geological changes on species diversification patterns was investigated in a widely distributed North American desert spider, Agelenopsis aperta (Araneae: Agelenidae), with particular reference to Pleistocene glacial cycles and earlier patterns of mountain building. Levels of sequence divergence obtained from the mitochondrial gene, cytochrome oxidase I, dated to the Pleistocene, eliminating Rocky Mountain orogeny as a cause of diversification, as orogeny ended 4 million years ago. The results of phylogenetic and network analyses showed the presence of three geographically defined clades, which were consistent with the presence of at least three glacial refugia: (i) east of the Rocky Mountains; (ii) between the Rocky Mountains and Sierra Nevadas; and (iii) west of the Sierra Nevadas. In addition, populations within the Rocky Mountains exhibited significantly lower genetic diversity than populations east of the Rocky Mountains and the haplotypes found within the Rockies were a subset of eastern haplotypes. These patterns suggest that a post-Pleistocene range expansion occurred out of an eastern glacial refugium into the Rocky Mountains. Examination of phylogeographical studies of other North American desert taxa indicated that mountain building explained diversification patterns more effectively for some taxa but Pleistocene climate change was more important for others, including A. aperta.     

Comparative phylogeography of salmonid fishes (Salmonidae) reveals late to post-Pleistocene exchange between three now-disjunct river basins in Siberia

Abstract. We use a comparative phylogeographical framework to evaluate the hypothesis of hydrological exchange during the Pleistocene among the now disjunct Lena, Amur, and Enisei basins in Siberia, and to provide evidence on the causal mechanism of their present day faunal dissimilarities. Approximately 600 bases of the mitochondrial control region were sequenced in five distinct lineages among three genera of salmonid fishes, Hucho, Brachymystax and Thymallus. All three basins were fixed for divergent (2–5.4%) lineages of Thymallus whereas a single shared haplotype was present in all three basins for Hucho taimen ( Pallas, 1773 ) and one shared haplotype between the Lena and Amur basins out of a total of five for blunt‐snouted and one out of five for sharp‐snouted Brachymystax lenok ( Pallas, 1773 ). For both blunt‐ and sharp‐snouted lenok the haplotypes found within each basin did not form clades, so no relationship between genotypes and geographical occurrence was found. Our data support relatively recent hydrological mixing of the major river drainage systems in eastern and far‐eastern Siberia, congruent with the hypothesis of large‐scale palaeo‐hydrological exchange stemming from glacial advance, retreat and melting during Pleistocene climate fluctuations. Furthermore, these results in conjunction with a comparison of overall faunal composition suggest that environmental differences rather than historical contingency may be responsible for the faunal dissimilarities of the Amur, Lena, and Enisei river basins.     

Deep phylogeographical structure and parallel host range evolution in the leaf beetle Agelasa nigriceps

To understand the mechanisms behind the diversification of herbivorous insects through insect–plant interactions, it is important to know how the insects change their diet breadth in response to environmental changes. In this study, we investigated the phylogeographical pattern of the leaf beetle Agelasa nigriceps to infer the evolutionary history of its host range. While this beetle commonly uses Actinidia arguta (Actinidiaceae) as a host plant, it has been recorded recently on Pterostyrax hispidus (Styracaceae), which is now increasing in abundance at some localities in Japan due to the indirect effects of high population size of a mammalian herbivore. Considerable variation among populations in the ability of Ag. nigriceps to use P. hispidus suggests that P. hispidus is a newly acquired host plant for this beetle. Phylogenetic analyses using mitochondrial DNA sequences and amplified fragment length polymorphism (AFLP) revealed a high degree of phylogeographical structure in Ag. nigriceps throughout Japan, which is consistent with the hypothesis that several glacial refugia existed in the Japanese archipelago. In contrast, no genetic structure associated with the host plants was detected. Both the mitochondrial DNA and AFLP analyses showed that populations that can use P. hispidus are polyphyletic. These results and geographical variation in host use suggest that the host range expansion to a novel host, P. hispidus, is a very recent and possibly ongoing phenomenon and has occurred independently in several regions. Our study illustrates that the host range of herbivorous insects can evolve repeatedly in response to similar environmental changes.     

Multiple causation of phylogeographical pattern as revealed by nested clade analysis of the bamboo viper (Trimeresurus stejnegeri) within Taiwan

In order to assess the utility of nested clade analysis, both standard phylogenetic algorithms and nested clade analysis were performed on a geographically widespread survey of mitochondrial DNA haplotypes of the bamboo viper, Trimeresurus stejnegeri, within Taiwan. Gross tree topologies were congruent for all analyses and indicated the presence of two geographically overlapping clades within Taiwan. The smaller lineage was restricted to the north and east coasts, whereas the larger lineage occupied all but the northern range of the species within Taiwan including the Pacific offshore populations of Green and Orchid Islands. The phylogeographical pattern supports the existence of at least one colonization event from the continent since the initial isolation of Taiwan from the mainland in the Pliocene. However, determining the exact number of colonization events was not possible due to the simultaneous vicariant forces of hypothesized continental landbridge connections and the occurrence of dramatic in situ orogenesis throughout the Pleistocene. Nested clade analysis provided multiple temporal and spatial population historical inferences that are not possible with standard analyses and therefore should become widely applied to future phylogeographical studies.