GEOLOGICAL AND CLIMATIC FORCES DRIVING SPECIATION IN THE CONTINENTALLY DISTRIBUTED TRILLING CHORUS FROGS (PSEUDACRIS)

Biogeography and the evolution of acoustic communication in the polyploid North American grey treefrog complex

After polyploid species are formed, interactions between diploid and polyploid lineages may generate additional diversity in novel cytotypes and phenotypes. In anurans, mate choice by acoustic communication is the primary method by which individuals identify their own species and assess suitable mates. As such, the evolution of acoustic signals is an important mechanism for contributing to reproductive isolation and diversification in this group. Here, we estimate the biogeographical history of the North American grey treefrog complex, consisting of the diploid Hyla chrysoscelis and the tetraploid Hyla versicolor, focusing specifically on the geographical origin of whole genome duplication and the expansion of lineages out of glacial refugia. We then test for lineage-specific differences in mating signals by applying comparative methods to a large acoustic data set collected over 52 years that includes >1500 individual frogs. Along with describing the overall biogeographical history and call diversity, we found evidence that the geographical origin of H. versicolor and the formation of the midwestern polyploid lineage are both associated with glacial limits, and that the southwestern polyploid lineage is associated with a shift in acoustic phenotype relative to the diploid lineage with which they share a mitochondrial lineage. In H. chrysoscelis, we see that acoustic signals are largely split by Eastern and Western lineages, but that northward expansion along either side of the Appalachian Mountains is associated with further acoustic diversification. Overall, results of this study provide substantial clarity on the evolution of grey treefrogs as it relates to their biogeography and acoustic communication.

Signatures of north-eastern expansion and multiple refugia: genomic phylogeography of the Pine Barrens tree frog, Hyla andersonii (Anura: Hylidae)

Range fragmentation poses challenges for species persistence over time and can be caused by both historical and contemporary processes. We combined genomic data, phylogeographical model testing and palaeoclimatic niche modelling to infer the evolutionary history of the Pine Barrens tree frog (Hyla andersonii), a seepage bog specialist, in eastern North America to gain a better understanding of the historical context of its fragmented distribution. We sampled H. andersonii populations across the three disjunct regions of the species range: Alabama/Florida (AF), the Carolinas (CL) and New Jersey (NJ). Phylogenetic relationships within H. andersonii were consistent between the nuclear species tree and mitochondrial analyses, indicating divergence between AF and CL/NJ (Atlantic clade) ~0.9 Mya and divergence of the NJ clade ~0.15 Mya. Several predictions of north-eastern expansion along the Atlantic coast were supported by phylogeographical analyses. Model testing using genome-wide single nucleotide polymorphism data and species distribution models both provided evidence for multiple disjunct refugia. This comprehensive phylogeographical study of H. andersonii demonstrates a long history of range fragmentation within an endemic coastal plain species and highlights the influence of historical climate change on the current distribution of species and their genetic diversity.

Hierarchical Hybrid Enrichment: Multitiered Genomic Data Collection Across Evolutionary Scales, With Application to Chorus Frogs (Pseudacris)

Determining the optimal targets of genomic subsampling for phylogenomics, phylogeography, and population genomics remains a challenge for evolutionary biologists. Of the available methods for subsampling the genome, hybrid enrichment (sequence capture) has become one of the primary means of data collection for systematics, due to the flexibility and cost efficiency of this approach. Despite the utility of this method, information is lacking as to what genomic targets are most appropriate for addressing questions at different evolutionary scales. In this study, first, we compare the benefits of target loci developed for deep- and shallow scales by comparing these loci at each of three taxonomic levels: within a genus (phylogenetics), within a species (phylogeography), and within a hybrid zone (population genomics). Specifically, we target evolutionarily conserved loci that are appropriate for deeper phylogenetic scales and more rapidly evolving loci that are informative for phylogeographic and population genomic scales. Second, we assess the efficacy of targeting multiple-locus sets for different taxonomic levels in the same hybrid enrichment reaction, an approach we term hierarchical hybrid enrichment. Third, we apply this approach to the North American chorus frog genus Pseudacris to answer key evolutionary questions across taxonomic and temporal scales. We demonstrate that in this system the type of genomic target that produces the most resolved gene trees differs depending on the taxonomic level, although the potential for error is substantially lower for the deep-scale loci at all levels. We successfully recover data for the two different locus sets with high efficiency. Using hierarchical data targeting deep and shallow levels: we 1) resolve the phylogeny of the genus Pseudacris and introduce a novel visual and hypothesis testing method that uses nodal heat maps to examine the robustness of branch support values to the removal of sites and loci; 2) estimate the phylogeographic history of Pseudacris feriarum, which reveals up to five independent invasions leading to sympatry with congener Pseudacris nigrita to form replicated reinforcement contact zones with ongoing gene flow into sympatry; and 3) quantify with high confidence the frequency of hybridization in one of these zones between P. feriarum and P. nigrita, which is lower than microsatellite-based estimates. We find that the hierarchical hybrid enrichment approach offers an efficient, multitiered data collection method for simultaneously addressing questions spanning multiple evolutionary scales. [Anchored hybrid enrichment; heat map; hybridization; phylogenetics; phylogeography; population genomics; reinforcement; reproductive character displacement.].

