New evidence for parallel evolution of colour patterns in Malagasy poison frogs (Mantella)

Cytogenetics of the Hybridogenetic Frog Pelophylax grafi and Its Parental Species Pelophylax perezi

Hybrid taxa from the genus Pelophylax can propagate themselves in a modified way of sexual reproduction called hybridogenesis ensuring the formation of clonal gametes containing the genome of only one parental (host) species. Pelophylax grafi from South-Western Europe is a hybrid composed of P. ridibundus and P. perezi genomes and it lives with a host species P. perezi (P-G system). Yet it is unknown, whether non-Mendelian inheritance is fully maintained in such populations. In this study, we characterize P. perezi and P. grafi somatic karyotypes by using comparative genomic hybridization, genomic in situ hybridization, fluorescent in situ hybridization, and actinomycin D-DAPI. Here, we show the homeology of P. perezi and P. grafi somatic karyotypes to other Pelophylax taxa with 2n = 26 and equal contribution of ridibundus and perezi chromosomes in P. grafi which supports F1 hybrid genome constitution as well as a hemiclonal genome inheritance. We show that ridibundus chromosomes have larger regions of interstitial (TTAGGG)n repeats flanking the nucleolus organizing region on chromosome no. 10 and a high quantity of AT pairs in the centromeric regions. In P. perezi, we found species-specific sequences in metaphase chromosomes and marker structures in lampbrush chromosomes. Pericentromeric RrS1 repeat sequence was present in perezi and ridibundus chromosomes, but the blocks were stronger in ridibundus. Various cytogenetic techniques applied to the P-G system provide genome discrimination between ridibundus and perezi chromosomal sets. They could be used in studies of germ-line cells to explain patterns of clonal gametogenesis in P. grafi and broaden the knowledge about reproductive strategies in hybrid animals.

An integrative taxonomic revision of slug-eating snakes (Squamata: Pareidae: Pareineae) reveals unprecedented diversity in Indochina

Slug-eating snakes of the subfamily Pareinae are an insufficiently studied group of snakes specialized in feeding on terrestrial mollusks. Currently Pareinae encompass three genera with 34 species distributed across the Oriental biogeographic region. Despite the recent significant progress in understanding of Pareinae diversity, the subfamily remains taxonomically challenging. Here we present an updated phylogeny of the subfamily with a comprehensive taxon sampling including 30 currently recognized Pareinae species and several previously unknown candidate species and lineages. Phylogenetic analyses of mtDNA and nuDNA data supported the monophyly of the three genera Asthenodipsas, Aplopeltura, and Pareas. Within both Asthenodipsas and Pareas our analyses recovered deep differentiation with each genus being represented by two morphologically diagnosable clades, which we treat as subgenera. We further apply an integrative taxonomic approach, including analyses of molecular and morphological data, along with examination of available type materials, to address the longstanding taxonomic questions of the subgenus Pareas, and reveal the high level of hidden diversity of these snakes in Indochina. We restrict the distribution of P. carinatus to southern Southeast Asia, and recognize two subspecies within it, including one new subspecies proposed for the populations from Thailand and Myanmar. We further revalidate P. berdmorei, synonymize P. menglaensis with P. berdmorei, and recognize three subspecies within this taxon, including the new subspecies erected for the populations from Laos and Vietnam. Furthermore, we describe two new species of Pareas from Vietnam: one belonging to the P. carinatus group from southern Vietnam, and a new member of the P. nuchalis group from the central Vietnam. We provide new data on P. temporalis, and report on a significant range extension for P. nuchalis. Our phylogeny, along with molecular clock and ancestral area analyses, reveal a complex diversification pattern of Pareinae involving a high degree of sympatry of widespread and endemic species. Our analyses support the "upstream" colonization hypothesis and, thus, the Pareinae appears to have originated in Sundaland during the middle Eocene and then colonized mainland Asia in early Oligocene. Sundaland and Eastern Indochina appear to have played the key roles as the centers of Pareinae diversification. Our results reveal that both vicariance and dispersal are responsible for current distribution patterns of Pareinae, with tectonic movements, orogeny and paleoclimatic shifts being the probable drivers of diversification. Our study brings the total number of Pareidae species to 41 and further highlights the importance of comprehensive taxonomic revisions not only for the better understanding of biodiversity and its evolution, but also for the elaboration of adequate conservation actions.

