Jacob A. Esselstyn

DO GEOLOGICAL OR CLIMATIC PROCESSES DRIVE SPECIATION IN DYNAMIC ARCHIPELAGOS? THE TEMPO AND MODE OF DIVERSIFICATION IN SOUTHEAST ASIAN SHREWS

Geological and climatic processes potentially alter speciation rates by generating and modifying barriers to dispersal. In Southeast Asia, two processes have substantially altered the distribution of land. Volcanic uplift produced many new islands during the Miocene—Pliocene and repeated sea level fluctuations during the Pleistocene resulted in intermittent land connections among islands. Each process represents a potential driver of diversification. We use a phylogenetic analysis of a group of Southeast Asian shrews (Crocidura) to examine geographic and temporal processes of diversification. In general, diversification has taken place in allopatry following the colonization of new areas. Sulawesi provides an exception, where we cannot reject within-island speciation for a clade of eight sympatric and syntopic species. We find only weak support for temporally declining diversification rates, implying that neither volcanic uplift nor sea level fluctuations had a strong effect on diversification rates. We suggest that dynamic archipelagos continually offer new opportunities for allopatric diversification, thereby sustaining high speciation rates over long periods of time, or Southeast Asian shrews represent an immature radiation on a density-dependent trajectory that has yet to fill geographic and ecological space.

Single-locus species delimitation: a test of the mixed Yule–coalescent model, with an empirical application to Philippine round-leaf bats

Prospects for a comprehensive inventory of global biodiversity would be greatly improved by automating methods of species delimitation. The general mixed Yule–coalescent (GMYC) was recently proposed as a potential means of increasing the rate of biodiversity exploration. We tested this method with simulated data and applied it to a group of poorly known bats (Hipposideros) from the Philippines. We then used echolocation call characteristics to evaluate the plausibility of species boundaries suggested by GMYC. In our simulations, GMYC performed relatively well (errors in estimated species diversity less than 25%) when the product of the haploid effective population size (Ne) and speciation rate (SR; per lineage per million years) was less than or equal to 105, while interspecific variation in Ne was twofold or less. However, at higher but also biologically relevant values of Ne × SR and when Ne varied tenfold among species, performance was very poor. GMYC analyses of mitochondrial DNA sequences from Philippine Hipposideros suggest actual diversity may be approximately twice the current estimate, and available echolocation call data are mostly consistent with GMYC delimitations. In conclusion, we consider the GMYC model useful under some conditions, but additional information on Ne, SR and/or corroboration from independent character data are needed to allow meaningful interpretation of results.

Phylogeny and biogeography of Philippine bent-toed geckos (Gekkonidae: Cyrtodactylus) contradict a prevailing model of Pleistocene diversification.

In the Philippines, Pleistocene sea level oscillations repeatedly connected and isolated neighboring islands. Hence, an understanding of the island platforms adjoined during periods of low sea level has led biologists to a suite of expectations that, taken together, represent a paradigm for the process of recent diversification in southeast Asia. We employ statistical tests of phylogenetic topology and population genetic analyses of widespread species of bent-toed geckos (Cyrtodactylus) to ascertain whether patterns of inter- and intra-specific diversity can be explained by a Pleistocene aggregate island model of diversification. Contrary to many classic studies of Philippine vertebrates, we find complex patterns that are only partially explained by past island connectivity. In particular, we determine that some populations inhabiting previously united island groups show substantial genetic divergence and are inferred to be polyphyletic. Additionally, greater genetic diversity is found within islands, than between them. Among the topological patterns inconsistent with the Pleistocene model, we note some similarities with other lineages, but no obviously shared causal mechanisms are apparent. Finally, we infer well-supported discordance between the gene trees inferred from mitochondrial and nuclear DNA sequences of two species, which we suspect is the result of incomplete lineage sorting. This study contributes to a nascent body of literature suggesting that the current paradigm for Philippine biogeography is an oversimplification requiring revision.

The role of repeated sea-level fluctuations in the generation of shrew (Soricidae: Crocidura) diversity in the Philippine Archipelago