Cytonuclear discordance, reticulation and cryptic diversity in one of North America's most common frogs

Complicated phylogenetic histories benefit from diverse sources of inference. Pseudacris crucifer (spring peeper) spans most of eastern North America and comprises six mtDNA lineages that form multiple contact zones. The putative Miocene or early Pliocene origins of the oldest lineages within Pseudacris crucifer imply sufficient time for species-level divergence. To understand why this species appears unified while congeners have radiated, we analyze and compare male advertisement calls, mitochondrial, and nuclear markers and speak to the complex processes that have potentially influenced its contemporary patterns. We find extensive geographic and topological mitonuclear discordance, with three nuclear lineages containing 6 more-structured mtDNA lineages, and nuclear introgression at some contact zones. Male advertisement call differentiation is incongruent with the genetic structure as only one lineage appears differentiated. Occupying the Interior Highlands of the central United States, this Western lineage also has the most concordant mitochondrial and nuclear geographic patterns. Based on our findings we suggest that the antiquity of common ancestors was not as important as the maintenance of allopatry in the divergence in P. crucifer genetic lineages. We use multiple lines of evidence to generate hypotheses of isolation, reticulation, and discordance within this species and to expand our understanding of the early stages of speciation.

Signature of the Paleo-Course Changes in the São Francisco River as Source of Genetic Structure in Neotropical Pithecopus nordestinus (Phyllomedusinae, Anura) Treefrog

Historical processes that have interrupted gene flow between distinct evolutionary lineages have played a fundamental role in the evolution of the enormous diversity of species found in the Neotropical region. Numerous studies have discussed the role of geographic barriers and Pleistocene forest refugia in the diversification of the region's biodiversity. In the present study, we investigated the relative contribution of these different factors to the evolutionary history of Pithecopus nordestinus, a Neotropical tree frog, which is amply distributed in the Brazilian Atlantic Forest and adjacent areas of the Caatinga biome. We used an extensive sample and multilocus DNA sequences to provide an overview of the intraspecific genetic diversity of P. nordestinus, characterize historical diversification patterns, and identify possible phylogenetic splits. We tested different scenarios of diversification based on Pleistocene Refugia and river barrier models using approximate Bayesian computation (ABC) and ecological niche modeling (ENM). The phylogenetic approach indicate the occurrence of processes of phylogeographic divergence in both time and space, related to historical shifts in the course of the São Francisco River during Plio-Pleistocene period, resulting in two principal, highly divergent clades. The ABC model provided strong statistical support for this scenario, confirming the hypothesis that the São Francisco River acted as an effective geographical barrier during vicariant events in the evolutionary history of P. nordestinus. We believe that the climatic changes that occurred during the Pleistocene also played a secondary role in the genetic signatures identified, reinforcing the divergence of populations isolated by physical barriers. These findings reinforce the conclusion that the two models of diversification (geographic barriers and refugia) are not mutually exclusive in the Neotropical domain but may interact extensively during the diversification of species on a regional scale.

Geographic separation and genetic differentiation of populations are not coupled with niche differentiation in threatened Kaiser's spotted newt (Neurergus kaiseri)