The evolutionary game of interspecific mutualism in the multi-species model

Mutualistic interspecific interactions, including Müllerian mimicry and division of labor, are common in nature. In contrast to antagonistic interactions, where faster evolution is favored, mutualism can favor slower evolution under certain conditions. This is called the Red King effect. Since Bergstrom and Lachmann (2003) proposed the Red King effect, it has been investigated only in two-species models. However, biological examples suggest that mutualism can include three or more species. Here, I modeled the evolutionary dynamics of mutualism in communities where involving two or more species, and in which all species mutually interact. Regardless of the number of species in the community, it is possible to derive conditions for stable equilibria. Although nonlinear relationships exist between the evolutionary rates and the evolutionary fate of each species in the multi-species communities, the model suggests that it is possible to predict whether faster evolution is favored or disfavored for the relatively rapidly evolving species; however, it is difficult to predict the evolutionary fate of species that evolve relatively slowly because their evolutionary dynamics are affected by the evolutionary fate of species evolving rapidly.

Mitochondrial Introgression, Color Pattern Variation, and Severe Demographic Bottlenecks in Three Species of Malagasy Poison Frogs, Genus Mantella

Madagascar is a biodiversity hotspot particularly rich in amphibian diversity and only a few charismatic Malagasy amphibians have been investigated for their population-level differentiation. The Mantellamadagascariensis group is composed of two rainforest and three swamp forest species of poison frogs. We first confirm the monophyly of this clade using DNA sequences of three nuclear and four mitochondrial genes, and subsequently investigate the population genetic differentiation and demography of the swamp forest species using one mitochondrial, two nuclear and a set of nine microsatellite markers. Our results confirm the occurrence of two main mitochondrial lineages, one dominated by Mantellaaurantiaca (a grouping supported also by our microsatellite-based tree) and the other by Mantellacrocea + Mantellamilotympanum. These two main lineages probably reflect an older divergence in swamp Mantella. Widespread mitochondrial introgression suggests a fairly common occurrence of inter-lineage gene flow. However, nuclear admixture seems to play only a limited role in this group, and the analyses of the RAG-1 marker points to a predominant incomplete lineage sorting scenario between all five species of the group, which probably diverged relatively recently. Our demographic analyses show a common, severe and recent demographic contraction, inferred to be in temporal coincidence with the massive deforestation events that took place in the past 1000 years. Current data do not allow to conclusively delimit independent evolutionary units in these frogs, and we therefore refrain to suggest any taxonomic changes.

Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard

Reptiles use pterin and carotenoid pigments to produce yellow, orange, and red colors. These conspicuous colors serve a diversity of signaling functions, but their molecular basis remains unresolved. Here, we show that the genomes of sympatric color morphs of the European common wall lizard (Podarcis muralis), which differ in orange and yellow pigmentation and in their ecology and behavior, are virtually undifferentiated. Genetic differences are restricted to two small regulatory regions near genes associated with pterin [sepiapterin reductase (SPR)] and carotenoid [beta-carotene oxygenase 2 (BCO2)] metabolism, demonstrating that a core gene in the housekeeping pathway of pterin biosynthesis has been coopted for bright coloration in reptiles and indicating that these loci exert pleiotropic effects on other aspects of physiology. Pigmentation differences are explained by extremely divergent alleles, and haplotype analysis revealed abundant transspecific allele sharing with other lacertids exhibiting color polymorphisms. The evolution of these conspicuous color ornaments is the result of ancient genetic variation and cross-species hybridization.

Phenotypic and genetic diversity in aposematic Malagasy poison frogs (genus Mantella)

