Origin of Madagascar's extant fauna: A perspective from amphibians, reptiles and other non‐flying vertebrates

Analysis of intrinsic evolutionary factors leading to microendemic distributions in New Caledonian leaf beetles

Microendemicity, or the condition of some species having local ranges, is a relatively common pattern in nature. However, the factors that lead to this pattern are still largely unknown. Most studies addressing this issue tend to focus on extrinsic factors associated with microendemic distributions, such as environmental conditions, hypothesising a posteriori about underlying potential speciation mechanisms, linked or not to these conditions. Here, we use a multi-faceted approach mostly focusing on intrinsic factors instead, namely diversification dynamics and speciation modes in two endemic sibling genera of leaf beetles with microendemic distributions, Taophila and Tricholapita, in a microendemicity hotspot, New Caledonia. Results suggest that the diversification rate in this lineage slowed down through most of the Neogene and consistently with a protracted speciation model possibly combined with several ecological and environmental factors potentially adding rate-slowing effects through time. In turn, species accumulated following successive allopatric speciation cycles, possibly powered by marked geological and climatic changes in the region in the last 25 million years, with daughter species ranges uncorrelated with the time of speciation. In this case, microendemicity seems to reflect a mature state for the system, rather than a temporary condition for recent species, as suggested for many microendemic organisms.

Ancestral transoceanic colonization and recent population reduction in a nonannual killifish from the Seychelles archipelago

Whether freshwater fish colonize remote islands following tectonic or transoceanic dispersal remains an evolutionary puzzle. Integrating dating of known tectonic events with phylogenomics and current species distribution, we find that killifish species distribution is not explained by species dispersal by tectonic drift only. Investigating the colonization of a nonannual killifish (golden panchax, Pachypanchax playfairii) on the Seychelle islands, we found genetic support for transoceanic dispersal and experimentally discovered an adaptation to complete tolerance to seawater. At the macroevolutionary scale, despite their long-lasting isolation, nonannual golden panchax show stronger genome-wide purifying selection than annual killifishes from continental Africa. However, progressive decline in effective population size over a more recent timescale has probably led to the segregation of slightly deleterious mutations across golden panchax populations, which represents a potential threat for species preservation in the long term.

Sympatric lineages in the Mantidactylus ambreensis complex of Malagasy frogs originated allopatrically rather than by in-situ speciation

Madagascar's biota is characterized by a high degree of microendemism at different taxonomic levels, but how colonization and in-situ speciation contribute to the assembly of local species communities has rarely been studied on this island. Here we analyze the phylogenetic relationships of riparian frogs of the Mantidactylus ambreensis species complex, which is distributed in the north of Madagascar and was originally described from Montagne d'Ambre, an isolated mountain of volcanic origin, currently protected within Montagne d'Ambre National Park (MANP). Data from mitochondrial DNA, and phylogenomic data from FrogCap, a sequence capture method, independently confirm that this species complex is monophyletic within the subgenus Ochthomantis, and identify two main clades within it. These two clades are separated by 5.6-6.8% pairwise distance in the mitochondrial 16S rRNA gene and co-occur in MANP, with one distributed at high elevations (940-1375 m a.s.l.) and the other at lower elevations (535-1010 m a.s.l.), but show almost no haplotype sharing in the nuclear RAG1 gene. This occurrence in syntopy without admixture confirms them as independent evolutionary lineages that merit recognition as separate species, and we here refer to them as high-elevation (HE) and low-elevation (LE) lineage; they will warrant taxonomic assessment to confidently assign the name ambreensis to one or the other. Populations of the M. ambreensis complex from elsewhere in northern Madagascar all belong to the LE lineage, although they do occur over a larger elevational range than in Montagne d'Ambre (285-1040 m a.s.l.). Within LE there are several phylogroups (LE1-LE4) of moderately deep divergence (1.5-2.8% in 16S), but phylogroup LE4 that occurs in MANP has a deeply nested phylogenetic position, as recovered separately by mitochondrial and sequence capture datasets. This suggests that HE and LE did not diverge by a local fission of lower and upper populations, but instead arose through a more complex biogeographic scenario. The branching pattern of phylogroups LE1-LE4 shows a clear south-to-north phylogeographic pattern. We derive from these results a testable hypothesis of vicariant speciation that restricted the HE lineage to MANP and the LE candidate species to a climatic refugium further south, with subsequent northwards range expansion and secondary colonization of MANP by LE. These results provide an example for complex assembly of local microendemic amphibian faunas on Madagascar.

Biogeography of worm lizards (Amphisbaenia) driven by end-Cretaceous mass extinction