The combination of niche modelling and landscape genetics (genomics) helps to disentangle processes that have shaped population structure in the evolutionary past and presence of species. Herein, we integrate a comprehensive genomic dataset with ecological parameters and niche modelling for the threatened Kaiser's newt, a newt species adapted to mountain spring-ponds in Iran. Genomic analysis suggests the existence of two highly differentiated clades North and South of the Dez River. Genetic variation between the two clades (76.62%) was much greater than within clades (16.25%), suggesting that the Dez River prevented gene flow. River disconnectivity, followed by geographic distance, contributed mostly to genetic differentiation between populations. Environmental niche and landscape resistance had no significant influence. Though a significant difference between climatic niches occupied by each clade at the landscape-scale, habitat niches at the local-scale were equivalent. 'Niche similarity analysis' supported niche conservatism between the two clades despite the southward shift in the climatic niche of the Southern clade. Accordingly, populations of different clades may occupy different climatic niches within their ancestral niche. Our results indicate that the change of climatic conditions of geographically and genetically separated populations does not necessarily result in the shift of an ecological niche.

Diversification of the widespread neotropical frog Physalaemus cuvieri in response to Neogene-Quaternary geological events and climate dynamics

Here we reconstructed the demographical history and the dispersal dynamics of Physalemus cuvieri through the Neogene-Quaternary periods by coupling DNA regions with different mutation rates, ecological niche modelling, reconstruction of spatio-temporal lineage dispersal and coalescent simulations. Still, to test alternative diversification scenarios we used approximate Bayesian computation. Molecular phylogenetic analysis recovered four deep and strongly supported clades, which we interpret as population lineages. The ancestral location reconstruction placed the root in southcentral Amazonia, and the dispersal events indicate that spatial displacement was widespread early in the diversification of this species. The demographical scenario of "Multiple Refugia" with recent lineage admixture was the most likely hypothesis to predict the observed genetic parameters of P. cuvieri. Our results revealed that Neogene orogenic events might have played a prominent role in the early diversification of P. cuvieri. The species shows deep divergences with strong regional population structure, despite its widespread distribution. Final uplift of the central Brazilian Plateau and formation of the river basins in Central South America played an important role in the origin, diversification and the maintenance of P. cuvieri lineages.

Geographic variation in hybridization across a reinforcement contact zone of chorus frogs (Pseudacris)

Reinforcement contact zones, which are secondary contact zones where species are diverging in reproductive behaviors due to selection against hybridization, represent natural laboratories for studying speciation‐in‐action. Here, we examined replicate localities across the entire reinforcement contact zone between North American chorus frogs Pseudacris feriarum and P. nigrita to investigate geographic variation in hybridization frequencies and to assess whether reinforcement may have contributed to increased genetic divergence within species. Previous work indicated these species have undergone reproductive character displacement (RCD) in male acoustic signals and female preferences due to reinforcement. We also examined acoustic signal variation across the contact zone to assess whether signal characteristics reliably predict hybrid index and to elucidate whether the degree of RCD predicts hybridization rate. Using microsatellites, mitochondrial sequences, and acoustic signal information from >1,000 individuals across >50 localities and ten sympatric focal regions, we demonstrate: (1) hybridization occurs and (2) varies substantially across the geographic range of the contact zone, (3) hybridization is asymmetric and in the direction predicted from observed patterns of asymmetric RCD, (4) in one species, genetic distance is higher between conspecific localities where one or both have been reinforced than between nonreinforced localities, after controlling for geographic distance, (5) acoustic signal characters strongly predict hybrid index, and (6) the degree of RCD does not strongly predict admixture levels. By showing that hybridization occurs in all sympatric localities, this study provides the fifth and final line of evidence that reproductive character displacement is due to reinforcement in the chorus frog contact zone. Furthermore, this work suggests that the dual action of cascade reinforcement and partial geographic isolation is promoting genetic diversification within one of the reinforced species.

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.

First estimates of the probability of survival in a small-bodied, high-elevation frog (Boreal Chorus Frog, Pseudacris maculata), or how historical data can be useful

In an era of shrinking budgets yet increasing demands for conservation, the value of existing (i.e., historical) data are elevated. Lengthy time series on common, or previously common, species are particularly valuable and may be available only through the use of historical information. We provide first estimates of the probability of survival and longevity (0.67–0.79 and 5–7 years, respectively) for a subalpine population of a small-bodied, ostensibly common amphibian, the Boreal Chorus Frog (Pseudacris maculata (Agassiz, 1850)), using historical data and contemporary, hypothesis-driven information–theoretic analyses. We also test a priori hypotheses about the effects of color morph (as suggested by early reports) and of drought (as suggested by recent climate predictions) on survival. Using robust mark–recapture models, we find some support for early hypotheses regarding the effect of color on survival, but we find no effect of drought. The congruence between early findings and our analyses highlights the usefulness of historical information in providing raw data for contemporary analyses and context for conservation and management decisions.