Intraspecific color variation has long fascinated evolutionary biologists. In species with bright warning coloration, phenotypic diversity is particularly compelling because many factors, including natural and sexual selection, contribute to intraspecific variation. To better understand the causes of dramatic phenotypic variation in Malagasy poison frogs, we quantified genetic structure and color and pattern variation across three closely related species, Mantella aurantiaca, Mantella crocea, and Mantella milotympanum. Although our restriction site-associated DNA (RAD) sequencing approach identified clear genetic clusters, they do not align with current species designations, which has important conservation implications for these imperiled frogs. Moreover, our results suggest that levels of intraspecific color variation within this group have been overestimated, while species diversity has been underestimated. Within major genetic clusters, we observed distinct patterns of variation including: populations that are phenotypically similar yet genetically distinct, populations where phenotypic and genetic breaks coincide, and populations that are genetically similar but have high levels of within-population phenotypic variation. We also detected admixture between two of the major genetic clusters. Our study suggests that several mechanisms-including hybridization, selection, and drift-are contributing to phenotypic diversity. Ultimately, our work underscores the need for a reevaluation of how polymorphic and polytypic populations and species are classified, especially in aposematic organisms.

Diversity in warning coloration: selective paradox or the norm?

Aposematic theory has historically predicted that predators should select for warning signals to converge on a single form, as a result of frequency‐dependent learning. However, widespread variation in warning signals is observed across closely related species, populations and, most problematically for evolutionary biologists, among individuals in the same population. Recent research has yielded an increased awareness of this diversity, challenging the paradigm of signal monomorphy in aposematic animals. Here we provide a comprehensive synthesis of these disparate lines of investigation, identifying within them three broad classes of explanation for variation in aposematic warning signals: genetic mechanisms, differences among predators and predator behaviour, and alternative selection pressures upon the signal. The mechanisms producing warning coloration are also important. Detailed studies of the genetic basis of warning signals in some species, most notably Heliconius butterflies, are beginning to shed light on the genetic architecture facilitating or limiting key processes such as the evolution and maintenance of polymorphisms, hybridisation, and speciation. Work on predator behaviour is changing our perception of the predator community as a single homogenous selective agent, emphasising the dynamic nature of predator–prey interactions. Predator variability in a range of factors (e.g. perceptual abilities, tolerance to chemical defences, and individual motivation), suggests that the role of predators is more complicated than previously appreciated. With complex selection regimes at work, polytypisms and polymorphisms may even occur in Müllerian mimicry systems. Meanwhile, phenotypes are often multifunctional, and thus subject to additional biotic and abiotic selection pressures. Some of these selective pressures, primarily sexual selection and thermoregulation, have received considerable attention, while others, such as disease risk and parental effects, offer promising avenues to explore. As well as reviewing the existing evidence from both empirical studies and theoretical modelling, we highlight hypotheses that could benefit from further investigation in aposematic species. Finally by collating known instances of variation in warning signals, we provide a valuable resource for understanding the taxonomic spread of diversity in aposematic signalling and with which to direct future research. A greater appreciation of the extent of variation in aposematic species, and of the selective pressures and constraints which contribute to this once‐paradoxical phenomenon, yields a new perspective for the field of aposematic signalling.

Evolutionary history of the river frog genus Amietia (Anura: Pyxicephalidae) reveals extensive diversification in Central African highlands

The African river frog genus Amietia is found near rivers and other lentic water sources throughout central, eastern, and southern Africa. Because the genus includes multiple morphologically conservative species, taxonomic studies of river frogs have been relatively limited. We sampled 79 individuals of Amietia from multiple localities in and near the Albertine Rift (AR) of Burundi, Democratic Republic of the Congo, and Uganda. We utilized single-gene (16S) and concatenated (12S, 16S, cyt b and RAG1) gene-tree analyses and coalescent species-tree analyses to construct phylogenetic trees. Two divergence dating approaches were used in BEAST, including secondary calibration points with 12S, 16S, cyt b and RAG1, and a molecular clock with the 12S, 16S, and cyt b genes. All analyses recovered Amietia as monophyletic with strong support, and revealed several well-supported cryptic lineages, which is consistent with other recent phylogeography studies of AR amphibians. Dating estimates were similar, and Amietia diversification is coincident with global cooling and aridification events in the Miocene and Pliocene, respectively. Our results suggest additional taxonomic work is needed to describe multiple new species of AR Amietia, some of which have limited geographic distributions that are likely to be of conservation concern.

Ecological and evolutionary influences on body size and shape in the Galápagos marine iguana (Amblyrhynchus cristatus)