Worm lizards (Amphisbaenia) are burrowing squamates that live as subterranean predators. Their underground existence should limit dispersal, yet they are widespread throughout the Americas, Europe and Africa. This pattern was traditionally explained by continental drift, but molecular clocks suggest a Cenozoic diversification, long after the break-up of Pangaea, implying dispersal. Here, we describe primitive amphisbaenians from the North American Palaeocene, including the oldest known amphisbaenian, and provide new and older molecular divergence estimates for the clade, showing that worm lizards originated in North America, then radiated and dispersed in the Palaeogene following the Cretaceous-Palaeogene (K-Pg) extinction. This scenario implies at least three trans-oceanic dispersals: from North America to Europe, from North America to Africa and from Africa to South America. Amphisbaenians provide a striking case study in biogeography, suggesting that the role of continental drift in biogeography may be overstated. Instead, these patterns support Darwin and Wallace's hypothesis that the geographical ranges of modern clades result from dispersal, including oceanic rafting. Mass extinctions may facilitate dispersal events by eliminating competitors and predators that would otherwise hinder establishment of dispersing populations, removing biotic barriers to dispersal.

New material of Beelzebufo, a hyperossified frog (Amphibia: Anura) from the late cretaceous of Madagascar

The extant anuran fauna of Madagascar is exceptionally rich and almost completely endemic. In recent years, many new species have been described and understanding of the history and relationships of this fauna has been greatly advanced by molecular studies, but very little is known of the fossil history of frogs on the island. Beelzebufo ampinga, the first named pre-Holocene frog from Madagascar, was described in 2008 on the basis of numerous disarticulated cranial and postcranial elements from the Upper Cretaceous (Maastrichtian) Maevarano Formation of Madagascar. These specimens documented the presence of a hyperossified taxon that differed strikingly from extant Malagasy frogs in its large size and heavy coarse cranial exostosis. Here we describe and analyse new, articulated, and more complete material of the skull, vertebral column, and hind limb, as well as additional isolated elements discovered since 2008. μCT scans allow a detailed understanding of both internal and external morphology and permit a more accurate reconstruction. The new material shows Beelzebufo to have been even more bizarre than originally interpreted, with large posterolateral skull flanges and sculptured vertebral spine tables. The apparent absence of a tympanic membrane, the strong cranial exostosis, and vertebral morphology suggest it may have burrowed during seasonally arid conditions, which have been interpreted for the Maevarano Formation from independent sedimentological and taphonomic evidence. New phylogenetic analyses, incorporating both morphological and molecular data, continue to place Beelzebufo with hyloid rather than ranoid frogs. Within Hyloidea, Beelzebufo still groups with the South American Ceratophryidae thus continuing to pose difficulties with both biogeographic interpretations and prior molecular divergence dates.

First large-scale DNA barcoding assessment of reptiles in the biodiversity hotspot of Madagascar, based on newly designed COI primers

BACKGROUND

DNA barcoding of non-avian reptiles based on the cytochrome oxidase subunit I (COI) gene is still in a very early stage, mainly due to technical problems. Using a newly developed set of reptile-specific primers for COI we present the first comprehensive study targeting the entire reptile fauna of the fourth-largest island in the world, the biodiversity hotspot of Madagascar.

METHODOLOGY/PRINCIPAL FINDINGS

Representatives of the majority of Madagascan non-avian reptile species (including Squamata and Testudines) were sampled and successfully DNA barcoded. The new primer pair achieved a constantly high success rate (72.7-100%) for most squamates. More than 250 species of reptiles (out of the 393 described ones; representing around 64% of the known diversity of species) were barcoded. The average interspecific genetic distance within families ranged from a low of 13.4% in the Boidae to a high of 29.8% in the Gekkonidae. Using the average genetic divergence between sister species as a threshold, 41-48 new candidate (undescribed) species were identified. Simulations were used to evaluate the performance of DNA barcoding as a function of completeness of taxon sampling and fragment length. Compared with available multi-gene phylogenies, DNA barcoding correctly assigned most samples to species, genus and family with high confidence and the analysis of fewer taxa resulted in an increased number of well supported lineages. Shorter marker-lengths generally decreased the number of well supported nodes, but even mini-barcodes of 100 bp correctly assigned many samples to genus and family.

CONCLUSIONS/SIGNIFICANCE

The new protocols might help to promote DNA barcoding of reptiles and the established library of reference DNA barcodes will facilitate the molecular identification of Madagascan reptiles. Our results might be useful to easily recognize undescribed diversity (i.e. novel taxa), to resolve taxonomic problems, and to monitor the international pet trade without specialized expert knowledge.

Mountain refugia play a role in soil arthropod speciation on Madagascar: a case study of the endemic giant fire-millipede genus Aphistogoniulus

