This article isThis article by protected copyright. All rightsreserved. 10.1111/mec.15189doi: differences to lead between the of thisversion c and Version Record.Please copyediting, paginationbeen throughthe andproofreadingtypesetting, process,may which This article acceptedhas been for publication andundergone fullpeer review buthasnot Phone: +34 954 232340 E Vespucio26,E Avda. Américo Biológica deDoñana,Estación EBD Vanina Tonzo Author forcorrespondence: CyprusNicosia, 2 26 Américo Vespucio 1 Vanina Tonzo agrasshopper in taxonomic designation current supportfor no but revealstructure deepgenetic data Genomic Article :Original type Article DR.JOAQUIN (Orcid ORTEGO ID 0000 : MRS. VANINA (OrcidID TONZO 0000 :
Department ofBiological Sciences, ofUniversity Cyprus; Department Biológica ofIntegrativeEcology,Estación deDoñana (EBD -
Acceptedmail: [email protected] Article 1 , AnnaPapadopoulou
–
41092; Seville, Spain 41092; Seville,
- 41092 Seville, Spain 41092 Seville,
- 2 CSIC, , Joaquín , Joaquín Ortego
- 0002
- 0003
species complex species - 5062
- 2709 1
- 1070)
- 1 Panepistimiou Av. 429X)
- CSIC); Av ite thisarticle as –
2109; da . This article isThis article by protected copyright. All rightsreserved. approache characterizing studyregion.Overall, thisstudyhighlights the theimportance ofintegrative geographic distance variation incompatible withecological/environmental speciation.Our supportlittle phenotypic results distribution variationof thatdonotmatch genetic current with taxonomic designations andis landscape genetic analysesrevealed popu that integrativephylogenetic reconstructions, Bayesian model from coveringof entireknown the the taxa. populations distribution two Coalescent phenotypic obtained (geometric morphometric boundaries, and understand factorshave the that contributedto Mediterranean quantitative phenotypic differences and distinct habitat assoc frommontane grasshoppers Pyreneanregion the wereoriginallydescribedbased that on sourcesother ofevidence.Specifically, wefocuson data by sequencing nor throughput the quantitative integration ofgenomic taxonomytraditional species to that hasled designations are genomic data developments inthe fieldofmolecular ecology diversity.classify Overthe biological decade,technological last T Abstract
Accepted Articleaxonomy
among popu s toidentify taxonomic unitselucidate and the evolutionary history ofspecies. has traditionallyreliedon has morphological andecological traits to interpret and
and changed - montane To current habitats). validate taxonomic species test designations, s
and lations limited connectivity across population abrupt the landscapes the paradigmthe analysis. biodiversity Herewe how of illustrate and a marked genetic and amarked structure lations assignedtothe two taxashowaspatial s ) and and systematics Omocestus antigai genome
- based species species based delimitation, and that is supported neither iations (alpineiations - wide SNPdata(ddRADSeq) genetic genetic have eased genera the advances mostly
divergence, wedivergence, and
and explained by information with vs O. navasi .
conceptual conceptual by high - based based , two tion oftion
This article isThis article by protected copyright. All rightsreserved. evidence could resul preferences morphological or diagnosable changes andadhering tothese two sources of therein2018 andreferences Platnick quantitative phenotypic differences (diagnosable species or phenetic concept; mostpopulations, descriptions species beentraditionally have onqualitative based or species concept;biological Queiroz species alternative concepts andthe different properties/criteria todefine used them unitfundamental that biological conforms remains contro aspecies formally described taxa establishment ofefficient conservationpolicies oftenputfocusonthat of the lists red loss rapidaccelerated ofspecies andthe degradation ofecosystems Wiens, 2007 formationspecies new taxa of result as a newtechnological and conceptual advancesthat biologistsallow to discover Dobzhansky, 1937; generatethat it is at the of core Describing 1. phylogenomicgenetics, inference Keywords:
Accepted ArticleIntroduction , ,
2007) 1981; Sokal at at resolution an unprecedented
biological diversitybiological and understanding the ddRADseq, ) . Theseadvances aretimely, currentbiodiversitycrisis asthe . Giventheinherent todifficulties forreproductive test isolation
( Jackson, Carstens, Morale Jackson, Carstens, Simpson t in an underestimation number species. ofthe biological & geometric morphometrics, (Agapow etal. Crovello Dobzhansky, 1970; ). However,speciationmay necessarily involve not shifts inniche ,
1944 molecular molecular ecology , , species delimitation
1970) ) . decade,this Inthe last ,
2004; Mace
and or distinct ecological traits ( gain newinsights into s, &O'Meara, M ayr and allied and allied disciplines integrative delimitation, species landscape ,
, 1942 2004) ecological and evolutionary processes
. Yet, the definition. Yet,the of the ) field 2017 , especially , especially in allopatric ; Sites ; Sites has experienced arevolution
versial due to the the –
e.g., requires the the requires processes & (Avise etal.
Jones Jones
Marshall –
i.e e.g., ., (i.e., the &
the many leading to leading to
Weisrock Nelson ,
( , 2003; see
1987; d e
& ,
This article isThis article by protected copyright. All rightsreserved. & (Edwards Knowles, Fujitaet 2012;Solís 2014; al., objectively identify independently evolving lineages into an analyticalintegrative framework combining (quantitative diverse sources ofinformation traits and genomic data) to Yang multilocus data Yang approach to alternativetest hypothesesdivergence oflineage 2014) ( methodsof inference newphylogenetic andmodel organisms, removingdata previous exponentially 2000 populations as distinct taxa ecological, phenotypic, genetic, tomovehelped awayfrom decisionsdegreeof subjective onthe differentiation( taxonomic statistical increased the delimitation framework rigorofspecies andhas has Fitzpatrick &Ruane,Guiher, Bryson, Yao,Ingrasci, &McCormack, Klicka, 2019;Burbrink, Bryson, Linkem, Lazcano, discover Leaché phylogenetic/coalescent DeWaard, AcceptedEkblom 2011; &Galindo, Article ) , . In particular, the . Inparticular,the coalescentmodelconsidered excellentmultispecies is an (MCM) . &
2013 In recent years,theadventofIn recent highthroughputDNA sequencing , Burbrink, McGuire, & Moritz, T
y Rannala he wide use ofDNA data he wideuse sequence in taxonomy (
of 2003; , 2012 ; Yang,2015 cryptic species cryptic increased increased our capability scalegenomic toobtain large from data non ( , 2010 Domingos, Lemmon, Colli, &Beheregaray, Lemmon, ) Tautz, Tautz, Arctander, Thomas,Minelli, Vogler, & . In particular, the . Inparticular,the incorp 2011; )
- and identify boundaries amongand identify boundaries recentlydiverged usingspecies ) Leavitt, Johnson,Leavitt, Clair, Goward, &St 2011; based based . More . More recent developments have incorporated thepossibility Fujita et Fujita et al. (Fujita et(Fujita al. and resolutionthe of species delimitation
etc.) thatforis needed different considering lineagesor - related constrain ,
2012 , 2012;
2012; Jones oration ofcoalescenttheory under anintegrative ; Yang &Rannala,2010 van Dijk,Auger,Jaszczyszyn, &Thermes, man - based approachesofspecies delimitation approaches ( t - s, while developmentpromoting further & Lemus,Knowles, &Ané, y taxonomic uncertainties uncertainties taxonomic Weisrock
e.g., ( Knowles &Carstens
Hebert, Cywinska, Ball, & Hebert, & Cywinska,Ball, 2003) e.g., ,
2018
) Pons al. et Niemiller, Near,
ha and applicationof the
technology has 2017; Rannala ; Wiens Servedio, & ve
allowed allowed
, 2015)
( e.g., 2006; e.g., - , 2007; model the
. Blair, Blair, &
& Fujita, of
This article isThis article by protected copyright. All rightsreserved. environmentabetween differentiation the and phenotypicdifferences often have data supported etc.)salinity, divergent naturalselection contrasting under presumed represent to ( description of practices and, ultimately, leadto misguided conservation polic compromise ourability to solid reach conclusions ontheoriginanddynamics ofbiodiversity can in turn, have anegative impact onsubsequent ecological and evolutionary studies, species populationto detect stru Sukumaran Knowles & with species boundaries ofresolution genomic can data lead 2018) delimitationimpact onspecies inferences withsexuallydeal dimorphic traits ( However, integrative speciesdelimitation approaches not are exempt oflimitations Acceptede.g., Article Sukumaran &Knowles recognized
Dowle, Morgan . More importantly, differentiation finestbrought detectionofgenetic the bythe highest
As mentionedAs above,a ( Fujita et 2012;al., Huang, ( Carstens et Carstens et al., 2013;O’Meara,2010; leads toreproductiveleads isolation
species
( taxonomic units are e.g.,
Dowle etal., -
based onthe concurrence ofecologicalandphen Richards l and cases ,
cturing may drive cturing mayartdrive to an , ( 201
Carstens, Pelletier, Reid,Carstens, Pelletier, &Satler, 201
classic classic
( phenotypic 7 e.g., of e 7 )
; . That is, in some. That is, in power situations, the higher ofgenomic data common oftaxonomic is the source controversy , Huang
cological cological speciation &
Lamichhaney al. et taxonomic studies 2014; Jones
Trewick (Solís 2018;
not consistent ofwith patterns to structure apotential confusion ofpopulation genetic ,
2018 divergence aconsequence is ofplasticity or divergent - Isaac, Isaac, Mallet, (Chan et al. Lemus ,
( 2014; ) environmental conditions Rundle &Nosil,2005; . Oneof limitations these is the & Sukumaran
Weisrock et al., et i Jones Jones factual increase in the number ofinferred
, separating , ies and ies inadequate management
, i.e., 2015)
2017;
2015) & &
the the formation of
Mace ,
. However Weisrock
2018 2013; O’Meara, 2010; 2013; O’Meara,2010;
& Noguerales, Cordero, , which can , whichcan have aconsiderable Knowles species basedspecies on , 2004 ) , genome suggest Schluter ,
, in many , in other ) 2018
otypic . Taxonomic inflation, ( , 2017 e.g., - ing new wide ).
elevation, , 2001 challenge These are often
distinctiveness that that the ) . species historical
ecological & ). G cases Ortego link link enetic when to
,
This article isThis article by protected copyright. All rightsreserved. populations mostlyare Accepted ofbothspecies followthewesttoeastorientation ofranges the Romera, Presa,Arnaldos, & 1999 quantitative andphenotypic traits asbodysize forewing such shape speciesThese two were describedoriginally based ondistinct habitat associa and Catalan Pre distributedthrough1908 are Pyrenees,including some the sisterputative speciespair elevationsoften isolatedathigh in different mountain systems Cordero,Navas, Noguerales, &Ortego, species) and complexis a diversified ofnine recently distributed species Peninsula in Iberian the (6 continuum Article allows phenomena and ecology populations experiencing contrasting & Bradburd manyfound in studies association 2014 not led to at selection a fewgenesor incipient speciation ) .
An to identify The grasshopper
(IBE) the reduction of geneflow
intermediate situationalongintermediate the that that between , is not only is not
n 2014) orthwestern species) Africa (3
and is and is might - Coastal ranges important . This
environmental dissimilarity and phenotypic generally deserve tobe protecteddeserve (
( relevant allopatric Sexton, Hangartner &Hoffmann,Sexton, Hangartner 2014;
(Shafer & (Shafer & Wolf subgenus pattern pattern Omocestus genomic ecological interpreted interpreted as ( in the context ofresolving taxonomic uncertainties Poniatowski, Defaut, Llucia . of
; Olmo Omocestus antigai Dreuxius
progressive islands involved local islands in processesadaptation ( environmental e.g.,
antigai
and evolutionary , 2017)
- 2013) Vidal
Mason 2015; &Taylor, speciation
Defaut, 1988(genus ( a signal
Cigliano, Braun,Eades, & 2019 Otte, . Most taxahave. Most (Bolívar Moritz . , genetic and phenotypic
Thus, the studyThus, the of 2002; Puissant
conditions
is distributedhigh at elevations (1450 , continuum of
, 2002
1897) and ongoing localadaptationprocesses processes
populations populations ). -
Pomares, & Fartmann, (Cigliano et(Cigliano al
, allopatric distributions andare has been has been Pyrenees, but Pyrenees, their Shafer &Wolf
2008) and genetic differentiationand genetic is the Omocestus
Omocestus navasi Soria these along thespeciation ( Clemente, García, . The known current
frequently reported in the pre differentiation - evolutionary termed isolation Carrasco , tions and certain ., 2019 Bolívar, 1878) ,
2013; Wang2013; - Pyrenees
that have et et al. , ;
) but
García . The 2009) Bolívar among , it
- also also by
. -
- - ,
This article isThis article by protected copyright. All rightsreserved. phenotypic sequencing (ddRADseq; design determine the of conservation status differentthe taxonomic within units divergenceecological boundaries and thefactors understand 2016b). progressive areas the sizebeing small andhighdegree of fragmentation ofitspopulations, overgrazing human and impacts high theandqualityofcontinuous decline in (Braud extend habitats al.,2016a).Incontrast, et the “ different as the taxonomy com ofthe bellmanninavasi speciesdistinct status(Clemente al.,1999)and et even describe subspecific anew taxon ( synonymize twotaxa the differences morphologicaland subtle variation Pomares is m) with2350 or grasslands associated alpine subalpine above tree the line, whereas endangered” dueto
Accepted Article circumscribed toMediterranean scrubby habita - , and the of decline elevation by grasslands occupied two proper In this study, genomicIn this wegenerated usingrestriction data Overall, t , assessment categories 2002; putative putative information abandonment oftraditional on Olmo - ground Puissant 2008) in Puissant et 2008) the decades(Cigliano last al. his complex makes anexcellent this model system totest species species are the included in this this taxon its extremely estimated small ofoccupancy area (60 plex also has has also importantplex implications for conservation - , Vidal, 2002; Vidal, which will through geometric morphometricanalyses the the Peterson, Weber, Kay, Fisher, Weber, Kay, Fisher, & Hoekstra, Peterson, manag ( Ragge & quality of its
was assessed as was assessedas “vulnerable”, main identifiedwith the threats
(Braud al.,2016a,b) et ement in turn
Poniatowski Poniatowski et al. Reynolds, 1998; that havethat contributedto genomic,phenotypic and practices O. antigai grassland grassland provide necessary provide the specific IUCN red list ofIUCN list red threatened species
have motivated differentauthors to .
ts at elevationsts at lower (700 are supposedles are tobe management practices management et practices (Braud al., alpine alpine habitatsand subalpine to due the Reynolds ,
. 2009) Omocestus navasi
. These strong ecological ,
198 baseline informationbaseline - site , (Bookstein 6
201 ) 2012) - , re associated DNA s affected bydirect 9 - - ). the the validate their 200 km
is classified as and management and Resolving the Resolving - 1600 m) complex and ,
19 obtained 2 92) but ) and the ) and O. navasi
to in some (Llucia to with shed O.
, -
This article isThis article by protected copyright. All rightsreserved. localities navasi specimensBetween 2012and2015,wecollected of 2.1. 2. study the area. populationsdistance among andresistancedistances defined bythe contraryon the shapingfactors adaptationto and processesoflocal different elevational andclimatic gradients are the mai system. we Inparticular, main Rannala, 2010 lineages them integrating while and inferred ofspatial geneticpatterns structure compared taxonomicofthe the status stud phylogenetic relationships among polymorphisms and data (SNP) genomic andphenotypicdata. hypothesisthe thatthe current taxonomic designation on light
AcceptedM Article Population sampling ethods
process und process we
navasi
the taxonomythe and (
(hereafter referred to as populations) Leach employed a
( ) , 8 . Second, we é genetic differentiationgenetic and genetic and phenotypicstructure is
, Fujita, &Bouckaert, 2014; Minin, populations) genomic erlying patterns spatial ofgenomic variation andphenotypic in the study suite of analyz data andquantitativetraitsdata evolutionary history of
tested tested the hypothesis Bayesian modelBayesian and In particular,weusedgenome ed alternative
all all O. navasi whether sampled ied populationsied withphylogenomic reconstructions phenotypic variationphenotypic trait coalescent environmental dissimilarity among populations bellmanni - populations populations based delimitation species approach (Table 1 (Dayrat that O. antigai mostly Omocestus antigai , Solís
-
to identify independently of complex the ecological based method (1 population) (1 population) ; Figure1 ,
2005; Jones
- of twotaxa. the Then,we Lemus explained by explained the geographical
- and wide single nucleotide - ). Samplesrepresent O. navasi
driven divergence et al., et
complex topography of s in the studysystem
is supported by is supported & to infer the the to infer from 16
(7 populations) 2015 Weisrock . First, . First, ; Yang & ; Yang& sampling sampling - evolving , we tested
es, 2018) is
one of the ,
O. or if, or if, ,
n This article isThis article by protected copyright. All rightsreserved. d function the using effects of oflocation, size, and rotation of relative positions the landmarksamong specimens performed Core Team,2017) using 1999) correspond to originally usedto those species thetwoputative distinguish ( pronotum, and forewing clusteringBayesian analyses (seeResults section). ( P most the one or two five phenotypic variation if available) Wea maximum selected of10 2.2. for s stored ( populations entire knownacross the distribution ofthe range S i.e., Clemente al., et istances betweenistances landmarks samples across opulations opulations with availablegeometric AcceptedTE Article geometric morphometric REO Geometric morphometric analyses populations
all species/subspecies) and specimens presenteddamaged( structures . Thecoordinate TPSDIG pecimensml vials with 2 96% ethanol in and preserved
D ISCOVERY generalized Procrustes analyses generalized and used landmarkand used
version sex
(MON, SAR, ERR,QUE,andTUR) 1999 es gpagen package
v
(< in thecomplex 2.2 ersion ersion ; s of landmarks(9
Llucia 3 ovipositor valves(infemales) ZEISSstereomicroscope witha ) could be ) could (Rohlf . Thissuperimpositionmethod minimizes sum the Geomorph 8; Carl MicroscopyZeiss GmbH,Germany). and - Pomares adult adult - , based geometric analyses morphometric
were representative of 2015) genomic
sampled . individuals from each W morphometric data
(Adams , and analyzed - e excluded e excluded
14 landmarks per mappedtrait) were images onthe 2002; Poniatowski etal.
analyses (GPA)
due to low population densities lowpopulation and/or due to (Rohlf & e.g., Otárola separately for eachsex to remove the . from morphometricanalyses geometric We imagesdigital ofhead, took
broken forewings or ovipositor forewingsor broken valves) & as implemented the in for which only afew only for which individuals Slice all all
included - Castillo genetic clusters p , opulation
three three 1990) them ,
2009 , all all tax
taxa basedtaxa onthe literature . Toexamine shape 2013) at −20° C until −20°Cuntil at (5 males and5females(5 males ;
Puissant, 2008Puissant, a to characterize . Specifically, T within thewithin complex (Clemente al. et inferred by inferred hese traitshese R - of
version 3.3.2( version - squared squared because needed needed ) f . we or We ,
. R
This article isThis article by protected copyright. All rightsreserved. Accepted Postlethwait,Cresko & 2011; Catchen, Bassham,Amores Hohenlohe, 201 &Cresko, demultiplexed ( libraries delimitationspecies (seebelow analyses Dreuxius individuals of population populations We obtainedgenomic forasubs data 2.3. matrices ofphenotypic differentiation ( Article speciesdelimitationfor integrative analysesusing phenotypic trait ( assess significance 0.05 asa threshold. significance modelsof determined MANOVA Wilk’s was using as the λteststatisticandα= populations if investigate trait shapesdifferent weresignificantly bet (PC1 andPC2) function analys differences taxabetween andamong populations, we Genomic library preparation and processing e differences betweenand taxa among populations sampled s (PCA) on plotTangentSpace
Peterson et al., (Cigliano etal. , males3 and 3 females whenavailable
of ( e.g., O and pre O. for . e.g.,
femoralis Adams antigai each andused trait the covariance covariance the matrix of
Friedline al.et - processed
, 2012) ,
& 201 O. navasi
. We thescores from retained twofirstprincipal components the Bolívar Otárola 9 ) are
and thatwasused Then, w using
presented inMethodsS1. 1908 -
navasi , Castillo et of94et specimens from representative the16sampl
2019). multivariate analyses of( variance P STACKS
ST e used on e used for details (Table (Table 1) ) ( , and aligned Procrustes shape Procrustes aligned , The first two PCsThe firsttwo
see for below details 2013
version 1.35 1.35 version O. navasibellmanni ; Table 1 ; Table e IBPP , aspecies ; as ). - Friedline Friedline al. et
way analys Details on preparation of the an perform
(Solís outgroup insome phylogenetic and ween andamongtaxa
). Additionally, w ( Raw sequences were Raw sequences - also also ed Catchen, Amores, Hohenlohe, Catchen, Amores,Hohenlohe, Lemus for eachtrait e
s ofvariance (ANOVA) to principal component principal separately for each separately foreach ). belonging subgenus tothe ,
2019
et al., et (~6 individuals per coordinates using the ) MANOVA . Statistical Statistical
2015) e genotyped six six e genotyped were sampled
also and to build ddRADseq )
3 to ) used
and ed
This article isThis article by protected copyright. All rightsreserved. values. independent default two Weruns and conducted ratesmutation setto becalculated by were distributions for (MCMC) generations polymorphic sites. sixteen individual analyses (Chifman implemented Kubatko,2014)as & in al. We 2.5. (PCA) by distance inourstudy 3.6) system Section (see sampled (Bradburd,Coop,&Ralph, populations 2018) ofgeneticstructureafter patterns controlling for thegeographicdistanceseparating the recently the spatialmethoddeveloped implemented the in &Donnelly,Stephens, 2000 structure genetic studyin our system We employed four 2.4. assemblingsequence filtering. anddata assembled
Accepted, Article
2012) estimated Population genetic Population genetic structure Phylogenetic Phylogenetic inference (Jombart, 2008) (Jombart, sampled
are are computationally highlydemanding
as implemented in (the onewith highestnumberthe ofretainedreads using using
species species tree alpha population PYRAD We complementary , sampling every
and . r
a Further details these Further on analyses are inMethodsS3. presented version version n s
SNAPP beta beta s using two using . )
B ; T EAST (ii) he (2, 200; 2,2 3.0.66
analyses analyses resulting datasetresulting contained FASTSTRUCTURE FASTSTRUCTURE
2 approaches to
, including , version version
coalescent 1, (Eaton
000 for 1 S NAPP , 000; 2,20 steps 2.4.3 ,
2014) : ,
P
000 - and (i) ; based AUP
exhaustively exhaustively and , and using , andusing different gamma prior and ( (Bouckaert et al. (Bouckaert et ( classic classic Raj, Stephens, Pritchard, 2014 & , ,
000 . * v for this reason, for this
we Methods S2providesdetailson all
, , giventhestrongsignal ofisolation iv) iv) a principal component analysis methods, 000). Theforward( . evaluated
Markov chain Markov chain
4.0a152 (Swofford, 2002). left left the R STRUCTURE
program ; 16 Figure S1 Figure explore explore
individuals and S re NAPP convergence with maining parameters at ,
2014) we CONSTRUCT
analyses
spatial Monte Carlo version version ) only from ofthe each
and u ) and reverse) and ( selected selected one patterns of S 1.3 (Pritchard, retained VD
to infer to infer Q (Bryant et (Bryant et S UARTETS ) T ; NAPP RACER
(iii) 858
v )
This article isThis article by protected copyright. All rightsreserved. Accepteddelimitation comparing marginalestimates likelihood model. foreach (MLE) of delimitation alternativespecies scenarios inan MSC explicit frameworkand bycalculating phylog alternative competing delimitation hypotheses,including species and phenotypicdata Yang &Rannala number ofindependently lineages evolving We applied 2.6. nonparametric bootstrapping Articleexamined) trees individuals and populations andsetting species tree undertaxa themultispecies coalescent modelandthen (MSC) assemblesquartets a these into (Bouckaert werecontained that in thehighestposterior 95% density (HPD), and maximumobtain trees, credibility different runsusing version Species Species delimitation by exhaustively evaluating all enetic clustering andBayesian analyses ( We the used S 1.6. We removed 10%oftreesburn as VD , multi from the Q
three (Chifman Kubatk & 2010) hypotheses UARTETS - , species
Bayesian coalescent 2010
to visualize theposteriordistributionto visualize of trees.
L entire genomicentire BFD ( OG I
uses SNPuses to data BPP ;
* C BFD included models
OMBINER
; method implemented in Solís * replicates assess support to branch ; Leach - o, 2014) Lemus informed bythe : i) T
version version REE data quartetspossible the hypothesis of asingle species ( é -
et al.,et based based species delimitation approa S infer phylogenetic relationships between of quartets
O. femoralis set et al., et ET . We performed A , two of them based 2.4.1.We used (63,860 NALYSER 2014)
- 2015) in and merged and log tree filesfrom the Table 2 Table
current taxonomycurrent SNAPP
and unlinked unlinked .
as (i.e. every combination of four tipswas
version version a single the the third ). an
S (Leaché etal., to (Leaché 2014)
This thecomparison method allows VD outgroup. T RE Q SNPs) - 2.4.1 toidentifyspecies trees UARTETS UARTETS E species modelspecies andfour only A
one (Chifman &Kubatko,2014).
NNOTATOR Specifically,
as well asthe as well and on moleculardata integrating integrating We constructed s D ENSITREE H analyses analyses including all ches determine to the used 1 ); ii
version version ) the ) the current 100 o
ur molecular version contrast species 1.8.3 to pecies pecies
( BPP 2.2.1
five ;
This article isThis article by protected copyright. All rightsreserved. delimitation analyses than locimoreselected withfiveor sites,asthese variable provide species more powerin dataset population analyses topolog Resultssection). (see g A11, unguidedanalysis) specified estimatingtree user(optionA10,guidedanalysis) or bythe species tree (option a to able test from sampled each see Thesescenario). analyses wereperformed difference is the inMLE twobetween competing (base models scenario ESS obtain statistic (2× > 200.TheBayesfactor(BF)test steps MCMC of 100,000 STRUCTURE corresponding clusteringby phylogenetic and Bayesian analyses ( model considering taxonomy Accepteduide Article Leaché, Leaché, and 2018McElroy, Trinh , d
and usingan
considering asdistinctconsidering speciesthetwomain
BPP analyses we ( ies had e.g., computationally tractable, of two species species of two (same usedfor as version version for
( the subspeciesthe
H Huang 5 to the the on )
(Table the the our single
3.3 canperform delimitationusing species analyses afixed input guide α population plus iterations apre (after iterations , fitted
2018; Noguerales etal. 2018; Noguerales six clusters genetic inferred by value settings of0.3. These weresufficient toensure convergenceand This ustoassessthe allowed species delimitation inferences 2 - (Yang
species model.species ) ( as input trees the as input topologies yielded by . We conduct H less variable less O. navasi bellmanni 2 : BFD O. antigai ,
2015; Yang *
analyses) analyses) two we for asimilar approach only ed apath of8 sampling analysis loci ( individuals fromindividuals outgroup the
- and burn
using amatrixusing withnomissing ( data , included twoindividualsfrom eachsampled & Huang, 2018 and O. navasi 2018; - in of Rannala
as aseparate species used a subsetof fromused a loci ourtotal genomic H potential potential genetic genetic Rancilhac al. et 10,000 Bayesian clustering Bayesian clustering analyses in 4 ); (Leaché ) , ; iii) ; iii) v
2014) ) a ) a
). T ln clusters clusters iterations) model ) and the taxonomy current impact
o evaluate theo evaluate impact BF) was calculated, whereBF) BF was calculated, & . We. Fujita
S including twoindividuals NAPP considering six species,
, (east used both options used both and for that 2019) ( , sampling each 100 H O. f ,
3 2010 - steps each consisting steps each and – alternative : west split)
three speciesthree altern emoralis . Specifically, S ). VD
To make 393 ative Q
UARTETS UARTETS but , that
tree
revealed SNPs to be
); the the
BPP we iv ; ) a ) a This article isThis article by protected copyright. All rightsreserved. values, parameter other S the females (see as geometricmorphometric data(PC1 andPC2)oftraitsmales three for and four traits for al., et offor continuous evolution quantitative underaBrownian traits motion process BPP 10 generations (10 0.30 and0.70 automatic adjustment ofthefinetune parameters, swapping allowing ratesbetweento range and divergence times, We respectively. for fourreplicates combination. ran each We used an and G distribution as in different demographic to scenarioscorresponding d https://github.com/airbugs/Dynastes_delimitation respectively sitesvariable (≥5) randomly thedifferent with selected ofloci the subsets specified minimum number of etal. Noguerales loci subset number ofemployed loci VD Accepted Article and(1, 10) topolog Q version version θ provide delimitationconsiderable powerforspecies analyses UARTETS = 2015) s
The integrative The integrative
of loci ( G (2, 2 ies .We 2.1 ) S τ written by J
(Yang , ection , yielded by = .2 000) and 000) and used fourprior combinations for gamma and distribution same forall settings 200, G ,
ran ran Huang and Knowles
(Rannala (Rannala 20 , us (1, 10), 000 samples), burn a after ,
18) 5 ing IBPP
2.2 2015) 00 IBPP . and ran and ran
). As in τ = two and 1,000 We created
- had BPP analyses analyses separately for sex, each incorporating geneticdata as wel P. Huangandavailable θ &
. We analysisfor each ran 100
approach wasbuiltupontheapproach architectureof oftheearly version G = R
(obtained usingoptionA01; (obtained Yang on (2, 2
G script four for independent replicates scenario. each (1, 10) and (1, 10) BPP our , , 000), where 000), where loci ,
2013; Yang s we BPP delimitationspecies inferences
( ( bpp_convert_Ama_sp.r 2016) ) .
ran threeran independent analysesusing input file input According toprevious t studies, τ
. We considered foursettings: parameter
= in of 100
G θ & (2, 2
and at at (Huang s
Rannala ifferent prior combinations ofgamma from outputfilefrom .loci the , τ 000), ,
000 generations. refer refer to ancestralpopulation the sizes
, ,
000 generations, sampling every Yang, 2015 2018) , θ
2010) = and
G . (2, 2 ( We impactthe of analyzed
e.g., numvarfrombppfile.r and incorporates models ,
we used ) ,
, 000) and 000) and Huang
S h NAPP ese
three random as guide tree These analyses , number
and 2018; (Solís τ
= PYRAD
G (1, 10), (1, 10), s - θ
Lemus of = ,
and s
l
This article isThis article by protected copyright. All rightsreserved. populationeach PCscores the of thetwoPCs, first whichexplained 80the 16 %and % ofthe 2.4 1km) (c. sec resolution. “ Then, weusedthe of 19 the dissimilarity eachof between pair ( using (Silva & Silva,2018) calculating pairwise Similarly, weestimated phenotypic differentiation betweeneachpairofsampled pop pairwise estimated differentiation genetic between eachpairofsampled populations calculating heterogeneity; and (3) isolation ( We potential analyzedthree thatcould factors have sh 2.7. both anathese (25 populations BPP analyses) sites SAR,QUE, ERR, andTUR were performed populations available for with phenotyp only http://biodiversityinformatics.amnh.org/open_source/gdmg P
AcceptedST Article - , 4 after a000 samples), conservative burn ) Landscape genetic Landscape genetic andphenotypic analyses
, BPP and differentiation among populations:geographical distance; (1) (2) environmental (Oksanen GEOGRAPHIC DISTANCE DISTANCE GEOGRAPHIC
F and lyses considering onlyphenotypicdata for sex. probability each Posterior (PP) of present including ST .
We ran values in values IBPP in orderto
et al.,et - day bioclimatic variables available in
model five P IBPP . 201 ST We geographical calculated the distance between sampled populations ARLEQUIN
s individuals/population (rather than two individuals/population asin values forthefirst twoPCsvalues usingthe ofeachtrait
a 7 increase the ofincrease accuracy estimates of was considered wellwas considered when supported PP>0.95. nalyses fornalyses 250 ) )
to perform component analysis aprincipal (PCA) andobtai , MATRIX GENERATOR MATRIX using a using - by
version 3.5 3.5 version sampled - resistance defined bytopographic complexity. random , 000 generations, sampling every10generations000 generations,
populations, weextractedfor each in rda” (Excoffier, Schneider, & Laval, 2005)
subset of500 loci subset of 300
function in function in the version version
, 000 generations. Finally, werepeated000 generations. aped genetic WORLDCLIM
1.2.3 ) . Toestimate environmental ic data( ic the the
with atleastfivevariable R
phenotypic variation of ( package , downloaded at30 arc , downloaded F i.e., ST ) and phenotypic
excluding excluding
R vegan
package package one
We the values version .
MON, ulations ulations Pstat n for -
This article isThis article by protected copyright. All rightsreserved. in multiple statistical tests andcalculate we (eight significant explanatory most adequatemodel for thevariability explaining inthe response only variable, where the model toensurethat model significantthe were significant usingabackward were selected were initially constructedexplanatory withall terms fitted geographical, resistance and env phenotypic differentiation ( with geographical, resistanceandenvironmental and pairwisedistances, (ii) population in cell connection scheme. matrix of ofresistancedistances betweenpairs all populations considering aneight weusedlayer, transformed 30arc to Topographic Mission (SRTM Digital Elevation for fromslope eachm cell a90 digital resolution elevationmodel from NASA Shuttle Radar the for eachpairofpopulationsusingEuclideandistances the between environmental variance,respectively. Afterward, wecalculated environmentaldissimilarity R R Accepted Article applied applied “dist” . (Wang traits We
function in variables , us a
for females for six and males) against 2013) Benjamini ed CIRCUITSCAPE
M
ultiple ultiple to analyze: to analyze: (i) pairwise p
(starting with(starting thesignificant ones)until least all retained terms within no additional - R sec (c. 1 km) resolution usingsec (c.1km) resolution . To investigate the . Toinvestigate the ofrole topographic complexity, wecalculated the terms - Hochberg false discovery rate discoveryHochberg false rate
M . Then, we tested . Then,wetested P atrix ST
version version were ) in relationwithgeneticdifferentiation) in ( - ironmental distances stepwise R
variable egressions with egressions retained.
4.0 the corresponding the ( procedure, byprogressively removing non McRae, McRae 2006; Beier,2007 & reached significance. opulation genetic differentiation (
the the As significance we tested we tested Data) and Data) the raster was final layer the samethe ofindependent subset variables, R QGIS andomization (MMRR) . Models correction q
- version version a high number of values ( i.e., of the for all dependent variablesfor all
full models)
The resul
using theusing 2.8. Then, onthis based to adjust
rejected obtained PCscoresusing F ST t
p
p.adjust ) and ) and terms phenotypic traits phenotypic was - and final and final models values for as implemented F ST )
the minimal to calculatea ) in relation ) in against this -
neighbor function -
This article isThis article by protected copyright. All rightsreserved. included). and SNPs number retained ofunlinked in final dataset the obtained with anaveragecontrol, of2.11 ×10 Illumina sequencingreturned anaverage of2.19×10 3.2. Genomicdata 0.0 F F by shapevariationofforewing (PC2: One Wilk’s= 0. λ = significant differences among sex 65% in variation forpronotumshapes ( The twofirstprincipal component scores,PC1andPC2,trait explainedmost foreach 3. 3.
10 10 1 0.0
Accepted Article Results . 13 , ,
are - minCov 40 42 Geometric morphometric analyses ; way ANOVAs multivariatethat revealed differentiation among was driven populations 0 ; PC2: ♂
2 = 3.08; = 2.25;
shown inFigure ; : 77%) ♂
:
Wilk’s= 0. λ F = 25% was = 25% 64,365 SNPs(63,860 940 10 p p and ovipositor valvesand ovipositor (
, = 0.005; PC2: = 0.03) in females and
36 , F
=
8 , 2
3 . 8 33
= s
047 S2 andS3 ; 0.305 p
= ,
sampled F F 0.0 , 54 6 10 p
reads retainedpersample were (Figure S1). Thetotal , 31 ,
♀: 87%
= 0. 142 40 . ) in males.) in
M = 2.47; = 2.47; . 2 F ♀: 95%) 7 959 ANOVAs based onPCANOVAs based head (PC2: populations 10
= 2 ,
; 40 ; ♂ . ♂
16 = 4.86; p : 76%) :
, SNPs whentheoutgroup .
= 0.021)and Wilk’s= λ p Princip
<
F 0.0 ( , head( p ♀: 10
6 < 0.001) and ovipositor valves< 0.001) ovipositor (PC1: and ,
01
reads sample.per After quality al componental plots for trait each and 37 Wilk’s λ=0.
)
= 2.22; but thetwo notbetween 0. ♀: 90% forewing (PC1 7 45
scores for scores , p F
; = 0.039 6 , ♂ PYRAD
0 35 35 : 91%)
= O. femoralis , all all 2 F ), pronotum (PC1: . :
80 00 for all individuals individuals for all F trait , forewing ( ,
10 , 192 p ,
s
36 . = 0. 50
showed
= taxa ( = 2 092 was 1 . . 75 674 ♀: ). ♀: ; p ,
p =
This article isThis article by protected copyright. All rightsreserved. FASTSTRUCTURE FASTSTRUCTURE exception of CRE more geographical sense populations same they five replicates. explainused to structure data inthe towas equal replicates four intwo and tofive replicate, equal inone and numberthe ofmodel components complexity that maximizes marginal waslikelihood 22 equal to threein replicates, equal to according to the to individuals and populations assignedtoa were singlegen groups populations for western analyses showed ahierarchical genetic structure, firstwith a split populations withtheirtaxonomicbut not (Figure assignment S butcriterion, probabilities ofthe log (LnPr(X|data STRUCTURE 3.3. 4 populations Accepted Articlea ). Theinferredgenetic Population genetic Population genetic structure - also ERR genetic genetic M into into genetic clusters other at smaller from geographical scales presented aximum - QUEfrom TUR
analyses yieldedanalyses an‘optimal’ for value clustering w TOS cluster the ere assigned todifferent assigned ere clusters genetic in putativezones contact between FASTSTRUCTURE
Δ for (Figure analyses , no other population showedanysignature, noother of admixture in genetic ‐ K likelihood scores likelihood some
analyses wi analyses criterion the ‘optimal’ clustering K populations populations BOR ; (ii) ; (ii) differences
= 2 and subsequent subdivision subdivision of = 2andsubsequent these two main groups weregroups consistent distribution geographic with the of the - SET FASTSTRUCTURE
th results mostlyresults mirroredthose - BAT asdoneby BAT S6
K from
). Assignment values ). Assignmentvalues
> 5 were extremely cases> 5were lowinall ( (Figure (Figure FASTSTRUCTURE - MON and BOIMON S6
analyses for analyses ) : (i) For : For (i) genetic genetic four in18replicatesand to equal seven in K STRUCTURE )) steadily)) increased up to solution by clusters K
etic cluster - FASTSTRUCTURE
to LLE analyses analyses
K = 3 additional geneticclustersadditional in
2 > obtained withobtained K
; was for ,
. However, 5 did not split populations5 didnot Figure Figure
STRUCTURE = 2 according tothe= 2according .
between and eastern
K K peaked for
, = 6
= 2( with admixture limited K S
; = 3 to = 3to 5 (iii) ). q Figure S4 Figure , whichmakes these these
STRUCTURE < 0.0001) included STRUCTURE population With the K K K
two pairs of two pairs
= 6. = 6(Figure = 3 b).
Δ and
M in the and,
K Model Model
ost , but
,
This article isThis article by protected copyright. All rightsreserved. external TUR sister placed asa clade of CRE 3 STRUCTURE which inferred same the clusters (Figure genetic that Population trees 3.4. admixedthe population of whereas PC1 location geographic ( of clusters presence twomain the genetic revealed by 2019) al., S S (Figure non similar result, seeBaiz, &CortésTucker, their thus, respective Figure 2 Figure ). For the east geneticgroup, the ). For 7
Accepted Articleb - ).
Phylogenetic Phylogenetic inference only spatial spatial model (Figure S separated separated CONSTRUCT were In agreement with The spatial model of
different forthe relationships group PC2 separated ; 7 differed in the differed the in
and assignments onlygenetic for Figure a
), butlayers beyond and grouped within the same genetic cluster at the withinthesame clusteratthe grouped genetic lower hierarchical level by the westernmostthe populations with respect to to fivewith respect the CRE easternmost populations FASTSTRUCTURE reconstructed S
analyses showed stronggeneticanalyses showed admixture 8 ) and main clusters genetic . bar virtuallybar were plots identical to forthose obtained
the rest ofwestern the rest inferred
previous TOS at an intermediate position. TOS atan 7
CONSTRUCT a ). Spatial statistical models supportupto hadstrong SNAPP in K relationships amongrelationships
s analyses BPP
= 3 1a analyses, analyses, -
, and ERR and
SNAPP populations NER contributed relatively little to littleto contributed relatively total covariance(Figure
showed better showed better model fitateveryvalue of ‐ BPP 4 Ortiz, 2019) - (Figure populations from ). QUE, whereas whereas QUE, K BOR ,
For group westgenetic the three PCA
analyses produced the sameanalyses produced the and topology and = rather than by than rather described species(Figure 3
- - SVDQUARTETS
4 MON rest fromthe of 2 grouped populationsgrouped according totheir
) were and some STRUCTURE .
-
SAR
clump partially SVDQUARTETS SVDQUARTETS nearby nearby
eastern eastern for all - CHI ed
yi
and compatible elded similar topologies - together inPCA (Figure populations belonging to ESP K populations andplaced - ERR
values FASTSTRUCTURE - PER placed TUR as - QUE K populations
(e.g. Whelan(e.g. et
(Figure (Figure 4) = 5 with - SET
methods (for a K K the the
than the the than = 3 - BAT.
3 , , ). - 5
This article isThis article by protected copyright. All rightsreserved. randomizingtest dataacrossindividuals phenotypic yielded support inconsistent phenotypic data IBPP few clustersgenetic (CHI the using (Figure SNAPP 2 ( scenario clustering analyses in ( delimitationSpecies with analyses the 3.5. no support hierarchic among inferred by populations populations tobelonging same the clusters (Figure genetic considerable fuzziness for in terms of As expected H ).
Accepted5 Article BPP ), which considered populations), whichconsidered assigned to each SVDQUARTETS
prior combinations prior and specific g Species Species delimitation analyses support analyses
and
analyses 4 the
). All sampled populations were consistentlysupported asdistinct specieswhen subset of1,000 andonlysome loci populations belonging nearby same tothe H al al structure yieldedgenetic by SVDQUARTETS
( low 1 split i.e. , ) was the one) wasthe results results
the relationshipthe amongth
polyphyly) node
( fully supported the ofacrosspresence 13species alltestedtopologies into differentspeciesinto Figure ) (Figure ) (Figure depending on - support support ed ESP andCRE STRUCTURE in probablydue some flowamong tree togene ofthe parts nearby
all populationsall different as usingeither species femalemale or 5 )
) , prior combinationsfor gammaand number distribution, ofloci considering . 4 for theseparationof the named (Figure taxa I BPP ) . (i.e. posterior probabilities for This evidencedis also all
analyses based onlydataanalyses based onmorphological
sex and as phylogenetic - ERR a uide for one species distinct species (Table BFD e se popula STRUCTURE - the the the the
QUE) were not supported as distinct notsupportedas underQUE) were taxa a * trees
method strongly a method strongly supported sixspecies model 3 tested -
4 analyses for subsets the of200and500loci
most external nodes most external tions alsotions showed considerable unc ,
but received much lower support but received muchsupport (Table lower demographic scenario , genetic cluster by FASTSTRUCTURE
the are in goodagreement with the
1a
SNAPP SNAPP almost result identical an ). 2 ). The second most). Thesecond In general, In general, t and bootstrapping values and bootstrapping analyses
identified by , (Figure
and s
1 ). ,
, whichshowed PCA he relationships but tended to but tended
showed 5 ).
and A validation Bayesian Bayesian -
supported ( Figure 4
ertainty indicate indicate
(Figure (Figure ( BPP , ).
This article isThis article by protected copyright. All rightsreserved. taxa groups ecological/environmental speciation. showOur results oftwo thepresence doesnot (i) match with populations assigned to ouranalysesCollectively, the rejected 4. multiple testing ( differentiation applying genetic after defined distances resistance bytopographic complexity, environmental dissimilarity or was showed that no trait significantly in correlated anyse final model ( withsignificant relationship geneticdifferentiation variablewas andthis from excluded the complexity ( distancesboth geographic ( genetic indicatedthat differentiationMMRR analyses ( 3.6. and across tested topologies different speciesfor morphological dataset highamountby the ofgenomic data. S 9 Discussion Accepted Article ) , indicating influence that the ofphenotypic traits in Landscape genetic Landscape genetic andphenotypic ( Wallis,Waters, Upton corresponding toaneast p p
= 0.004) (Table 3 = 0.004) = 0.374; = 0.374; 3 Table all q some external - s values
(i.e. without g (i.e. without their respectivetheir O. antigai p , & Craw, , &Craw,
and demographic priors (Figure S
> 0.05) > 0.05) (Tab = 0.001) resistancedistances by and defined topographic ). MMRR analyses for phenotypic differentiation). MMRRforphenotypic ( analyses - ). However, environmental dissimilari west split west split analogousthatfound to in
and nodes
Benjamini Remarkably enomic data) also enomic also data) 2016) current current O. navasi taxonomic designation isincompatible and(ii) with , although resultsvaried the strongly analyses le 4 le
and a markedand a differentiation genetic at local ). hypothesis oftwo and indicated species that
- Hochberg false discovery rate discoveryHochberg false rate
show a distribution of geneticvariationshow adistribution that ,
IBPP
F analyses analyses IBPP tended tosupport ST
x with geographical distances, ) was significantly correlated) wassignificantly with
analyses analyses is probablyoverridden 9 ). only
many based onrandomized ty didshow not a
the into split main other Pyreneanother between sexes corrections for
genetic genetic P ST ) This article isThis article by protected copyright. All rightsreserved. do notmatch currenttaxonomic designations and, accordingly, populations sampled inthe structure analyses thatindicate populations split two in our indicate nosupportfor analyses, results the anda phenotypic data comprehensive suiteofphylogenomic andBayesian clustering have agenetic basisof scenario in line witha theirand determine whether specific associationsphenotypic and habitat slight differences now, nomolecular conservation assessment programs et 2016 (Braud al., aretwo taxa validspecies now(Cigliano considered etal. and biometrical analyses by Clemente al.(1999) et resurr and songs), recommended synonymy the of population (1998)Ragge &Reynolds attributedsubtle their morphological diversity and Bolívar ( of differentcriteriaconsidered the to establish the boundaries thetwoputative between taxa and As inmanyspeciesother complexes 4.1. characterizing studyregion the scalesspatial limited reflecting connectivity across population abruptthe landscapes Accepted ArticleClemente al.,1999; et Clemente al. et (1990) O. navasi Integrative species Integrative de
in the transition the between XIX andXXcenturies
geographic isolation,denieddifferences inacoustic communication (
has changed through time withthe conjunction discussionthe in about validity - based analysis hadbeenperformed themonophyly test ofthetwo to taxa Ragge
confirmed their limitation & ( e.g., Reynolds
( Noguerales, Cordero,Noguerales,
e.g.,
, Leaché al., et 2018
taxonomic 1998 O. navasi ecological current current taxonomy ) . Afterwere thedescr twospecies
a to distinction, Reynolds (1986) anddistinction, Reynolds(1986) , main groups genetic b speciat O. antigai ected species the
; ,
Hochkirch Hochkirch al.et 201 & ( Bolívar ) Ortego
9 the the taxonomy of ion. Using genomic and ) and treated as such in treatedassuch ) and . Phylogenomic andgenetic . However, newacoustic , 18 ,
201 9 to effect the of 7; 6 O. navasi , a
2016) ) Bolívar (east/west) .
O. antigai i.e., . Upto ibed by ,
courtship 1908) and the that that
This article isThis article by protected copyright. All rightsreserved. Acceptedpop high powerof numbermarkers alarge of (hundredsof tens to thousands) to fine detect overestimation beparticularlyexacerbated when can vast using genome Leach populations thancapturingrather speciesdivergence fromisolation truespecies boundaries andtend to over such methods not are abletos coalescent model(MSC) (Rannala limitationsthe delimitation ofBayesianspecies approachesbased onthe multi alsoFig.(see 2bin considering the Accordingly FASTSTRUCTURE assigned to thesame groups genetic clusteringby PCABayesian analyses ( and Articleremarkable that many populati nearby assupported distinct km and assupported combinations forgammadistribution fullyanalyses presence supported the of presencethe of nominated taxa, speciesarephylogeneticallyputative inte two distinct habitats (alpine ulation structureulation ( é , Zhu, Ranna Despite phylogenomicDespite recover not reconstructions belong , the species , the delimitati a ing to theing to same cluster genetic distinct distinct species several several highest number ofspecies,
molecular and Leaché la, &la, Yang taxa CONSTRUCT Hime Noguerales etal., 2018; al.,2016; et plausible species.
under under -
et al.,et based vs tatistically distinguish structuredue genetic to population . Mediterranean/montane) supposedly occupied bythetwo and o ,
a few prior combinationsa fewprior and 201 ) were ) were 2018) species delimitation analyses using on bymodel mostsupported & 9 nly . ; Yang, 2013;Yang
ons (BOR I Sukumaran n most cases, . consistently some Recent Recent empirical and have theoretical studies shown 13 rspersed (Figure Independently ofthe
one per geneticcluster inferred by species species s nearby populationsnearby
(CHI - MON, NER & across across - delimited asdistinct species (Figure
ESP and CRE each ( Knowles Huang, 2018; - split genetically differentiated &
ing 1 of the 17 all Rannala, 2010),pointing outthat ).
-
the monophyly ofthetwo tested topologies SAR, BOI specific specific ,
number used, ofloci 2017) BFD Sukumaran separated by less separated byless - ERR * sampled Leaché al., et
. BPP analyses wastheone analyses
The problem of species The problem ofspecies topologies. - - - wide data due the wide data to
LLE and SET QUE) and & population BFD S -
and Knowles species
TRUCTURE STRUCTURE were not were not * 2018; It It is also than prior
identified identified BPP - -
s grain grain BAT) , 20
w 4
). as ,
This article isThis article by protected copyright. All rightsreserved. requirements considered tobe di intersecting ofevidence suggestlines multiple that the delineated entitiesdo notmeet the coalescent diagnosesmodel frequently structure ratherthanspecies, different genetic ameliorate this problem by of emp loci Knowles, 2017 framework prevent cannot overestimation species al. et (Noguerales of delimitationspecies approach implemented in seems little or very individuals yielded populations as distinct species and a test validation randomizing dataacross phenotypic between only ( help could to genomic other with data sour Pyron, Hsieh,Lemmon, Hendry,2016 Lemmon, & distance can bedelimited asdistinctwhen hundredsspecies using ofloci( intraspecific populations presentinggenotypicdifferences simply 2017 Pyron et al.,Pyron et 2016;
Accepteddelineated Article BPP ). pres
to be overriddenbyto be In this line, and Beyond thetaxonomic inflation studiesRecent have currently recognized ented verysubtledifferences loyed are are loyed likely tohave BFD influence
species, species, mitigate ). *
and integrative and integrative In sum,t r an ecent empirical studies Sukumaran 2017 Knowles, & indicating
the problem ofconfounding almost on the .
he stronggenetic structure inour highnumberstudy systemand the the highamount data.the ofgenetic suggested that ces ofinformationces morphological,ethological,etc.) (ecological,
outcome of identical result. identical result. taxa stinct taxaaccording toalternative contemporary species that integrating phenotypic and genetic intodata this analyses
( extraordinarily see among of some pairs Figure issue
IBPP implemented in species delimitation approaches integrating have found
linked withlinked thefactthatmultispecies
Thus, incorporating s analyses analyses IBPP
S2 andS3 ).
inflated inflated the numberidentified ofspecies population ) does Although most traitsanalyzed phenotypic .
that geographically continuous geographicallythat continuous and the not
) IBPP , Our our populations populations notdifferand did really really
structure withstructure speci does not does not seem to tohelp impact ofphenotypic traits
IBPP analyses resulted morphological help to the reduce number ,
analyses supported most 2018 suggest ; Sukumaran & Huang, 2018; from isolation by isolation that the traits es
limits has no This article isThis article by protected copyright. All rightsreserved. do not 3) probably populations. hierarchically structured in congruence with clusteringBayesian analyses supported thatgeneticvariation thewithin complex is Phylogenomic reconstructions andspatial patterns of from structure inferred genetic 4.2. conservation intraspecific with by now 1999; Slatkin, ouropinion,thestudiedshould populations 1993).Thus,in betreated 2010), such mostpattern reflects plausibly with a 2007 of indicative eastern and diverging ar low elevation populations ( environmental studied traits overlappingof largerly populations, phenotypesacross taxa,andfor clusters all genetic (see 1in concepts Table de Queiroz, 2007). Theselines ofevidence (i)theinclude: presence
Accepted Articlee always Factors structuring
) the the . Although rapid even
– possibility of
as a single as a allopatric and explains whythe clearly clearly stepping
However, the strong signal strong However, ofisolation the evolutionary significantevolutionary lack of lack and management (Figure S (Figure - driven divergence and the strong forthe signal isolation recover intrinsic speci - stone speciation process stone speciation process (e.g. VanderWerf,Young, Yeung,&Carlon, western desi 2
and S admixture at contact zoneshappensevenbetweenadmixture twomost atcontact the genomic andgenomic phenotypic variation es gnating
ecological the west/eastthe
taxon with genetic strong structure, whichis
reproductive isolation spatial 3 genetic groups strategies strategies ; phenetic species concept the analyses analyses in units units
the presencethe of + main
phyloge (Braud (Braud et al.,2016a, b). split (Moritz, 2002)that should populatio genetic clusters genetic the geographical distribution ofthe
( revealed byphylogenetic reconstructions - genotypic cluster species concept by ne CONSTRUCT CONSTRUCT - distance might bepotentiallycompatible tic tic species concepts - ( n connectivity by biological species
phylogenetically interspersedhigh and - distance (
);
(ii) Noevidence for (e.g. blur
east most structure genetic and r
= 0.71, (Hutchison & Templeton,(Hutchison & - west be ) ; concept
(iii) Genetic (iii) clusters considered not in p genetic groupsgenetic
< 0.0001; < 0.0001; Table ) compat ), which
sampled ( de Queiroz, in future classic classic –
at least at i
ble ble and is )
as
This article isThis article by protected copyright. All rightsreserved. Wallis al. et highrecolonize elevationsinterglacial during periods & Hugues, andtherange presence likely ofice reconstructionsgeological documenting heavily areas theglaciated in central of part the 2010; Valbuena 2016; ( havebeensplits many reported for plant andanimal other from species P the west ofanimaginary line around located Segre (~1ºE). the river Similar east a equilibrium) artifactual consequence geneticof differentiation gradual and fromgenuine resulted historical processes systemstudy evidence suggest that results from made clustering analyses byclassic ( clustering analyses ( a remarkable congruencebetweenthestructure hierarchical genetic inferred bynon by clusters inferred F other analysesall e.g., SNAPP nalyses nalyses showed abasal intwo split clades corresponding topopulations eastlocated and Accepted Article ST
= 0.35) and extraordinarily high andextraordinarilyhigh = 0.35) drift estimates ofgenetic divergence g after forall ,
Alvarez 2018 Charrier, Dupont, Escaravage, Pornon, & , et al In agreement with In agreement with SVDQUARTETS ) 2011) . ( . merely driven bygeographical distance (Bradburd al.,2018). et Overall, Overall, - and 2016) Presas, Mateos, Vila - Ureña et al.
. Ourstudy ahighdegreeof systempresents differentiation genetic (mean indicate indicate that
STRUCTURE for numerous alpine STRUCTURE these might allowed have montane/alpine species tosurvive glaciations and , and STRUCTURE that
results results BPP ,
2018) the
) and PCA,beinga ) and thelatter methodof free , analyses FASTSTRUCTURE strong genetic revealed by structure all other analyses is CONSTRUCT are - . -
Farre, Sluys,Farre, &Riutort,
(Jombart, Devillard, & 2010)Balloux, free refugia at longitudinal its extremes This
and consistent consistent breakswith the transverse proposed taxa
( east F FASTSTRUCTURE - value >
- and thesupport ideathat
west might split bewell explained bydetailed ( make i.e. 2014; ) and ) and strong strong
0. little biological sense our little in
( 5 Mila, Carranza,Mila, Guillaume, Clobert, & Charrier 2014; et al., both ; Pritchardet al., 2000). Also,
(i.e. migration(i.e. and gene isolation
results for results 2012; phylog Bidegaray ) enetic reconstructions
rather than K the the fragmentation of
= 2, . All these linesof the the Valbuena ( -
- yrenean yrenean region Ehlers, Gibbard, all phylogenetic all west genetic Batista et al., assumptions assumptions being specific tic -
spatial drift
by - there there is an Ureña enetic enetic This article isThis article by protected copyright. All rightsreserved. Mediterranean for eitherdriven ecologically divergence the or taxonomic separationbetween supposedly the a populations local/regionalat scalesfollowed alongitudinal cline rather than asegregation of Cordero, landscapes on structur obtainedto those byprevious studies showing the impactand ofsteepslopes complex studiedtaxon,malesthe brachypterous and being females micropterous, comparableand are topographic (Table3 roughness by boththe geographical among distance populations and resistance distances definedb geneticAccordingly, landscape analysesindicatedthat differentiationgenetic wasexplained geneless increased flowthroughthe abrupt landscape admixture ( ata located similar longitudebut higher elevations showednoevidence of genetic ongoing geneflow between the populations TOS (andinto a lesserextent andPER in SAR),contact the located zonebetween populations,nearby suggest astrong effect ofgeographic isolation been assigned cluster,with toaunique genetic onlyafew casesofgenetic admixture between geneticgroupswest andeast fine spatialscales,at with of also mountain ranges distributions duringancestral glacial and presenceperiods the of along refugia Pleistocene Accepted lpine and Mediterranean Article west and east geneticwest andeast groups montane and Genomic alsorevealed data &
e.g., Ortego O. navasi
BOI, whichsuggestsahigher LLE), isolation populations ofalpine and has ,
2016 alpine organisms. played a key role on the diversificationplayed akeyroleonthe ing variation genetic inmontane/alpine grasshoppers
a and the alpineand the - ) montane populations ( up to . Genetic structure analyses showed split that ofthedifferentthe them (Figure (Figure ). T three genetic clusters comprisedthree clusters genetic each within main ofthetwo
presented signatures ofgenetic admixture
(Figure hese results might the low dispersal reflect capabilityof strong genetic differentiation and 2 ). The fact that almost). Thefactthat all
O. antigai s 1and
e.g., 2). It 2). It isnotew . This was corroborated . Thiswascorroborated by
SET nosupport andBAT), indicating s characterizingthe foothills. Pyrenean
not o orthy that orthy populations other sampled nly (Wa a deep genetic a deep structure of ng lowland speciesbut
populations populations ,
, 2013)
(Noguerales, (Noguerales, evidencing our landscape our landscape . Only have y This article isThis article by protected copyright. All rightsreserved. Accepted of Europeanconservation status grasshoppers has assigned the twotaxato different IUCN for can have species conservation andmanagement. Therecent assessment ofthe S1) Table synonymized unique a into taxon taxonomic ofthe status complex that andindicate delimitationspecies and analyses, distinctivenessecological ofinformationsources to delimit speciesthathave beendescribed solelymorphology by or divergedallopatrictaxaandrecently highlightsthe exemplifiesThis study problems the associated with 4.3. 2016b& Ortego, Article Billman,Ingley, Johnson, Belk, & 2014; Laiolo risk,predation microhabitat structure or adaptations by ecologicalandevolutionary aspects consideredin not Makinen, & 2006;Merila, Leinonen, O'Hara, Cano,&Merila, 2008 selectiondriven differences foundamong populations not are a consequenceof gen orenvironmentaldistances dissimilarity. Overall, suggest this found by populations among differentiation, notexplained were genetic geographical variationgenetic inthecomplex ( analyses,which genetic no revealed effect ofenvironmental dissimilarity on structuring Conclusions O ur geometric morphometric analyses showedthatthe subtle pheno .
Our study also Our study ( )
Keller, Alexander, Holderegger, &Edwards,2013;Leinonen, Cano, .
illustrates implications the taxonomic that incorrect designations ( e.g.,
Jones : i.e., phenotypic and did ecologicaldata notsupport the current O. antigai
isolation &
Weisrock
(Bolívar, 1897) - by ,
Illera, Illera, & Obeso, - O. navasi , environment;
importance ofintegrating different 2018)
to different resources feeding species validationtests involving this this study selection, assexual such
. Phylogenetic inferences,. Phylogenetic Bayesian
and
(for list ofsynonyms a s
that that the Shafer O. antigai 2013; etic driftetic orenvironmental )
andmight beexplained weak & Noguerales
Wolf
typic typic differences must phenotypic ,
2013) be
( e.g., , Cordero, , Cordero, ,
see .
- This article isThis article by protected copyright. All rightsreserved. this of version earlier an on comments constructive and helpful for referees anonymous three sequencing. Illumina for Genomics) Applied for Centre (The Pereira Sergio Rodríguez Francisco libraries genomic Ampa to grateful are We ACKNOWLEDGEMENTS Knowles, 2016). distribution variationof genetic could 2011 analysesthat genetic studies genetic structure, strong phenotypiclittle variation, environmentalniche. andawide Future 2002). significant units populations bealsoshould considered conservationpotential infuture as plans evolutionary b ofthe status shouldbeprobablydownlistedto“ taxon ecologicalthe andhabitat requirements assessmentsRed List theresults, conservation ofthe status complex totalreconsider requires et (Braud al. categories Red List ) .
Accepted Article However, the presenceHowever, the of two ) greatly help to greatly helpto and Overall, our analyses performing reconstructions ands demographic inferred ,
2016a)
that
, Víctor Noguerales and Pedro J. Cordero for their help during fieldwork, during help their for Cordero J. Pedro and Noguerales Víctor , distributional shifts (Hochkirch al.(Hochkirch et
probably . Given that . Giventhat the andpopulationsizesof range and further integrat - áce fr rvdn valuable providing for Sánchez O. antigai o Hidalgo ro
understand point point to e
d available eserve beprotected to within among and marked genetic groups
as “vulnerable” , the presencethe of bya asingle species characterized
2016) of species the - much Galiana the the information on , considering historical wider than previously than thought,theconservationwider
o hr upr drn te rprto of preparation the during support her for
the the (Braud et(Braud al. linked linked to Pleistocene glacia processes near threatened” patiotemporally explicit and managed independently current past past advice
and many genetically
O. navasi glacier glacier populations
that have shaped in spatial data analyses and and analyses data spatial in ,
2016b)
O. antigai as a as a “ extent (
ation ( Braud et al. . view ofIn the our endangered
(e.g., Massatti(e.g., & Ehlers et al.
in future IUCNin future
are larger and landscape We also thank also We sub l cycles the the - (Moritz , structured
” species 2016a, spatial ,
,
This article isThis article by protected copyright. All rightsreserved. Accepted for download fr BioProject Raw DATA ACCESSIBILITY manuscriptthe samples. AP VT, AUTHORCONTRIBU Article fellowship,research respectively. (SEV Ochoa Severo a by supported were JO and AP work, this During Competitividad. y Economía de Ministerio from 73159) CGL2017 and Fund Development Regional European the and Competitiveness This resources. computer to (CESGA) Doñana. de Biológica Estación (LAST Teledetección y Geográfica Información de Sistemas de Laboratorio and article. I llumina have reads beendepositedat NCBIthe
Logistical support was provided by La by provided was support Logistical
n J cnevd n dsge te td ad nlss Al uhr cletd the collected authors All analyses. and study the designed and conceived JO and VT and Doñana's Singular Scientific Singular Doñana's and PRJNA543714
performed the lab work and analysed the data guided by guided data the analysed and work lab the performed - with 83433 om Repository DryadDigital the help ofJO and - P TIONS . T a spotd y an by supported was VT ).
. Morphometric dataandMorphometric files input for work
e lo hn t Cnr d Sproptcó d Galicia de Supercomputación de Centro to thank also We
inputs frominputs was funded by the Spanish Ministry of Economy and and Economy of Ministry Spanish the by funded was - 2012 - 0262) and a Ramón y Cajal (RYC Cajal y Ramón a and 0262) - Te boratorio de Ecología Molecular (LEM Molecular Ecología de boratorio chnical Infrastructure (ICTS Infrastructure chnical AP ( doi:10.5061/dryad.249v096 .
P peotrl elwhp (BES fellowship predoctoral FPI Sequence Read ArchiveSequence (
(ERDF) all all analyses are available AP and AP (CGL2014 - RBD) ). SRA
- JO. 2013 – EBD) from from EBD) for access for ) VT -
- 54671 under under
12501) 12501) - -
2015 EBD) EBD) wrote wrote - P -
This article isThis article by protected copyright. All rightsreserved. AcceptedBidegaray D.,Tucker,P.K.,& Baiz, M. C.,Arnold, J., Avise, J. R. M., Ball, Bermingham, E., Lamb,T.,Neigel, J. E.,.Saunders, ArticleAlvarez Agapow,M., Bininda P. Adams,C., &Otarola D. REFERENCES 4 analysis ofgeometric morphometric shapedata. Molecular Ecology, 25 by revealed andpalaeodistributionmodellingphylogeography of spider. anendemic P., &Arnedo,M.A.(2016). 1069. doi:10.1111/mec.14966 isolation reproductive in a primatezone. hybrid 489 geneticspopulation andsystematics. N. C. Intraspecific phylogeography (1987). 64 Peninsula.northern section oftheIberian (Platyhelminthes, Tricladida)inmicrorefugia during theLast Glacial Maximum inthe persistencethe species oftheland planarian studies. J. Marshall, C., &Purvis, concept A.(2004).Theimpactonbiodiversity ofspecies - (4), 393 (4), Presas, M.,Mateos, E.,Vila (3), 491 (3), - Batista, L.,Sanchez - 522. doi:10.1146/annurev.ecolsys.18.1.489
Quarterly Review ofBiology,79 - 399. doi:10.1111/2041 - 499. doi:10.1016/j.ympev.2012.05.001 - - Emonds, R.P.O. Castillo, E. (2013).
Cortes (9), 2046 - Gracia, A.,Santulli, G.,Maiorano, L.,Guisan, A.,Vogler, A. Imprintsmultiple of refugia glacial in Pyreneesthe - -
Ortiz, L. Ortiz, (2019). Multipleforms selectionshape of Farre, Farre, M., S - - 2064. doi:10.1111/mec.13585 210x.12035 , Crandall, K.A.,Gittleman, , Crandall, L.,Mace, G. M., J. Annual ReviewofEcology18 andSystematics, Geomorph (2), 161 (2), Molecular andEvolution, Phylogenetics luys, R., & Riutort, R., &Riutort, M. luys, (2012). - Microplana terrestris
The mitochondrial DNA bridge
Molecular Ecology, 28 Methods Methods in Ecology andEvolution, - : an R package for the collection forthe : anRpackage and 179. doi:10.1086/383542
(Muller, 1774)(Muller, (5), 1056 Evidence for
between between - ,
This article isThis article by protected copyright. All rightsreserved. Accepted Braud, Y.,Hochkirch, A., Fartmann, T.,Presa,J.J.,Monnerat, C.,Roesti, C.,… Inferring G.S.,G. M.,&Ralph,P.L.(2018). Bradburd, continuous Coop, Bouckaert, R. (2010). Wu,Bouckaert, R., Heled,D.,Vaughan,T., J.,Kuhnert, H., Xie,. C.D., Article deBolívar, I. (1908).Algunosortopterosnuevos Marruecos España, yCanarias. Bolívar, I. (1897).Catálogo sinóptico delosortópteros la fauna ibérica. F.L.Bookstein, (1992). C.,Bryson,Blair, W., R. Linkem, McCormack, J.,& C.W.,Lazcano,D.,Klicka, E. J. Univer doi:10.1111/1755 rattlesnakes (Viperidae: Crotalus). phylogeneticilluminate relation Cryptic diversity(2019). Mexican inthe highlands: Thousands help ofUCE loci https://www.iucnredlist.org/species/16084614/70647955 3.RLTS.T16084614A70647955.en (Publication no. (2016 doi:10.1534/genetics.118.301333 geneticpopulation structure across space. Bioinformatics, 26 Computational 10 Biology, (2014). Hi Real SociedadEspañolade Naturaes Sciencias a ). sity Press. BEAST Omocestus antigai 2: ASoftwareevolutionary platform analysis. for Bayesian Morphometric tools for toolsMorphometric for landmark data e.T16084614A70647955.
- DENSIT 0998.12970 (10), 1372(10), , 4 (4): 203 REE (4), e1003537.(4), : TheIUCN of RedList : making senseofphylogenetic ofsets trees. - - 1373. doi:10.1093/bioinformatics/btq110
storia 8 storia Natural, 232. ships, species limits anddivergence timesmontaneof
). Retrieved from
Molecular Ecology Resources,19 Genetics, 210 doi:10.1371/journal.pcbi.1003537 http://dx.doi.org/10.2305/IUCN.UK.2016 , 317
Threatened Species - 334. . Cambridge, UK: Cambridge (1), 33(1),
- 52. and discrete Annaes de Drummond, A.J. (2), 349
2016 Dusoulier, F. Boletin de la PLoS
- 365. -
This article isThis article by protected copyright. All rightsreserved. AcceptedClemente, M.E.,Garcí Catchen, J. W.,M., Amores, A.,Hohenlohe,J.H.(2011). P.,Cresko, &Postlethwait, W.Catchen, J., A.,&Cresko, Hohenlohe, P.A.,Bassham,S., Amores, (2013). A. B. Carstens, T.C., Pelletier, A.,Reid,J.D.(2013). N.M.,&Satler, Howtofailatspecies Article F.Burbrink, M., T.,Yao,Ingrasci, Bryson,R.W., H., T.J.,&Ruane,S.(2011). Guiher, D.,Bouckaert, Bryant, J.,Rosenberg,R., Felsenstein, N.A., &RoyChoudhury,A.(2012). Braud, Y.,Hochkirch, A., Presa,J.J., C., Roesti, C.,… Rutschmann, F.,Monnerat, Dusoulier, Inferring species species treesInferring directlyfrom geneticm biallelic https://www.iucnredlist.org/species/69641425/69641429 3.RLTS.T69641425A69641429.en (Publication no.e.T69641425A69641429. F. (2016b). taxonómicas especí navasi Corroboration status of specific of the 1 Genomes Genetics, STACKS doi:10.1111/mec.12354 analysis tool forset population genomics. delimitatio doi:10.1016/j.ympev.2011.05.009 pyromelana at MogollonSpeciation the Ri doi:10.1093/molbev/mss086 a fullcoalescentanalysis.
Bolívar : Building andgenotypingloci : Building n. Omocestus navasi ). Molecular Ecology, 22 Molecular Phylogenetics 60 andEvolution, , 1908(Orthoptera,Acrididae). a, M.D. ficas de (3), 171
, Arnaldos, M.I.,, Arnaldos, Romera,J. E.,&Presa, J.(1999). Molecular Biologyand Evolution, 29 Omocestus antiga
- 182. doi:10.1534/g3.111.000240 : The IUCN Red m Kingsnake Mountain ( Arizona inthe
). Retrieved from (17), 4369 Omocestus antigai Omocestus antigai de novo Molecular Ecology, 22 http://dx.doi.org/10.2305/IUCN.UK.2016 i [Confirmació
- List ofThreatened SpeciesList 2016 ( from short 4383. doi:10.1111/mec.12413 Bolívar
arkers: Bypassing in arkers: Bypassing gene trees , 1897) y
(3), 445 - ( read sequences. Bolívar n del
(8), 1917 (8), O. navasi (11), 3124(11), as posiciones - , 1897) and 454. Lampropeltis - 1932.
G3 S Bolívar - TACKS 3140.
- O. Genes , : an - This article isThis article by protected copyright. All rightsreserved. natureDobzhansky, T.(1937). Genetic ofspecies differences. d B.(2005).Dayrat, integrativetaxonomy. Towards M.Cigliano, M.,Braun,H.,Eades,D.C.,&Otte, (2019). Chifman, &Kubatko,L. (2014). J., O.,Dupont,Charrier, P.,Porn A.Chan, K. Alexander, O., M.,Grismer, L.L.,Su,Y.C.,Grismer,H., & Quah, E.S. J.L., Clemente, M.E.,Garcí Accepted concepts species delimitation.e Queiroz, K. (2007).Species and Article Brown, R. M. (2017). Brown, R.M.(2017). Catalanad'Historia Natural,58 Institucio (In Secci Natural Acrididae).].1908 (Orthoptera, 420. doi:10.1086/280726 879 85 5.0/5.0.RetrievedVersion from model. doi:10.1371/journal.pone.0092976 inthePyrenees. (Ericaceae) reveals the complex phylogeographic of history doi:10.1111/mec.14296 Frogs.in MalaysianTorrent comparison population and genomics approach to cryptic elucidate species boundaries (3), 407 (3), secta: yprepirineo Orthoptera) del Pirineo Catalano - 886. doi:10.1080/10635150701701083 Bioinformatics, 30 Bioinformatics, - 415. doi:10.1111/j.1095 ó n Biol a, M. D., & Presa, J. J. J. (1990).a, M.D.,&Presa, sobrelosAcridoidea Nuevosdatos ó gica, 95 Species delimitation withgeneflow:methodological A on, A.,&Escaravage, Microsatellite marker N.(2014). analysis
(23), 3317(23),
Quartet fromQuartet inference the SNP under coalescent data PLoS Molecular Ecology,26 (3 Boletin de la RealBoletin dela EspaSociedad - http://orthoptera.speciesfile.org/ 4), 27
One, 9
- 8312.2005.00503.x - 3324. doi:10.1093/bioinformatics/btu530 - 50. (3), e92976.(3), , 37 Biological Journal of Biological Journal Linnean the Society,
- 44. Rhododendron ferrugineum
(20), 5435(20), Orthoptera SpeciesFile Orthoptera - American Naturalist, 71 American Naturalist,
Aragonés. Systematic Biology,Systematic 56 - 5450.
ñ Butlleti de la ola de ola de Historia .
, 404 (6), (6), - This article isThis article by protected copyright. All rightsreserved. AcceptedFriedline Excoffier, G.,&Schneider, S. L.,Laval, (2005). Ekblom,R., & J.,Gibbard,Ehlers, P.L., &Hugues,P.D.(2011). L.(2014).SpeciesdetectionandindividualEdwards, D.L.,&Knowles, assignment in ArticleEaton, D.A.R. (2014). Dowle, E.J.,Morgan Domingos, F.,Colli,G.R.,Lemmon, A.,Lemmon, E.M.,&Beheregaray, L.B.In (2017). Dobzhansky Phylogenetics 107 andEvolution, in diversity endemica Cerrado (Squamata: lizard Tropidurus). shadows:the Phylogenomics delimitation species andcoalescent unveil cryptic Press. University doi.org/10.1111/mec.15069 gradient.latitudinal A. J.(2019)Evolutionary genomicsmoth ofgypsy populations sampled along 47 packag ofecology non Acloser chronology: look B Society delimitation:species can integrative increase data efficacy? Bioinformatics, 30 doi:10.1186/s12862 gene grasshopper.despite flow endangered inan , C. J., - 50. , T. (1970). , T.(1970). e for geneticspopulation data analysis.
Galindo, Galindo, ApplicationsJ. (2011). of generation sequencing next inmolecular
Faske - B iological 281 Sciences, - , T.M., Richards, M.,& - model organisms. PYRAD Genetics oftheGenetics Evolutionary Process
(13), 1844(13), - Molecular Ecology, 014
Lind : assemblyof - 0216
(1st ed.). (1st Amsterdam, Elsevier. TheNetherlands: , B.M.,
- 1849. doi:10.1093/bioinformatics/btu121 -
x Trewick, differentiationTrewick, S. A.(2014).Morphological
, 455 Heredity, 107 Heredity, , 20132765 Hobson de novo - 465. doi:10.1016/j.ympev.2016.12.009 28 A RLEQUIN (9), 2206(9), , E.M., Quaternary glaciations
RADseq loci forphylogenetic analyses. RADseq loci . doi:10.1098/rspb.2013.2765 Evolutionary Bioinformatics 1 Online, (1), 1 BMC
Parry
ver. 3.0: an 3.0:an integrated softwarever. - 2223 - . NY: NewYork, Columbia 15. doi:10.1038/hdy.2010.152
Evolutionary Biology,Evolutionary 14 , D., . Proceedings of the Royal ofthe Proceedings Dyer Molecular , R.J.… –
Extent andExtent
Eckert
a
, 216. ,
,
This article isThis article by protected copyright. All rightsreserved. AcceptedspeciesHuang, J.P.,&Knowles, L.(2016).The versus subspecies conundrum: beenandwhatHuang, J.P.(2018).Whatcan have bedelimited asspecies usingmolecular Hochkirch, A., Nieto, A.,García ArticleHime, S., P.M.,Hotaling, E.Voss,S.R., Grewelle, R. E.,Shaffer, O'Neill, H.B.,& P.D.N., A.,Ball,Hebert, Cywinska, S. L.,&DeWaard,Biological J.R.(2003). Garcí M.Fujita, K.,Leaché a - doi:10.1111/jeb.13131 (Orthoptera: Acrididae). evolution body size and sizedimorphismsexual inshort Navas, V.,Cordero, V., Noguerales, P.J.,J.Ecological drivers& Ortego, (2017). of & Evolution, 27 Coalescent doi:10.1093/sysbio/syv119 ap Bayesian delimitationQuantitative frommultiple integrating single datatypeswithin a doi:10.1093/isd/ixy003 (Dynast underthedata multi Luxembourg, Luxembourg: Publications Office Union. oftheEuropean Tumbrinck, J.(2016). Molecular Ecology, 25 delimitation:species an empirical test case in D.W.Weisrock, Theinfluenceof number (2016). andinformation locus on content 270Sciences, identifications through DNA barcodes. es: Dynastidae). - based speciesbased delimitation inanintegrative taxonomy. proach in hercules in proach beetles. (1512), , A. D., Burbrink, F. T., McGuire, J. A., & Moritz, J.A.,&Moritz, , A.D., F.T.,McGuire, C.(2012). Burbrink, (9), 480 (9), - species coalescent model? A casestudy using 313 European andbush ofgrasshoppers,crickets red list
(23), 5959
- Insect SystematicsInsect and 2 Diversity, 488. doi:10.1016/j.tree.2012.04.012 Jo ‐ - Criado, M.,Cálix,Braud, F.M.,… Y.,Buzzetti, 321. doi:10.1098/rspb.2002.2218
urnal ofEvolutionaryurnal 30 Biology, - 5974. doi:10.1111/mec.13883 Systematic Biology,65 Proceedings of theProceedings of Royal SocietyB
an salamander.endangered Mexican
- horned grasshoppers
(2)
(8), 1592 (4), 685 (4), , 1
- 10 Trends in Ecology Trends in Hercules . - - 699. 1608.
- ‐ Biological crickets beetles . This article isThis article by protected copyright. All rightsreserved. AcceptedKnowles, L.L.,& Keller, I., Alex K.S.,&Weisrock,Jones, rejectthehypothesis D.W.Genomicdata ofsympatric (2018). Jombart, S., T.,Devillard, F. &Balloux, (2010). of Discriminant analysis principal Jombart, T.(2008). Article N.D.,B.Jackson, A.E.,&O'Meara, Carstens, B.C. C.,Morales, Species (2017). N.J. Isaac, B., Mallet, J.,&Mace, itsinfluenceG. M.(2004).Taxonomicinflation: on S.J.,Ingley, Billman,divergenceJ., Belk,M. Morphological C., &Johnson,J.B.(2014). E. Hutchison, W., geneticD. A.R.(1999).Correlation ofpairwise &Templeton, and PL bypredation environmentdriven andbetweenspecies of within distributionon the of geneticvariability. geographic Inferring therelativedistance measures: influencesof geneflowanddrif trees. trees. Evolutionary 26 Biology, phenotypic and genetic divergence alongaltitudinal gradients inanimals. 2378 speciationof inaclade ecological BMC Genetics methodcomponents: new ofgenetically a foranalysis the populations.structured Bioinformatics, 24 doi:10.1093/sysbio/syw117 delimitation flow.with gene doi:10.1016/j.tree.2004.06.004 macroecology andconservation. oS -
2393. doi:10.1111/evo.13606 Systematic Biology,56 One, 9 ander, J. ander, R.,& M., Holderegger, WidespreadEdwards, P.J.(2013).
Carstens, B.C.(2007).Delimiting specieswithout monophyleticgene Adegenet (2), , 11
e90274. , 94.doi:10.1186/1471 (11), 1403(11), : a Rpackageforthemultivariate analysis of genetic (12), 2527 (12), doi:10.1371/journal.pone.0090274
(6), 887(6), Systematic Biology, Systematic 66 - 1405. doi:10.1093/bioinformatics/btn129
Trends Ecology 19 in &Evolution, Desmognathus
- 2543. doi:10.1111/jeb.12255 - 895. doi:10.1080/10635150701701091 Evolution, 53 - 2156 - 11
salamanders. - 94 (5), 799 (5),
(6), 1898(6), -
812. - Evolution, 72 Brachyrhaphis 1914.
(9), 464
Journal of Journal markers. markers. -
469. (11), (11),
fishes. t This article isThis article by protected copyright. All rightsreserved. AcceptedLeaché Leaché Leaché Leaché Leinonen, T., O'Hara, R.B., Merila, Cano,J.M.,& (2008). Comparative studiesof ArticleLeinonen, T.,J.(2006). ofbody Cano,J.M.,Makinen,H.,&Merila, Contrasting patterns Lamichhaney, J.,Almen, M.S., S.,Berglund, P.,Illera, J.C., Laiolo, & Obeso, R.(2013).Local climate J. determines intra , A. D., & Fujita, M. K. (2010). Bayesian species delimitation, A.D., M.K.(2010).Bayesianspecies &Fujita, in West African forest , A.D., Z , A.D., M.T., A.(2018). McElroy, &Trinh, Agenomic oftaxonomic evaluation , A.D., Fujita, M. K.,Minin,V.N., delimitationR. &Bouckaert, Species (2014). by genome sequencing. by genome sequencing. L.(2015).A., .Andersson, E doi:10.1111/jeb.12213 communities. interspecific variation in species. species. 27 trends through (genus time horned lizards coast in doi:10.1093/sysbio/syu018 genomeusing 277Sciences, geckos ( Evolutionary 21 Biology, quantitative andn trait 9101.2006.01182.x sticklebacks. andneutralshape genetic divergence i (13), 2884(13), Hemidactylus fasciatus Systematic Biology, 68 hu, T. Q., Rannala, B., &Yang,Z.H.Thespectreofhu, T. (2019). toomany - Journal of Journal Evolutionary 19 Biology, (1697), 3071(1697), 2895. doi:10.1111/mec.14715 - Journal ofJournal Evolutionary 26 Biology, wide SNP data. data. wide SNP
eutral marker divergence: a meta
Nature, 518 sexual sizedimorphismsexual mountain in grasshopper (1), 1 (1), -
3077. doi:10.1098/rspb.2010.0662 - Systematic Biology, 63 (1), 168(1), 17. doi:10.1111/j.1420 volution ofDarwin's beaks their revealed finchesand ). Proceedings of RoyalSociety B the (7539) n marine populations ofthreespine andlake Maqbool, K.,Grabherr, M., Maqbool, Martí - 181. doi:10.1093/sysbio/syy051 , 371
- 375 (6), 1803 (6),
(10), 2171 (10), Phrynosoma . doi:10.1038/nature14181 (4), 534(4), - 9101.2007.01445.x - analysis. analysis. - 1812. doi:10.1111/j.1420 - 2183. -
542. ). Molecular Ecology, Journal of - - Biological
and
nez
- Barrio, Barrio,
- This article isThis article by protected copyright. All rightsreserved. AcceptedMilá McRae, B.H.,&Beier, P.(2007). Circuit theory predictsgene flow animal inplantand McRae, B.H.(200 E.(1942).Mayr, R., Massatti, Contrasting foralternative &Knowles,L.(2016). support modelsofgenomic Article N.A.,Mason, &Taylor, S.A.(2015).Differentially expressed match genes billmorpholog G.M.Mace, The ofta (2004). role Llucia S.D.,Leavitt, St Johnson,L.A.,Goward,T.,& (2011).Clair, Speciesdelimitation in , B.,Carranza, S.,Guillaume, &Clobert,J. (2010). O., - Pomares, D.(2002).Revision(Insecta) oftheOrthoptera ofCatalonia (Spain). de laSociedad Aragonesa, Entomológica 7 [Revisión de losortópteros (Insecta: Orthoptera) de Cataluña (España).]. Phylogenetics 60 andEvolution, chemically diverse taxonomically difficu extreme dispersal limitation inthePyrenean brook newt America, 104 populations. doi:10.1111/j.0014 Press. sedges.in alpine basedonmicrvariation Molecular Ecology, 24 and plumagelargely despite undifferentiatedgenomes a songbird. in Holarctic 711 Royal B the ofLondonSeries Transactions of Society - 719. doi:10.1098/rstb.2003.1454
Systematics and the ofSpecies Origin 6). Isolation byresistance.6). Isolation Proceedings of the National the Proceedings of Academyofthe Sciences Uni of (50), 19885(50), Molecular Ecology, 25 - Xanthoparmelia 3820.2006.tb00500.x lt lt lichen (12), 3009 ohabitat preference:ohabitat species - 19890. doi:10.1073/pnas.0706568104 xonomy conservation. inspecies - forming Anexample fungi: morphologically and from (3), 317 - 3025. doi:10.1111/mec.13140 (Parmeliaceae) in North(Parmeliaceae) in America.
Evolution, 60 (16), 3974(16), -
332. doi:10.1016/j.ympev.2011.05.012 , 1 - 226. . New York, NY: ColumbiaUniversity
- - 3986. doi:10.1111/mec.13735
(8), 1551 specific effects ofclimate change Marked genetic structuring and - Biological Sciences,359 Calotriton asperCalotriton Philosophical - 1561.
Molecular Monografias Statesof ted (Amphibia: (1444),
y This article isThis article by protected copyright. All rightsreserved. Accepted F.Oksanen, J.,Blanchet, G.,F V.,Noguerales, P.J.,&Ortego,J.(2018). Cordero, V.,Noguerales, García Article V.,Noguerales, P.J.,&Ortego,J.(2016 Cordero, Niemiller, M.L.,Fitzpatrick,B.(2012).Delimiting using Near,T.J.,& species Nelson, G.,&Platnick, N. (1981). C.(2002).Moritz, Strategies toprotectbiological diversity evolutionaryand the processes University Press. University sustainthat it. 19 Salamandridae) revealed by genome from H. (2017). doi:10.1111/mec.14504 radiationofevolutionary grasshoppers. alternativeto evaluate phylogenetic anddelimitation hypotheses species in arecent doi:10.1111/jeb.12915 montane grasshopper. environment andcore 16 top doi:10.1111/j.1558 subterraneus multilocus data:Diagnosing cry (1), 108 (1), ,
ography ography in anarrow 96. http://cran.r doi:10.1186/s12862 - vegan 120. doi:10.1111/j.1365
(Teleostei: Amblyopsidae). Sys - Navas, V.,Cordero,P.J.,& : Community EcologyPackage. 2.4 Rpackageversion
- tematic 51 Biology, project.org/package=vegan - 5646.2011.01480.x - Journal ofEvolutionaryBiology, 29 margin effects onrange riendly, M., Kindt, riendly, M.,Kindt, R., Legendre,P.,McGlinn,D.,…Wagner, -
endemicmontane grasshopper. - Systematics andSystematics Biogeography 016 ptic diversity inthe southern cavefish, - 0663 - 294X.2009.04441.x - wide AFLP but not mtDNA.wide AFLPbutnot (2), 238 (2), - 7
Molecular Ecology, 27
Evolution, 66 a ).
Integrating genomicand phenotypicdata Ortego, of J.(2016b). The role
- Hierarchical geneticHierarchical shapedby structure 254. - wide phenotypicvariationina
(3), 846(3), B
MC . NewYork, NY:Columbia (11), 2129
Evolutionary Biology, (5), 1229 (5), - 866. Molecular Ecology, Molecular - 2142. Typhlichthys Typhlichthys ‐ - 4. Retrieved 1244. This article isThis article by protected copyright. All rightsreserved. Accepted F. W.,Pyron, R.A.,Hsieh, A.R., Lemmon, &Hendry,C. Lemmon, E.M., R.(2016). S. The Puissant, (2008). ofpresence J.P.(2000).Pritchard, &Donnelly, InferenceK., Stephens,M., ofpopulationstructure using T. Pons, J.,Barraclough, G.,Gó Poniatowski, D.,Defaut Article B. Peterson, Weber,K., J methodsO'Meara, delimitation Newheuristic jointspecies B.C.(2010). for andspeciestree Olmo - Vidal, J.M. inference. inference. Habitatge. Spain:Barcelona, 2002Libro Rojo oftheOtthopteraCataloniaandRead /Atlas DataBook2002 Saltamartins, SomeretesLlagostes, /Atlas Grills, deCataluña y de LosOrtópteros Integrating Entomocenotiques,Orthopteriques et 13 nouvelle sous en presence d' France description (Caelifera,of Acrididae, anewsubspecies Gomphocerinae). multilocus genotype data. undescribed insects. A.P.(2006). Vogler, fauna ofregion Pyrenean the doi:10.1371/journal.pone.0037135 model andnon RADseq:Aninexpensive methoddigest for
Systematic Biology, 59 (2002). phylogenomic andmorphologicalphylogenomic datatoassess speciescandidate - espece (Caelifera,Acrididae, Gomphocerinae).]. - model species. , B.,Llucia , . N., Kay, E. H., Fisher, H. S., & Hoekstra, H.E. . N.,Kay,H.,Fisher, Double H. S., E. (2012). &Hoekstra, Generalitat Catalunya, de Departament de AmbientMedi i Atlas Dels Ortòpters de Dels i VermellAtlas Ortòpters Catalunya Llibre 2002/ Systematic Biology, 55
Omocestus navasi Omocestus Sequence mez Genetics, 155Genetics, - Pomares, D.,&Fartmann, T.(2009). The Orthoptera - Omocestus navasi - Zurita, A.,Duran,Zurita, S., J., Cardoso, D.P., Hazell, .
a field guide. PLoS - based delimitation species for DNA the taxonomy of (1), 59
One, 7 - (2) Bolívar , 69 73. doi:10.1093/sysbio/syp077 , 945 (4), 595(4), (5), e37135.(5), - de novo Articulata Beiheft, 14 73.
, description1908, avec d'une Bolívar - 959.
-
609. SNP discovery and SNPgenotyping discovery in
, 1908inFrance with
Materiaux Materiaux , 1 - 143.
[Sur la
- . This article isThis article by protected copyright. All rightsreserved. AcceptedRundle, H.D., speciation. Ecological &Nosil, P.(2005). & Slice,D.( Rohlf, F.J., Rohlf, The F.J.(2015). Reynolds, W.J.(1986).Adescriptionofthesong Yang,Z Rannala, B.,& Article L., Rancilhac, Goudarzi, Hemami,F., Gehara,M., R., Elmer, M. K.R.,Vences, M.,& M., Raj, A.,Stephens, & Pritchard,J. W. Reynolds,Ragge, D.R.,& J.(1998). R Core Team(2018). western Europewestern for Computing.Austria: Foundation Statistical Linnean 177 the Society, delimitation brownmodels in andred doi:10.1111/j.1461 doi:10.2307/2992207 ofsuperimposition landmarks. 9 111 onAcrididae) with notes delimitation.species Phylogenetics 133 andEvolution, ongenomenewts (Salamandridae) based S.Steinfarz, Phylogeny delimitation (2019). andspecies of near Eastern doi:10.1534/genetics.114.164350 structurepopulation SNPin large sets. data - 12. doi:10.4404/hystrix - 116. doi:10.1080/00222938600770101 R . Colchester,UK: : A language : Alanguage and environm . H.(2013).Improved reversible jump algorithms Bayesian for TPS 1990). Extensionsofthe procrustes method optimal for the - 0248
series ofseries software. Genetics, 194
(4), 937 (4), .2004.00715.x its taxonomicits position. - 26.1
K. (2014). - Systematic Zoology,39 Systematic 11264 -
949. doi:10.1111/zoj.12392 The songsgrasshoppers the andof of crickets , 189 Harley Books. (1), 245 -
bellied snakes ( -
197. doi:10.1016/j.ympev.2019.01.003 Hystrix - FAST wide RADseq data analysis. wide RADseq analysis. data - ent for statistical computing. Vienna, forent statistical
253. doi:10.1534/genetics.112.149039 Genetics, 197 Omocestus broelemanniOmocestus STRUCTURE
- Italian of Journal Mammalogy,26
Journal of History,20 Journal Natural
Ecology Letters, 8 Storeria (1), 40 (1), (2), 573 : Variational inference: Variational of -
59. ). Zoological Journal ofZoological Journal - 589 (3), 336 (3),
. (Orthoptera: Molecular Neurergus - 352.
(1), (1),
(1), (1), This article isThis article by protected copyright. All rightsreserved. Soria Solis Sokal, R.R.,&Crovello, T.J.(1970).Biological species concept Slatkin, M. (1993).Isolation bydistance equilibrium in andnonequilibriumpopulations W., J. Sites, & Simpson, G.(1944). G. Da Wolf,W. evidenceShafer, A.B.A.,& J.B.Widespread ecological forincipient (2013). P.,Hangartner, isolationSexton, J. S.A. (2014).Genetic B.,&Hoffmann,A. byenvironment EcologySchluter, ofspecies. D.(2001). andthe origin
Accepted Article Silva, S.B.,&DaSilva,A. Silva, (2018). - - Lemus, delimit L.L.,&Ane,C.(2015).Bayesian species C.,Knowles, Carrasco, Z.,Comeault, V.,Gompert, A., Farkas,T.E.,Parchman,L., Johnston, A. T. speciation: ametaspeciation: doi:10.1111/evo.12258 or distance: Which ofgeneflowism patterns 372 parallel parallel speciation. P.(2014). J. S.,.Nosil, Stick insect genom doi:10.1111/evo.12582 genesmultiple traits and inaunified framework. 104American Naturalist, 47 Evolution, 5347(03)00184 biology. Press. differentiation traits. inphenotypic doi:10.1111/ele.12120 - 380.
Marshall, J. C.(2003).DelimitingaRenaissance systematic species: issuein Trends Ecology in &Evolution, 18
(1), 264 (1), - 8 Tempo and Mode in EvolutionTempo andModein
- analysis ofisolation Science, 344 - 279.
(936), 127 (936),
Pstat (6185), 738 The RJournal - : An R Package to assess population: AnRPackage toassess 153 - by . doi:10.1086/282646 - ecology. - (9), 462(9), 742. doi:10.1126/science.1252136 es reveal natural selection's es reveal role in ost common? ost common? . New York, NY: ColumbiaUniversity Trends in Ecology & Evolution, 16 Evolution, 69 Evolution, ,
10 - Ecology Letters, 16 470. doi:10.1016/s0169 (1), 447 (1), -
Evolution, 68 A critical evaluation. -
454. (2), 492 (2),
ation combiningation - 507. (7), 940 (1), 1 - -
15. - . 950. (7), (7), This article isThis article by protected copyright. All rightsreserved. Accepted effects oflandscapeWang, Examining heterogeneity full J.(2013). the I. genetic onspatial Waters, G.P., Wallis, J. alpine D. Transverse M., Upton,P.,&Craw, (2016). speciation E. van Dijk, L., Auger, H.,Jaszczyszyn, Y.,& VanderWerf, N.W., E.A.,Young,L.C.,Yeung, &Carlon, D. B. Step (2010). ArticleValbuena D.,Arctander,Tautz, P., Min Swofford, L.(2002). D. coalescentSukumaran, Knowles, L.(2017).Multispecies J.,& delimits structure, not taxonomy. Version Sunderland,Methods). 4. MA: Sinauer Associates. America, 114 species. ecological ecological isolation. Amultiple matrixvariation: regression for approach quantifyingand geographic doi:10.1016/j.tree.2016.08.009 byglaciation.driven doi:10.1016/j.tig.2014.07.001 sequencinggeneration technology. doi:10.1007/s10592 endemics the within elepaio. inspeciation molecular Hawaii'sdivergence flycatchers: supportsnewisland doi:10.1371/journal.pone.0200214 endemic across newtspecies the Pyrenees. . 5347(02)00041 - Steinfa Ure ñ a, E.,Oromi, N., Soler Proceedings of theNational Academy of rtz, S.(2018).Jailed diversityin the mountains: Genetic and structureofan Trends inEcology &Evolution, 18 Trends (7), 1607 (7), - 1 PAUP
- Evolution, 67 009 Trends elli, A., A.P.(2003).ApleaforDNAThomas,H., &Vogler, R. *. Phylogenetic Analysis Using Analysis *. Phylogenetic Parsimony Other (*and - 1612. doi:10.1073/pnas.1607921114 - 9958
in Ecology &Evolution, 31 Conservation Genetics, 11
- -
1 Membrives, F.,Camarasa, A.,Carranza,S.,Amat, S.,.
(12), 3403
Trends in Genetics,Trends in 30
Thermes, Tenyearsofnext C.(2014). P - LoS 3411. doi:10.1111/evo.12134 (2), 70
One, 13 Sciences of the United States of oftheUnitedStates Sciences
- 74. doi:10.1016/s0169 (4), 1283(4), (12), (12), 916 (8) (9), 418
, e0200214
- - 1298. - 926. 426. . ping stone
- - This article isThis article by protected copyright. All rightsreserved. Accepted speciesYang, Z.H.,&Rannala, delimitationB. (2010).Bayesian multilocus using sequence The Yang, Z.H.(2015). ArticleWiens, delimitation& Servedio, Species M.R.(2000). J.J., in systematics: inferring SpeciesWiens, delimitation:for(2007). Newapproaches discovering J. diversity. Whelan, N.V., M.P.,Sipley, B.N.,Weber, Galaska, J.M.,Johnson,P.D., K.M., Halanych, Wang, I. Ecology, morphologic variat& HelmsRiverscape genetic S. (2019). B. 5649 America, 107 data. Zoology,61 Current 267Sciences, differencesdiagnostic between species. Systemat J., & Bradburd,J., & G.S.(2014).Isolationby environment. - Proceedings ofthe Academy National of ofSciences the United States of 5662. doi:10.1111/mec.12938 ic ic 56 Biology, 28 (7), 1593 al variation al of threatenedr the round (20), 9264 (20), (1444), 631(1444), BPP (5), (5), 854
- program fortreeprogram estimation species delimitation. andspecies 1610
(6), 875 - 9269. doi:10.1073/pnas.0913022107 - 636. doi:10.1098/rspb.2000.1049 . doi - 865. - : 878. doi:10.1080/10635150701748506 10.1111/mec.15032
Proceedings the of Royal Society B ion, migrationion, and patterns, ocksnail, ocksnail,
Molecular 23 Ecology, Leptoxis ampla.Molecular
- Biological (23), (23), This article isThis article by protected copyright. All rightsreserved. (legendpriors bottom). at represent the meanprobability (PP)for posterior differentcombinations of demographic (gamma, (200, 500,and1000loci three alternative topologies (a: Figure asin 1. codes Table which at and(triangles) Figure 3 1. Table blackvertical separate individuals lines from populations. different segmentscolored probability ofbelonging representing the to cluster the with that color.Thin 2 Figure STRU map antigai and 6) Figure ( Figure 1 legendsFigure
Accepted Articleto K indicate the main the cluster populations genetic which at indicate were assignedaccording to CTURE
= 6) and= 6)
4 2 (triangles) and
individual α, β α, Results ofspecies delimitationResults analyses in Principal component analysis (PCAs) ofgenetic Genetic assignment ofindividuals (a) Bayesian
analyses for analyses )
for ancestral population (θ)andage size root (τ)
b) map of distribution showingthe sampled the populations of O. navasi ( b ) CONSTRUCT s
were assigned were assigned phylogenetic O. navasi ; two individuals/population); two K (dots). (dots). Node support in Node support termspro ofposterior
= 6. Population codes asinTable = 6.Population 1.
BPP ( C from from (dots) the Pyrenees. in Colorsof trianglesanddots olors ; b: according to tree tree reconstructed with K SNAPP
of indicate trianglesanddots = 2 to = 2to
based of on theresults ; c: ; c: K
SVDQUARTETS STRUCTURE = 4). Individuals Individuals are = 4). BPP , and . Analyses were performed were . Analyses
variation for four gamma com prior SNAPP
analyses for analyses ), three different subsets), three of loci . Coloredboxesateach node
babilities (
(node support support (node
( Population codes asin a
) the the partitioned into STRUCTURE Omocestus antigai main geneticcluster K
for = 6. Omocestus SNAPP Population Population presented binations binations
( using using from on the
tree K
K )
in = This article isThis article by protected copyright. All rightsreserved. Accepted Article asin 1. codes Table demog of combinations (PP) fordifferent (θ)androotage size (τ). Colored boxesateachnode mean represent the probability posterior SVDQUARTETS were performedanalyses three using alternativ phenotypic forseparately sex each both combining genomic (500loci and5individuals/population) and Figure Population codes inTable 1. as and bootstrapping values ( 5
Results ofspecies delimitationResults analyses in data (treesondata left) the and onlyusing ) and four) and gamma combinations prior (gamma,for populationα, β) ancestral
for
SVDQUARTETS
raphic priors (legend priors bottom at raphic right). Population
tree e topologies(a: ) is phenotypic indicated IBPP . Analyses were performed were . Analyses
data (treesondata right).All the
in parentheses node. foreach BPP ; b: ; b: SNAPP ; c: ; c: This article isThis article by protected copyright. All rightsreserved. sampled Table 1 AcceptedO. femoralis O. navasi O. navasi O. navasi O. navasi O. navasi ArticleO. navasi O. navasi O. navasi O. navasi O. antigai O. antigai O. O. antigai O. antigai O. antigai O. antigai Species antigai
Geographical location, Geographical andnumberelevation ofgenotypedindividuals ( population
bellmanni navasi navasi navasi navasi navasi navasi navasi navasi
(locality) Sierra de Espuña Sierra de Batère de Tour Tost de Montan Perves Espés Chía Saravillo Nerín de Monrepós Puerto Borau Setcases l'Home de Turó Queralbs Err laCreueta de Coll Llessui Taüll Boí Locality
of the
three
putative taxa ESP CHI SAR NER MON BOR SET TUR QUE ERR CRE LLE BOI Code ESP BAT TOS PER
42.44307 42.55418 42.56089 42.59047 42.35067 42.67388 42.42756 41.77167 42.36627 42.38830 42.30131 42.44554 42.47945 37.86515 42.50785 42.23790 42.35250 Latitude
and the outgroup
Longitude - - - 0.58389 0.43703 0.23046 0.01062 2.26630 2.44335 2.14902 2.08812 1.99350 1.03668 0.86732 2.57641 1.38052 0.83709 0.39762 0.57932 1.57115
.
Elevation (m)Elevation 1361 1691 1313 1623 1267 1357 2145 1622 2070 2150 1928 1979 2046 1514 1430 1182 1370 n ) for each ) foreach
6 7 2 7 5 7 7 7 6 3 6 6 5 n 6 7 7 6
This article isThis article by protected copyright. All rightsreserved. (calculated factors as 2 defining eachdelimitation alternative species ( hypothesis Table 2 H H H H H Species delimitation hypothesis ( 5 4 3 2 1
Accepted : (BAT+SET) (TOS) (CRE+ERR+QUE+TUR) (PER+SAR+CHI+NER+ESP) (MON+BOR) (BOI+LLE) : (BAT+TOS+CRE+ERR+QUE+TUR+SET) Article: (BAT) (PER+SAR+CHI+TOS+NER+ESP+MON+BOR) (CRE+ERR+BOI+LLE+QUE+TUR+SET) : (BOI+LLE+CRE+ERR+QUE+TUR+SET) (PER+SAR+CHI+NER+ESP+MON+BOR+BAT+TOS) : (BAT+PER+SAR+CHI+TOS+NER+ESP+MON+BOR+CRE+ERR+BOI+LLE+QUE+TUR+SET)
Results of Results BFD * analyses analyses testing thesupportofcompeting* species delimitation Thehypotheses. table shows scheme theclustering ×
ln BF) and their rank. codesln BF)andtheir Population inTable1. as H i )
(PER+SAR+CHI+NER+ESP+MON+BOR+BOI+LLE)
H i ). Forwe hypothesis, marginal each show likelihood estimates (ML
6 2 3 2 1 Species
- - - - - MLE 4,384.48 5,221.06 5,384.07 5,601.22 4,794.85
- 1.67 × 10 2.00 × 10 2.43 × 10 8.21 × 10 2
× lnBF
3 3 3 2
1 3 4 5 2 Rank
E), their E), their Bayes
This article isThis article by protected copyright. All rightsreserved. Accepted Article distances. differentiationgenetic ( Table 3
Multiple Matrix Multiple Matrix with RandomizationRegression (MMRR) for
term Rejected terms Retained Variable
Environmental distance Environmental (slope) distance Resistance distance Geographical Intercept
F
ST ) in relation with geographical, resistance relation withgeographical, ) in and environmental
Coefficient - 0.007 0.255 0.108 -
- - 5.278 6.841 1.061 3.079 t
pairwise population pairwise 0.374 0.004 0.001 0.991 p
This article isThis article by protected copyright. All rightsreserved. Accepted ArticleHochberg significant and geographical, resistance and environmental distances. phenotypic differentiation ( Table 4
Multiple Matrix Multiple Matrix Regressions with Randomization (MMRR) for false discovery rate Ovopositor valve Ovopositor valve Forewing Forewing Pronotum Pronotum Head Head
Geographical distance Geographical distance Geographical distance Geographical distance Geographical distance Geographical distance Geographical distance Geographical distance F Environmental distance Resistance distance (slope) F Environmental distance Resistance distance (slope) F Environmental distance Resistance distance (slope) F Environmental distance Resistance distance (slope) F Environmental distance Resistance distance (slope) F Environmental distance Resistance distance (slope) F Environmental distance Resistance distance (slope) F Environmental distance Resistance distance (slope) retained into retained ST ST ST ST ST ST ST ST
– –
PC2 PC1 – – – –
PC2 PC1 PC2 PC1
P
final – – ST
PC2 PC1 corrections ) in relationwithgeneticdifferentiation) in (
models
.
of Table shows p
------values toadjust formultipletests. statistical 1.882 0.125 0.989 0.448 0.617 1.401 0.637 1.614 0.210 1.110 3.159 0.388 2.695 2.775 0.868 0.963 2.781 1.423 1.047 1.699 0.257 0.408 0.282 0.222 2.197 0.190 0.735 1.240 0.267 0.960 0.416 1.264 t
Females
No independent variablewas
q 0.791 0.791 0.686 0.950 0.791 0.791 0.919 0.919 0.791 0.919 0.919 0.919 0.686 0.919 0.919 0.791 0.400 0.919 0.919 0.520 0.686 0.919 0.405 0.919 0.791 0.919 0.919 0.791 0.791 0.919 0.791 0.400
q -
values
after after
------0.135 1.223 2.045 0.782 2.185 0.695 1.152 0.318 0.715 0.587 1.713 0.376 1.735 0.429 0.052 2.993 1.106 0.048 1.471 0.031 2.455 0.386 0.054 0.858 – – – – – – – – t
applying
Males F
pairwise population pairwise
ST
) and ) and 0.980 0.973 0.864 0.980 0.288 0.853 0.980 0.720 0.973 0.496 0.840 0.980 0.288 0.973 0.973 0.853 0.980 0.973 0.973 0.973 0.853 0.973 0.973 0.840
q – – – – – – – –
Benjamini
-
This article isThis article by protected copyright. All rightsreserved. Accepted Article 1 Figure (b) (a) BOR
MON NER SAR CHI ESP PER BOI LLE TOS
CRE ERR QUE TUR SET BAT
This article isThis article by protected copyright. All rightsreserved. Accepted Article 2 Figure (b) (a) CONSTRUCT STRUCTURE
K = 4 K = 3 K = 2 K = 6 K = 5 K = 4 K = 3 K = 2
This article isThis article by protected copyright. All rightsreserved. Accepted Article 3 Figure
This article isThis article by protected copyright. All rightsreserved. Accepted Article 4 Figure
This article isThis article by protected copyright. All rightsreserved. Accepted Article 5 Figure