Quantifying the spatiotemporal dynamics in a chorus frog (Pseudacris) hybrid zone over 30 years

Although theory suggests that hybrid zones can move or change structure over time, studies supported by direct empirical evidence for these changes are relatively limited. We present a spatiotemporal genetic study of a hybrid zone between Pseudacris nigrita and P. fouquettei across the Pearl River between Louisiana and Mississippi. This hybrid zone was initially characterized in 1980 as a narrow and steep “tension zone,” in which hybrid populations were inferior to parentals and were maintained through a balance between selection and dispersal. We reanalyzed historical tissue samples and compared them to samples of recently collected individuals using microsatellites. Clinal analyses indicate that the cline has not shifted in roughly 30 years but has widened significantly. Anthropogenic and natural changes may have affected selective pressure or dispersal, and our results suggest that the zone may no longer best be described as a tension zone. To the best of our knowledge, this study provides the first evidence of significant widening of a hybrid cline but stasis of its center. Continued empirical study of dynamic hybrid zones will provide insight into the forces shaping their structure and the evolutionary potential they possess for the elimination or generation of species.

Integrative Taxonomic Approach for Describing a New Cryptic Species of Bush Frog (Raorchestes: Anura: Rhacophoridae) from the Western Ghats, India

A new cryptic species of bush frog Raorchestes honnametti sp. nov. is described from the south-eastern part of the Western Ghats, India. This newly described species belongs to the Charius clade and is morphologically similar to other clade members--R. charius and R. griet. Therefore, an integrative taxonomic approach based on molecular and bioacoustic analysis along with morphology was used to delimit the new species. Raorchestes honnametti sp. nov., is currently known only from Biligiri Rangaswamy Temple Tiger Reserve, a part of Biligiri Rangaswamy horst mountain range (a mountain formed due movement of two faults) formed during the Late Quaternary period (1.8-2.58 Ma). Discovery of cryptic species from a highly speciose and well-studied genus Raorchestes hints at the possible existence of several more cryptic species in this genus. We discuss the possible reasons for crypsis and emphasize the need for continued systematic surveys of amphibians across the Western Ghats.

Systematics of the Rhinellamargaritifera complex (Anura, Bufonidae) from western Ecuador and Panama with insights in the biogeography of Rhinellaalata

The Rhinellamargaritifera species group consists of 17 species of toads distributed in tropical and subtropical South America and eastern Central America. The identity of some of its species is poorly understood and there are numerous undescribed cryptic species. Among them, the status of Rhinellamargaritifera is one of the most problematic. Its range includes lowland rainforests separated by the Andes, the Chocoan rainforest to the west and the Amazonian rainforest to the east. This distribution is puzzling because the Andes are an old and formidable barrier to gene flow and therefore should generate vicariant speciation between disjunct lowland populations. Herein we clarify the taxonomy of populations of the Rhinellamargaritifera complex from Central America and the Chocó region of South America. The morphological and genetic variation of Rhinellamargaritifera was examined from 39 populations from Chocó, 24 from the upper Amazon region of Ecuador, and 37 from Panama, including the holotype of the Panamanian Rhinellaalata. Phylogenetic analyses were performed based on mitochondrial genes 12S rRNA, 16S rRNA, and cytochrome c oxidase I (COI) and the nuclear gene Tyrosinase (Tyr). The genetic and morphological data show that Panamanian and Chocoan populations are conspecific. In the phylogeny, populations from Chocó and Panama form a well-supported clade. The morphology of the holotype of Rhinellaalata falls within the variation range of Panamanian and Chocoan populations. Based on all this evidence, we assign the populations from western Ecuador and Panama to Rhinellaalata and demonstrate that the unusual distribution pattern of "Rhinellamargaritifera" on both sides of the Andes was an artifact of incorrectly defined species boundaries.

Speciation within Columnea section Angustiflora (Gesneriaceae): islands, pollinators and climate