Oceanic islands are often inhabited by endemic species that have undergone substantial morphological evolutionary change due to processes of multiple colonizations from various source populations, dispersal, and local adaptation. Galápagos marine iguanas are an example of an island endemic exhibiting high morphological diversity, including substantial body size variation among populations and sexes, but the causes and magnitude of this variation are not well understood. We obtained morphological measurements from marine iguanas throughout their distribution range. These data were combined with genetic and local environmental data from each population to investigate the effects of evolutionary history and environmental conditions on body size and shape variation and sexual dimorphism. Our results indicate that body size and shape are highly variable among populations. Sea surface temperature and island perimeter, but not evolutionary history as depicted by phylogeographic patterns in this species, explain variation in body size among populations. Conversely, evolutionary history, but not environmental parameters or island size, was found to influence variation in body shape among populations. Finally, in all populations except one, we found strong sexual dimorphism in body size and shape in which males are larger, with higher heads than females, while females have longer heads than males. Differences among populations suggest that plasticity and/or genetic adaptation may shape body size and shape variation in marine iguanas. This study will help target future investigations to address the contribution of plasticity versus genetic adaptation on size and shape variation in marine iguanas.

Individual and Geographic Variation of Skin Alkaloids in Three Swamp-Forest Species of Madagascan Poison Frogs (Mantella)

Seventy skins of three mantellid frog species from Madagascan swamp-forest habitats, Mantella aurantiaca, M. crocea, and M. milotympanum, were individually examined for skin alkaloids using GC/MS. These poison frogs were found to differ significantly in their alkaloid composition from species of Mantella originating from non-flooded rainforest in eastern Madagascar, which were examined in earlier work. Only 16 of the previously detected 106 alkaloids were represented among the 60 alkaloids from the swamp-forest frogs of the present study. We hypothesize this difference is related mainly to habitat but cannot exclude a phylogenetic component as the three swamp-forest species are a closely related monophyletic group. The paucity of alkaloids with unbranched-carbon skeletons (ant-derived) and the commonness of alkaloids with branched-carbon skeletons (mite-derived) indicate that oribatid mites are a major source of alkaloids in these species of mantellids. Furthermore, most of the alkaloids have an oxygen atom in their formulae. Differences in alkaloids were observed among species, populations of the same species, and habitats. In M. aurantiaca, small geographic distances among populations were associated with differences in alkaloid profiles, with a remote third site illustrating even greater differences. The present study and an earlier study of three other mantellid species suggest that oribatid mites, and not ants, are the major source of alkaloids in the species of mantellids examined thus far.

Consuming viscous prey: a novel protein-secreting delivery system in neotropical snail-eating snakes

BACKGROUND

Efficient venom delivery systems are known to occur only in varanoid lizards and advanced colubroidean snakes among squamate reptiles. Although components of these venomous systems might have been present in a common ancestor, the two lineages independently evolved strikingly different venom gland systems. In snakes, venom is produced exclusively by serous glands in the upper jaw. Within the colubroidean radiation, lower jaw seromucous infralabial glands are known only in two distinct lineages-the basal pareatids and the more advanced Neotropical dipsadines known as "goo-eating snakes". Goo-eaters are a highly diversified, ecologically specialized clade that feeds exclusively on invertebrates (e.g., gastropod molluscs and annelids). Their evolutionary success has been attributed to their peculiar feeding strategies, which remain surprisingly poorly understood. More specifically, it has long been thought that the more derived Dipsadini genera Dipsas and Sibynomorphus use glandular toxins secreted by their infralabial glands to extract snails from their shells.

RESULTS

Here, we report the presence in the tribe Dipsadini of a novel lower jaw protein-secreting delivery system effected by a gland that is not functionally related to adjacent teeth, but rather opens loosely on the oral epithelium near the tip of the mandible, suggesting that its secretion is not injected into the prey as a form of envenomation but rather helps control the mucus and assists in the ingestion of their highly viscous preys. A similar protein-secreting system is also present in the goo-eating genus Geophis and may share the same adaptive purpose as that hypothesized for Dipsadini. Our phylogenetic hypothesis suggests that the acquisition of a seromucous infralabial gland represents a uniquely derived trait of the goo-eating clade that evolved independently twice within the group as a functionally complex protein-secreting delivery system.

CONCLUSIONS

The acquisition by snail-eating snakes of such a complex protein-secreting system suggests that the secretion from the hypertrophied infralabial glands of goo-eating snakes may have a fundamental role in mucus control and prey transport rather than envenomation of prey. Evolution of a functional secretory system that combines a solution for mucus control and transport of viscous preys is here thought to underlie the successful radiation of goo-eating snakes.

