recent extinctions disturb path to equilibrium diversity in Caribbean bats

Similar documents
Zoogeography of Antillean Bats

Review of genetic diversification of bats in the Caribbean and biogeographic relationships to Neotropical species based on DNA barcodes

Bat Management Plan for Jamaica

MAMMALS OF THE WEST INDIES.

LIZARD EVOLUTION VIRTUAL LAB

PHYLOGENY OF PHYLLOSTOMID BATS (MAMMALIA: CHIROPTERA): DATA FROM DIVERSE MORPHOLOGICAL SYSTEMS, SEX CHROMOSOMES, AND RESTRICTION SITES

The Origin of Species: Lizards in an Evolutionary Tree

COMPARATIVE MORPHOLOGY OF THE TONGUE AND PALATE IN NOCTILIONIDAE AND MORMOOPIDAE (CHIROPTERA: NOCTILIONOIDEA).

A Revision of Extant Greater Antillean Bats of the Genus Natalus

Carlos E. L. Esbérard 1 & Davor Vrcibradic 2

Foraging Behavior, Home Range, Movements and Activity Patterns of Epicrates inornatus (Boidae) at Mata de Plátano Reserve in Arecibo, Puerto Rico

How to Tell the Difference Between Native Rock Iguanas and Invasive Green Iguanas. By Elaine A. Powers Illustrated by Anderson Atlas

Lecture 11 Wednesday, September 19, 2012

Evolution of Biodiversity

West Indian Mammals from the Albert Schwartz Collection: Biological and Historical Information

Dynamic evolution of venom proteins in squamate reptiles. Nicholas R. Casewell, Gavin A. Huttley and Wolfgang Wüster

Inferring Ancestor-Descendant Relationships in the Fossil Record

Fossil Vertebrates from the Bahamas

Biodiversity and Distributions. Lecture 2: Biodiversity. The process of natural selection

The Origin of Species: Lizards in an Evolutionary Tree

Morphological diagnoses of higher-level phyllostomid taxa (Chiroptera: Phyllostomidae)

Bio 1B Lecture Outline (please print and bring along) Fall, 2006

The Making of the Fittest: LESSON STUDENT MATERIALS USING DNA TO EXPLORE LIZARD PHYLOGENY

INQUIRY & INVESTIGATION

LABORATORY EXERCISE: CLADISTICS III. In fact, cladistics is becoming increasingly applied in a wide range of fields. Here s a sampling:

Evolution as Fact. The figure below shows transitional fossils in the whale lineage.

Chapter 22 Darwin and Evolution by Natural Selection

Cyprus biodiversity at risk

HISTORICAL BIOGEOGRAPHY OF WEST INDIAN VERTEBRATES

CLADISTICS Student Packet SUMMARY Phylogeny Phylogenetic trees/cladograms

Evolution of Birds. Summary:

What are taxonomy, classification, and systematics?

Taxonomic Review of the Wingmite Genus Cameronieta (Acari: Spinturnicidae) on Neotropical Bats, with a New Species from Northeastern Brazil

Postilla PEABODY MUSEUM OF NATURAL HISTORY YALE UNIVERSITY NEW HAVEN, CONNECTICUT, U.S.A.

Cladistics (reading and making of cladograms)

Genetic divergence of early song discrimination between two young songbird species

AP Biology. AP Biology

Title: Phylogenetic Methods and Vertebrate Phylogeny

September Population analysis of the Spaniel (English Springer) breed

Evolution by Natural Selection

STAT170 Exam Preparation Workshop Semester

Five albino bats from Guerrero and Colima, Mexico

BIMONTHLY BULLETIN of the CAYMAN ISLANDS DEPARTMENT of ENVIRONMENT S TERRESTRIAL RESOURCES UNIT

Molecular Phylogeny and Biogeography of West Indian Teiid Lizards of the Genus Ameiva

Biodiversity and Extinction. Lecture 9

A Mitochondrial DNA Phylogeny of Extant Species of the Genus Trachemys with Resulting Taxonomic Implications

Evolution by Natural Selection

Modern Evolutionary Classification. Lesson Overview. Lesson Overview Modern Evolutionary Classification

USING DNA TO EXPLORE LIZARD PHYLOGENY

Nomination of Populations of Dingo (Canis lupus dingo) for Schedule 1 Part 2 of the Threatened Species Conservation Act, 1995

UNIT III A. Descent with Modification(Ch19) B. Phylogeny (Ch20) C. Evolution of Populations (Ch21) D. Origin of Species or Speciation (Ch22)