Despite many advances in evolutionary biology, understanding the proximate mechanisms that lead to speciation for many taxonomic groups remains elusive. Phylogenetic analyses provide a means to generate well-supported estimates of species relationships. Understanding how genetic isolation (restricted gene flow) occurred in the past requires not only a well-supported molecular phylogenetic analysis, but also an understanding of when character states that define species may have changed. In this study, phylogenetic trees resolve species level relationships for fourteen of the fifteen species within Columnea section Angustiflorae (Gesneriaceae). The distributions of sister species pairs are compared and ancestral character states are reconstructed using Bayesian stochastic mapping. Climate variables were also assessed and shifts in ancestral climate conditions were mapped using SEEVA. The relationships between morphological character states and climate variables were assessed with correlation analyses. These results indicate that species in section Angustiflorae have likely diverged as a result of allopatric, parapatric, and sympatric speciation, with both biotic and abiotic forces driving morphological and phenological divergence.

High levels of diversity uncovered in a widespread nominal taxon: continental phylogeography of the neotropical tree frog Dendropsophus minutus

Species distributed across vast continental areas and across major biomes provide unique model systems for studies of biotic diversification, yet also constitute daunting financial, logistic and political challenges for data collection across such regions. The tree frog Dendropsophus minutus (Anura: Hylidae) is a nominal species, continentally distributed in South America, that may represent a complex of multiple species, each with a more limited distribution. To understand the spatial pattern of molecular diversity throughout the range of this species complex, we obtained DNA sequence data from two mitochondrial genes, cytochrome oxidase I (COI) and the 16S rhibosomal gene (16S) for 407 samples of D. minutus and closely related species distributed across eleven countries, effectively comprising the entire range of the group. We performed phylogenetic and spatially explicit phylogeographic analyses to assess the genetic structure of lineages and infer ancestral areas. We found 43 statistically supported, deep mitochondrial lineages, several of which may represent currently unrecognized distinct species. One major clade, containing 25 divergent lineages, includes samples from the type locality of D. minutus. We defined that clade as the D. minutus complex. The remaining lineages together with the D. minutus complex constitute the D. minutus species group. Historical analyses support an Amazonian origin for the D. minutus species group with a subsequent dispersal to eastern Brazil where the D. minutus complex originated. According to our dataset, a total of eight mtDNA lineages have ranges >100,000 km². One of them occupies an area of almost one million km² encompassing multiple biomes. Our results, at a spatial scale and resolution unprecedented for a Neotropical vertebrate, confirm that widespread amphibian species occur in lowland South America, yet at the same time a large proportion of cryptic diversity still remains to be discovered.

The phylogeography of Eurasian Fraxinus species reveals ancient transcontinental reticulation

To investigate the biogeographical history of ashes species of the Eurasian section Fraxinus and to test the hypothesis of ancient reticulations, we sequenced nuclear DNA (nETS and nITS, 1075 bp) for 533 samples and scored AFLP for 63 samples of Eurasian ashes within the section Fraxinus. The nITS phylogeny retrieved the classical view of the evolution of the section, whereas nETS phylogeny indicated an unexpected separation of F. angustifolia in two paraphyletic groups, respectively found in southeastern Europe and in the other parts of the Mediterranean basin. In the nETS phylogeny, the former group was closely related to F. excelsior, whereas the later was closely related to F. mandshurica, a species which is restricted nowadays to northeastern Asia. This topological incongruence between the two loci indicated the occurrence of an ancient reticulation between European and Asian ash species. Several other ancient reticulation events between the two European species and the other species of the section were supported by the posterior predictive checking method. Some of these reticulation events would have occurred during the Miocene, when climatic variations may have lead these species to expand their distribution range and come into contact. The recurrent reticulations observed among Eurasian ash species indicate that they should be considered as conspecific taxa, with subspecific status for some groups. Altogether, the results of the present study provide a rare documented evidence for the occurrence of multiple ancient reticulations within a group of temperate tree taxa with modern disjunct distributions in Eurasia. These ancient reticulation events indicate that the speciation process is slow in ashes, necessitating long periods of geographical isolation. The implications for speciation processes in temperate trees with similar life history and reproductive biology are discussed.

De novo transcriptome analysis of Liriodendron chinense petals and leaves by Illumina sequencing