Cryptic genetic diversity is paramount in small-bodied amphibians of the genus Euparkerella (Anura: Craugastoridae) endemic to the Brazilian Atlantic forest

Morphological similarity associated to restricted distributions and low dispersal abilities make the direct developing “Terrarana” frogs of the genus Euparkerella a good model for examining diversification processes. We here infer phylogenetic relationships within the genus Euparkerella, using DNA sequence data from one mitochondrial and four nuclear genes coupled with traditional Bayesian phylogenetic reconstruction approaches and more recent coalescent methods of species tree inference. We also used Bayesian clustering analysis and a recent Bayesian coalescent-based approach specifically to infer species delimitation. The analysis of 39 individuals from the four known Euparkerella species uncovered high levels of genetic diversity, especially within the two previously morphologically-defined E. cochranae and E. brasiliensis. Within these species, the gene trees at five independent loci and trees from combined data (concatenated dataset and the species tree) uncovered six deeply diverged and geographically coherent evolutionary units, which may have diverged between the Miocene and the Pleistocene. These six units were also uncovered in the Bayesian clustering analysis, and supported by the Bayesian coalescent-based species delimitation (BPP), and Genealogical Sorting Index (GSI), providing thus strong evidence for underestimation of the current levels of diversity within Euparkerella. The cryptic diversity now uncovered opens new opportunities to examine the origins and maintenance of microendemism in the context of spatial heterogeneity and/or human induced fragmentation of the highly threatened Brazilian Atlantic forest hotspot.

Ongoing speciation in the Tibetan plateau Gymnocypris species complex

Local adaptation towards divergent ecological conditions often results in genetic differentiation and adaptive phenotypic divergence. To illuminate the ecological distinctiveness of the schizothoracine fish, we studied a Gymnocypris species complex consisting of three morphs distributed across four bodies of water (the Yellow River, Lake Qinghai, the Ganzi River and Lake Keluke) in the Northeast Tibetan Plateau. We used a combination of mitochondrial (16S rRNA and Cyt b) and nuclear (RAG-2) genetic sequences to investigate the phylogeography of these morphs based on a sample of 277 specimens. Analysis of gill rakers allowed for mapping of phenotypic trajectories along the phylogeny. The phylogenetic and morphological analyses showed that the three sparsely rakered morphs were present at two extremes of the phylogenetic tree: the Yellow River morphs were located at the basal phylogenetic split, and the Lake Keluke and Ganzi River morphs at the peak, with the densely rakered Lake Qinghai morphs located between these two extremes. Age estimation further indicated that the sparsely rakered morphs constituted the oldest and youngest lineages, whereas the densely rakered morph was assigned to an intermediate-age lineage. These results are most compatible with the process of evolutionary convergence or reversal. Disruptive natural selection due to divergent habitats and dietary preferences is likely the driving force behind the formation of new morphs, and the similarities between their phenotypes may be attributable to the similarities between their forms of niche tracking associated with food acquisition. This study provides the first genetic evidence for the occurrence of convergence or reversal in the schizothoracine fish of the Tibetan Plateau at small temporal scales.

Postglacial colonization of the Qinling Mountains: phylogeography of the swelled vent frog (Feirana quadranus)

BACKGROUND

The influence of Pleistocene climatic fluctuations on intraspecific diversification in the Qinling-Daba Mountains of East Asia remains poorly investigated. We tested hypotheses concerning refugia during the last glacial maximum (LGM) in this region by examining the phylogeography of the swelled vent frog (Feirana quadranus; Dicroglossidae, Anura, Amphibia).

METHODOLOGY/PRINCIPAL FINDINGS

