Page 1of50 between this versionandtheVersionofrecord. Pleasecitethis articleasdoi:10.1111/geb.12716. through thecopyediting, typesetting, paginationandproofreadingprocess, whichmayleadtodifferences This istheauthormanuscript acceptedforpublicationandhasundergone full peerreviewbuthasnotbeen authors all Itescu of behalf Yuval on Sincerely, Biogeography Ecology and Global in for and at level publication ready now is our manuscript find you will hope We comment. each with of we how dealt by description point a find point Below please accordingly. manuscript the of these and revised most followed We suggestions. important their for and referee the editor the and thanks review supportive were see We to happy the comments. added her own Santos, Dr. Ana Editor, the and referee anonymous a by was new reviewed manuscript suggestion). The (Ref. GEB 2017 0100; title modified according referee’s to modified title reptiles’ (Ref. GEB 2017 0100; island co-occurring of our manuscript version attach aWe revised McGill, Dear Prof. Biogeography and EcologyGlobal Chief in Editor McGill, Professor Brian 07.12.2017
Author Manuscript This article isprotected by copyright. All rights reserved.
Global EcologyandBiogeography . . ‘Inconsistent patterns of body size evolution in in size evolution body of patterns ‘Inconsistent
suggestion and it would strengthen your results, but I will leave it up to you to decide if youyou decide to if I to up your it but leave will results, strengthen would andsuggestion it a good is I this think selection. based model you AIC that use and new reviewer suggests This your submission. previous greatly since has improved Iyour it and think paper, share view of his I of reviews. round first the in involved not was a that by referee evaluated was Your manuscript your manuscript. happened as with it to, like would we takes that longer it process.editorial Sometimes for the delay the in I apologize First, want to Author: the to Comments Ana Santos, Editor: TO AUTHORS COMMENTS EDITOR'S last paragraph, as part of your take home message. your message. of take home as paragraph, part last the you in theory former that rejecting you are stress to In need little. particular, although very discussion, the in also rewording needed is of this theory”. Some former rejects/contradicts “this that also reptiles.” And Aegean on evolution size body driversuniversal of are no “There like saying something sentence a straightforward need start to conclusions abstract the main the in Similarly, “Divergent”. than rather coherent” or “Non “Heterogeneous” I where title, use the heterogeneity such e.g. in put to front, in text the more, tweaking a bit it you I highlight should why think is That timeliness. and its of paper the importance the with line I’m in 100% so yougreat job, a correlations make on each based species separately explain to trying rather than species among differ responses that By highlighting rules”. “island of general lack the front and about make a in paperyou this put choseI’m to glad really so result, important most the is your species all predictors over of power predictive general lack the me To of ones. and a few small comments, got significant I’ve two more data. species with the analyses on additional the eliminate good to was decision a it me to so straightforward, an robust seems version this but round, first get the I review a to nice is work. didn’t This Author the to Comments 2 Referee: . Author the to Comments 1 Referee: TO AUTHORS COMMENTS REVIEWER Editor Dr. Ana Santos, islands. on evolution generality of theory size the contrast our results that discussion of the paragraph concluding the and in abstract the in and highlighted appendix scores a new in AICc included we below, As we specify comments. supporting for editor the her thank We universal. not rule is island the that generalyour conclusions more ofsuggestion highlighting his follow to you is do I to advise What strongly manuscript. the these analyses in include should
Author Manuscript This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography
Page 2of50 Page 3of50 predictor across all species? Here the total weight of a predictor for a species it would be the sum be sum the would it a species for Hereof weight a predictor total the species? all predictor across file. supplementary a new in for species each ( AICc 2) ≤ models best paper in cite this now We modest. too (2) is parameters additional adding to applies AIC penalty I which concur, results, perhaps raw the not although supplement, the in AICcprovide results after but readingit, I with certainly live could a biostatistician I’m not unconvincing. Since bit
low sample size, as some of yours are. Here, the pick would be to compare models according to models compare be to would yours pick the Here, are. of as some size, sample low with for datasets suited particularly are AICc And scores model. the in of variables amount the way, penalizes for it you correct the in used AIC I whether wonder Diniz Filho responseyour to a use AIC to not scores out, stand reasoning I the found will conclusions the so general drivers, lack the of on insights sound provides scores based R2 on of models your selection you that in I agree with would is Although selection. model remaining text the in problem main the me To it. to according changed text We the suggestion. forThanks this results from a different approach. approach. from a different results similar includes it if be more trustworthy yourwill But work informative. is often it how about you just tells it of relationships, the signs homogeneity the a on pick is provide certainly not does approach this What paradigms. analytical are across results that consistent show to useful I quite find So still AICs”. used they if had general patterns be would there they Surely use R2s! “but readers saying I other but hear also can meta analysis”, the with that showed we“but have youII saying can hear think of species. the regardless informative be a to general is predictor, it closer the 16, to closer the across species, all And is. a variable more the informative 1, closer to for the are AICw, models tested all Given that present. where is models it all of of AICw the per AIC weights all summing why just not so tested, are parameters weights same the if are AIC but datasets, based different on among comparable models not AIC are scores weights? based AIC on similar something why doing your not So meta analysis. is here, approach sound is frequentist a why reason actually real the your and results, compelling of most the me To the AICc of scores present and now get themselves, our impression readers the let to suggestion referee’s the accept we Nevertheless, ours. as such for studies selection model criteria based over information selection based model prefer for p value to good why example a serve as our results think as we point regarding this added also a of paragraph discussion We text. the and second, the want it), you not do still diseases of infectious best the (like models informative group a of non of best the AIC may lowest the be simply with why: model the first, explains compellingly and clearly http://onlinelibrary.wiley.com/doi/10.1111/1365 2664.13060/abstract) anygood?” model ‘best’ the criteria, is but usingselection information “Model 2017, Ecology, of Applied (Journal paper et al.’s recent Nally and Mac informative, most the is model which scores determine AICc) to (or AICon even relying comfortable with not are weof revision), round previous the in referees reply the to (and in text the weAs in stated be confusing. would quickly, and so approach the information abandon youI shouldn’t think species. Anyway, for several model best as the competing probability say,you’llwith, or 30% 20 havesure models I’m informative. most the is model a that probability the particular is, that Akaiketheir weight;
Author Manuscript This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography
Now added. added. Now as a valuable contribution for GEB in my opinion, and these other analyses can be published can be published analyses and these other my in for opinion, GEB contribution as a valuable out and stands and straightforward, clear is work The current sound. is few species aanalyses on more the detailed out of leaving decision I the that say. really to think I little Besides have that weights analysis. AIC the include we to 5, prefer not Appendix in AICc based models the show do although we predictor. Therefore, of effect the of the direction the information no about provides it mentions, rightly as referee the especially since, confusing, somewhat is meta analysis the in shown what we show have to using weights AIC that think but suggestion this for referee the thank We comment). (see first my result your important most highlight to similar, or “Strikingly” with paragraph I’d this L366. start text. the from sentence the we omitted interspecific differences, rather but the right), own its in studyquestion an interesting study perhaps (but focus of this are the not evolution size body differences in and since inter sexual the readers, byconfusion our any avoid want to we or only Since males. females only we if test even maintained species, is different in size affect body factorsdifferent that conclusion, general our but studythe species, some in of exist do single to predictors response and their females in males differences between misleading: a bit be may it phrased as we sentence originally the that realized now We incorrect) grammatically is (it sentence this and rephrase material), (as supplementary results these show L242 243. Please book 1967 inclusive rather more the than and Wilson, of MacArthur paper for 1963 the citation the just and retain (1957) for Darlington a citation we so added mechanism, (often contentious) any to point than rather islands, nature only of depauperate the emphasize want to actuallyWe ibt. just ETIB than rather mention rephrase and just theories, need cite No these latter to anagenesis. than be more important to starts wherezone cladogenesis radiation determine a area and increasing or that isolation with, familiar you one be may mention richness, to species increases also island the heterogeneity within habitat as,such e.g., that topics, many other includes (ibt) theory biogeography island (ETIB). Currently Biogeography Island Theory of Equilibrium you & Wilson’s MacArthur refer to here certain, Line74. be To Sentence rephrased. reptiles…” all data for the analysing informative, are more predictors which determine to and meta analyses linear models use “We analyse” for I’d “We substitute In 52 line you use. of analyses kind which you tell don’t In Abstract. methods the follow: comments minor Some elsewhere.
Author Manuscript
This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography
Page 4of50 Page 5of50 continuously to the caption to understand the figure. figure. the understand to caption the to continuously have jump to not reader the does so Y axes, variable the in names put Figures and Please 2 3. Done. for sentences. or separate 2 3 out canIprobably make read. to and difficult fragmentary sentence, is it this L483 486. Rephrase prediction. Heaney’s failed supporting that first the we not are as wording, other strong such use “strikingly” or to appropriate not is it specifically here that we think but suggestion, this for referee the thank We The names of the variables on the Y axes are now clearly marked in figures 2 & 3. & 2 3. figures clearly in marked now are Y the axes on The of names variables the
Author Manuscript This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography
Research paper Research Ann Arbor, MI 48109, USA48109, MI Arbor, Ann 4NS,UK E1 London snakes. Key words: 7 Greece 84, 6 5 4 USA 3 02138,MAUSA Cambridge, 2 1 Pafilis Panayiotis Valakos D. Sagonas Kostas ([email protected]), Slavenko Alex Donihue Colin ([email protected]), Itescu Yuval reptiles. island co-occurring in size evolution body of patterns Inconsistent School of Biological and Chemical Sciences, Queen Mary University of London,of University Mary Queen ChemicalSciences, andof Biological School USA TX 79409, Lubbock, Tech TexasSciences,University, Biological of Department 76203,Denton, TX ofTexas,North University Sciences, Biological of Department Israel6997801, TelAviv University, TelofAviv Zoology, School Department of Biology, National and Kapodistrian University of Athens, Athens 157- Athens of Athens, University andKapodistrian Biology, ofNational Department School of Natural Resources and Environment, University of Michigan, Ann Arbor, Arbor, Ann Michigan,of Universityand Environment, Resources of Natural School Department of Organismic and Evolutionary Biology, Harvard University, HarvardBiology, andofEvolutionary Organismic Department
Author Manuscript
6 ([email protected]), Johannes Foufopoulos Johannes ([email protected]),
1 body size, competition, island area, isolation, lizards, predation, predation, lizards, isolation, area, island competition, size, body
([email protected]), Rachel SchwarzRachel ([email protected]), 1 ([email protected]), Stephanos A. Roussos A.Stephanos ([email protected]), 6 This article isprotected by copyright. All rights reserved. ([email protected]), ShaiMeiri ([email protected]),
Global EcologyandBiogeography
5 ([email protected]), Efstratios Efstratios ([email protected]), 2 ([email protected]) 1 ([email protected]) 1 1 7 ([email protected]), ([email protected]), 3,4
,
1
22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 Page 6of50 Page 7of50
Author Manuscript 73 of references: Number maintext:in the of words Number 300abstract: inthe of words Number [email protected] 6409403; +972 +972640phone:9811, fax: Israel; Aviv6997801, Tel Author: Corresponding title running Short This article isprotected by copyright. All rights reserved. : Body size evolution in insular reptiles reptiles insular in evolution size Body
Yuval Itescu, Department of Zoology, Tel AvivUniversity, Zoology, of Department Itescu,Yuval Global EcologyandBiogeography
5515 5515
2 30 29 28 27 26 25 24 23
predictors to directly assess the effects of island area, isolation (both spatial and isolation islandofarea, effects assess tothedirectly predictors species. single separately,foreach and snakes lizards and for reptiles, all for informative arewhichpredictors determine to meta analyses and models linear used We bodysize evolution. on competition andpredation inter specific temporal), largeof aarray quantifiedalso We 273 islands. species across eightsnake and representingeightlizard individuals 10,000 nearlydata on surveycollect to literature Methods studiedtaxa Major period Time Location conditions.insular similarly to respond species co occurring whether and patterns are, the common whatevolution, body size primarily drive characteristics island which tested andthesame archipelago in co occurring species reptile 16populationsof insular multiple across bodyvariation studied size We evolution. body size promote these conditions thatand islands, remote onsmall, relaxed pressures are predation interspecificand competition that predicts Theory evolution. ofbody sizedrivers patterns andboth poorly understand stillwe research much Aim ABSTRACT Animal body sizes are often remarkably variable across islands, but despitebut islands,variable across often remarkably sizes are body Animal Author Manuscript
We combined field work, museum measurements, and a comprehensive comprehensive aand measurements, work, museum field We combined Aegean Sea islands. SeaAegean
1984 2016. 1984 2016. This article isprotected by copyright. All rights reserved. Reptiles. Reptiles. Global EcologyandBiogeography
3 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 Page 8of50 Page 9of50 particularities interact to drive the course of sizeofevolution. interact the course to drive particularities eachpredictor snakes.While andlizardsbetween andwithindifferfamilies patterns reptile body size evolution on islands. Instead, species’ autecology and island andspecies’ autecology Instead, islands. on evolution body size reptile ofprocesses and patterns over simplify generalizations thatoverarching conclude onislands. We trajectories body sizeevolutionary general of theories contradicts This conclusions Main isolation. with decreases richness,and predator and competitor area, with increases generally size snake we whereas tested, predictors the byofany hardly areaffected lizards group, aAsrecovered. trend was nogeneral leastspecies,one at sizein body influenced Results
Author Manuscript Body size varies with different predictors across the species we studied, andwe speciesstudied, the differentacross predictors varies with sizeBody This article isprotected by copyright. All rights reserved. No factor emerges as a predominant driver of Aegean reptile sizes. reptileof driverAegean predominant aemerges as No factor Global EcologyandBiogeography
4
66 65 64 63 62 61 60 59 58 57 56 55
become smaller and more isolated (e.g., Darlington, 1957; MacArthur & Wilson, &MacArthur Wilson, 1957;Darlington, (e.g.,moreisolated andsmaller become oftenperceivedas islandsis on Trait evolution contradicted.and proposed been havedrivers patterns and and multiple of attention, has received decades populations body sizes. Additionally, he hypothesized that food limitation is the most important important most isthe limitation thatfood he hypothesized Additionally, sizes. body differentattrends different produces but all sizes, important at isequally predation
Authoronbody factors theseecological selection indirectof suggesting islands, on variation bodysize for explanations Alternative themainland. onislandslargeand competition interspecific and islands, on medium sized predation islands,on small agent selection oftheeffect that Heones.hypothesized to small also than animals to large important more isones limitation thatto food and large animals smallthan to important more iscompetition interspecificthat He hypothesized animal. and theislandfocal of thesize with changes theseoffactors of each thatthe effect suggested (1978) Heaney 2005).1982;Lomolino, Melton, 1978; (Case, onislands sizeevolution drive body to Manuscript thought processes commonly are these ecological resourcelimitation, with Together 1982). (Melton, competition intraspecific stronger consequently and densities population promotehigher in turn, which, pressures, predation and competition interspecificrelaxed to experience thought animals are insular Therefore, 1963). islands as species poor becoming depauperate, generally are faunas 2015Insular & Winkler, 2014;ItescuLeisler Meiri, al.,2007; et 2016 cf. Svenning, & Faubry Lomolino, 1973;2005; size (Van body Valen, a totendmedium evolve animals rule”, insular suggesting “island isthe pattern debated The 2008).most et al., Köhler 2005; Lomolino, 1973; Valen, (Vantaxa across consistent and predictable strong, variation insular among bodysize of evolution in pathwaysthe general Identifying INTRODUCTION This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography
5 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 Page 10of50 Page 11of50 becoming smaller on both smaller and larger islands (Heaney, 1978). However, MeiriHowever,1978). (Heaney, larger islands andsmaller on both smaller becoming large mammals thatwhileof decreases mammals small of size thebody predicted smalleron increase size small in that suggested species (2005)Lomolino size. body be completed (Anderson & Handley, 2002). However, accelerated trait evolution on traitevolution accelerated However, 2002). Handley, & (Anderson completed be
Authorto optimum anchangetowards unidirectional aallow toyetpassed not has isolation since time wheresufficient systems bodysize in with associated be to thought is Temporalisolation differ.wouldnot islands farand on sizesclose body 1999), & Ives, SeaAegeanreptiles,Foufopoulos (asofthe case in unlikely is dispersal over water where that, suggested (2002) Handley &Anderson 1979). Arnold, 1978; (Heaney, pressures and competitionpredation as such factors reflecting indirect, possiblysize ison body isolation spatial of The2000).effect (Heaney, likely more makes rates gene (i.e., and flow) immigration reduces al.,2010), et Handley,Meik &2002, (e.g., Anderson nearest mainland thefromdistance island as Manuscript isolation,usuallycalculated Spatialtime. inspacein and defined be both can Isolation island area. to sizeof linearresponse or a for prediction, Heaney's nofor support (2005) al. found et islands, largestareintermediate sized sized on mediumanimals that predicted further He1998). Taper, Marquet & 1982;also (seeMelton, withareaincreasing increases (1978) Heaney (Meiri,2007). lizards (Boback, and 2003), snakes 2006), 2005, al., (Meirietmammals in relationships such for support found however,minor studies, Other islands. on such becomesmaller species largewhile islands isolated more and to animal affect thought characteristics island thearetwomain isolationand Area Palkovacs,2003). 1994; Levins, (Adler& proposed have beentraits,alsohistory life on directselection via size, This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography in situ in adaptations adaptations
6 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91
population densities, and consequently stronger intraspecific competition and strongercompetition intraspecific densities,consequently and population aggressiveness (Pafilis et al., 2009; Donihue et al., 2016), which in turn favors large large in favorsturn which 2016), et al., Donihue et2009; (Pafilis al., aggressiveness Authorhigher promote thought to iscompetitors) also predators interspecific and from (both releaseEcological 1978). etal., seeDunham but 2003; Hasegawa, Case,1978; 1970; Schoener, 1966; (Soulé, sizechanges promotes nicheand shifts allows competition interspecific Relaxedal.,2007). et 1978;Vervust (Heaney, adaptations anti predatory onselectionsize related byrelaxingdirect specieslarge in decrease size and species small in increase drive issize thoughtpredationto Release from2016). et al., Donihue 2010; Meik Meiri,et al., 2007; 2003; Boback, 1989; &Moriguchi, Hasegawa 1966; (cf.Soulé, results providedinconsistent unclear.studies Previousis evolution size reptilebody affectisolation and area island How2008). 2007,Meiri, 2003; Guyer, Manuscript&1978;Boback e.g., Case, (see characteristic clade specific abe to seems mainland theto on islandscompared diminish or tosizesgrow tendbody reptile Whether (Caribbean snakes and geckos ( lizards smallesttheworld's oneand on hand, dragons Komodo giantand tortoises insular forms: extreme sized their for well known are Reptiles species.small stronger in effectshowing a possibly with the later competition, interspecific predation toand apply should area for suggestedisland pattern the therefore and interspecific competition andpredation reflectsarea island suggested, Heaneyas Furthermore, isolation). of highdegrees lowand larger at and isolation, at intermediate smallest arespecies intermediate sized (i.e., that well as isolation possiblyfor is true area prediction island for (1978) Levins,1994),& Heaney's (Adler area island decreasing as changessame phenotypic theto drive isexpected isolation As2015).increased (Aubret, suggested been also islandshas isolated recently
Brookesia This article isprotected by copyright. All rights reserved. Tetracheilostoma chameleons; Hedges & Thomas, 2001; Glaw et al., 2012) and2012) etal., Glaw Thomas,2001; & Hedges chameleons; Global EcologyandBiogeography threadsnakes; Hedges, 2008) on theother. 2008) onHedges, threadsnakes; Sphaerodactylus
7 139 138 137 136 135 134 133 132 131 130 129 128 127 126 125 124 123 122 121 120 119 118 117 116 115 Page 12of50 Page 13of50 potentially enables us to distinguish between patterns driven by the island conditions conditions bytheisland driven betweenpatterns tous enables distinguish potentially
2016), making this system ideal to study size evolution on islands. We aimed to test to aimed Weevolution onislands. size study to ideal thissystem making 2016), Author etal., Itescu size, documented maximalreaches the even some species in (and greatly alsovaries thesepopulationsislands on in reptile Body size composition. faunal isolation inand area, widelyThese islandsvary theSea. Aegean in 273islands from snake species), (eight eight lizardandspecies individuals16reptile of10,000 nearly fordata bodysize encompassing coverage databaseunprecedented ofa assembled We species specific. those thatstudiedare and we design Thisstudy (Meiri al.,2008). species et co occurring such among uniform relatively butremain etc., primaryvegetation,productivity, climate, latitude, as Manuscriptsuchregions different across vary that factors island specific eliminatetoone allows samearchipelago the within co occurringspecies Comparing archipelago. single a within reptile species multiple ofpopulations multiple island agents across selection potentialof multiple quantifyingthe effect questionsdirectly by these approached we Therefore,2011). Pringle, &Losos 2011;2010, Cox, & (cf.Calsbeek species single fordebated even a sometimes influentialmostis is factor Which conditions. insular to respondsimilarly species co occurringwhether and patterns are, common the whatevolution,body size driveprimarily characteristics island which unclear remains It1978).al.,(Dunham genuset or 2010) etal., Meik1966; Soulé, (e.g., species single of a populations across ofpredictors some the effects or 2014), examined et al., Itescu 2007; 2003; Meiri, Guyer,& Boback 1978; Case, (e.g., populations or species island and mainland studieseither contrasted evolution reptilebodysize date, insular To 2003). Palkovacs, 1982; (Melton, isadvantageous whereearlysizematurity is expected smaller although sizes, This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography
8 163 162 161 160 159 158 157 156 155 154 153 152 151 150 149 148 147 146 145 144 143 142 141 140
Mediodactylus kotschyi Mediodactylus
Koenig in Bonn, Natural History Museum of Crete, Goulandris Natural HistoryNatural Goulandris Crete,of Museum Natural HistoryBonn, in Koenig AuthorAlexander Forschungsmuseum Zoologische München, Staatssammlung (Zoologische collections museum in eight specimens measured further We (1984 2016). years 33 periodicallyover spring summer duringandthe field specimens in measured We collection Data 2008). etal., (Valakos islets verysmall even inhabiting thegeckowithSea islands,inhabitAegean species reptile Fifty islands. varyacross greatly resource availability composition and faunal Consequently, 2005). al., et(Fielding areas maquis agricultural and sclerophyllous evergreen ‘phrygana’), Manuscript called (locally scrub Mediterranean patchworkdwarfof aare landscapes the and 2010) Poulakakis, (Lymberakis & diverse are geologicalhistories Their area. in ofsix magnitude orders varyingacross several islands thousand has Sea Aegean The system Study METHODS 2008).etKöhler al., 2005; (Lomolino, taxa across consistent are islands on patterns 4. size Body 1978); (Heaney, size dependent ispredation of not effect the while species large thanstronglymore species small competition affects Interspecific (Heaney,3. 1978); and competition predation of isolation, degrees intermediate and islands largestintermediate sized specieson are Medium 2.sized 2005); Lomolino, 1978; (Heaney, thefocal species for predators competitors and fewer harbor space),and time and(inisolation in increase in area, decrease as islands smaller species become largewhile size in species increase Small 1.hypotheses: several This article isprotected by copyright. All rights reserved. and thelacertid and Global EcologyandBiogeography Podarcis erhardii being most common, common, mostbeing
9 186 185 184 183 182 181 180 179 178 177 176 175 174 173 172 171 170 169 168 167 166 165 164 Page 14of50 Page 15of50 p=0.23; p=0.23;
of isolation for islands isolated during the past 20,000 years (since the end of thelast theend (since pastof years 20,000 during the isolated islandsisolation for of AuthorPeriods tools.Earth Google usingwere calculated Distances sinceisolation. time the index, isolationtemporal island alargerand the closest fromthedistance studied also we Ives,Therefore 1999).& (Foufopoulos archipelago AegeaninSea the especially (Itescu,2017), isolation effective quantify notadequately maythis index systems island land bridge isolation.for spatial However, ofindex as mainlandan nearest the testedfrom thedistance Meirial.,2005),we 2005;et Lomolino, 2003; Boback, included; e.g., bodysize of (studies biogeographyislandstudies most Following analyzes. further inall field data literatureand museum,pooled therefore We Manuscripterhardii ( sourcesbetween differencesdata found we no withsufficient species two the For source.data each and females for males five of five at least bodysize data recorded wewhich islands from only compared islands. We several forin thefield measured specimens museum collections and in measured specimens bodysizeof mean the we compared nottheresults,bias did data sources of multiple our use ensure To the literature. fromand thein recorded field datawe wasfrom calculated species each bodymassfor Mean snakes. for length total and (SVL)lizards, for length snout vent reportedsizeindices: commonly themost used 2). We Appendix 1, islands(Table 273 species originating from snake lizard of individualseightand eight 9951adult for size bodysex recordedand We 1 seefora list). sources;Appendix 97 extracted from were (data possible as reptiles islandAegean as many forbodysize information recorded and surveyedliterature we Finally, comprehensively History). ofNatural Museum Yale Peabody and University, HarvardZoologyat of Comparative Museum History, Natural of Museum National French Museum, HistoryNaturalBritish Museum, : field mean SVL=61.05 mm, museum=60.43 mm, n=38 islands, t=1.22, n=38islands, mm, museum=60.43 mm, meanSVL=61.05 : field M.kotschyi This article isprotected by copyright. All rights reserved. : field=43.35 mm, museum=43.68 mm, n=25, t= 1.05, p=0.30). p=0.30). t= 1.05, n=25, mm, museum=43.68 mm, : field=43.35 Global EcologyandBiogeography P. 10
211 210 209 208 207 206 205 204 203 202 201 200 199 198 197 196 195 194 193 192 191 190 189 188 187
population. Predators were defined as all mammals and reptiles likely to prey upon prey to reptiles mammalslikely and as defined all were Predators population. Service Navyof Hydrographic the maps Hellenic bathymetric withregional wasconnected itwhich to landmass thelast andaisland focal between potentially compete with adults of smaller ofspecies. withadults compete potentially ofpredatory counts valuesthecombinedand richness birdpredatory between thecorrelation tested we richnessvalues predatorour affect significantlydid not birds
Author species oflargethat assumingjuveniles snakes),snakes (forother or lizards) (for otherlizards defined as rare.were tooCompetitors themare listfor faunal specific and byisland reptilespredationarthropods, on either on reports since centipedes) scorpions,spiders,(e.g., arthropods predatorypotential quantify not could We our fromdatabase. counts bird to confidentexclude feel therefore (p<0.01). We 0.90 (r)was coefficient correlation 2017).The (Itescu et al., lists birdpredatory reliable reasonablywe have did islandsforwhich 41 across reptiles predatoryand mammal ManuscriptthatTo excluding ensure lists. island specific reliable to itcreate impossible making sites, breeding their hunt awaythemtofrom well allows migration seasonal their and islands across mobility since their excluded wereBirds (Appendix species4). focal of speciesfocal each competitor and predator potential the count to observations our own field and theliteraturebased onlists faunal island specific We assembled 3). (Appendixliterature the from assembled wereLGM the than isolatedearlier islands of isolation Temporal finesuchresolution.a at isolation time estimate accurately can method presumeour donot and acknowledge wevalues didalthough round not we resolution, yearand one to werea calculated Estimations (www.hnhs.gr/geoindex/). resolution fine were depthsdrawn from Maximum 2011). et al., Foufopoulos 1999; Ives,&Foufopoulos (seeLGM the changesince sea levelof charts specific depth themaximum fordata by crossing calculated –were maximum"LGM") glacial This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography 11
236 235 234 233 232 231 230 229 228 227 226 225 224 223 222 221 220 219 218 217 216 215 214 213 212 Page 16of50 Page 17of50 positive relationship where Heaney's prediction holds, since it asserts small species sincespecies itsmallholds, asserts Heaney's where prediction relationship positive
transformed body mass. Here we expected to find a significant negative trend for for negative trend tosignificant a wefind bodymass. expected Here transformed Author against log coefficients correlationof the the value absolute regressedwe classes, size all important isacross predation equallywhile species largethan for species important small for more iscompetition thatinterspecific prediction Heaney'stest To 2005). et Meirial., (seepositiveshowslopes species and zero, large to equal slopes have would species slopes,medium sized bodysize predictor show negative would tofind a expected ofspecies. bodymassWe each log transformed against variables, sixpredictorof andthe bodylength) each body (i.e. betweensize relationship ofthe coefficient the correlation regressed we bodysize thespecies by affected Manuscriptvariablesis predictor of ofthe species and each sizeevolution between relationship thethat prediction (1978) Heaney's testTo curvilinear. islinearor richness) competitor andpredator richness, isolation,time since island, larger closest the fromdistance the distancemainland, area, (islandfrom islandcharacteristics six the ofbodysize ofwith each therelationship whether species each for explored first We difficult. isextremely manner invasive non the in ina sexfield their determiningreliably ofsince sex regardless individuals means. Onlyfor femaleand male bysize mean averaging bodypopulation calculated we size dimorphism, sexualto due sizebiases avoidTopredictor variables.in the variation maximizehence, and islands ofthe rangeaswell as populations of species numberand us maximize theto allowed This approach islands. at least eight sampledacross population per individuals at least ofthree hadwe measurements which reptile species for 16 examined We Analyses This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography Ablepharus kitaibelii Ablepharus kitaibelii weof amean all used 12
260 259 258 257 256 255 254 253 252 251 250 249 248 247 246 245 244 243 242 241 240 239 238 237
package (Laliberté, 2011). Finally, we examined for each species which of the sixof the which species examined weforeach Finally, 2011).(Laliberté, package it iscomparable across since these analyses in predictorvariable the mass as body being the best of a collection of poor models; Arnold, 2010; Mac Nally et al., 2017). MacetNallyal., 2010;Arnold, poorofmodels; collection aof thebest being levelswere significance when pvalues>0.05 (i.e. hadsize with associated poorly
that our general conclusions are robust for using different model selection approaches. approaches. model selection differentusing robustfor are ourconclusions general that Author5toAppendixhighlight in species for models each bestAICc based present the we Nevertheless, conservative. moreproved p value approach maligned themuch Thus, merelyscores thelowest(models with them non informative making estimated), werewhich predictors had oftenscores with the lowest models because the scores) orAICc (AIC selection model criteria for information theAkaike using avoided We quadratic terms. linearand both using(α<0.05), on p values based procedure selection model backward stepwise bya followed species, eachfor regressiontest multiple Manuscripta performedwespecies. thisendTo models bestacross selected the compared and size onislands body predicts its that best themodel comprise characteristics island Rthe 'metacor' wereusing performed Themeta analyses we analyzed. valuesr the standardizeto all species, predictors) for tested (against body size log10 transformed of regressions from thecorrelationcoefficients used separately. group We each for sizes, effect as all species in variable predictor eachsize ofagainst body regressions linear of the (r)of coefficients correlation meta analysis random effect Laird DerSimonian We conducted separately. lizards)(snakes andsuborder each across generalorin studied we reptiles all bodysize across significantly affects characteristic islandwhetherany to explore approach meta analytic took thena We 2016). al.,etFeldmanMeiri, 2013; (Feldman& length body than comparisons specific inter for suitableismore andtherefore bodyshape, driven bybiases eliminating taxa, usedholds. We predictionpredationif the no trend for andcompetition interspecific This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography 13
285 284 283 282 281 280 279 278 277 276 275 274 273 272 271 270 269 268 267 266 265 264 263 262 261 Page 18of50 Page 19of50 predictor body size relationships were non significant, and no predictor was importantwas predictor no non significant, and wererelationships size predictor body butmost species, leastone of bodysizeat with wasthe wecorrelated tested predictors theanalyzedvalues. standardize to species all for predictors) Following model selection, selection, model Following 6).statisticsAppendix in model univariatesee (Fig. 1;modelsfull univariate inwastested area when size between and area, twonegative) positive and (three linearrelationship asignificant showed species only five and significant, area island snake (the reptilespecies 16 ofoutone2).In only (Table within families separately,evenand lizards and snakes within 16 species examined, the across greatlyvaried predictorof each sizes effect modestus and snakes Onlyspecies.the the16studied signs) across trend and variablesincluded ofterms (in models 15 different found 1). We (Table all species for ofEachspecies. the highlyacross inconsistent bodysizewere for models best The RESULTS tested(against bodysize log10 transformed ofregressions usedvaluesfrom ther we variabledependent the(r)as coefficients we correlation analyzed where allcases In all analyses. normalityintest) Shapiro Wilk transformation; normalbefore not were 16species ofsix of the distributions(residual size bodyneeded, where and isolation, timesinceislandarea, log10 transformed weparametric oftests assumptions the meet Tosamemodel. in 5)the ≥ factor (variance inflation others with co linear highly werewhich predictors discarded We predictor. tested per points islands) (i.e., data threeof amaximum include to modelsby limiting over parameterization avoided We
Authorammodytes Vipera Manuscript and and
Natrix natrix Natrix This article isprotected by copyright. All rights reserved. shared a similar model. For two snake species ( species snake twoFor model. sharedsimilar a ) no predictors were significant. Explanatory powerand Explanatory significant. were )predictors no Global EcologyandBiogeography T.fallax Telescopus fallax Telescopus still showed the same quadratic pattern, pattern, and quadratic thesame stillshowed ) was a quadratic model of bodyofmodelfor quadratic size was )a Elaphe quatuorlineata Elaphe Eirenis Eirenis 14
309 308 307 306 305 304 303 302 301 300 299 298 297 296 295 294 293 292 291 290 289 288 287 286
predictor, negatively correlated with body size in three species and positively so in in so positively species threeinand bodysize with correlated negatively predictor, significantfrequently most thewas distance the mainland snakes, from For predictor. p=0.23. R p=0.23.
DISCUSSION DISCUSSION Authorpattern).thelizard of opposite the (i.e., island larger theclosest was thedistance from size body snake affect significantly did not that variable only predictor The sinceisolation. with time and mainland the fromthedistancedeclined with and richness, predator and competitor with well as island as area, withincreased significantly however,bodysize, Snake that). weakatone (and a lizards on (positive) effect generalhave a to seemed island larger theclosest fromdistanceoverall. the reptiles bodyOnly in on size effect significant a wehad of tested thepredictors none revealedthat4) (Table meta analyses The ManuscriptR p=0.40. (slope=0.06±0.07, predation bothof the importance showed that (n=16) bodymass against regressions richness length competitor body and richness length predator bodyfrom coefficients valuesof the correlation the absolute 2). Regressing 3, (Table Fig. non significant were inthisanalysis richness competitor andrichnesspredator area,island for results The since isolation). time largerandtheisland closest from distance mainland, the (distance from indices threeisolation thefornegative relationship asignificant found we mass,body its log against species study each in variable eachpredictor against bodyregressionof size (r)fromthe coefficient thecorrelation regressedwe When one. significant most wasthefrequently island area lizards for fact, In in four). negative species and(positivetwo onlyin in lizards emerged relationships linear significant
2 =0.10) for body size variation is not size dependent (Fig.3). not issize dependent variation bodysize for =0.10) This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography 2 =0.05) and interspecific competition (slope=0.08±0.06, (slope=0.08±0.06, competition interspecific =0.05) and 15
333 332 331 330 329 328 327 326 325 324 323 322 321 320 319 318 317 316 315 314 313 312 311 310 Page 20of50 Page 21of50 patterns emerged, except that most predictors were uninformative for most species (a (a species most for wereuninformative predictorsmost that except emerged, patterns between body size and island area, Boback’s (2003) meta analysis revealed no revealed meta analysis (2003) area, Boback’s island bodysize and between
snakes in general support the findings of Meik et al. (2010, 2012), but for most most for but2012),(2010, ofet al. Meik the supportfindings ingeneral snakes Author( rattlesnakes inspeckled correlation positive a2012)strong(2010, found et al. Meik them,between and correlation negative correlation found (1989)a Moriguchi & snakes, Hasegawa For inconsistent. often are resultstested, directly arereptile bodysize on and ofisolation area effects the when However,2005).2004;Lomolino, Ziv, & (Melton,1982; Filin intensify inhabitants their traitsofon the environment theinsularof effects the isolated, moreandgetsmaller islandsthatas biogeographyis inisland perception common A increases Manuscriptisolationtemporal geographicasand smaller become totend species large while islands,isolated more largeron become totend species reptile said,small That snakes. in and lizards those in conditions to environmental on size ofislands body response thein strikingdifferences revealedalso results Our reptile body onsize. insular effect nogenerallikely have richness competitor and richness, predator area,Island species). across inconsistentare (i.e., patterns direction particular a in reptiles studied the ofbodysize drives insular significantly indices isolation ofnonethethree that showed sizesof effect the thismeta analysis withofInfinding, sorts).line consistency fewconsistent ofVery factors. different combination bya were affected species allnearly showedthat species across best ofmodels thea andcomparison species, of the majority even affected variablesconsistently of thepredictor None significance. and (linear, curvilinear) shape sign,varyingin variables predictor the of effects with studied,we thespecies diversityamong pattern great found studied. We we factors theto respond sizes theirway in idiosyncrasy the greatshow reptiles island Aegean This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography Crotalus mitchellii Crotalus ). Our meta analysis results ).Ourfor meta analysis 16
358 357 356 355 354 353 352 351 350 349 348 347 346 345 344 343 342 341 340 339 338 337 336 335 334
patterns in the predicted way. Moreover, isolation, regardless of the index tested,theindex of isolation,Moreover,regardless predicted way. in the patterns surprising somewhat iseightspecies) of (for sevenout lizards ofmodels best proved false. Predator and competitor richness did not affect body size of most of ourofmostof size notbodyaffectdid richness competitor and Predator false. proved biased.
the biotic interactions directly would result in stronger patterns. This, however, however, This, stronger patterns. in wouldresult directly interactions biotic the Authorof the testing effect thatwe2005), expected (Lomolino, size withbody associations significant found andbiotic effects for proxies as isolation area islandand used haveothersAs interactions. biotic of weak effects only wefound few, Surprisingly, size isnot island on variation forsizepredation of thattheimportance prediction his butspecies,largesupports strongly than species smallmore influences competition thatinterspecific prediction(1978) refuteHeaney's results also Our archipelago. theSea Aegean leastin on islands, at reptilebodyofvariation patterns size on overallhas noeffect 2005), Lomolino, 1978; theory(Heaney, to in contrast however, Manuscript area, Island 2005). Lomolino, 1973;(VanValen, size evolution body reptile driving in ofisolation islandrole highlight aresults These patterns. theopposite shows bodysize reptile drives examined we theisland ofcharacteristics that none clear Itis of inreptiles. case the predictions (1978) Heaney's ofmost support cannot We trend.nogeneral to species points the relationship across changes of direction thethatthe fact However,2016). see butet al., 2007,Donihue Meiri, 2004; al., 1978; etet al.,LososDunham 1966; (Soulè, studies previousof results considering theinofisland area inclusion The frequent relationship. negativewe a found whereas isolation or geographic temporal size either and betweenno association found (2010) al. and et (2003) Meik sizeBoback ishump shaped). body with relationship butthe is, (for to model selection bodyofaccording size predictor important particularly ais notarea island species)oftheouteight (seven studiedwe species This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography Telescopus fallax Telescopus it 17
383 382 381 380 379 378 377 376 375 374 373 372 371 370 369 368 367 366 365 364 363 362 361 360 359 Page 22of50 Page 23of50 predation pressure index in the literature (e.g., Pérez Mellado et al., 1997; Cooper et et 1997; Cooper et Pérez Melladoal., literature(e.g., the in pressureindex predation body size evolution on islands are intraspecific competition and resourcelimitation and competition intraspecific onislandsare evolution size body
Author2008).2007,Meiri, 2003; 2003; Boback, Hasegawa, Moriguchi, &1989; Hasegawa 1987; Shine, (Schwaner,1985; snakes forwell, asespecially onislands sizes body withincreased reptiles associated in usually are preyavailability and diversity, preyof inpreysize, terms resources, Richer 1993). &Schwaner, see Case but 2013, etal.,1978; Keehn (Case,islands sizeson withlarger toassociated be thought also is ofstudy someourspecies, characterizes and conspecifics resources against defending involveswhich Territoriality, 2009).al.,(Pafilis et pressureislow predation where competition, intra specific ofstronger large under conditions asize to evolveneed the byhas explainedbeen on islands gigantism example,For 1982). 1978; Melton, (Case, Manuscript to shape oftenthought this arenot intest but didstudy thatwe importantfactors Two islands. thanpredator poor, ratherpredator free, on only apparent are changes evolutionary issignificant that possibility Another 2014). (Meirietal., intensity competition ofindex richness abe weak may Similarly,competitor 2017). et al., Itescu 2005; (Meiri et al., intensitypredation poorlyrichnessreflects predator 2004), al., use aas suspectits common that despite we Therefore, competitors\predators. twodominant orone pressurethan selection astrongerimpose necessarily donot predators and competitors possibly inefficient, many, thatmaybe implies richness predatorof and thecompetitor weak effect think We isgreater. richness competitor andpredator largersizeswhere tendencytowards showeda meta analysis agroup),the generalas (i.e., inreptiles or for notlizards for but snakes, For model). bestintheir included factors hadthesespecies of (less the thanquarter a species study This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography 18
407 406 405 404 403 402 401 400 399 398 397 396 395 394 393 392 391 390 389 388 387 386 385 384
predators are rare, foraging and basking times may increase, thereby allowing allowing thereby mayincrease, times basking foragingand rare, are predators founder effects. orresponse) genetic adaptive (in to contrast responses phenotypic populations is that where adult mortality rates are low (e.g., where predation is low) wherepredation are rates(e.g., mortality low where adult isthat populations se per islands (Kolbe et al., 2012). Thus body size variation is not necessarily orsolely necessarily not variation is body size Thus2012). et(Kolbeal., islands youngandremote small,in patterns, especially bodysize in shaping role have a also may effects Founder &2008). HasegawaMori, 1989;(King,lives their throughout growreptiles longerandsurvive individuals because attained are adultsizes larger variation across for size body non adaptive explanation alternative Another underselection. truetrait gape size with the being spandrel, an evolutionary than more no inis variationsnakes size that body proposed (2009) al. et Vincent 2012). (Aubret, bodysize variation population level drive may pressure) also predation size biased ratheronthanbodysize growth rates on directselection Ofgrowth.course, enhanced where Additionally, 1993).Shine, Madsen 1991; & 1976; Forsman, (Case, differences size body non adaptive mayproduceresource availability in variation from resulting populations insular variability growth rateacross plastic that indications Thereare growth).faster permits food more (e.g., changesor plastic adaptation) (i.e., genetic changes of varyresult as a ratescan growth individual example, For patterns. reptilebodysize on have agents selection suggested commonly effect minor thesystem, ourconsideringinstudy relevant seem especially possibilities These plastichabitat drivenas such possibilities out alternative rule cannot one adaptations, underof a framework iscommonly onvariationexamined islands bodysize While 2009). etal.,(Pafilissizes largerselectingbody for consequently reptiles, insular drive in cannibalism may resource limitation Furthermore,
Author or intensity limitation, differencescompetition in food ontogenetic to (e.g., due Manuscript This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography 19
431 430 429 428 427 426 425 424 423 422 421 420 419 418 417 416 415 414 413 412 411 410 409 408 Page 24of50 Page 25of50 patterns found for snakes follow the common predictions regarding insular evolution evolution insular regarding predictions thecommon snakes follow found for patterns island. larger theclosest distancefrom from apart factorsmost respondsto size body
eggs), or that the variation of one group is adaptive, while that of the other is led by theotheris led of that groupwhile isadaptive, oneof thatthevariation or eggs), Authortheir lizards and uponpossibly prey snakes,but quality food offer for which islands, of onratsremote theabsence (e.g., differently snakes lizards and that affect factors ofor combination reflects factor ofanother isolation theeffect thateither speculate we However, other.the thanstrongly moregroupone indices affect isolation different reasons discern the pointwe thiscannot At islands.remote on more larger grow to species) smaller (i.e., lizards tendencyofgenerally by driven probablythe is largerisland theclosest thedistancefrom trendfor the negative In contrast, increase. factorsas dwarfism these largertowards species) (i.e., generally of snakes tendency bythe driven andare isolation time since themainland distanceof from the cases thetrendin thenegative Likely, complex. somewhat bodyis size onreptile isolation Manuscriptof thattheeffectoverall suggest theseparately, groupsof two each found for we part. Thethemostpatterns 1978), in species (Heaney, predictedsmallfor patterns the notlizards do follow of however,those Interestingly, 1978). (Heaney, specieslarge in Generally, sincethe withisolation. andtime themainland fromdistance with decline richness,predator to and richnessand competitor withincreaseisland area, totends It snake In contrast, closestisland. thewithfromlarger declinetodistance tendency weakfor except a island toconditions, general responseno shows bodysize Lizard each other. topicture opposite anexactlyshow they we found Instudied. fact, we respond to the factors sizes howtheir in differ level, markedly thesub order at snakes, is lizards resultsand that fromour emerged that pattern important and interesting An near patterns future. the in non adaptive and between adaptive discrimination 2017)better may Raia,allow & (Diniz Filho approaches noveland adaptive, This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography 20
456 455 454 453 452 451 450 449 448 447 446 445 444 443 442 441 440 439 438 437 436 435 434 433 432
preys upon its prey species wherever they co exist, and that we can correctly identify we identify correctlycan that co exist, theyand species wherever prey upon its preys species predator thata(e.g., someassumptions general richness)necessitates predator not allow did species thestudied available for currently phylogenies population level body sizes because for some islands we only had data from few individuals. Weindividuals. few onlyfromdata wehad islands some for because sizes body levendis Podarcis strauchi, Blanus schweizeri, Macrovipera
Authormean population level assessing of accurately theexpensecome at may however, This, values). in predictor variance the maximizing also species (thus each for islands thenumberof maximize to aimedwe therefore population and the was analysis ofmain unit island.Our per sampled individuals islandsand sampledof numbers the between trade off exists inherentthe rest,an For beproperly analyzed. to islands speciesSeveral(e.g., species. certainfor some populations and infor samples short quitewe came still reptiles examined 10,000 adult samplesizes.Withalmost small is to face, challenging mostdrawback, and theimportant mostPerhaps results. our biasednot havethey we islands, areconfident species and regarding all assumptions Manuscriptthesekept we Sinceconsistently predators). andimportant competitors all richness, (e.g., competitor ourofpredictors of some Thenature patterns. evolutionary of investigation generallyfurther facilitate will which studies, molecular level island importance of weTherefore,the highlight thisaspect. robustlyexamine to us the However, non phylogenetic models. and phylogenetic bycomparing addressed beto needs onvariationislands size body shaping in condition sharedancestral of noted.beThe role shouldseveral andcaveats methodologicalchallenges and logistic alwayswillonecreate thiswideas scope ofa such study a comparative Conducting 2007lizards). for Meiri, andsnakes for 2003 Guyer, & Boback rule,island(e.g.,cf. the beenbefore haveshown islands onevolution body size snake oflizard Inconsistentpatterns and isnon adaptive. signalsor founder effect strong This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography ) simply occur on too occur simply few ) on 21
481 480 479 478 477 476 475 474 473 472 471 470 469 468 467 466 465 464 463 462 461 460 459 458 457 Page 26of50 Page 27of50 because of any theoretical views about its merit, but because we prefer to err on on the errprefertobecause webutmerit, its viewsabout theoretical of any because 5). (Appendix We uninformative as rejected method based thatthep value predictors pOften, <(at0.05). predictors onlysignificant with thethanbest model predictors
side of caution. side Authorvaluesnot p use thus We islands, important. actuallyare on sizeevolution animal with relationship tenuous mosthave the that thatvariables, many madehaveinfer us would liberal more thatiseven methodselection model aUsing2017). et al., Benjamin (e.g.,2013; Johnson, enough isnot 0.05conservative p= anything, if that, think –including predictors AICcmore thelowestcontained with models the however, fewer lowestcontain AICc the with models, groupof undistinguishable or model, the didno 5).In Appendix case in thoseTable2to in models one (cf.based p value the Manuscript than lessconservative far simplyprovesmethod theAICcIn cases many cases). these in bad(i.e. equally goodequally thatare predictors thatmultiple arethere show only ofsize body in variation the weexplains tested the of none predictors that inferusing p values we example, For inference. biological problematic forbe can alonescores onAIC relying thus and pooroverall often are AICc) models AIC(orprefer. toThebest approach selection ofmodel which issue thestatistical for example a providecompelling results Our ourresults. from conclusions valid robustlydraw to usthisenables think reptiles.Wefor leastscope, at unprecedented bothof speciesare of these sampledwithin each populationswe insular of number the weand species numberofexamined the Moreover, valid.are conclusions ourconfidence have wespecies,thus the best sampled few even across apparent speciesare found we across inconsistencies patternsandgeneral the Yet, cases. some inresults our havemay affected powerstatistical that low possibility the acknowledge This article isprotected by copyright. All rights reserved. Natrix natrix Natrix Global EcologyandBiogeography and Eirenis modestus Eirenis . Using AICc, we could could AICc, we Using . 22
506 505 504 503 502 501 500 499 498 497 496 495 494 493 492 491 490 489 488 487 486 485 484 483 482
predictors of animal trait evolution, and into the need to quantify relevant ecological ecological quantifyneedrelevant tothe into of and trait animalevolution, predictors isprobablyisland bodyon size evolution that thus conclude We reached. being study inthis least at island byconditions, identity than speciesby more affected be
collections they manage. We thank the editor Dr. Ana Santos, the reviewer Dr. J.A.F.Dr.reviewer Santos,the editorAna Dr. thethank manage.they We collections intheworkingwhile Crete), MuseumHistory (Natural Papadimitrakis Manolis AuthorLymberakisandPetros Koenig, Bonn), Alexander Forschungsmuseum (Zoologische Böhme and Butt, UrsulaWolfgang Dennis Rödder, München), Staatssammlung (Zoologische Glaw Franzen,and Frank Michael of the assistance appreciate deeply We work.help with field all their for Werheimer Keren Or and Vezyrakis AlexandrosTamar, KarinTallowin, Sviri, Oliver Sharoni, Achiad Eylon Savion, GayaMosses, MichaelErez Lewin,Maza, Amir Levinkind, Shani Kurzits, Menachem LamaKhoury, SimonJamison, Hartmann, Anat Goshen, Golan,Eden Efrat Gainsbury, AlisonFreud, NirTamar Einav, Brown, Stav AvivBejerano, thank We ACKNOWLEDGMENTS Manuscript systems. study differentfor effects monotonous as island characteristics from paradigm, away evolutiontheinsular of majorre thinking afor theneed elucidates study This size evolution. of processes and patterns the complex over simplify generalizations and island specific, and species generalizations from did,prevent tounjustified weas simultaneously, mechanisms drivingpotential severaltestingof the importance highlightsspecies, across greatly differedmodelsthe best that and thespecies,of majority aaffected consistently tested wethe predictors of that none Thefact taxa. different across general are islands bodyofon size patterns thatevolutionary theoryasserting the This contradicts system. to seemsvariation bodysize reptileThus, environment. insulartothe idiosyncratically respond speciesdifferentlarge aextentthat to indicate thisof study results The This article isprotected by copyright. All rights reserved. Global EcologyandBiogeography 23
531 530 529 528 527 526 525 524 523 522 521 520 519 518 517 516 515 514 513 512 511 510 509 508 507 Page 28of50 Page 29of50 processes driving the evolution of morphological and life history trait variation ofvariation traitlifehistory and of morphological evolution drivingthe processes REFERENCES REFERENCES Authorfauna.insularon of theecology environment insulartheof implications bythe fascinated alsoHe is animalsspecifically.and general insularin animals theand patterns theinbiogeographic interested candidate isaPhD Itescu Yuval BIOSKETCH 2.in inAppendix isavailable thisstudy analyzedthedata All ACCESSIBILITY DATA Manuscript6. Appendix 5. Appendix 4. Appendix 3. Appendix 2. Appendix 1. Appendix INFORMATION SUPPORTING declare. interests to ofhave noconflict The authors #1005/12. (ISF) grant Foundation Science Israel byan study funded Thewas (Permit 20305/824). Change Climate andEnergy, Environment, of Greek Ministrythe from withpermission conducted was work All of thispaper. improveearlier versions to contributedgreatly that reviews, constructive and providinginsightful for tworeviewers andanonymous Diniz Filho,
A list of potential predator species included in our analyses.ourincluded in species predatorApotentiallist of LGM. before islandsisolated time of isolation Literaturefor sources Full statistics of single predictor linear and curvilinear models. curvilinear models. linear single predictor and Fullof statistics scores. AICc species basedon each for Best models usedanalyses.Data for data. size body listliterature containing forReference This article isprotected by copyright. All rights reserved.
Global EcologyandBiogeography
24
551 550 549 548 547 546 545 544 543 542 541 540 539 538 537 536 535 534 533 532
Naturalist Machery, E., Maxwell, S.E., McCarthy, M., Moore, D., Morgan, S.L., Munafó, M., Munafó, S.L., Morgan, D., M., Moore, McCarthy, S.E., Maxwell, E., Machery, A., R., Little,Lupia, J.,D.,List, Laibson, M.,Kirchler,Jones,M., Jeon,J.H., J.P.A., G.,Ioannidis, Imbens, K.,Imai, D.J., H., Hruschka, Hoijtink,T.H., Ho, L., Held, L.V.,Hedges, J.D., Hadfield, A., Greenwald, D.P.,Green,E., Green, S., Goodman, R., E.I.,Gonzalez, George, M.,Forster,A.P.,Field, Fidler,F., E.,Fehr, C., Efferson, K., A.,Easwaran, Dreber, Z., Dienes, P.,Boeck,De T.D., Cook,Clyde,M., C.D., Chambers, D.,Cesarini,C., Camerer,L., Brown,B., K.A.,Brembs, Bollen, R., Berk, E. J., B.A.,Wagenmakers, Nosek,M., J.O.,Johannesson, D.J.,Berger, Benjamin, snakes. neonate insularsize ofinrates body evolutionary theandIsland F.(2015) colonisation Aubret, size? onselectionbirth of consequence anonadaptive snakes tigerin variations size on islands: are adult evolutionBody size F.(2012) Aubret, Criterion. Information using Akaike's modelselection andparameters Uninformative (2010) T.W. Arnold, 286 adaptations. islandand systematicstortoises: giant (1979)their IndianOcean E.N. Arnold, rate. evolutionaryand selection, sloths:biogeography, in Dwarfisminsular (2002) Handley,R.P.&C.O. Anderson, Biology of Review Quarterly syndromerodentpopulations. in Theisland(1994) R.Levins, & G.H. Adler, Author Manuscript127 145. , , 179
Philosophical Transactions of the Royal Society B: Biological Sciences SocietytheRoyal of Transactions Philosophical
, 756 767., Heredity This article isprotected by copyright. All rights reserved. Journal of WildlifeManagement of Journal , 115 Global EcologyandBiogeography , , 69 , 349 356., , 473 490. 473 490. , Evolution , , 56 , 1045 1058. 1045 1058. , , , 74 , 1175 1178. 1175 1178. , The American American The The , , 25
573 572 571 570 569 568 567 566 565 564 563 562 561 560 559 558 557 556 555 554 553 552 Page 30of50 Page 31of50 predation and competition as agents of selection. selection. of agents as competition and predation populations. Behaviour Rouder,J.,M., Perugini, Pericchi, T.H.,L., Parker,B.,Nyhan, S., Nakagawa, 262 populations. lizard among intensity predation varywith autotomy costlyof effectiveness EaseL.J. (2004)and V.Vitt, & Pérez Mellado, W.E., Cooper, varanid lizards. Australian differences in bodysize Island/mainland(1993) T.D. Schwaner, &T.J. Case, vertebrates. terrestrial intrends bodysize insular for general explanation A (1978) T.J. Case, obesus Sauromalus chuckwalla,of the betweenpopulations size differences Body (1976) T.J. Case, reply.Cox & Calsbeek (2011) R.M. Cox,& R. Calsbeek, ofimportance therelative assessing Experimentally (2010) R.M. Cox,& R. Calsbeek, snakes. sizein body optimal an evidencefor Empirical (2003) & Guyer, C. S.M. Boback, island infrom snakes: evolutionevidence (2003) sizeBodyS.M. Boback, significance.statistical (2017) V.E.Redefine Johnson, & J. Zinman, Young,C., Y.,Xie, R.L., Wolpert, C., S.,D.J., Winship, Watts, Vazire,T., Zandt, Van D.,Tingley,B., Sinclair, T., Sellke, F.D., V.,Schönbrodt, Savalei, J., Rousseau, , 243 255. , Author Manuscript Evolution
, DOI:10.1038/s41562 017 0189 z. DOI:10.1038/s41562 017 0189 z. ,
Ecology Copeia , This article isprotected by copyright. All rights reserved. 57 . , , , Ecology , 345 451., 2003 59 , 1 18. 1 18. , , 81 94., Oecologia Global EcologyandBiogeography , 57 , 313 323. , , , 94 , 102 109. , Nature , 465 Nature , 613 616., , , 475 Journal of Zoology Journal of Nature HumanNature , E3 E3. , 26
, , 594 593 592 591 590 589 588 587 586 585 584 583 582 581 580 579 578 577 576 575 574
Biological Sciences Biological in evolution size body Ecology andBiogeography Ecology Gardens, Kew, London, Surrey, UK. London, Surrey, Kew, Gardens, (2005) Mathew, B. &N.J. Turland, J., Fielding, thetuatara. and snakes,lizards,amphisbaenians ratesof diversification and sizesBody S.(2016) Meiri,I. R.A., Mayrose, & N., Pyron, A.,Sabath, Feldman, Society theLinnean of snakes. in allometry Meiri,Length–mass S.(2013) & A. Feldman, 187 EvolutionaryPattern. aObscured Have Basic Extinctions Human Driven Rule: theIsland ofResurrection Svenning,(2016) J. C.& S. Faurby, Tale. Cautionary A Island in SizeLizards:Body J.W.(1978) D.W.Gibbons, & Tinkle, A.E., Dunham, lizard. island of an ecology onthefunctional intraspecific aggression and ofavailability impact food the testing or Feedfight: (2016) A. Herrel, J. & Foufopoulos, Brock,K.M., C.M., Donihue, brain– andgenetics quantitativeRule, (2017)Island P.Raia, J.A.F. & Diniz Filho, York.New Sons, & Wiley (1957)P.J. Darlington,
Author Manuscript812 820. , Ecology This article isprotected by copyright. All rights reserved. , 284 , Homo floresiensis Homo 108 Zoogeography: the geographical distribution of animalsof distribution thegeographical Zoogeography: , 20171065. , , 59 , 161 172., , , Global EcologyandBiogeography Functional Ecology Functional 25 , 1230 1238. 1230 1238. , , 187 197. 187 197. , . Proceedings of the Royal SocietyB:Royal of the Proceedings , Flowers of FlowersCrete 30 , 566 575., The American Naturalist TheAmerican , Royal Botanic Royal , Biological Journal Journal Biological Global Global . John . , , 27
614 613 612 611 610 609 608 607 606 605 604 603 602 601 600 599 598 597 596 595 Page 32of50 Page 33of50 Research AuthorIzuIslands.on the thesnakes of historytraitslife and size habits,body in food Geographicvariation H.(1989) Moriguchi, & M. Hasegawa, Herpetology of study. long term field a from gigantic? Evidence be to longer orlive snakegrowfaster insular gigantic Doesa (2008) Mori,A.& M. Hasegawa, IzuIslands. the of snakes on thelizard and the studies of areviewreptiles: insularterrestrial of diversification Ecological (2003) M. Hasegawa, Manuscript( chameleons leafof species new microendemic andminiaturized discoveryreptiles: of smallest world'sthe Rivaling (2012) Vences,T.M.& Townsend,J.,M. F., Köhler, Glaw, 119 129. islands. of Mediterranean ratesreptiles from extinction elevated and Climate (2011)change A.R. Ives, A.M. & Kilpatrick,J., Foufopoulos, to extinction. vulnerability andattributes history islands:life land bridge on extinctions Reptile (1999) J. Ives, A.R. & Foufopoulos, 253 267. preysize.of effects populations: adder among body maximumsize dimorphismand sexual size in Variation (1991) A. Forsman, change. body mass and dispersability of loss unifying the theory of evolution: insular New (2004) Ziv,Y.I.& Filin, , 7
, 59 68., Brookesia
, 3 , 145 154., This article isprotected by copyright. All rights reserved. ) from northern Madagascar. Madagascar. ) northern from Global EcologyandBiogeography Evolutionary Ecology Research Ecology Evolutionary The Journal of Animal Ecology Animal ofThe Journal The American Naturalist TheAmerican PlosOne Current Herpetology in East East Herpetology in Current The American Naturalist TheAmerican Global Environmental Environmental Global South American Journal Journal SouthAmerican , , 7 , e31314. e31314. , , , 6 , 115 124. 115 124. , , 153 , , , 1 25., 60 , , 177 , , 28
, , 635 634 633 632 631 630 629 628 627 626 625 624 623 622 621 620 619 618 617 616 615
Zootaxa Asia predation, drives lizard tail loss on islands.ontail loss lizard drives predation, Itescu, Y., Schwarz, R., Meiri, S. & Pafilis, P. (2017) Intra specific competition, not competition, Intra specific P.(2017)Pafilis, & S. Meiri, Y.,R.,Schwarz, Itescu, disagree. Turtles general? rule island Isthe S.(2014) Meiri, & P.C.H. Pritchard, P.,Raia, N.E., Y.,Karraker, Itescu, characteristics island of Y.(2017) Itescu, 173. the Indies. West diminutivelizardfrom new a vertebrates: limitamniote size in the lower(2001) At Thomas, R. S.B.& Hedges, Leptotyphlopidae: (Squamata: threadsnakes of species new two insnakes: size limit Atthelower S.B.(2008) Hedges, 59 74. biogeography. ofmodel island equilibrium the onperspective large scale long term, a disequilibrium: Dynamic(2000) L.R. Heaney, ( squirrel tri colored the mammals: from evidence sizeof insular body Island and (1978) area L.R. Heaney, ofJapan. Society
Author ManuscriptHerpetological The414 432. pp. Goris), and R.C.T. Hikida M.(ed. Matui, by , , 1841
, 1 30. , Evolution on islands: testing the extent and consistency of the of effects and consistency the extent testing islands: Evolutionon This article isprotected by copyright. All rights reserved. Callosciurus prevosti Callosciurus . Ph.D. thesis, Tel Aviv University, Israel. AvivTel thesis, University, Ph.D. . Global EcologyandBiogeography Global Ecology and Biogeography and Ecology Global ) of Southeast Asia. SoutheastAsia. ) of Journal of AnimalEcology of Journal Leptotyphlops Caribbean journal of Scienceof journal Caribbean Global Ecology andEcologyBiogeography Global ) from the Lesser Antilles. theLesserAntilles. from) Evolution , ,
23 , 689 700. , , , 86 , 32 , 66 74. 66 74. , , , 29 44., 37 , 168 , , , 9 29 , ,
654 653 652 651 650 649 648 647 646 645 644 643 642 641 640 639 638 637 636 Page 34of50 Page 35of50 Herpetologica ofSciences Academy National Bulletin Herpetological Author masses. bills and beyond ofwarblers: Evolutionisland H. (2015) Winkler, & B. Leisler, http://CRAN.R project.org/package=metacor (2011) E. Laliberté, Science with lizards. aexperiment field differentiation: despiteadaptive persist effects Founder J.B.(2012) Losos, Spiller,& D.A. T.W., Schoener,M., Leal, J.J., Kolbe, Manuscript&Ecologyin Evolution Trends broken. cannot rulesbe (2008)Island R.W. &Moyà