Liriodendron chinense (Hemsl.) Sarg is an endangered species and occupies a pivotal position in phylogenetic studies of flowering plants, while its genomic resources are limited. In this study, we performed transcriptome sequencing for L. chinense petals and leaves using the Illumina paired-end sequencing technique. Approximately 17.02-Gb clean reads were obtained, and de novo assembly generated 87,841 unigenes, with an average length of 778 bp. Of these, there were 65,535 (74.61%) unigenes with significant similarity to publically available plant protein sequences. There were 3386 genes identified as significant differentially expressed between petals and leaves, among them 2969 (87.68%) were up-regulated and 417 (12.31%) down-regulated in petals. Metabolic pathway analysis revealed that 25 unigenes were predicted to be responsible for the biosynthesis of carotenoids, with 7 genes differentially expressed between these two tissues. This report is the first to identify genes associated with carotenoid biosynthesis in Liriodendron and represents a valuable resource for future genomic studies on the endangered species L. chinense.

Phylogeographic analyses strongly suggest cryptic speciation in the giant spiny frog (Dicroglossidae: Paa spinosa) and interspecies hybridization in Paa

Species identification is one of the most basic yet crucial issues in biology with potentially far-reaching implications for fields such as conservation, population ecology, and epidemiology. The widely distributed but threatened frog Paa spinosa has been speculated to represent a complex of multiple species. In this study, 254 individuals representing species of the genus Paa were investigated along the entire range of P. spinosa: 196 specimens of P. spinosa, 8 specimens of P. jiulongensis, 5 specimens of P. boulengeri, 20 specimens of P. exilispinosa, and 25 specimens of P. shini. Approximately 1333 bp of mtDNA sequence data (genes 12S rRNA and 16S rRNA) were used. Phylogenetic analyses were conducted using maximum parsimony, maximum likelihood and Bayesian inference. BEAST was used to estimate divergence dates of major clades. Results suggest that P. spinosa can be divided into three distinct major lineages. Each major lineage totally corresponds to geographical regions, revealing the presence of three candidate cryptic species. Isolation and differentiation among lineages are further supported by the great genetic distances between the lineages. The bifurcating phylogenetic pattern also suggests an east-west dispersal trend during historic cryptic speciation. Dating analysis estimates that P. spinosa from Western China split from the remaining P. spinosa populations in the Miocene and that P. spinosa from Eastern China diverged from Central China in the Pliocene. We also found that P. exilispinosa from Mainland China and Hong Kong might have a complex of multiple species. After identifying cryptic lineages, we then determine the discrepancy between the mtDNA and the morphotypes in several individuals. This discrepancy may have been caused by introgressive hybridization between P. spinosa and P. shini.

Integrated fossil and molecular data reveal the biogeographic diversification of the eastern Asian-eastern North American disjunct hickory genus (Carya Nutt.)

The hickory genus (Carya) contains ca. 17 species distributed in subtropical and tropical regions of eastern Asia and subtropical to temperate regions of eastern North America. Previously, the phylogenetic relationships between eastern Asian and eastern North American species of Carya were not fully confirmed even with an extensive sampling, biogeographic and diversification patterns had thus never been investigated in a phylogenetic context. We sampled 17 species of Carya and 15 species representing all other genera of the Juglandaceae as outgroups, with eight nuclear and plastid loci to reconstruct the phylogeny of Carya. The phylogenetic positions of seven extinct genera of the Juglandaceae were inferred using morphological characters and the molecular phylogeny as a backbone constraint. Divergence times within Carya were estimated with relaxed Bayesian dating. Biogeographic analyses were performed in DIVA and LAGRANGE. Diversification rates were inferred by LASER and APE packages. Our results support two major clades within Carya, corresponding to the lineages of eastern Asia and eastern North America. The split between the two disjunct clades is estimated to be 21.58 (95% HPD 11.07-35.51) Ma. Genus-level DIVA and LAGRANGE analyses incorporating both extant and extinct genera of the Juglandaceae suggested that Carya originated in North America, and migrated to Eurasia during the early Tertiary via the North Atlantic land bridge. Fragmentation of the distribution caused by global cooling in the late Tertiary resulted in the current disjunction. The diversification rate of hickories in eastern North America appeared to be higher than that in eastern Asia, which is ascribed to greater ecological opportunities, key morphological innovations, and polyploidy.

Marine incursion into East Asia: a forgotten driving force of biodiversity