We obtained complete mitochondrial ND2 gene sequences of 224 individuals from 34 populations of Feirana quadranus for phylogeographic analyses. Additionally, we obtained nuclear tyrosinase gene sequences of 68 F. quadranus, one F. kangxianensis and three F. taihangnica samples to test for mitochondrial introgression among them. Phylogenetic analyses based on all genes revealed no introgression among them. Phylogenetic analyses based on ND2 datasets revealed that F. quadranus was comprised of six lineages which were separated by deep valleys; the sole exception is that the Main Qinling and Micang-Western Qinling lineages overlap in distribution. Analyses of population structure indicated restricted gene flow among lineages. Coalescent simulations and divergence dating indicated that the basal diversification within F. quadranus may be associated with the dramatic uplifts of the Tibetan Plateau during the Pliocene. Coalescent simulations indicated that Wuling, Daba, and Western Qinling-Micang-Longmen Mountains were refugia for F. quadranus during the LGM. Demographic analyses indicated that the Daba lineage experienced population size increase prior to the LGM but the Main Qinling and the Micang-Western Qinling lineages expanded in population size and range after the LGM, and the other lineages almost have stable population size or slight slow population size decline.

CONCLUSIONS/SIGNIFICANCE

The Qinling-Daba Mountains hosted three refugia for F. quadranus during the LGM. Populations that originated in the Daba Mountains colonized the Main Qinling Mountains after the LGM. Recent sharp expansion of the Micang-Western Qinling and Main Qinling lineages probably contribute to their present-day secondary contact.

Molecular phylogeny supports the paraphyletic nature of the genus Trogoderma (Coleoptera: Dermestidae) collected in the Australasian ecozone

To date, a molecular phylogenetic approach has not been used to investigate the evolutionary structure of Trogoderma and closely related genera. Using two mitochondrial genes, Cytochrome Oxidase I and Cytochrome B, and the nuclear gene, 18S, the reported polyphyletic positioning of Trogoderma was examined. Paraphyly in Trogoderma was observed, with one Australian Trogoderma species reconciled as sister to all Dermestidae and the Anthrenocerus genus deeply nested within the Australian Trogoderma clade. In addition, time to most recent common ancestor for a number of Dermestidae was calculated. Based on these estimations, the Dermestidae origin exceeded 175 million years, placing the origins of this family in Pangaea.

Evolutionary history of the genus Tarentola (Gekkota: Phyllodactylidae) from the Mediterranean Basin, estimated using multilocus sequence data

Background

The pronounced morphological conservatism within Tarentola geckos contrasted with a high genetic variation in North Africa, has led to the hypothesis that this group could represent a cryptic species complex, a challenging system to study especially when trying to define distinct evolutionary entities and address biogeographic hypotheses. In the present work we have re-examined the phylogenetic and phylogeographic relationships between and within all Mediterranean species of Tarentola, placing the genealogies obtained into a temporal framework. In order to do this, we have investigated the sequence variation of two mitochondrial (12S rRNA and 16S rRNA), and four nuclear markers (ACM4, PDC, MC1R, and RAG2) for 384 individuals of all known Mediterranean Tarentola species, so that their evolutionary history could be assessed.

Results

Of all three generated genealogies (combined mtDNA, combined nDNA, and mtDNA+nDNA) we prefer the phylogenetic relationships obtained when all genetic markers are combined. A total of 133 individuals, and 2,901 bp of sequence length, were used in this analysis. The phylogeny obtained for Tarentola presents deep branches, with T. annularis, T. ephippiata and T. chazaliae occupying a basal position and splitting from the remaining species around 15.38 Mya. Tarentola boehmei is sister to all other Mediterranean species, from which it split around 11.38 Mya. There are also two other major groups: 1) the T. mauritanica complex present in North Africa and Europe; and 2) the clade formed by the T. fascicularis/deserti complex, T. neglecta and T. mindiae, occurring only in North Africa. The cladogenesis between these two groups occurred around 8.69 Mya, coincident with the late Miocene. Contrary to what was initially proposed, T. neglecta and T. mindiae are sister taxa to both T. fascicularis and T. deserti.

Conclusions

At least in the Iberian Peninsula and Northwest Africa, the lineages obtained have some geographic coherency, whilst the evolutionary history of the forms from Northeast Africa remains unclear, with a paraphyletic T. fascicularis with respect to T. deserti. The separation between the T. mauritanica complex and the clade formed by the T. fascicularis/deserti complex, T. neglecta and T. mindiae is coincident with the uplift of the Atlas Mountain chain, and the establishment of two distinct bioclimatic regions on each side of the barrier.

Description of a new Malagasy treefrog (Boophis) occurring syntopically with its sister species, and a plea for studies on non-allopatric speciation in tropical amphibians