Evolution by Natural Selection

Chapter 16: Evolution Lizard Evolution Virtual Lab Honors Biology. Name: Block: Introduction

Supplementary Fig. 1: Comparison of chase parameters for focal pack (a-f, n=1119) and for 4 dogs from 3 other packs (g-m, n=107).

Living Planet Report 2018

RWO 166. Final Report to. Florida Cooperative Fish and Wildlife Research Unit University of Florida Research Work Order 166.

Trait-based diversification shifts reflect differential extinction among fossil taxa

VERTEBRATA PALASIATICA

Amphibian Ark Conservation Needs Assessment - Cuba, March 2011 Page 1

Quiz Flip side of tree creation: EXTINCTION. Knock-on effects (Crooks & Soule, '99)

Soleglad, Fet & Lowe: Hadrurus spadix Subgroup

Interpreting Evolutionary Trees Honors Integrated Science 4 Name Per.

What defines an adaptive radiation? Macroevolutionary diversification dynamics of an exceptionally species-rich continental lizard radiation

Global comparisons of beta diversity among mammals, birds, reptiles, and amphibians across spatial scales and taxonomic ranks

Yr 11 Evolution of Australian Biota Workshop Students Notes. Welcome to the Australian Biota Workshop!! Some of the main points to have in mind are:

Accepted Manuscript. News & Views. Primary feather vane asymmetry should not be used to predict the flight capabilities of feathered fossils

Systematics, Taxonomy and Conservation. Part I: Build a phylogenetic tree Part II: Apply a phylogenetic tree to a conservation problem

Turtles (Testudines) Abstract

Caecilians (Gymnophiona)

of Nebraska - Lincoln

AMITY. Biodiversity & Its Conservation. Lecture 23. Categorization of Biodiversity - IUCN. By Prof. S. P. Bajpai. Department of Environmental Studies

Phylogeny of the Sciaroidea (Diptera): the implication of additional taxa and character data

Drivers of Extinction Risk in Terrestrial Vertebrates

ESTIMATING NEST SUCCESS: WHEN MAYFIELD WINS DOUGLAS H. JOHNSON AND TERRY L. SHAFFER

What is the evidence for evolution?

Malcolm L. McCallum Biodiversity and Conservation

Required and Recommended Supporting Information for IUCN Red List Assessments

Turks and Caicos Islands Turtle Project (TCITP)

Use of Agent Based Modeling in an Ecological Conservation Context

Supporting Online Material for

COULD YOU HAVE RIDDEN A HORSE MILLIONS OF YEARS AGO? Horse evolution goes back more than 55 million years

Natural Sciences 360 Legacy of Life Lecture 3 Dr. Stuart S. Sumida. Phylogeny (and Its Rules) Biogeography

September Population analysis of the Akita breed

Biol 160: Lab 7. Modeling Evolution

Rosenberger et al.: Capercaillie eggshell pigmentation, maculation and thickness

September Population analysis of the Beagle breed

September Population analysis of the Whippet breed

September Population analysis of the Schnauzer breed

New records of bats from Saint Barthélemy, French West Indies

Are node-based and stem-based clades equivalent? Insights from graph theory

Differential human impact on the survival of genetically distinct avian lineages

Supporting Online Material

September Population analysis of the Bearded Collie breed

Introduction to phylogenetic trees and tree-thinking Copyright 2005, D. A. Baum (Free use for non-commercial educational pruposes)

September Population analysis of the Fox Terrier (Wire) breed

ANIMAL RABIES IN NEPAL AND RACCOON RABIES IN ALBANY COUNTY, NEW YORK

These small issues are easily addressed by small changes in wording, and should in no way delay publication of this first- rate paper.

September Population analysis of the Giant Schnauzer breed

Predic'ng propaga'on of dengue with human mobility:

September Population analysis of the French Bulldog breed

Transcription:

Log-likelihood In the format provided by the authors and unedited. recent extinctions disturb path to equilibrium diversity in Caribbean bats Luis Valente, 2, rampal S. etienne 3 and Liliana M. Dávalos 4, 5 * SUPPLEMENTRY INFORMTION VOLUME: RTICLE NUMBER: 0026 Rate 0. 0.2 0.3 0.4-45 -40-35 -30 λ c µ γ λ a Likelihood 20 25 30 35 40 45 ssumed archipelago age (Ma) Supplementary Figure Sensitivity analysis of the DI model with archipelago age. Effect of varying assumed archipelago age on maximum likelihood values and parameter estimates of the DI model, M = 00. γ - per lineage rate of immigration; λ a - per lineage rate of anagenesis; λ c - per lineage rate of cladogenesis; µ - per lineage rate of extinction. NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

a 0 20 40 60 80 00 b 0 20 40 60 80 00-2 - 0 2 3 4 5 LLR 0 2 4 6 8 0 LLR Supplementary Fig. 2 - Bootstrap likelihood ratio test. a) distribution of logarithms of the likelihood ratio of diversity-dependent (DD) and diversity-independent (DI) model for islands simulated under DI. b) same as a) but for data simulated under the DI model. Black arrow - logarithm of the likelihood ratio for the real data; red arrow - logarithm of the likelihood ratio for a significance level of a = 0 05. NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 2 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