Episodic marine incursion has been a major driving force in the formation of present-day diversity. Marine incursion is considered to be one of the most productive 'species pumps' particularly because of its division and coalescence effects. Marine incursion events and their impacts on diversity are well documented from South America, North America and Africa; however, their history and impacts in continental East Asia largely remain unknown. Here, we propose a marine incursion scenario occurring in East Asia during the Miocene epoch, 10-17 Ma. Our molecular phylogenetic analysis of Platorchestia talitrids revealed that continental terrestrial populations (Platorchestia japonica) form a monophyletic group that is the sister group to the Northwest Pacific coastal species Platorchestia pacifica. The divergence time between the two species coincides with Middle Miocene high global sea levels. We suggest that the inland form arose as a consequence of a marine incursion event. This is the first solid case documenting the impact of marine incursion on extant biodiversity in continental East Asia. We believe that such incursion event has had major impacts on other organisms and has played an important role in the formation of biodiversity patterns in the region.

Ancient tepui summits harbor young rather than old lineages of endemic frogs

The flattop mountains (tepuis) of South America are ancient remnants of the Precambrian Guiana Shield plateau. The tepui summits, isolated by their surrounding cliffs that can be up to 1000 m tall, are thought of as "islands in the sky," harboring relict flora and fauna that underwent vicariant speciation due to plateau fragmentation. High endemicity atop tepui summits support the idea of an ancient "Lost World" biota. However, recent work suggests that dispersal between lowlands and summits has occurred long after tepui formation indicating that tepui summits may not be as isolated from the lowlands as researchers have long suggested. Neither view of the origin of the tepui biota (i.e., ancient vicariance vs. recent dispersal) has strong empirical support owing to a lack of studies. We test diversification hypotheses of the Guiana Shield highlands by estimating divergence times of an endemic group of treefrogs, Tepuihyla. We find that diversification of this group does not support an ancient origin for this taxon; instead, divergence times among the highland species are 2-5 Ma. Our data indicate that most highland speciation occurred during the Pliocene. Thus, this unparalleled landscape known as "The Lost World" is inhabited, in part, not by Early Tertiary relicts but neoendemics.

Machine learning identifies specific habitats associated with genetic connectivity in Hyla squirella

The goal of this study was to identify and differentiate the influence of multiple habitat types that span a spectrum of suitability for Hyla squirella, a widespread frog species that occurs in a broad range of habitat types. We collected microsatellite data from 675 samples representing 20 localities from the southeastern USA and used machine-learning methodologies to identify significant habitat features associated with genetic structure. In simulation, we confirm that our machine-learning algorithm can successfully identify landscape features responsible for generating between-population genetic differentiation, suggesting that it can be a useful hypothesis-generating tool for landscape genetics. In our study system, we found that H. squirella were spatially structured and models including specific habitat types (i.e. upland oak forest and urbanization) consistently explained more variation in genetic distance (median pR(2)  = 47.78) than spatial distance alone (median pR(2)  = 23.81). Moreover, we estimate the relative importance that spatial distance, upland oak and urbanized habitat have in explaining genetic structure of H. squirella. We discuss how these habitat types may mechanistically facilitate dispersal in H. squirella. This study provides empirical support for the hypothesis that habitat-use can be an informative correlate of genetic differentiation, even for species that occur in a wide range of habitats.

Different processes lead to similar patterns: a test of codivergence and the role of sea level and climate changes in shaping a southern temperate freshwater assemblage

Background

Understanding how freshwater assemblages have been formed and maintained is a fundamental goal in evolutionary and ecological disciplines. Here we use a historical approach to test the hypothesis of codivergence in three clades of the Chilean freshwater species assemblage. Molecular studies of freshwater crabs (Aegla: Aeglidae: Anomura) and catfish (Trichomycterus arealatus: Trichomycteridae: Teleostei) exhibited similar levels of genetic divergences of mitochondrial lineages between species of crabs and phylogroups of the catfish, suggesting a shared evolutionary history among the three clades in this species assemblage.

Results

A phylogeny was constructed for Trichomycterus areolatus under the following best-fit molecular models of evolution GTR + I + R, HKY + I, and HKY for cytochrome b, growth hormone, and rag 1 respectively. A GTR + I + R model provided the best fit for both 28S and mitochondrial loci and was used to construct both Aegla phylogenies. Three different diversification models were observed and the three groups arose during different time periods, from 2.25 to 5.05 million years ago (Ma). Cladogenesis within Trichomycterus areolatus was initiated roughly 2.25 Ma (Late Pliocene - Early Pleistocene) some 1.7 - 2.8 million years after the basal divergences observed in both Aegla clades. These results reject the hypothesis of codivergence.