Repeated sea-level fluctuations during and prior to the Pleistocene had a dramatic effect on the distribution of land and connection among islands between the Asian and Australian continents. These cycles of connection and isolation have long been recognized as an important factor determining the distribution and organization of biodiversity in the Philippines. However, surprisingly few studies have tested for predicted patterns of genetic diversity derived from a Pleistocene sea-level model of geography. Here, we examine evidence for fit to such a model in a widely distributed lineage of shrews (Crocidura). The topology of relationships among Crocidura from the Philippines is concordant with a Pleistocene sea-level model, but (1) AMOVAs reveal that genetic diversity is explained at least as well by modern islands as by Pleistocene islands; (2) Mantel tests reveal a significant influence of isolation by distance; and (3) the degree of genetic divergence between some populations connected by dry land during the last glacial maximum reveals isolation that almost certainly predates the most recent glacial activity. We further employ multiple strategies for inferring time-calibrated phylogenies, but these result in widely varying time estimates for the invasion of SE Asian islands by shrews. Overall, our results suggest Pleistocene sea-level fluctuations have been an important, but not dominant factor shaping shrew diversity.

A NEW SPECIES OF DESMALOPEX (PTEROPODIDAE) FROM THE PHILIPPINES, WITH A PHYLOGENETIC ANALYSIS OF THE PTEROPODINI

Abstract We describe a new species of flying fox of the genus Desmalopex from Mindoro Island, Philippines. Discrete and mensural morphological characters distinguish the new species from other flying foxes in Southeast Asia. The new species shares several probable morphological synapomorphies with Desmalopex leucopterus, including features of the pelage, patagia, dentition, and cranium, suggesting that the 2 species are closely related. We present phylogenetic analyses of mitochondrial DNA sequences, which support the taxonomic status of the new species and the recently revalidated genus Desmalopex. Together, D. leucopterus and the new species form a well-supported clade that may be sister to Pteropus + Acerodon, or perhaps more distantly related to these genera. Discovery of the new species highlights the need for continued biodiversity inventories in the Philippines, where new taxa are being discovered at a remarkable rate.

Phylogeography and historical demography of Polypedates leucomystax in the islands of Indonesia and the Philippines: evidence for recent human-mediated range expansion?

Southeast Asias widespread species offer unique opportunities to explore the effects of geographical barriers to dispersal on patterns of vertebrate lineage diversification. We analyzed mitochondrial gene sequences (16S rDNA) from a geographically widespread sample of 266 Southeast Asian tree frogs, including 244 individuals of Polypedates leucomystax and its close relatives. Our expectation was that lineages on island archipelagos would exhibit more substantial geographic structure, corresponding to the geological history of terrestrial connectivity in this region, compared to the Asian mainland. Contrary to predictions, we found evidence of numerous highly divergent lineages from a limited area on the Asian mainland, but fewer lineages with shallower divergences throughout oceanic islands of the Philippines and Indonesia. Surprisingly and in numerous instances, lineages in the archipelagos span distinct biogeographical provinces. Phylogeographic analyses identified four major haplotype clades; summary statistics, mismatch distributions, and Bayesian coalescent inference of demography provide support for recent range expansion, population growth, and/or admixture in the Philippine and some Sulawesi populations. We speculate that the current range of P. leucomystax in Southeast Asia is much larger now than in the recent past. Conversion of forested areas to monoculture agriculture and transportation of agricultural products between islands may have facilitated unprecedented population and range expansion in P. leucomystax throughout thousands of islands in the Philippine and Indonesian archipelagos.

The Challenges of Resolving a Rapid, Recent Radiation: Empirical and Simulated Phylogenomics of Philippine Shrews

Phylogenetic relationships in recent, rapid radiations can be difficult to resolve due to incomplete lineage sorting and reliance on genetic markers that evolve slowly relative to the rate of speciation. By incorporating hundreds to thousands of unlinked loci, phylogenomic analyses have the potential to mitigate these difficulties. Here, we attempt to resolve phylogenetic relationships among eight shrew species (genus Crocidura) from the Philippines, a phylogenetic problem that has proven intractable with small (< 10 loci) data sets. We sequenced hundreds of ultraconserved elements and whole mitochondrial genomes in these species and estimated phylogenies using concatenation, summary coalescent, and hierarchical coalescent methods. The concatenated approach recovered a maximally supported and fully resolved tree. In contrast, the coalescent-based approaches produced similar topologies, but each had several poorly supported nodes. Using simulations, we demonstrate that the concatenated tree could be positively misleading. Our simulations also show that the tree shape we tend to infer, which involves a series of short internal branches, is difficult to resolve, even if substitution models are known and multiple individuals per species are sampled. As such, the low support we obtained for backbone relationships in our coalescent-based inferences reflects a real and appropriate lack of certainty. Our results illuminate the challenges of estimating a bifurcating tree in a rapid and recent radiation, providing a rare empirical example of a nearly simultaneous series of speciation events in a terrestrial animal lineage as it spreads across an oceanic archipelago.