Based on concordant differences in male advertisement call, tadpole morphology, and absence of haplotype sharing in mitochondrial and nuclear DNA markers, we describe a new species of treefrog from Ranomafana National Park in the southern central east of Madagascar. In its adult stage Boophis narinsi sp. n. is highly similar to its sister species, Boophis majori, but appears to differ in having longer hindlimbs. The genetic divergences between these two species (2.5-3.3% in a fragment of the 16S rRNA gene, depending on fragment length and individual haplotype analyzed) are below the threshold typically characterizing distinct species of anurans. Together with their relatively small and largely overlapping ranges and their sympatric occurrence in Ranomafana National Park, this indicates that they potentially could have originated rather recently by adaptive speciation under parapatric or sympatric conditions. Most studies on amphibian speciation have so far by default assumed vicariant speciation. We suggest that alternative speciation scenarios should be considered in future works and characterize settings in which more reliable assessments of adaptive parapatric or sympatric speciation could be carried out.

Multiple lines of evidence for demographic and range expansion of a temperate species (Hyla sarda) during the last glaciation

Many temperate species experienced demographic and range contractions in response to climatic changes during Pleistocene glaciations. In this study, we investigate the evolutionary history of the Tyrrhenian tree frog Hyla sarda, a species inhabiting the Corsica-Sardinia island system (Western Mediterranean basin). We used sequence analysis of two mitochondrial (overall 1229 bp) and three nuclear (overall 1692 bp) gene fragments to assess the phylogeography and demographic history of this species, and species distribution modelling (SDM) to predict its range variation over time. Phylogeographic, historical demographic and SDM analyses consistently indicate that H. sarda does not conform to the scenario generally expected for temperate species but rather underwent demographic and range expansion mostly during the last glacial phase. Palaeogeographic data and SDM analyses suggest that such expansion was driven by the glaciation-induced increase in lowland areas during marine regression. This unusual scenario suggests that at least some temperate species may not have suffered the adverse effects of glacial climate on their population size and range extent, owing to the mitigating effects of other glaciations-induced palaeoenvironmental changes. We discuss previous clues for the occurrence of such a scenario in other species and some possible challenges with its identification. Early phylogeographic literature suggested that responses to the Pleistocene glacial-interglacial cycles were expected to vary among species and regions. Our results point out that such variation may have been greater than previously thought.

Spiny frogs (Paini) illuminate the history of the Himalayan region and Southeast Asia

Asian frogs of the tribe Paini (Anura: Dicroglossidae) range across several first-order tectono-morphological domains of the Cenozoic Indo-Asian collision that include the Tibetan Plateau, the Himalayas, and Indochina. We show how the tectonic events induced by the Indo-Asian collision affected the regional biota and, in turn, how the geological history of the earth can be viewed from a biological perspective. Our analysis of a concatenated dataset comprising four nuclear gene sequences of Paini revealed two main radiations, corresponding to the genera Nanorana (I) and Quasipaa (II). Five distinct clades are recognized: Tibetan plateau clade (I-1), Himalaya clade (I-2), environs of Himalaya-Tibetan plateau clade (I-3), South China clade (II-1), and Indochina clade (II-2). This pattern of relationships highlights the significance of geography in shaping evolutionary history. Building on our molecular dating, ancestral region reconstruction, and distributional patterns, we hypothesize a distinct geographic and climatic transition in Asia beginning in the Oligocene and intensifying in the Miocene; this stimulated rapid diversification of Paini. Vicariance explains species formation among major lineages within Nanorana. Dispersal, in contrast, plays an important role among Quasipaa, with the southern Chinese taxa originating from Indochina. Our results support the tectonic hypothesis that an uplift in the Himalaya-Tibetan plateau region resulting from crustal thickening and lateral extrusion of Indochina occurred synchronously during the transition between Oligocene and Miocene in reaction to the Indo-Asian collision. The phylogenetic history of Paini illuminates critical aspects of the timing of geological events responsible for the current geography of Southeast Asia.

Base composition, selection, and phylogenetic significance of indels in the recombination activating gene-1 in vertebrates

BACKGROUND

