Evolutionary history of the Afro-Madagascan Ixora species (Rubiaceae): species diversification and

20 Pages·2013·0.91 MB·English

AnnalsofBotany112:1723–1742,2013
doi:10.1093/aob/mct222,availableonlineatwww.aob.oxfordjournals.org
Evolutionary history of the Afro-Madagascan Ixora species (Rubiaceae): species
diversification and distribution of key morphological traits inferred from dated
molecular phylogenetic trees
J.Tosh1,2,*,S.Dessein3,S.Buerki4,I.Groeninckx1,A.Mouly5,6,B.Bremer6,E.F.Smets1,7,8andP.DeBlock3
1LaboratoryofPlantSystematics,KULeuven,KasteelparkArenberg31,POBox2437,BE-3001Leuven,Belgium,2AshdownHouse
School,ForestRow,EastSussexRH185JY,UK,3NationalBotanicGardenofBelgium,DomeinvanBouchout,BE-1860Meise, D
Belgium,4JodrellLaboratory,RoyalBotanicGardens,Kew,Richmond,SurreyTW93AE,UK,5LaboratoireChrono-environnement, ow
UMRCNRS6249,Universite´deFranche-Comte´,16RoutedeGray,F-25030Besanc¸oncedex,France,6BergiusFoundation,Royal nlo
SwedishAcademyofSciencesandBotanyDepartment,StockholmUniversity,SE-10691,Stockholm,Sweden,7NationalHerbarium ad
e
ofTheNetherlands,LeidenBUiondiivveerrssiittyyBNraatnucrha,liPs,OPBOoBxo9x591541,7N,LN-2L3-20300R0ARLAeLideeind,enT,hTehNeeNtheetrhlearnladnsdasnd8NetherlandsCentrefor d from
*Forcorrespondence.E-mail:[email protected] h
ttp
s
Received:9April2013 Returnedforrevision:20May2013 Accepted:6August2013 Publishedelectronically:18October2013 ://a
c
a
d
e
†BackgroundandAimsPreviousworkonthepantropicalgenusIxorahasrevealedanAfro-Madagascanclade,butas m
yetnostudyhasfocusedindetailontheevolutionaryhistoryandmorphologicaltrendsinthisgroup.Heretheevo- ic.o
lutionaryhistoryofAfro-MadagascanIxoraspp.(acladeofapprox.80taxa)isinvestigatedandthephylogenetictrees up
comparedwithseveralkeymorphologicaltraitsintaxaoccurringinMadagascar. .c
o
†MethodsPhylogeneticrelationshipsofAfro-MadagascanIxoraareassessedusingsequencedatafromfourplastid m
regions(petD,rps16,rpoB-trnCandtrnL-trnF)andnuclearribosomalexternaltranscribedspacer(ETS)andinternal /ao
b
transcribed spacer (ITS) regions. The phylogenetic distribution of key morphological characters is assessed. /a
Bayesian inference (implemented in BEAST) is used to estimate the temporal origin of Ixora based on fossil rtic
evidence. le
†KeyResultsTwoseparatelineagesofMadagascantaxaarerecovered,oneofwhichisnestedinagroupofEast -ab
s
Africantaxa.DivergenceinIxoraisestimatedtohavecommencedduringthemidMiocene,withextensiveclado- tra
†geCnoensicsluoscicounrsriBnogthinlitnheeaAgefsroo-fMMaaddaaggaassccaanncIlxaodreadeuxrhinibgitthmeoPrplihoocleongeicoanlwinanrodvs.ationsthatarerarethroughouttherest ct/11
2
ofthegenus,includingatrendtowardspauciflorousinflorescencesandatrendtowardsextremecorollatubelength, /9
suggestingthatthesameecologicalandselectivepressuresareactingupontaxafrombothMadagascanlineages. /1
7
Novelecologicalopportunitiesresultingfromclimate-inducedhabitatfragmentationandcorollatubelengthdiver- 23
sificationarelikelytohavefacilitatedspeciesradiationonMadagascar. /2
7
6
8
Keywords:Rubiaceae,Ixora,Afro-Madagascan,molecularphylogenetics,moleculardating,biogeography,ETS, 91
ITS,petD,rps16,rpoB-trnC,trnL-trnF. 4 b
y
g
u
e
s
t o
INTRODUCTION 1998). In continental Africa, Ixora mainly occurs in the n 0
Guineo-Congolian Regional Centre of Endemism (RCE) (fol- 9 A
The pantropical genus Ixora is one of the largest genera in lowingWhite,1983),butalsointheAfromontanearchipelago- p
Rubiaceae,withapprox.530speciesofshrubsandsmalltrees likeRCE,andextendsintotheZambezianRCE,theSwahilian ril 2
0
that typically grow in humid rain forest (Davis et al., 2009). RCE and the Swahilian/Maputaland regional transition zone 1
9
The centre of species diversity for the genus is in South-East (RTZ)(DeBlock,1998).InMadagascar,Rubiaceaearemostnu-
Asia,inparticularBorneo(Lorenceetal.,2007).Althoughno merousandspeciesrichintheevergreenhumidforests(Davis
modern monograph of Ixora exists, there have been a number and Bridson, 2003). Ixora is no exception to this, occurring
ofrevisionsfocusingonspecificgeographicalregions(e.g.De mostfrequentlyinthehumidevergreenforest(littoral,lowland
Block, 1998, revision of continental African Ixora spp.; De and montane) on the eastern coast of Madagascar, although
Block,2013,revisionofMadagascanIxoraspp.).Phylogenetic Ixora spp. also occur in the semi-deciduous forest of
studiesofIxorahaveprimarilyfocusedonthetribalplacement Madagascar(DeBlock,2003,2013).
and circumscription of the genus (Andreasen and Bremer, Ixora is one of the most easily recognizable genera in
1996,2000;Mouly,2007;Moulyetal.,2009a).Mostrecently, Rubiaceae,inpartduetotheoftenstrikinginflorescencesandtet-
Moulyetal.(2009b)identifiedsomewell-supported,geograph- ramerous flowers (Fig. 1). Diagnostic features for the genus
icallydefinedlineages,includingan‘Afro-Madagascan’clade. (adaptedfromDeBlock,2007)includearticulatedpetioles,nar-
Thereareapprox.80Afro-MadagascanIxoraspp.distributed rowly tubular tetramerous flowers, bilobed stigmas, bilocular
equallybetweencontinentalAfricaandMadagascar(DeBlock, ovaries and fruits (or, rarely, with more than two locules),
#TheAuthor2013.PublishedbyOxfordUniversityPressonbehalfoftheAnnalsofBotanyCompany.Allrightsreserved.
ForPermissions,pleaseemail:[email protected]
1724 Toshetal.—EvolutionaryhistoryoftheAfro-MadagascanIxoraspecies(Rubiaceae)
A B C
D E F
D
o
w
n
lo
a
d
e
d
fro
m
h
ttp
s
G H ://a
c
a
d
e
m
ic
.o
u
p
.c
o
m
/a
o
b
/a
rtic
le
-a
I J K bs
tra
c
t/1
1
2
/9
/1
7
2
3
/2
7
6
8
9
1
4
b
y
L M gu
e
s
t o
n
0
9
A
p
ril 2
0
1
9
FIG.1. ExampleofmorphologicalvariationinIxora.(A)InflorescenceofIxoraregalis.(B)InflorescenceofIxoraelliotii.(C)InflorescenceofIxoraemirnensis.(D)
PendulousinflorescenceofIxoramangabensis.(E)InflorescenceofIxoradensithyrsa.(F)InflorescenceofIxorasiphonantha.(G)FrontviewofaflowerofIxora
guillotii.(H)FloweringbranchofIxorarakotonasoloi.(I)MaturefruitsofIxoraguillotii.(J)FruitsofIxoraquadrilocularis.(K)FloweringnodeofIxorahomolleae.
(L)ArticulatepetiolesofIxorafinlaysoniana.(M)ArticulatepetiolesofIxorahomolleae.
uniovulateloculesandseedswithalargeadaxialhilarcavity.In flowers (De Block, 1998, 2003). This is particularly the case
contrast, identification at the species level is difficult, with for the African representatives of the genus, which De Block
speciesdistinguishedonthebasisofminorandoftencontinuous (1998) described as ‘extremely homogeneous’ in theircharac-
characters,typicallyinvolvingfeaturesoftheinflorescenceand ters. On Madagascar, there are several morphological traits
Toshetal.—EvolutionaryhistoryoftheAfro-MadagascanIxoraspecies(Rubiaceae) 1725
occurringinIxorathatareabsentinthecontinentalAfricantaxa 25mL, and contained 1mL of each primer (100ng mL21),
andrareinthegenusasawhole.Theseinclude:(1)reductionof 0.35mLofBiotaqDNApolymerase,2.5mLof10×NH reac-
4
thenumberofflowersperinflorescencetowardssolitaryflowers; tionbuffer,1.5mLof50mMMgCl2,2.5mLof10mMdNTPs,
(2)increasefromtwo-tofour-locularovaries;and(3)increase 1mLofbovineserumalbumin(BSA;0.4%)and2mLoftotal
towards large flowers (corolla tubes .15cm long) and fruits genomic DNA. Amplification of petD, rps16 and trnL-trnF
(DeBlock,2007,2008,2013). used the following temperature profile: 948C for 3min; 32
Inthepresentstudy,wefurtherinvestigatethephylogenetic cycles of 948C for 1min, 508C for 1min, 728C for 1.5min;
relationships of Madagascan and continental African Ixora finalextensionof728Cfor7min.Theamplificationprofilefor
spp. using molecular sequence data from four plastid regions rpoB-trnCwas:948Cfor3min;32cyclesof948Cfor1min,
(petD,rps16,rpoB-trnCandtrnL-trnF)andnuclearribosomal 538C for 1min, 728C for 2min; final extension of 728C for
external transcribed spacer (ETS) and internal transcribed 7min.
spacer (ITS) regions. The purpose of this study is to improve PCR mixesfornuclear regions were the same as forplastid D
o
taxon sampling of both African and Madagascan species in regions, except that 1mL of dimethylsulfoxide (DMSO) was w
n
orderto:(1)testexistinghypothesesconcerningtheevolutionary addedper25mL.TheITSamplificationprofilewas:948Cfor lo
a
affinitieswithinandbetweenAfricanandMadagascanspecies; 3min;32cyclesof948Cfor1min,528Cfor1min,728Cfor de
d
(2)assessthedistributionofkeymorphologicalinnovationsof 1min;finalextensionof728Cfor7min.TheETSamplification fro
the Madagascan species on our molecular phylogenetic trees; profilewas:978Cfor1min;40cyclesof978Cfor10s,558Cfor m
and(3)investigatetheageofspeciesdiversificationanddispersal 30s,728Cfor30s;finalextensionof728Cfor7min.Allamp- http
usingmoleculardatingtechniques. lification products were purified using Nucleospin purification s
columns and sent to Macrogen Inc. (Seoul, South Korea) for ://a
c
sequencing. ad
MATERIALS AND METHODS em
ic
TaxonsamplingandDNApreparation Sequencealignmentandphylogeneticanalyses .ou
p
Extensive fieldwork was undertaken in eastern and northern Contiguous sequences were assembled and edited using the .co
m
Madagascarinordertocollectherbarium,alcoholandDNAma- Stadensoftwarepackage(Stadenetal.,1998).Sequenceswere /a
terial of Madagascan Ixora spp. This material was used in the manually aligned in MacClade v. 4.04 (Maddison and ob
macocleescsuiolanrs,arnedprtehseemntoinrgphaoplporgoixc.al50stuspdeyc.iWesethinatcloucdceudr6th7roIuxogrha- Mvaaridadtiiosonn.A,2l0l0v2a)riwabitlheonuutcdliefofitcidueltpyodsuietiotonlsowwelreevveelsriofifesdeqaguaeinncset /article
outtheglobaldistributionofthegenus(Table1).Inparticular, the original electropherograms. Gaps were treated as missing -ab
osepaucerhciteasxM.oWnadhsaaegmraespcplaionnsgsibisslpeef,cowiceuesseindtcoloundteeMdstamdusaplgteaipcsilceeasnacacmnedossniAoopnfrhsiycflaoynr. d‘(s2ai0tma0;p0l)pe.oTitonendmteiilanliclmyodiziinnefgtoh’remmtiaemttihveoeadnidnodcfoeSlmsimpwumteaorteinosncaoalndedefdfoOructshrieonqtgeuriertnehdae stract/11
2
Thirty-eight Madagascan accessions were included, represent- forphylogeneticanalyses,weexcludedduplicateaccessionsofa /9
ing24species.Weincluded16accessions(14species)fromcon- species if the sequences from each accession were identical /17
tinental Africa, and from west and east tropical Africa. The acrossallsixdatasets. 23
remaining 12 Ixora accessions are Asian (three species), Congruenceofthedatasetswasassessedusingthepartition /27
6
Mascarene (two species), Neotropical (four accessions from homogeneitytestimplementedinPAUP*v.4.0b10(Swofford, 89
three species) and Pacific Island (three species) taxa. 2003).Allconstantanduninformativecharacterswereexcluded. 14
Vangueria madagascariensis, representing the closely related Onethousandpermutationcycleswererun,eachconsistingofa by
tribe Vanguerieae (e.g. Kainulainen et al., 2013) was selected heuristic maximum parsimony (MP) search of ten random se- gu
e
asanoutgroup. quenceadditionreplicateswithTBR(treebisectionandrecon- st o
TotalgenomicDNAwasisolatedfromeithersilicagelcollec- struction) branch swapping, holding ten trees at each step and n
tions, fresh leaf material from the living collections of the savingnomorethanfivetreesperreplicate. 09
NationalBotanicGardenofBelgium(NBGB)orherbariumma- MaximumparsimonyanalyseswereperformedwithPAUP* Ap
terial(BR,MO,P)usingastandardcetyltrimethylammonium v. 4.0b10. We conducted equal weighted parsimony heuristic ril 2
bromide(CTAB)protocol(DoyleandDoyle,1987).Asreported tree searches on: (a) individual data sets; (b) a combined 01
9
elsewhere (e.g. Rajaseger et al., 1997), Ixora leaves may be plastiddataset;and(c)acombinedplastid–nuclearDNAdata
highlycoriaceousandcancontainhighlevelsofphenoliccom- set.Eachanalysisconsistedof1000randomsequenceaddition
poundsthatmayaffectthequalityofisolatedDNA.Therefore, replicates,holdingtentreesateachstep,withTBRbranchswap-
we purified isolated DNA using Nucleospin purification pingandMulTreesineffect,DELTRANoptimizationandsaving
columns (Macherey-Nagel), following the manufacturer’s nomorethantentreesperreplicate.Supportforcladeswaseval-
instructions. uated with 1000 full-heuristic bootstrap pseudo-replicates
(Felsenstein,1985),usingthesamesettingsasoutlinedabove.
Bayesian analyses were implemented in MrBayes 3.1
Amplificationandsequencing
(HuelsenbeckandRonquist,2001).ThemodelofDNAsubstitu-
Primers for amplification of plastid and nuclear ribosomal tion for each region was determined using Modeltest v. 3.06
DNA(nrDNA)regionsarelistedinTable2.PCRandcyclese- (PosadaandCrandall,1998)undertheAkaikeinformationcri-
quencing was performed using a Perkin Elmer GeneAMPw terion(AIC;SupplementaryDataTableS1).Fourindependent
9700 thermocycler. Plastid PCR mixes were made up to Bayesian analyseswithfourchainswere runforeachdataset,
1726 Toshetal.—EvolutionaryhistoryoftheAfro-MadagascanIxoraspecies(Rubiaceae)
TABLE 1. Taxonaccessiondata(onlyfirstcollectorlistedforvoucher)
Taxon Voucher/Herbarium/Countryoforigin ETS ITS petD rps16 rpoB-trnC trnL-trnF
Vangueriamadagascariensis VanCaekenberghe82(BR),Mozambique — — HG315109 HG315176 HG315244 HG315312
J.F.Gmel.
IxoraaluminicolaSteyerm. Pre´vost4160(P),FrenchGuiana HG315378 HG315441 — HG315177 HG315245 —
I.ambrensisDeBlock Gautier5006(BR),Madagascar HG315379 HG315442 HG315110 HG315178 HG315246 HG315313
I.amplidentataDeBlock* DeBlock976(BR),Madagascar HG315380 HG315443 HG315111 HG315179 HG315247 —
I.ankazobensisDeBlock* Tosh30(BR),Madagascar HG315381 HG315444 HG315112 HG315180 HG315248 HG315314
I.ankazobensisDeBlock* Tosh245(BR),Madagascar HG315382 HG315445 HG315113 HG315181 HG315249 HG315315
I.ankazobensisDeBlock* DeBlock943(BR),Madagascar HG315383 HG315446 HG315114 HG315182 HG315250 HG315316
I.ankazobensisDeBlock* DeBlock857(BR),Madagascar HG315384 HG315447 HG315115 HG315183 HG315251 HG315317
I.ankazobensisDeBlock* Tosh256(BR),Madagascar HG315385 HG315448 HG315116 HG315184 HG315252 HG315318 D
o
I.batesiiWernham Dessein1455(BR),Cameroon — HG315449 HG315117 HG315185 HG315253 HG315319 w
n
I.borboniaeMouly& VanCaekenberghe42(BR),Mauritius HG315386 HG315450 HG315118 HG315186 HG315254 HG315320 lo
B.Bremer ad
I.brachypodaDC. Bradley1022(MO),Gabon HG315387 HG315451 HG315119 HG315187 HG315255 HG315321 ed
I.brachypodaDC. Walters1437(MO),Gabon — HG315452 HG315120 HG315188 HG315256 HG315322 fro
I.brevifoliaBenth. Delpretes.n.(BR),Brazil HG315388 HG315453 HG315121 HG315189 HG315257 HG315323 m
I.capuronianaDeBlock* Tosh400(BR),Madagascar HG315389 HG315454 HG315122 HG315190 HG315258 HG315324 h
I.caulifloraMontrouz. Mouly267(P),NewCaledonia HG315390 HG315455 HG315123 HG315191 HG315259 HG315325 ttp
s
I.chinensisLam. VanCaekenberghe316(BR),China — — HG315124 HG315192 HG315260 HG315326 ://a
I.collina(Montrouz.) Mouly236(P),NewCaledonia HG315391 HG315456 HG315125 HG315193 HG315261 HG315327 c
a
Beauvis. d
e
I.crassipesBoivinexDe Groeninckx80(BR),Madagascar HG315392 HG315457 HG315126 HG315194 HG315262 HG315328 m
Block ic
.o
I.cremixoraDrake Mouly659(P),ComoroIslands HG315393 HG315458 HG315127 HG315195 HG315263 HG315329 u
p
I.cremixoraDrake DeBlock987(BR),Madagascar HG315394 HG315459 HG315128 HG315196 HG315264 HG315330 .c
I.densithyrsaDeBlock DeBlock1773(BR),Madagascar HG315395 HG315460 HG315129 HG315197 HG315265 HG315331 om
I.elliotiiDrakeexDeBlock DeBlock1977(BR),Madagascar HG315396 HG315461 HG315130 HG315198 HG315266 HG315332 /a
I.emirnensisBaker DeBlock1786(BR),Madagascar HG315397 HG315462 HG315131 HG315199 HG315267 HG315333 ob
I.emirnensisBaker DeBlock1788(BR),Madagascar HG315398 HG315463 HG315132 HG315200 HG315268 HG315334 /a
I.ferrea(Jacq.)Benth. Merello1716(MO),Commonwealthof HG315399 HG315464 HG315133 HG315201 HG315269 HG315335 rtic
Dominica(LesserAntilles) le
I.ferrea(Jacq.)Benth. Taylor11693(MO),PuertoRico HG315400 HG315465 HG315134 HG315202 HG315270 HG315336 -ab
II..ffoolliioicsaalyHxieGrune´de`s OTonsahna355266(B(PR)),,CMaamdeargoaosncar HHGG331155440012 HHGG331155446667 HHGG331155113356 HHGG331155220034 HHGG331155227712 HHGG331155333378 strac
I.foliicalyxGue´de`s DeBlock696(BR),Madagascar HG315403 HG315468 HG315137 HG315205 HG315273 HG315339 t/1
1
I.franciiSchltr. Mouly241(P),NewCaledonia HG315404 HG315469 HG315138 HG315206 HG315274 HG315340 2
I.guillotiiHoch. DeBlock2091(BR),Madagascar HG315405 HG315470 HG315139 HG315207 HG315275 HG315341 /9/1
I.guillotiiHoch. Tosh408B(BR),Madagascar HG315406 HG315471 HG315140 HG315208 HG315276 HG315342 7
2
I.guineensisBenth. Gereau5601(MO),Ghana HG315407 HG315472 HG315141 HG315209 HG315277 HG315343 3
I.hartianaDeBlock Bamps4320(BR),Angola HG315408 HG315473 HG315142 HG315210 HG315278 HG315344 /27
I.hierniiScott-Elliot Adam23101(P),SierraLeone HG315409 HG315474 HG315143 HG315211 HG315279 HG315345 68
I.hippoperiferaBremek. Dessein1669(BR),Cameroon HG315410 HG315475 HG315144 HG315212 HG315280 HG315346 91
I.hippoperiferaBremek. VanValkenburg3083(WAG),Gabon HG315411 HG315476 HG315145 HG315213 HG315281 HG315347 4 b
I.homolleaeDeBlock& Tosh107(BR),Madagascar HG315412 HG315477 HG315146 HG315214 HG315282 HG315348 y
Govaerts† gu
I.homolleaeDeBlock& Tosh207(BR),Madagascar HG315413 HG315478 HG315147 HG315215 HG315283 HG315349 es
Govaerts† t o
n
I.lagenifructaDeBlock* DeBlock2036(BR),Madagascar HG315414 HG315479 HG315148 HG315216 HG315284 HG315350 0
I.macilentaDeBlock Dessein1404(BR),Cameroon HG315415 HG315480 HG315149 HG315217 HG315285 HG315351 9 A
I.mangabensisDC. Tosh128(BR),Madagascar HG315416 HG315481 HG315150 HG315218 HG315286 HG315352 p
I.mangabensisDC. Tosh130(BR),Madagascar HG315417 HG315482 HG315151 HG315219 HG315287 HG315353 ril 2
I.mangabensisDC. DeBlock2040(BR),Madagascar HG315418 HG315483 HG315152 HG315220 HG315288 HG315354 0
1
I.mangabensisDC. DeBlock2053(BR),Madagascar HG315419 HG315484 HG315153 HG315221 HG315289 HG315355 9
I.masoalensisDeBlock* Razafimandimbison654(BR),Madagascar HG315420 HG315485 HG315154 HG315222 HG315290 HG315356
I.microphyllaDrake DeBlock985(BR),Madagascar HG315421 — HG315155 HG315223 HG315291 HG315357
I.minutifloraHiern Dessein1440(BR),Cameroon HG315422 HG315486 HG315156 HG315224 HG315292 HG315358
I.mocquerysiiDC. Malcomber2805(MO),Madagascar HG315423 HG315487 HG315157 HG315225 HG315293 HG315359
I.moramangensisDeBlock* Tosh255(BR),Madagascar HG315424 — HG315158 HG315226 HG315294 HG315360
I.moramangensisDeBlock* DeBlock837(BR),Madagascar HG315425 — HG315159 HG315227 HG315295 HG315361
I.narcissodoraK.Schum. DeBlock418(BR),Kenya HG315426 HG315488 HG315160 HG315228 HG315296 HG315362
I.nematopodaK.Schum. Dessein1449(BR),Cameroon HG315427 HG315489 HG315161 HG315229 HG315297 HG315363
I.nitens(Poir.)Mouly& Friedmann2631(P),Mauritius HG315428 HG315490 HG315162 HG315230 HG315298 HG315364
B.Bremer
I.perrieriDeBlock* DeBlock841(BR),Madagascar HG315429 HG315491 HG315163 HG315231 HG315299 HG315365
I.perrieriDeBlock* Tosh232(BR),Madagascar HG315430 HG315492 HG315164 HG315232 HG315300 HG315366
I.platythyrsaBaker DeBlock773(BR),Madagascar HG315431 HG315493 HG315165 HG315233 HG315301 HG315367
Continued
Toshetal.—EvolutionaryhistoryoftheAfro-MadagascanIxoraspecies(Rubiaceae) 1727
TABLE 1. Continued
Taxon Voucher/Herbarium/Countryoforigin ETS ITS petD rps16 rpoB-trnC trnL-trnF
I.praetermissaDeBlock Dessein1519(BR),Cameroon HG315432 HG315494 HG315166 HG315234 HG315302 HG315368
I.quadrilocularisDeBlock* Tosh85(BR),Madagascar HG315433 HG315495 HG315167 HG315235 HG315303 HG315369
I.rakotonasoloiDeBlock Tosh316(BR),Madagascar HG315434 HG315496 HG315168 HG315236 HG315304 HG315370
I.regalisDeBlock* DeBlock835(BR),Madagascar HG315435 HG315497 HG315169 HG315237 HG315305 HG315371
I.regalisDeBlock* DeBlock2083(BR),Madagascar HG315436 — HG315170 HG315238 HG315306 HG315372
I.scheffleriK.Schum. Luke9162(P),Tanzania HG315437 HG315498 HG315171 HG315239 HG315307 HG315373
I.siphonanthaOliv. Tosh389(BR),Madagascar HG315438 HG315499 HG315172 HG315240 HG315308 HG315374
I.sp.‘Brunei’ Malcomber2980(MO),Brunei — HG315500 HG315173 HG315241 HG315309 HG315375
I.sp.‘Malaysia’ Billiet7327(BR),Malaysia HG315439 HG315501 HG315174 HG315242 HG315310 HG315376
I.tanzaniensisBridson Luke9304(P),Tanzania HG315440 HG315502 HG315175 HG315243 HG315311 HG315377 D
o
w
n
*sp.nov.ined. lo
a
†nom.nov.ined. d
e
d
fro
m
h
TABLE 2. Amplificationprimersforplastidandnuclearregions ttp
s
Region Primer Primersequence(5’–3’) Reference ://a
c
a
d
e
petD PetB1365 TTGACYCGTTTTTATAGTTTAC Lo¨hneandBorsch(2005) m
PetD738 AATTTAGCYCTTAATACAGG ic
rpoB-trnC rpoB-F CACCCRGATTYGAACTGGGG Shawetal.(2005) .ou
trnC-R CKACAAAAYCCYTCRAATTG p.c
rps16 rps16-F AAACGATGTGGTARAAAGCAAC Shawetal.(2005) o
m
rps16-R AACATCWATTGCAASGATTCGATA /a
trnL–trnF trnL-c CGAAATCGGTAGACGCTACG Taberletetal.(1991) ob
trnF-f ATTTGAACTGGTGACACGAG /a
ETS 18S-ETS GCAGGATCAACCAGGTGACA Negro´n-OrtizandWatson(2002) rtic
ETS-ERIT CTTGTATGGGTTGGTTGGA BaldwinandMarkos(1998) le
ITS ITS1 TCCGTAGGTGAACCTGCGG Whiteetal.(1990) -ab
ITS4 TCCTCCGCTTATTGATATGC stra
c
t/1
1
2
/9
startingfromrandomtrees,for5milliongenerations,sampling assessed the distribution of these morphological characters by /17
trees every 1000 generations. TRACER v. 1.4 (Rambaut and mapping unambiguous character state changes onto our com- 23
Drummond,2007)wasusedtoassessifthesearchhadreached bined Bayesian MJ consensus tree using MacClade v. 4.04 /27
6
stationarity and to check each parameter had an effective (MaddisonandMaddison,2002). 89
1
sample size (ESS) .100. The initial 1250 (25%) trees were 4
b
discardedasaconservativeburn-in.Post-burn-intreesfromthe y
Divergencetimeestimation g
four independent analyses were pooled and summarized by a u
e
50% majority rule consensus tree using PAUP* v4.0b10 to ABayesianapproachwasappliedtoinferthetemporalframe- st o
obtainposteriorprobabilities(PPs). workoftheevolutionofIxora.Duetothelimitedfossilrecord n
0
that could be unequivocally assigned to the most recent 9
A
common ancestor of Ixora, an expanded family-level data set p
Morphologyandoptimizationofmorphologicalcharacters was constructed and divergence time estimation was inferred ril 2
0
Thematerialforthetaxonomicalandmorphologicalstudyof based on fossils (see the following section for more detail). 1
9
IxorainMadagascarconsistsofpreservedsamplesandherbar- This large-scale analysis allowed estimation of the temporal
iummaterialofthefollowinginstitutions:BM,BR,G,K,MO, originofIxoratogetherwitha95%confidenceinterval,which
P,S,TAN,TEF,UPS,W,WAGandZ(abbreviationsofinstitu- wassubsequentlyusedasapriortoperformasecondanalysisfo-
tionsfollowHolmgrenetal.,1990).Intotal,.1000herbarium cusing only on this genus. The family-level data set included
collections, each with several duplicates, were studied. representatives of all the major lineages in Rubiaceae (see
Morphological terminology generally follows Robbrecht Appendix 1) and was based on the group II plastid introns
(1988).FortheMadagascanspecies,morphologicalcharacters petD, rps16 and the trnL-trnF spacer (Groeninckx, 2009). In
were scored on herbarium material available for this species. the large-scale analysis, a sub-set of taxa representative of the
For the continental African species, morphological characters main clades in Ixora was included. These taxa were selected
were taken from the revision of African Ixora (De Block, basedonpreliminaryphylogeneticanalyses.Weoptedforthis
1998).Wechosemorphologicaltraitsofinterest(i.e.uniflorous approachratherthanpoolingallthedatasetsandtaxatogether
inflorescences,four-locularovaries,extremecorollatubelength) mainlytoavoidencounteringproblemsrelatedtomissingdata.
andcharactersofpotentialtaxonomicsignificance(Table3).We In addition to the group II plastid introns, the Ixora data set
1
7
2
8
TABLE 3. Selectedmorphologicalcharacters
Inflorescence Calyx Flower D
o
w
Taxon Sessile/pedunculate Lax/compact Tubelength(mm) Lobelength(mm) Number Corollatubelength(cm) Ovary T n
o lo
s a
Africantaxa I.batesii Sessile Lax 0.4–0.5 0.4–0.6 30–70 0.5–1.3 Bilocular h de
e d
II..fborlaicohsaypoda PPeedduunnccuullaattee CLaoxmpact 00..34––00..58 00..43––01..70 5045––29000 03..93––21.13 BBiillooccuullaarr tal. from
I.guineensis Sessile Lax 0.3–0.8 0.2–0.6 30–90 0.5–2.3 Bilocular — h
I.hartiana Pedunculate Lax 0.2–0.5 0.3–0.6 9–30 1.2–2.4 Bilocular E ttp
III...mhhiiapecrpnioliepinetraifera PSSeuedsbsu-isnleecsuslialeteorsessile LLCaaoxxmpact 000...222–––000...475 000...222–––000...456 30390–––1922000 110...649–––232...732 BBBiiilllooocccuuulllaaarrr volutio s://acad
I.minutiflora Sessile Compact 0.2–0.4 0.1–0.4 9–30 1.0–2.0 Bilocular na em
III...pnnreaamrecatiestosrompdoiodssraaa SPSeeedssssuiinlleeculate LLLaaaxxx 000...238–––001...565 000...382–––010...826 503200–––1750000 103..,104–––027..,955 BBBiiilllooocccuuulllaaarrr ryhis ic.oup
II..stacnhzeafflneiernisis SPeedssuinleculate CCoommppaacctt 00..34––00..76 00..26––01..40 2300––5900 20..28––32..22 BBiillooccuullaarr tory .com
o /a
Madagascantaxa(Clade2) II..admenpsliitdheynrtsaata PPeedduunnccuullaattee LCaoxmpact 00..25––01..40 02..55––41..025 509––16200 11.58––23.33 BBiillooccuullaarr fthe ob/a
I.foliicalyx Sub-sessile Compact 2.0–4.0 1.5–6.0 45–90 3.5–8.5 Bilocular A rtic
II..hgoumillooltlieiae (PSeudbu-n)csueslastiele LUanxiflorous 05..05––11.00 50..075––122.5 50–1510 25..60––48..70 FBoiluorc-ulolacrular fro-M le-ab
IIIIIII.......pmmlsrqaeilupoiaggcahcteardyqonoltriunihpsifealyrhornurycystclushaltialaaiaris PP(PPPPSeeeeeeuddddddbuuuuuu-nnnnnn)ccccccsuuuuuuesllllllaaaaaasttttttieeeeeele LCLLLLLaaaaaaoxxxxxxmpact 0050040.......2550550–––––––0128117.....(00004–1.5) (1320130––.....055095)5(–––––––61211911....0504)5 (1553001–(332––)–(3––11–1–48259590050)) 532114(1......5253040–––––––146582).....1105505(5––1262) BBFBBFBooiiiiillllluuooooorrccccc--uuuuulloolllllaaaaaccrrrrruullaarr Ixoadagascan stract/112/9/1723
Madagascantaxa(Clade3) I.ambrensis Pedunculate Lax 0.4–0.5 0.2–0.35 25–50 3.2–3.9 Bilocular ra /27
s 6
I.ankazobensis Sessile(+pedunculate) Lax 0.5–0.75 0.3–0.75 7–18 4.0–6.3 Bilocular p 8
e 9
I.capuroniana Sessile Laxorcompact 0.3–0.6 0.1–0.8 40–90 1.5–3.5 Bilocular c 1
I.crassipes Sub-sessile Lax 0.75–1.25 0 15–50 17–22.5 Bilocular ies 4 b
I.cremixora Sessile Lax 0.2–0.75 ,0.3 50–120 4.5–7.8 Bilocular (R y g
IIII....mmeemlaaliisnorogtniaaeilbneesninsssisis SPPSeeeeddssssuuiinnlleeccuullaattee LLLLaaaaxxxx 0000....5222––––2000....0534 0000....5533––––1100....5045 38990––––23595000 2211....2008––––4313....0251 BBBBiiiillllooooccccuuuullllaaaarrrr ubiacea uest on 0
I.moramangensis Pedunculate Lax 0.25–0.5 0.25–0.8 9–30 2.1–3.5 Bilocular e) 9 A
II..rpaekroriteorniasoloi SPeedssuinleculate ULanxiflorous 00..33––00..55 00..32––00..54 30–910 24..25––27..58 BBiillooccuullaarr pril 2
0
1
9
AdaptedfromDeBlock(1998,2013).
Toshetal.—EvolutionaryhistoryoftheAfro-MadagascanIxoraspecies(Rubiaceae) 1729
included one additional plastid (rpoB-trnC) and two nuclear basedonfossilrecordsofpollenfromtheOligocene(Graham,
(ETSandITS)regions.Theseregionsprovedusefultoresolve 2009);(d)thestemoftheGaliumandRubiacladewasassigned
phylogeneticrelationshipsinIxorafurther. anoffsetof5.3Ma,followingthefirstoccurrenceoffossilpollen
With the exception of the numberof runs and length of the duringthelateMioceneforthisgroup(Graham,2009);(e)the
Monte Carlo Markov chain (MCMC), the settings of the stem of the Chiococca clade was calibrated with an offset of
Bayesiananalyseswereidenticalforthetwodatasets.Aparti- 5.3MabasedonleafmaterialfromthelateMiocene(Graham,
tioned Bayesian MCMC analysis was conducted in BEAST 2009);(f )thestemoftheEmmenopteryscladewassetwithan
v.1.7witharelaxedlog-normalmolecularclockandaYulespe- offset of 48Ma, following the first record of Rubiaceae fruits
ciationmodel.Thepartitionswereunlinkedforthemodelofevo- in the mid Eocene (Graham, 2009); and (g) an offset of
lution,butlinkedfortheestimationofthemolecularclockand 14.5Ma was assigned to the Gardenia clade, based on pollen
thetreetopology.Theotherparametersweresetasdefault(see datafromthemidMiocene(Graham,2009).Basedontheposter-
Drummond and Rambaut, 2007). The best-fit models foreach iordistributionofthedatinguncertaintyonthestemofIxora,the D
o
DNAregionwerekeptidenticaltothoseintheMrBayesanalyses most recent common ancestor of the ingroup in the second w
n
(seeabove).Inthecaseofthefamily-leveldataset,fourrunsof BEASTanalysiswassetwithanormaldistribution,ameanof lo
a
10 million generations were performed, sampling a tree every 18Maandans.d.of2.5. de
d
1000 generations. To improve the tree topology research, the fro
MrBayesconsensustreewasprovidedasthestartingtreeinthe m
BEASTfamily-levelanalysis.InthecaseoftheIxoradataset RESULTS http
(which was rooted using Vangueria infausta), two runs of 5 s
million generations were performed, sampling a tree every The368novelsequencesgeneratedinthisstudywerecombined ://a
1000generations,andtheanalyseswereseededwitharandom with31sequencespreviouslygeneratedbyMouly(2007),result- cad
inginatotalcombineddatasetof399sequences,representing e
tree. For each parameter, convergence of runs was confirmed m
bTyRtAhCeEexRamv.i1n.a4tio(RnaomfbthaeuitrarnedspDecrtuimvempoonsdte,r2io0r0d4i)s.tIrnibaudtidointisoinn, asppepcroiexs.w5e0reIxgoernaersapllpy.lLoweveaclsroossf aglelnseixticrevgaiorinastioinnvebsetitgwaeteedn ic.oup
(Table4).Thetotalnumberofpotentiallyparsimonyinformative .c
we considered the MCMC sampling sufficient when the ESS o
charactersrangedfrom21inrps16andtrnL-trnFto58intheITS. m
was .200 using TRACER v. 1.4 (Rambaut and Drummond, In terms of percentage variability,the ETS proved tohave the /ao
2004). A maximum clade credibility tree with median branch b
highestproportionofpotentiallyparsimonyinformativecharac- /a
olTenrnegtehtAhesnsnaeonttdoaft9ot5rrev%e.s1ca.o5fnt.e4fird(bDeunrrucnem-iimnnto(efnrovdraalenaodcnhRnraoumdneb,saauwbta,us2rn0b-0ui7nil)tpbueasrsiinoeddg tpeertsD(1(32..43%%))..TOhuerltoawxoenstspaemrcpelnintaggiencvlaurdiaebdilmituylwtipalseoabcsceersvseidonins rticle-a
for 13 Ixora spp. We observed no intraspecific sequence vari- bs
of1millionwasapplied). ation in nine of these species. In contrast, varying amounts of tra
c
intraspecific sequence divergence were observed in the other t/1
1
Fossilcalibration five species; multiple accessions from these five species were 2
/9
includedinallsubsequentphylogeneticanalyses. /1
Toestimateabsoluteagesforlineagedivergencesinthecase 7
2
ofthefamily-leveldataset,wehaveusedsevenfossilcalibration 3
/2
points(seebelow)inRubiaceaetosetageconstraintsforseveral 7
Phylogeneticanalyses 68
nodes.AsinBuerkietal.(2011),foreachcalibrationpoint,the 9
1
oldest fossil record was selected and the upper (younger) Plastiddatasets.Phylogeneticanalysesofindividualplastiddata 4
b
bound of the geological interval (Gradstein et al., 2004) in sets(petD,rpoB-trnC,rps16,trnL-trnF)generatedlargelyunre- y
g
which the fossil was foundwasused torepresent theage con- solvedandpoorlysupportedphylogenetictrees(datanotshown). u
e
s
straint.Inthefamily-levelBayesiananalysis,allthecalibration Thepartitionhomogeneitytestdidnotrevealanysignificantin- t o
points were modelled as follows: log-normal distribution, congruencebetweenplastiddatasets(P¼0.69)andtherewas n
0
mean¼0.5,s.d.¼1,offset¼agefossil(seebelow). nosupportedincongruence(BayesianPP.0.95),sothesefour 9
A
Beforedetailingthefossilsusedhere,wewouldliketoreport regionswerecombinedforallsubsequentanalyses. p
thecurrentknowledgeonfossilrecordsthathavebeententatively Thecharacteristicsoftheindividualandcombinedplastiddata ril 2
0
assignedtoIxora.PalaeobotanicalremainsofIxoraarepoorly setsarelistedinTable4.InboththeMPandBayesiananalyses 1
9
known, and Graham (2009) reported the presence of Ixora (Fig.2),theAsianandPacificOceantaxaaresistertotherest
pollen dating from the Miocene on the Marshall Islands of the genus [bootstrap (BS) 70, PP 1.00]. The Mascarene
(Micronesia).Inhisstudy,Graham(2009)stressedthatthisevi- Ixora clade is in turn sister to the Neotropical and
dencecouldbechallengedandthatfurtherstudieshavetobecon- Afro-MadagascanIxoraspp.(BS77,PP1.00).Thesisterrela-
ducted to assign these pollen unequivocally to existing taxa. tionship between the three Neotropical Central–South
Therefore,wehavenotconsideredthisrecordinourdivergence AmericantaxaandtheAfro-Madagascantaxaisweaklytomod-
timeestimations.Thefollowingcalibrationpointswereselected: eratelysupportedintheMP(BS56)andBayesian(PP0.96)ana-
(a) the stem group of Rubiaceae was set with an offset of 54 lyses (Fig. 2). In the Afro-Madagascan group, we identified a
millionyears(Ma),followingthefirstoccurrenceofRubiaceae clade containing African and Madagascan taxa (BS 72, PP
fossils during the Eocene (Malcomber, 2002); (b) the stem of 1.00; Clade 1) and an exclusively Madagascan clade (BS 59,
theFarameacladewassetwithanoffsetof40Ma,basedondis- PP0.99;Clade2).Thereisalackofresolutionandweaknode
tinctivediporatepollenfromthelateEocene(Graham,2009);(c) supportinClades1and2duetoapaucityofpotentiallyparsi-
theCoprosmacladewasconstrainedwithanoffsetof23.8Ma mony informative characters (Table 4). In Clade 1 there is
1730 Toshetal.—EvolutionaryhistoryoftheAfro-MadagascanIxoraspecies(Rubiaceae)
TABLE 4. Characteristicsofindividualandcombineddatasetsandtreestatistics
petD trnL-trnF rps16 rpoB-trnC Plastid ETS ITS nrDNA Total
Numberoftaxa 56 57 57 55 57 52 52 49 57
Totallength(bp) 1023 828 722 1016 3589 409 631 1040 4629
Non-parsimonyinformative 49 41 46 75 211 51 71 122 333
characters
Potentiallyparsimonyinformative 24(2.3%) 21(2.5%) 21(2.9%) 33(3.2%) 99(2.7%) 55(13.4%) 58(9.2%) 108(10.4%) 212(4.6%)
characters(%oftotal)
Indels 3 2 1 3 9 0 0 0 9
Treelength 81 68 76 140 376 159 196 357 758
Consistencyindex 0.914 0.956 0.921 0.850 0.872 0.780 0.791 0.759 0.801
Retentionindex 0.955 0.969 0.917 0.811 0.889 0.898 0.819 0.839 0.849 D
o
Numberoftreessaved 100 68 511 4887 8722 1151 9580 3464 5362 w
n
lo
a
d
e
d
fro
strong support for a clade of Madagascan and tropical East 1 is supported in both the Bayesian and MP analyses (BS 73, m
Africantaxa(BS85,PP1.00;Clade3). PP1.00)followingtheexclusionofI.nematopoda(notshown). h
ttp
In Clade 1, the group of I. foliosa and I. hippoperifera is s
nrDNA data sets. We were unable to obtain ETS and ITS weaklysupported(BS57,PP0.92),buttwoothergroupsoftrop- ://a
c
sequences for Vangueria madagascariensis, due to amplifica- icalWest–CentralAfricantaxaarewellsupported:I.guineensis, ad
tion difficulties and ambiguous sequence reads, respectively. I. minutifolia, I. hiernii and I. batesii (BS 92, PP 1.00) and em
As our primary interests concern the Afro-Madagascan I.hartiana,I.macilentaandI.praetermissa(BS85,PP1.00). ic.o
elementofIxora,phylogeneticanalysesfromseparateandcom- Supportforthesisterrelationshipsbetween thesethreeWest– up
bined nuclear DNA analyses were rooted using Asian and/or Central African clades, and between the widespread .co
m
Pacific Ocean Ixora spp. Relationships in the Afro- I. brachypoda, is weak or lacking (Fig. 3). Branch lengths /a
Madagascan Clade were not affected by the root choice. The between these clades are short, which may be a contributing ob
chaTrhaectetoripsotilcosgoyfothfethneucElTeaSrDphNyAlogdeatnaestiectstraereel(iSstuepdpilnemTaebnltear4y. fCaecntotrralinAftrhiecapnhsypleocgieense(tSicupupnlecmeretanitnatryyDbeattwaeFeing.tSh3e)s.e West– /article
DataFig.S1)issimilartothatofthecombinedplastidtree,but NestedinClade1isthestronglysupportedClade3(BS100, -a
b
therearesomedifferences.Clades1and2arenotfullyrecovered PP 1.00), comprised of East African and Madagascan species stra
in the ETS topology, with the deepest nodes in the (Fig. 3). The branch length subtending Clade 3 is relatively c
Afro-Madagascan group being unresolved. Furthermore, there long, with eight synapomorphies (Supplementary Data Fig. t/11
2
isweaktomoderatesupport(BS78,PP0.93)forthesisterrela- S3).InClade3,theEastAfricanI.narcissodoraandI.tanzanien- /9
tionshipoftheGuineo-CongolianI.nematopodawithrespectto sis are sister to the Madagascan taxa (BS 66, PP 0.94). /17
2
Neotropical Ixora and the remaining Afro-Madagascan Ixora PhylogeneticrelationshipsbetweentheseMadagascantaxaare 3
/2
(discussedbyMoulyetal.,2009b).However,Clade3isstrongly poorly resolved (Fig. 3; Supplementary Data Fig. S3). 7
6
supportedintheETStopology(BS94,PP1.00).TheITStop- However,thereareafewspeciesrelationshipsthatarewellsup- 89
ologyispoorlyresolvedandweaklysupported,inparticularin ported,suchasI.elliotiiandI.rakotonasoloi (BS99,PP1.00) 14
the MP strict consensus (Supplementary Data Fig. S2). Few andI.crassipesandI.cremixora(BS98,PP1.00).Theacces- by
g
clades are well supported in both MP and Bayesian analyses, sionsofI.mangabensisdidnotgrouptogetherinouranalyses, u
e
andClades1–3werenotrecovered. due to the presence of two distinct plastid haplotypes in the st o
fouraccessionsofthisspecies. n
0
Totalcombineddataset.Despitesometopologicalinconsisten- The monophyly of Clade 2, comprised exclusively of 9 A
ciesbetweenthecombinedplastidandnuclearDNAdatasets, Madagascan taxa, is supported in the combined plastid– p
results from the simultaneous partition homogeneity test did nuclearDNAdataset(BS70,PP1.00).Therearetwomainsub- ril 2
0
not reveal significant incongruence between data sets (P¼ cladesrecoveredintheMPstrictconsensusandBayesianmajor- 1
9
0.10).Thecombinedplastid–nuclearDNAdatasetcomprised ityruleconsensustrees,althoughsupportforthesesub-cladesis
4629characters.Intotal,212ofthe545variablecharactersare lacking (Fig. 3). Although the partition homogeneity test
potentiallyparsimonyinformative(Table4). revealed no significant incongruence between individual data
The Mascarene clade is sister to the Neotropical and sets, the phylogenetic placement of certain taxa (e.g.
Afro-Madagascanspecies(BS84,PP1.00).Themonophylyof I. densithyrsa and I. regalis 2) differed between plastid and
theAfro-Madagascangroupissupported(BS65,PP1.00),and nucleardatasets.Intheplastidphylogenetictree,I.regalis2is
in this group there is support for Clades 1–3 (Fig. 3). In the nested in the group containing I. regalis 1, I. guillotii and
Bayesian analysis, I. nematopoda and I. scheffleri are sister to I. amplidentata (Fig. 2). In contrast, I. regalis 2 is unresolved
therestofClade1(PP1.00).Althoughthisrelationshipisrecov- in the ETS phylogenetic tree (Supplementary Data Fig. S1)
eredintheMPstrictconsensus,itisnotsupportedbythebootstrap andsistertoI.homolleaeandI.quadrilocularisintheITSphylo-
analysis. The phylogenetic position of I. nematopoda differs genetic tree (Supplementary Data Fig. S2). Similarly,
between the ETS and the plastid and ITS data sets (Fig. 2; I. densithyrsa andI. mocquerysii demonstratedifferingphylo-
SupplementaryDataFigsS1andS2).ThemonophylyofClade genetic affinities in the plastid and nuclear data sets (Fig. 2;
Toshetal.—EvolutionaryhistoryoftheAfro-MadagascanIxoraspecies(Rubiaceae) 1731
V. madagascariensis
0·73/* I. aluminicola
1·00/56 I. ferrea Neotropical lxora
I. brevifolia
I. amplidentata
0·62/* I. guillotti
1·00/78 I. regalis 1
I. regalis 2
I. lagenifructa
1·00/75 II.. dheonmsoitlhleyarsea 1 e 2 Madagascan lxora
0·91/* d
0·99/59 0·98/60 II.. hqoumadorlilleoaceu la2ris Cla
I. mocquerysii D
o
1·00/77 I. foliicalyx w
n
1·00/71 I. microphylla lo
I. platythyrsa ad
I. siphonantha ed
I. ankazobensis 1 fro
1·00/85 I. perrieri m
I. ankazobensis 2 h
I. ankazobensis 3 ttp
s
1·00/86 I. capuroniana ://a
1·00/64 I. masoalensis ca
I. emirnensis d
e
0·96/56 1.00/85 1·00/62 III... mccrraaesnmsgiixapoberesan 1sis 1 ade 3 MEaasdta Agfarisccaann l xaonrda mic.ou
I. cremixora 2 Cl p.c
1·00/96 I. elliotii om
I. rakotonasoloi /a
o
I. mangabensis 2 b
1·00/70 0·99/52 1·00/99 III... manmaorrbcarismesnaosndisgoeransis Clade 1 /article-a
I. tanzaniensis bs
I. batesii tra
1·00/67 1·00/66 II.. ghuieinrneieinsis ct/11
2
I. minutiflora /9
1·00/100 I. brachypoda 1 /1
7
I. brachypoda 2 West and Central 23
1·00/100 1·00/72 0·86/54 II.. fhoipliopsoaperifera African lxora /276
8
0·91/* I. hartiana 9
1·00/57 I. macilenta 14
I. praetermissa b
y
0·70/* I. nematopoda g
u
I. scheffleri e
s
1·00/100 I. nitens Mascarene lxora t o
I. borboniae n 0
I. cauliflora 9
1·00/100 I. collina Pacific lxora A
p
0·99/60 0·59/50 II.. fcrhainnceinisis ril 20
1·00/97 I. sp. Malaysia Asian lxora 19
I. sp. Brunei
FIG.2. CombinedplastidBayesianmajorityruleconsensustree.Bayesianposteriorprobabilities.0.5andbootstrapvalues.50%areindicatedabovebranches(PP/
BS).Asterisks(*)denotenodesthathavebootstrapsupport,50intheMPanalysis.Clade1,‘Afro-Madagascanclade’;Clade2,‘Madagascanclade’;Clade3,‘East
African–Madagascanclade’.
Supplementary Data Figs S1 and S2). Exclusion of I. regalis I. regalis 1, I. guillotii, I. amplidentata, I. lagenifructa,
2andI.densithyrsa(datanotshown)fromtheBayesiananalyses I.siphonanthaandI.mocquerysii(PP0.95;notshown).
results in increased support forthe twosub-cladesof Clade 2,
with both sub-clades supported by a PP of 1.00. In the
Distributionofkeymorphologicalcharacters
first subclade, I. foliicalyx is sister to I. homolleae and
I.quadrilocularis(PP1.00;notshown).Inthesecondsub-clade, The phylogenetic distribution of key morphological charac-
I. microphylla and I. platythyrsa are sister to the group of tersisillustratedinFig.4.Pedunculateandsessileinflorescences
1732 Toshetal.—EvolutionaryhistoryoftheAfro-MadagascanIxoraspecies(Rubiaceae)
V. madagascariensis
1·00/84 I. aluminicola
0·64/* I. brevifolia Neotropical lxora
I. ferrea
0·92/* I. amplidentata
0·88/* I. lagenifructa
0·56/* I. guillotii
0·66/* I. regalis 1
0·67/* I. densithyrsa 2
0·69/* I. mocquerysii e Madagascan lxora
I. siphonantha d
1·00/85 I. microphylla a
1·00/70 I. platythyrsa Cl
I. foliicalyx
0·72/* 1·00/95 I. homolleae 1
1·00/84 1·00/99 I. homolleae 2
0·54/* I. quadrilocularis D
I. regalis 2 o
1·00/54 0·94/* 1·00/83 III... aaannnkkkaaazzzooobbbeeennnsssiiisss 123 wnloa
d
I. perrieri e
0·99/* 1·00/98 II.. cmaapsuoraolneinasnias d fro
I. emirnensis 3 m
0·94/66 1·00/98 II.. mcraasnsaigpaebsensis 1 de MEaasdta Agfarisccaann l xaonrda http
1·00/86 I. cremixora 1 a s
1·00/65 I. cremixora 2 Cl ://a
1·00/99 I. elliotii c
1·00/100 I. rakotonasoloi ad
I. mangabensis 2 1 e
1·00/78 0·55/* 0·60/* 1·00/851·10·00/01/0808 IIIIII...... manthmamaaoanrrrbczctaiiarialsmeennsnniaaoestndanisgoserisansis Clade mic.oup.com
I. praetermissa /a
0·55/* I. batesii ob
1·00/100 1·00/890·81/53 1·00/92 0·54/* III... hmgiueinirnunetiieiflnosrais WAferiscta ann lxdo Craentral /article
1·00/100 I. brachypoda 1 -a
I. brachypoda 2 b
1·00/* 0·92/57 II.. fhoipliopsoaperifera stra
0·77/* II.. nscehmeafftloeprioda ct/11
1·00/100 1·00/100 III... nbcaioteurbnliofslonriaae Mascarene lxora 2/9/17
0·67/60 I. collina Pacific lxora 23
0·92/61 I. francii /2
0·89/59 I. chinensis 7
1·00/97 I. sp. Malaysia Asian lxora 68
I. sp. Brunei 9
1
4
FIG. 3. Combinedplastid–nuclearBayesianmajorityruleconsensus.Bayesianposteriorprobabilities.0.5andbootstrapvalues.50%areindicatedabove by
branches.Asterisks(*)denotenodesthathavebootstrapsupport,50intheMPanalysis.Clade1,‘Afro-Madagascanclade’;Clade2,‘Madagascanclade’;Clade g
3,‘EastAfrican–Madagascanclade’. ue
s
t o
n
0
occur in African taxa and both clades of Madagascan the Miocene–Pliocene boundary. Based on the three-gene 9
A
taxa (Fig. 4A). All taxa from Clade 2 and the African Rubiaceae-wide data set (Supplementary Data Fig. S4, Table p
I. nematopoda and the Madagascan I. masoalensis (Clade3) S2),thecrownageforIxoreae(node152)isestimatedat16.67 ril 2
possess calyx lobes .1mm long (Fig. 4B). Calyx tubes millionyearsold(9.67–27.55,95%highestposteriordensity; 01
9
.1mmlongarefoundinallbuttwospeciesfromClade2,in hereafter HPD). The subsequent age estimates (discussed
addition to the African I. narcissodora and the Madagascan below)arebasedontheresultsofthesecondarydatinganalysis
I.masoalensisandI.crassipes(Fig.4C).Specieswithuniflorous (Fig.5;SupplementaryDataTableS3).Theonsetofdivergence
inflorescences or corolla tubes .15cm long occur in both betweentheAsian-PacificOceanIxora(node59)andtherestof
Madagascan clades, but species with four-locular ovaries are the genus is estimated at 15.37Ma (7.39–22.89, 95% HPD).
onlyfoundinClade2(Fig.4D–F). The estimated age of divergence between the Neotropical and
Afro-Madagascan taxa (node 66) is 9.51Ma (4.47–14.94, 95
% HPD), with the crown age for the Afro-Madagascan group
Divergencetimeestimation (node 67) estimated at 7.95 million years old (3.71–12.52, 95
Theestimateddivergencetimesfornodesofinterestaresum- % HPD). The crown age of Clade 1 (node 68) is estimated at
marizedinFig.5andSupplementaryDataFig.S4,TablesS2and 7.22 million years old (3.36–11.46, 95% HPD). The crown
S3.Theonsetofdiversificationofthegenusstartedduringthe ageofClade2(node98)isestimatedat6.24millionyearsold
mid Miocene, with the emergence of most of the lineages at (2.88–10.03, 95% HPD), with two separate periods of