EVIDENCE FOR CLIMATE‐DRIVEN DIVERSIFICATION? A CAUTION FOR INTERPRETING ABC INFERENCES OF SIMULTANEOUS HISTORICAL EVENTS

Approximate Bayesian computation (ABC) is rapidly gaining popularity in population genetics. One example, msBayes, infers the distribution of divergence times among pairs of taxa, allowing phylogeographers to test hypotheses about historical causes of diversification in co‐distributed groups of organisms. Using msBayes, we infer the distribution of divergence times among 22 pairs of populations of vertebrates distributed across the Philippine Archipelago. Our objective was to test whether sea‐level oscillations during the Pleistocene caused diversification across the islands. To guide interpretation of our results, we perform a suite of simulation‐based power analyses. Our empirical results strongly support a recent simultaneous divergence event for all 22 taxon pairs, consistent with the prediction of the Pleistocene‐driven diversification hypothesis. However, our empirical estimates are sensitive to changes in prior distributions, and our simulations reveal low power of the method to detect random variation in divergence times and bias toward supporting clustered divergences. Our results demonstrate that analyses exploring power and prior sensitivity should accompany ABC model selection inferences. The problems we identify are potentially mitigable with uniform priors over divergence models (rather than classes of models) and more flexible prior distributions on demographic and divergence‐time parameters.

Evolutionary novelty in a rat with no molars

Rodents are important ecological components of virtually every terrestrial ecosystem. Their success is a result of their gnawing incisors, battery of grinding molars and diastema that spatially and functionally separates the incisors from the molars. Until now these traits defined all rodents. Here, we describe a new species and genus of shrew-rat from Sulawesi Island, Indonesia that is distinguished from all other rodents by the absence of cheek teeth. Moreover, rather than gnawing incisors, this animal has bicuspid upper incisors, also unique among the more than 2200 species of rodents. Stomach contents from a single specimen suggest that the species consumes only earthworms. We posit that by specializing on soft-bodied prey, this species has had no need to process food by chewing, allowing its dentition to evolve for the sole purpose of procuring food. Thus, the removal of functional constraints, often considered a source of evolutionary innovations, may also lead to the loss of the very same traits that fuelled evolutionary diversification in the past.

Species Interactions during Diversification and Community Assembly in an Island Radiation of Shrews

Background Closely related, ecologically similar species often have adjacent distributions, suggesting competitive exclusion may contribute to the structure of some natural communities. In systems such as island archipelagos, where speciation is often tightly associated with dispersal over oceanic barriers, competitive exclusion may prevent population establishment following inter-island dispersal and subsequent cladogenesis. Methodology/Principal Findings Using a combination of tools, we test the hypothesis that the distributions of shrew (Crocidura) species in the Philippines are the result of competitive exclusion preventing secondary invasion of occupied islands. We first compare ecological niche models between two widespread, allopatric species and find statistical support for their ecological similarity, implying that competition for habitat between these species is possible. We then examine dispersion patterns among sympatric species and find some signal for overdispersion of body size, but not for phylogenetic branch length. Finally, we simulate the process of inter-island colonization under a stochastic model of dispersal lacking ecological forces. Results are dependent on the geographic scope and colonization probability employed. However, some combinations suggest that the number of inter-island dispersal events necessary to populate the archipelago may be much higher than the minimum number of colonization events necessary to explain current estimates of species richness and phylogenetic relationships. If our model is appropriate, these results imply that alternative factors, such as competitive exclusion, may have influenced the process of inter-island colonization and subsequent cladogenesis. Conclusions/Significance We interpret the combined results as providing tenuous evidence that similarity in body size may prevent co-occurrence in Philippine shrews and that competitive exclusion among ecologically similar species, rather than an inability to disperse among islands, may have limited diversification in this group, and, possibly other clades endemic to island archipelagos.