The Recombination Activating Proteins, RAG1 and RAG2, play a crucial role in the immune response in vertebrates. Among the nuclear markers currently used for phylogenetic purposes, Rag1 has especially enjoyed enormous popularity, since it successfully contributed to elucidating the relationships among and within a large variety of vertebrate lineages. We here report on a comparative investigation of the genetic variation, base composition, presence of indels, and selection in Rag1 in different vertebrate lineages (Actinopterygii, Amphibia, Aves, Chondrichthyes, Crocodylia, Lepidosauria, Mammalia, and Testudines) through the analysis of 582 sequences obtained from Genbank. We also analyze possible differences between distinct parts of the gene with different type of protein functions.

RESULTS

In the vertebrate lineages studied, Rag1 is over 3 kb long. We observed a high level of heterogeneity in base composition at the 3(rd )codon position in some of the studied vertebrate lineages and in some specific taxa. This result is also paralleled by taxonomic differences in the GC content at the same codon position. Moreover, positive selection occurs at some sites in Aves, Lepidosauria and Testudines. Indels, which are often used as phylogenetic characters, are more informative across vertebrates in the 5' than in the 3'-end of the gene. When the entire gene is considered, the use of indels as phylogenetic character only recovers one major vertebrate clade, the Actinopterygii. However, in numerous cases insertions or deletions are specific to a monophyletic group.

CONCLUSIONS

Rag1 is a phylogenetic marker of undoubted quality. Our study points to the need of carrying out a preliminary investigation on the base composition and the possible existence of sites under selection of this gene within the groups studied to avoid misleading resolution. The gene shows highly heterogeneous base composition, which affects some taxa in particular and contains sites under positive selection in some vertebrate lineages in the 5'-end. The first part of the gene (5'-end) is more variable than the second (3'-end), and less affected by a heterogeneous base composition. However, in some vertebrate lineages the 5'-end of the gene is not yet widely used for phylogenetic studies.

Genome 10K: a proposal to obtain whole-genome sequence for 10,000 vertebrate species

The human genome project has been recently complemented by whole-genome assessment sequence of 32 mammals and 24 nonmammalian vertebrate species suitable for comparative genomic analyses. Here we anticipate a precipitous drop in costs and increase in sequencing efficiency, with concomitant development of improved annotation technology and, therefore, propose to create a collection of tissue and DNA specimens for 10,000 vertebrate species specifically designated for whole-genome sequencing in the very near future. For this purpose, we, the Genome 10K Community of Scientists (G10KCOS), will assemble and allocate a biospecimen collection of some 16,203 representative vertebrate species spanning evolutionary diversity across living mammals, birds, nonavian reptiles, amphibians, and fishes (ca. 60,000 living species). In this proposal, we present precise counts for these 16,203 individual species with specimens presently tagged and stipulated for DNA sequencing by the G10KCOS. DNA sequencing has ushered in a new era of investigation in the biological sciences, allowing us to embark for the first time on a truly comprehensive study of vertebrate evolution, the results of which will touch nearly every aspect of vertebrate biological enquiry.

The predictability of evolution: glimpses into a post-Darwinian world

The very success of the Darwinian explanation, in not only demonstrating evolution from multiple lines of evidence but also in providing some plausible explanations, paradoxically seems to have served to have stifled explorations into other areas of investigation. The fact of evolution is now almost universally yoked to the assumption that its outcomes are random, trends are little more than drunkard's walks, and most evolutionary products are masterpieces of improvisation and far from perfect. But is this correct? Let us consider some alternatives. Is there evidence that evolution could in anyway be predictable? Can we identify alternative forms of biological organizations and if so how viable are they? Why are some molecules so extraordinarily versatile, while others can be spoken of as "molecules of choice"? How fortuitous are the major transitions in the history of life? What implications might this have for the Tree of Life? To what extent is evolutionary diversification constrained or facilitated by prior states? Are evolutionary outcomes merely sufficient or alternatively are they highly efficient, even superb? Here I argue that in sharp contradistinction to an orthodox Darwinian view, not only is evolution much more predictable than generally assumed but also investigation of its organizational substrates, including those of sensory systems, which indicates that it is possible to identify a predictability to the process and outcomes of evolution. If correct, the implications may be of some significance, not least in separating the unexceptional Darwinian mechanisms from underlying organizational principles, which may indicate evolutionary inevitabilities.