Conclusions

The similar genetic distances between terminal sister-lineages observed in these select taxa from the freshwater Chilean species assemblage were formed by different processes occurring over the last ~5.0 Ma. Dramatic changes in historic sea levels documented in the region appear to have independently shaped the evolutionary history of each group. Our study illustrates the important role that history plays in shaping a species assemblage and argues against assuming similar patterns equal a shared evolutionary history.

Phylogeny, life history, and ecology contribute to differences in amphibian susceptibility to ranaviruses

Research that identifies the potential host range of generalist pathogens as well as variation in host susceptibility is critical for understanding and predicting the dynamics of infectious diseases within ecological communities. Ranaviruses have been linked to amphibian die-off events worldwide with the greatest number of reported mortality events occurring in the United States. While reports of ranavirus-associated mortality events continue to accumulate, few data exist comparing the relative susceptibility of different species. Using a series of laboratory exposure experiments and comparative phylogenetics, we compared the susceptibilities of 19 amphibian species from two salamander families and five anurans families for two ranavirus isolates: frog virus 3 (FV3) and an FV3-like isolate from an American bullfrog culture facility. We discovered that ranaviruses were capable of infecting 17 of the 19 larval amphibian species tested with mortality ranging from 0 to 100%. Phylogenetic comparative methods demonstrated that species within the anuran family Ranidae were generally more susceptible to ranavirus infection compared to species from the other five families. We also found that susceptibility to infection was associated with species that breed in semi-permanent ponds, develop rapidly as larvae, and have limited range sizes. Collectively, these results suggest that phylogeny, life history characteristics, and habitat associations of amphibians have the potential to impact susceptibility to ranaviruses.

Post-glacial expansion into the Paleozoic Plateau: evidence of an Ozarkian refugium for the Ozark minnow Notropis nubilus (Teleostei: Cypriniformes)

Genetic variation was examined within the Ozark minnow Notropis nubilus using complete mtDNA cytochrome b gene sequences from 160 individuals representing 30 localities to test hypotheses on the origin of the distribution. Phylogenetic analyses revealed three strongly supported clades of haplotypes consistent with geographic distributions: a clade from the Western Ozarks, a clade from the Southern Ozarks and a clade from the Northern Ozarks and upper Mississippi River basin. The estimated mean ages of these clades indicated that they diverged during pre-Illinoian glacial cycles extending from the late Pliocene into the early Pleistocene. Results of demographic analyses based on coalescent approaches supported the hypothesis that the Paleozoic Plateau was not a refugium for N. nubilus during periodic glacial advances. There is evidence of a genetic signature of northern expansion into the Paleozoic Plateau from a Southern Ozarkian refugium. Populations expanded out of drainages in the Northern Ozarks into the Paleozoic Plateau during the late Pleistocene. Subsequently, the two regions were isolated due to the recent extirpation of intervening populations caused by the loss of suitable habitat.

A likelihood framework for estimating phylogeographic history on a continuous landscape

Due to lack of an adequate statistical framework, biologists studying phylogeography are abandoning traditional methods of estimating phylogeographic history in favor of statistical methods designed to test a priori hypotheses. These new methods may, however, have limited descriptive utility. Here, we develop a new statistical framework that can be used to both test a priori hypotheses and estimate phylogeographic history of a gene (and the statistical confidence in that history) in the absence of such hypotheses. The statistical approach concentrates on estimation of geographic locations of the ancestors of a set of sampled organisms. Now we use (2) to derive the likelihood of the ancestral geographic coordinates and the value of the scaled dispersal parameter, given the observed geographic coordinates (assuming known topology and branch lengths). Using a maximum likelihood approach, which is implemented in the new program PhyloMapper, we apply this statistical framework to a 246-taxon mitochondrial genealogy of North American chorus frogs, focusing in detail on one of these species. We demonstrate three lines of evidence for recent northward expansion of the mitochondrion of the coastal clade of Pseudacris feriarum: higher per-generation dispersal distance in the recently colonized region, a noncentral ancestral location, and directional migration. After illustrating one method of accommodating phylogenetic uncertainty, we conclude by discussing how extensions of this framework could function to incorporate a priori ecological and geological information into phylogeographic analyses.