To elucidate the speciation mechanisms prevalent within hotspots of biodiversity, and the evolutionary processes behind the rise of their species-rich and endemic biota, we investigated the phylogeny of the giant fire-millipede genus Aphistogoniulus Silvestri, 1897, a Malagasy endemic. This study is the first comprehensive (molecular and morphological) phylogenetic study focusing on millipede (class Diplopoda) speciation on Madagascar. The morphological analysis is based on 35 morphological characters and incorporates ten described as well as two newly described species (A. rubrodorsalisn. sp. and A. jeekelin. sp.) of Aphistogoniulus. The molecular analysis is based on both mitochondrial (COI and 16S), and nuclear genes (complete 18S rDNA), together comprised of 3031 base pairs, which were successfully sequenced for 31 individual specimens and eight species of Aphistogoniulus. In addition to the null-model (speciation by distance), two diversification models, mountain refugia and ecotone shift, were discovered to play a role in the speciation of soil arthropods on Madagascar. Mountain refugia were important in the speciation of the A. cowani clade, with three species occurring in the Andringitra and Ranomafana Mountains in the southeast (A. cowani), the Ambohijanahary and Ambohitantely Mountains in the mid-west (A. sanguineus), and the Marojejy Mountain in the northeast (A. rubrodorsalisn. sp.). An ecotone shift from the eastern rainforest to the unique subarid spiny forest of Mahavelo was discovered in the A. vampyrus-A. aridus species-pair. In the monophyletic A. diabolicus clade, evidence for divergent evolution of sexual morphology was detected: species with greatly enlarged gonopods are sister-taxa to species with normal sized gonopods. Among the large-bodied Spirobolida genera of Madagascar, Colossobolus and Sanguinobolus were found to be close sister-genera to Aphistogoniulus. Forest destruction has caused forest corridors between populations to disappear, which might limit the possible resolution of biogeographic analyses on Madagascar.

The origins of the giant pill-millipedes from Madagascar (Diplopoda: Sphaerotheriida: Arthrosphaeridae)

Giant pill-millipedes (order Sphaerotheriida) are large-bodied millipedes without poison glands which can roll-up into a complete ball. Their disconnected area of distribution spanning South Africa, Madagascar, India, SE Asia, Australia and New Zealand makes them interesting model organisms for biogeographic studies. The here presented phylogeny is based on a molecular dataset covering all areas of distribution with a special focus on Madagascar, where some species of giant pill-millipedes show island gigantism, reaching the size of a baseball. For our study, two mitochondrial genes (partial 16S rRNA and COI) as well as the complete nuclear 18S rDNA were sequenced. While many recent vertebrate studies hint that the ancestors of the recent Malagasy fauna crossed the >350 km wide Mozambique Channel several times, no such crossing was discovered in the Sphaerotheriida. For the first time in a molecular phylogenetic study of soil arthropods, a Madagascar-India group, the family Arthrosphaeridae, is recovered, hinting to a Gondwanan origin of the Sphaerotheriida. The Malagasy-Indian family Arthrosphaeridae forms a monophyletic, statistically well-supported group in all obtained trees. The giant pill-millipedes from Madagascar are paraphyletic because the Malagasy genus Sphaeromimus is the sister-taxon of the Indian Arthrosphaera. In Sphaeromimus, an ecotone shift occurred only once: the spiny forest species Sphaeromimus musicus forms the sister-clade to the species collected in rainforests and littoral rainforests. The two species of the Malagasy genus Zoosphaerium which express island gigantism form a monophyletic group in some trees, but these trees lack good statistical support. Deeper nodes inside the Sphaerotheriida, like the position of the Australian genera Procyliosoma and Epicyliosoma, the Southeast Asian family Zephroniidae and the South African genus Sphaerotherium could not be resolved. This study is the first genetic study inside the order Sphaerotheriida and provides a proper basis for future molecular biogeographic studies in millipedes and soil organisms from Madagascar.

Phylogeny and Biogeography of a Cosmopolitan Frog Radiation: Late Cretaceous Diversification Resulted in Continent-Scale Endemism in the Family Ranidae

Ranidae is a large anuran group with a nearly cosmopolitan distribution. We investigated the phylogenetic relationships and early biogeographic history of ranid frogs, using 104 representatives of all subfamilies and families, sampled from throughout their distribution. Analyses of approximately 1570 bp of nuclear gene fragments (Rag-1, rhod, Tyr) and approximately 2100 bp of the mitochondrial genome (12S rRNA, tRNAVAL, 16S rRNA) indicate that the monophyly of several taxa can be rejected with high confidence. Our tree is characterized by a clear historical association of each major clade with one Gondwanan plate. This prevalence of continent-scale endemism suggests that plate tectonics has played a major role in the distribution of ranid frogs. We performed dispersal-vicariance analyses, as well as analyses constrained by paleogeographic data, to estimate ancestral distributions during early ranid diversification. Additionally, we used molecular clock analyses to evaluate whether these scenarios fit the temporal framework of continental breakup. Our analyses suggest that a scenario in which the ancestors of several clades (Rhacophorinae, Dicroglossinae, Raninae) reached Eurasia via the Indian subcontinent, and the ancestor of Ceratobatrachinae entered via the Australia-New Guinea plate, best fits the paleogeographic models and requires the fewest number of dispersal/vicariance events. However, several alternatives, in which part of the ranid fauna colonized Laurasia from Africa, are not significantly worse. Most importantly, all hypotheses make clear predictions as to where to expect key fossils and where to sample other living ranids, and thus constitute a strong basis for further research.