a b 0 0 20 30 40 50 60 0 0 20 30 40 50 60 70-0. 0.0 0. 0.2 0.3 0.4 μ - λ c 0.0 0. 0.2 0.3 0.4 0.5 μ - λ c Supplementary Figure 3 Distribution of difference between ML rate of extinction (µ) and ML rate of cladogenesis (λ c ). Based on datasets simulated with the parameters of the DI model for T = 20 (a) and T = 45 (b). When µ - λ c > 0, equilibrium dynamics are preferred (0 value marked with dotted line for visibility). NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 3 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

Proportion of equilibrium diversity 0.0 0.2 0.4 0.6 0.8.0.2 20 25 30 35 40 45 ssumed archipelago age (Ma) Supplementary Figure 4 Proportion of extant noctilionoid diversity compared with expected diversity at equilibrium. For all assumed archipelago ages T between 20 and 45 Ma. Black line high human impact scenario; grey line low human impact scenario; blue dotted line equilibrium value. NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 4 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

a b Number of lineages 0 5 0 5 20 25 Number of lineages 0 5 0 5 20 0 0 20 30 40 Time (Ma) 0 0 20 30 40 Time (Ma) Supplementary Figure 5 - Expected future number of independent colonist lineages under the DI model assuming high human impact, starting from current extant number of lineages. Scenario assuming all species went extinct as a result of human activity. a) ssuming archipelago age T = 20 Ma; b) assuming archipelago age T = 45 Ma. Black dotted line prehuman number of lineages. Light and dark shaded areas show the 2.5-97.5 and the 25-75 percentiles, respectively. NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 5 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

a b 0 200 400 600 0 200 400 600 20 40 60 80 Number of species 5 0 5 20 25 Number of colonization events c d 0 200 400 600 0 200 600 000 5 0 5 20 Number of endemic clades 0 2 4 6 8 0 Number of anagenetic species Supplementary Figure 6 Goodness of fit of the DI model, assuming archipelago age T = 20 Ma. Histograms show the total number of (a) species, (b) colonisation events, (c) endemic clades and (d) anagenetic species (endemic species with no relatives on the G) in 5,000 phylogenetic datasets simulated with the ML parameters of the DI model. Black vertical lines indicate median values across simulated datasets and arrows indicate the observed values in the G data. NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 6 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

a 0 00 300 500 b 0 00 300 500 30 40 50 60 70 80 90 Number of colonisation events 50 00 50 200 250 300 350 Number of cladogenesis events c 0 200 400 600 d 0 50 50 250 5 0 5 20 25 30 Number of anagenesis events 50 00 50 200 250 300 Number of extinction events Supplementary Figure 7 Number of events during the course of entire simulation, assuming archipelago age T = 20 Ma. Histograms show the total number of (a) colonisation, (b) cladogenesis, (c) anagenesis and (d) extinction events that took place during the course of 20 million years in each of 5,000 datasets simulated with the ML parameters of the DI model. Black vertical lines indicate median values across simulated datasets. NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 7 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

a b Number of endemic species 0 20 40 60 80 Number of endemic species 0 0 20 30 40 0 5 0 5 20 Time (Ma) 0 5 0 5 20 Time (Ma) Supplementary Figure 8 - Expected future noctilionoid bat endemic diversity in the Greater ntilles for low human impact scenario. Scenario assuming only three species went extinct as a result of human activity. Based on diversity-independent model, starting from current extant diversity, with archipelago ages of T = 20 Ma (a) and T = 45 Ma (b). Black dotted line - pre-human diversity levels; blue dotted line expected number of endemic species at equilibrium. Light and dark shaded areas show the 2.5-97.5 and the 25-75 percentiles, respectively. NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 8 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

Supplementary Table Noctilionoid bat species of the Greater ntilles. Including extinct and extirpated taxa. Taxon Brachyphylla cavernarum Brachyphylla nana nana Brachyphylla nana pumila Erophylla bombifrons Erophylla sezekorni Phyllonycteris aphylla Clade Code (Supplementary Table 2) Source of phylogenetic data Geographical distribution WI Endemic G + Bahamas G + TC G Endemic Status G + Bahamas + TC G Endemic Phyllonycteris major * G Endemic Extinct 2 Phyllonycteris poeyi obtusa Phyllonycteris poeyi poeyi Mormoops blainvillei B G + Bahamas G Endemic WI Endemic Mormoops magna B * G Endemic Extinct 3 Mormoops megalophylla Monophyllus plethodon C D SUPPLEMENTRY INFORMTION WI + Continent Extirpated 3 WI Endemic Extirpated 2 Monophyllus redmani D G + Bahamas + TC Desmodus puntajudensis E * G Endemic Extinct 4,5 Pteronotus macleayii F Pteronotus quadridens F Glossophaga soricina G G + Bahamas G + Bahamas WI + Continent Cubanycteris silvai H a G Endemic Extinct 6 riteus flavescens Phyllops falcatus I I G Endemic G Endemic Phyllops silvai I a G Endemic Extinct 7 Phyllops vetus I a G Endemic Extinct 8 Stenoderma rufum Tonatia saurophila I J G Endemic WI + Continent Extirpated 9 rtibeus anthonyi K a G Endemic Extinct 0 Pteronotus parnellii portoricensis Pteronotus parnellii pusillus L L G Endemic G Endemic Pteronotus pristinus L * WI + Continent Extinct 3 Pteronotus sp. nov. L * G Endemic Extinct Noctilio leporinus Pteronotus parnelli parnellii N M WI + Continent G + Bahamas NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 9 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

Macrotus waterhousii jamaicensis O 2 G Endemic Macrotus waterhousii minor O 2 G + Bahamas + TC Macrotus waterhousii O 2 G Endemic waterhousii (Hispaniola) Macrotus waterhousii O * WI Endemic Extinct 3 waterhousii (Puerto Rico & nguilla) rtibeus jamaicensis P 4 WI + Continent - Lineage extinct from the G archipelago a Source trees deposited in Dryad (doi:0.506/dryad.9mk20) * - Sampled as missing species in clade G - Greater ntilles WI - West Indies TC - Turks & Caicos NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 0 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

Supplementary Table 2 Independent colonisation events of noctilionoid bats in the Greater ntilles. Estimated times of colonisation in millions of years from present. nagenetic clades are those with a single endemic species whereas cladogenetic clades have two or more endemic species. Independent colonisation event Clade code Number of species Time of colonisation 95% CI Clade type Brachyphyllini + Phyllonycterini 9 7.97 5.46-20.45 Cladogenetic Mormoops blainvillei B 2 4.7 9.32-9.37 Cladogenetic Mormoops megalophylla C 4.7 9.32-9.37* Non-endemic Monophyllus D 2 3.94.63-6.22 Cladogenetic Desmodus puntajudensis E 2.39 9.50-6.08 nagenetic Pteronotus macleayii F 2 0.72 8.48-2.94 Cladogenetic Glossophaga soricina G 7.5 5.54-8.74* Non-endemic Cubanycteris H 6.83 5.08-8.78 nagenetic Sternodermatina I 5 6.6 5.20-7.24 Cladogenetic Tonatia saurophila J 4.5 2.84-6.5* Non-endemic rtibeus anthonyi K 3.32.64-5.20 nagenetic Pteronotus parnellii pusillus L 4 2.59 2.42-4.07 Cladogenetic Noctilio leporinus M 2.55.67-3.56* Non-endemic Pteronotus parnellii parnellii N.72.7-2.44 nagenetic Macrotus waterhousii O 4 0.28 0.06-0.58 Cladogenetic rtibeus jamaicensis P 0.03 0.0-0.04 Non-endemic Lineage extinct from the archipelago. * Upper bound for age of colonization. NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

Supplementary Table 3 DISIE maximum-likelihood analyses for M = 20. ML parameters of the diversity-independent (DI) and diversity-dependent models (DD) for the two different assumed archipelago ages. Results for consensus dataset. λ c - per lineage rate of cladogenesis; µ - per lineage rate of extinction; K carrying capacity; γ - per lineage rate of immigration; λ a - per lineage rate of anagenesis; LogLik Log likelihood. Dataset rchipelago age (Ma) Model λ c µ K γ λ a LogLik High human impact 20 DI 0.26 0.29-0.37 0.3-6.69 DD 0.4 0.33 5 0.46 0.24-5.44 45 DI 0.29 0.38-0.70 0.33-8.64 DD 0.42 0.40 8 0.76 0.25-7.60 Low human impact 20 DI 0.29 0.32-0.097 0.36-90.70 DD 0.43 0.34 8 0.095 0.30-89.92 45 DI 0.33 0.4-0.8 0.37-92.5 DD 0.43 0.42 25 0.5 0.32-9.64 Contemporary dataset 20 DI 0.22 0.24-0.082 0.34-88.45 DD 0.33 0.27 5 0.084 0.32-87.82 45 DI 0.25 0.33-0.03 0.37-90.05 DD 0.33 0.34 20 0.03 0.33-89.63 NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 2 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

Supplementary Table 4 95% confidence intervals of the ML parameter estimates for the high human impact dataset. Obtained by fitting DISIE models across 00 posterior datasets incorporating phylogenetic uncertainty, for mainland pool size M = 00 species. λ c - per lineage rate of cladogenesis; µ - per lineage rate of extinction; K carrying capacity; γ - per lineage rate of immigration; λ a - per lineage rate of anagenesis. rchipelago age (Ma) Model λ c µ K γ λ a 20 DI 0.24 0.32 0.26-0.38-0.025 0.033 0.5 0.26 DD 0.27 0.38 0.26 0.39 25-0.024 0.032 0. 0.25 45 DI 0.27 0.35 0.36 0.46-0.032 0.040 0.2 0.26 DD 0.29 0.4 0.36 0.47 33-0.03 0.039 0.04 0.24 NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 3 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.

Supporting information references. Rojas, D., Warsi, O. M. & Dávalos, L. M. Bats (Chiroptera: Noctilionoidea) challenge recent origin of neotropical diversity. Syst. Biol. 65, 432 448 (206). 2. nthony, H. E. Two new fossil bats from Porto Rico. Bull. m. Museum Nat. Hist. 37, 565 568 (97). 3. Silva-Taboada, G. Fossil Chiroptera from cave deposits in central Cuba, with description of two new species (Genera Pteronotus and Mormoops) and the first West Indian record of Mormoops megalophylla. cta Zool. Cracoviensia 9, 33 73 (974). 4. Suárez, W. Taxonomic status of the Cuban vampire bat (Chiroptera: Phyllostomidae: Desmodontinae: Desmodus). Caribb. J. Sci. 4, 76 767 (2005). 5. Woloszyn, B. & Mayo, N. Postglacial remnants of a vampire bat (Chiroptera: Desmodus) from Cuba. cta Zool. Cracoviensia 9, 253 265 (974). 6. Mancina, C.. & García-Rivera, L. New genus and species of fossil bat (Chiroptera: Phyllostomidae) from Cuba. Caribb. J. Sci. 4, 22 27 (2005). 7. Suárez, W. & Díaz-Franco, S. new fossil bat (Chiroptera: Phyllostomidae) from a quaternary Cave Deposit in Cuba. Caribb. J. Sci. 39, 37 377 (2003). 8. nthony, H. new rabbit and a new bat from Neotropical regions. Bull. m. Museum Nat. Hist. 35, 335 337 (97). 9. Koopman, K. & Williams, E. Fossil Chiroptera collected by H.E. nthony in Jamaica. m. Museum Novit. 99-920 59, 29 (95). 0. nthony, H.. Mammals collected in eastern Cuba in 97: with descriptions of two new species. Bull. m. Museum Nat. Hist. 4, 625 674 (99).. Morgan, G. in Biogeography of the West Indies: patterns and perspective (eds. Woods, C.. & Sergile, F. E.) 369 407 (CRC Press, 200). 2. Fleming, T. H., Murray, K. L. & Carstens, B. C. in Evolution, ecology, and conservation of island bats (eds. Fleming, T. H. & Racey, P..) 6 50 (University of Chicago Press, 200). 3. Choate, J. R. & Birney, E. C. Sub-recent Insectivora and Chiroptera from Puerto Rico, with the description of a new bat of the genus Stenoderma. J. Mammal. 49, 400 42 (968). 4. Dávalos, L. M., Wei, W. & Russell,. L. Ecological release and genetic signatures of population growth in ntillean rtibeus jamaicensis. Bat Res. News 53, 78 79 (202). NTURE ECOLOGY & EVOLUTION DOI: 0.038/s4559-06-0026 www.nature.com/natecolevol 4 207 Macmillan Publishers Limited, part of Springer Nature. ll rights reserved.