Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
Gourds afloat: a dated phylogeny reveals an Asian origin of
the gourd family (Cucurbitaceae) and numerous oversea
dispersal events
Hanno Schaefer, Christoph Heibl and Susanne S Renner
Proc. R. Soc. B 2009 276, 843-851
doi: 10.1098/rspb.2008.1447
Supplementary data
"Data Supplement"
http://rspb.royalsocietypublishing.org/content/suppl/2009/02/20/276.1658.843.DC1.ht
ml
References
This article cites 35 articles, 9 of which can be accessed free
Subject collections
Articles on similar topics can be found in the following collections
http://rspb.royalsocietypublishing.org/content/276/1658/843.full.html#ref-list-1
taxonomy and systematics (58 articles)
ecology (380 articles)
evolution (450 articles)
Email alerting service
Receive free email alerts when new articles cite this article - sign up in the box at the top
right-hand corner of the article or click here
To subscribe to Proc. R. Soc. B go to: http://rspb.royalsocietypublishing.org/subscriptions
This journal is © 2009 The Royal Society
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
Proc. R. Soc. B (2009) 276, 843–851
doi:10.1098/rspb.2008.1447
Published online 25 November 2008
Gourds afloat: a dated phylogeny reveals an
Asian origin of the gourd family (Cucurbitaceae)
and numerous oversea dispersal events
Hanno Schaefer*, Christoph Heibl and Susanne S. Renner
Systematic Botany, University of Munich, Menzinger Strasse 67, 80638 Munich, Germany
Knowing the geographical origin of economically important plants is important for genetic improvement
and conservation, but has been slowed by uneven geographical sampling where relatives occur in remote
areas of difficult access. Less biased species sampling can be achieved when herbarium collections are
included as DNA sources. Here, we address the history of Cucurbitaceae, one of the most economically
important families of plants, using a multigene phylogeny for 114 of the 115 genera and 25 per cent of the
960 species. Worldwide sampling was achieved by using specimens from 30 herbaria. Results reveal an
Asian origin of Cucurbitaceae in the Late Cretaceous, followed by the repeated spread of lineages into the
African, American and Australian continents via transoceanic long-distance dispersal (LDD). North
American cucurbits stem from at least seven range expansions of Central and South American lineages;
Madagascar was colonized 13 times, always from Africa; Australia was reached 12 times, apparently always
from Southeast Asia. Overall, Cucurbitaceae underwent at least 43 successful LDD events over the past
60 Myr, which would translate into an average of seven LDDs every 10 Myr. These and similar findings
from other angiosperms stress the need for an increased tapping of museum collections to achieve extensive
geographical sampling in plant phylogenetics.
Keywords: Bayesian molecular clock; biogeography; dispersal; economic plants; museomics;
test of monophyly
1. INTRODUCTION
Molecular clock analyses suggest that the majority of
lineages of legumes that occur on islands are younger than
30 Myr (Lavin & Beyra Matos 2008) and that plant
diaspores from source areas hundreds or thousands of
kilometres away regularly reach isolated Arctic islands
(Alsos et al. 2007), island-like mountains in Eastern Africa
and mountain ranges in the Northern Cape, South Africa
(Galley et al. 2007). Striking dispersal events have also been
documented for the flora of Hawaii ( Wagner et al. 1990),
the montane region of New Zealand ( Winkworth et al.
2005) and many other island systems. Such frequent longdistance dispersal (LDD) implies that long-established
views on the origin of economically important plants may
need to be re-evaluated based on drastically enlarged
geographical sampling. An example is the origin of
Cucumis sativus, the cucumber. Cucumber ranks among
the top 10 vegetables in world production (Chen et al.
2004). Until 2006, it was thought that the genus Cucumis
had 32 species and was essentially African. Only C. sativus
and C. hystrix were thought to occur naturally in India,
China, Burma and Thailand (Ghebretinsae & Barber
2006). However, broader geographical species sampling
revealed that C. sativus is closer to 13 species from
Australia, India, Yunnan and Indochina than to any
African species (Renner & Schaefer 2008).
Biogeographic inference for economically important
plants is complicated by human transport of seeds
between continents for at least 10 000 years (Smith
1997; Sanjur et al. 2002; Dillehay et al. 2007). The extent
of the anthropogenic transfer, however, is difficult to work
out without comprehensive phylogenetic frameworks,
which can be prohibitively expensive if worldwide
collecting of material is required. In the economically
important plant family Cucurbitaceae, these difficulties
have led to the geography of the closest relatives of
watermelon (Citrullus lanatus), cucumber (Cucumis sativus), loofah (Luffa acutangula), bitter gourd (Momordica
charantia), chayote (Sechium edule), ivy gourd (Coccinia
grandis), snake gourd (Trichosanthes cucumerina) and
creeping cucumber (Melothria pendula) remaining ambiguous. Natural LDD of cucurbit diaspores may be frequent
because many are adapted for transport by birds or wind,
or they can withstand long periods in water (Cayaponia,
Fevillea, Hodgsonia, Lagenaria, Luffa and Sicana; Ridley
1930; Whitaker & Carter 1954).
Here, we use worldwide sampling, based on museum
specimens, to infer the biogeographic history of Cucurbitaceae, a family consisting of climbers or trailers of tropical
and subtropical regions that are typically strongly
seasonal, lacking aboveground parts during part of the
year. These traits have caused cucurbits to be undercollected (Gentry 1991), resulting in dozens of species still
known from only one or two collections even in the world’s
leading herbaria (e.g. De Wilde & Duyfjes 2007). Of the
approximately 960 accepted species of Cucurbitaceae,
approximately 40 per cent are endemic in the American
* Author and address for correspondence: Imperial College London,
Silwood Park campus, Ecology and Evolutionary Biology, Buckhurst
Road, Ascot, Berkshire SL5 7PY, UK (hanno.schaefer@imperial.ac.uk).
Electronic supplementary material is available at http://dx.doi.org/10.
1098/rspb.2008.1447 or via http://journals.royalsociety.org.
Received 6 October 2008
Accepted 4 November 2008
843
This journal is q 2008 The Royal Society
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
844 H. Schaefer et al. Biogeography of the Cucurbitaceae
continent, and the remainder occur in Africa (28%),
Asia (26%), Australia (2%) and Europe (1%; Schaefer &
Renner in press).
Based on chloroplast sequences from all but one of
the 115 genera and 25 per cent of the 960 species, and
employing specimens from 30 herbaria and up to 172 years
old, we address here the following questions: did Cucurbitaceae initially diversify in Asia, in America, or in Africa and
Madagascar? What are the geographical sources of the
world’s major Cucurbitaceae floras? Are the transoceanic
geographical ranges in the genera Cayaponia, Lagenaria,
Luffa and Sicyos anthropogenic or the result of natural
LDD? We also use Cucurbitaceae to illustrate the still
barely tapped potential of museum collections to achieve
less biased geographical sampling than has traditionally
been employed in tropical plant phylogenetics.
2. MATERIAL AND METHODS
(a) Taxon sampling
We generated 126 sequences, representing 32 species from
seven genera not sampled in previous studies (Anangia,
Cucumeropsis, Gomphogyne, Hodgsonia, Papuasicyos, Pseudosicydium and Zanonia). GenBank accession numbers (EU436320–
EU436422) and vouchers for newly sequenced taxa are listed in
table 1 in the electronic supplementary material. Accession
numbers and voucher information for additional Cucurbitaceae sequences from our earlier studies are given in Zhang et al.
(2006), Kocyan et al. (2007), Schaefer et al. (2008a) and Nee
et al. (submitted). Fourteen sequences from Genbank were
included to represent Indomelothria (EF065456), Neoachmandra (EF065484–86), Urceodiscus (EF065464) and Zehneria
( EF065485, EF065489, EF065491–493, EF065497,
EF065499–500 and EF065502). This resulted in a sampling
of 114 of the 115 genera currently recognized in Cucurbitaceae
(Schaefer & Renner in press). The only genus of Cucurbitaceae
not yet sequenced is Khmeriosicyos W. J. de Wilde & B. Duyfjes,
which is only known from the Cambodian type collection.
Judging from morphology, it is expected to group with other
Asian Benincaseae. As outgroups, we used 15 species of the
Cucurbitales families Anisophyllaceae, Begoniaceae, Coriariaceae, Corynocarpaceae, Datiscaceae and Tetramelaceae,
based on Zhang et al. (2006). Of these, the Begoniaceae, Datiscaceae and Tetramelaceae, with the Cucurbitaceae, constitute
a morphologically and molecularly well-defined clade,
traditionally called the core Cucurbitales (Zhang et al. 2006).
(b) DNA extraction, amplification, sequencing and
alignments
Total genomic DNA was isolated from herbarium specimens
or, more rarely, silica-dried leaves with a commercial plant
DNA extraction kit (NucleoSpin, MACHEREY-NAGEL,
Düren, Germany), following the manufacturer’s manual. We
amplified the rbcL and matK genes, the trnL intron and the
trnL-F and rpl20–rps12 intergenic spacers. Polymerase chain
reactions (PCRs) were performed with the standard protocol
and primers described in Kocyan et al. (2007), and products
were purified with the Wizard SV PCR clean-up kit (Promega
GmbH, Mannheim, Germany). Cycle sequencing was
performed with BigDye Terminator cycle sequencing kits
on an ABI Prism 3100 Avant automated sequencer (Applied
Biosystems, Foster City, California, USA). Sequences were
edited with SEQUENCHER v. 4.6 (Gene Codes, Ann Arbor,
Proc. R. Soc. B (2009)
Michigan, USA) and aligned by eye, using MACCLADE v. 4.06
(Maddison & Maddison 2003).
The data matrices comprised 245 ingroup species plus 15
outgroup species. The lengths of the individual loci were
1356 aligned nucleotides for the rbcL gene, 1195 for the matK
gene, 667 for the tRNA-Leu (trnL) intron (after exclusion of
a poly A run and a highly variable microsatellite region), 803
for the tRNA-Leu–tRNA-Phe (trnL-F) intergenic spacer and
1010 for the rpl20–rps12 intergenic spacer. The combined
dataset comprised 5031 aligned nucleotides.
(c) Phylogenetic analysis
Maximum-likelihood (ML) tree searches and ML bootstrap
searches were performed using RAxML v. 7.0.3 (Stamatakis
et al. 2008; available at http://phylobench.vital-it.ch/raxml-bb/)
and GARLI v. 0.951 (Zwickl 2006; available at www.bio.utexas.
edu/faculty/antisense/garli/Garli.html). RAxML and GARLI
searches relied on the GTRCGCI model (six general timereversible substitution rates, assuming gamma rate heterogeneity and a proportion of invariable sites), with model
parameters estimated over the duration of specified runs.
Analyses in RAxML were run both with the combined
unpartitioned data and with a model that partitioned the rbcL
gene from the remaining non-coding regions. GARLI does
not allow data partitioning. The data matrix and trees have
been deposited in TREEBASE (www.treebase.org; study
number S2210).
(d) Molecular clock analysis
Estimation of divergence times relied either on a strict clock
or on a Bayesian relaxed clock with autocorrelated rates
( Thorne et al. 1998). Very short (‘zero-length’) branches are
known to cause problems for time estimation algorithms, and
we therefore reduced their number by using the best-scoring
ML tree for 147 taxa instead of the full 260-taxon tree. The
clock tree was rooted on Coriariaceae and Corynocarpaceae,
instead of Anisophyllaceae, because the latter are extremely
rich in autapomorphies, contributing towards rate heterogeneity near the base.
For the Bayesian approach, we used BASEML from the PAML
package ( Yang 2007) and MULTIDIVTIME ( Thorne et al. 1998;
Thorne & Kishino 2002) in LAGOPUS, an R package written
by Heibl & Cusimano (2008). LAGOPUS checks the input
data for consistency, automates the assignment of constraints
to nodes and connects the executables of the mentioned
software packages in a pipeline. Model parameters for the 147taxon matrix were estimated in BASEML, and branch lengths
and their variance then calculated in ESTBRANCHES, all under
the F84CG model (the only model implemented in MULTIDIVTIME). Priors for MULTIDIVTIME were as follows: based on
outgroup fossils (below), the prior on the mean age of the
root node was set to 84 Myr, with an equally large standard
deviation. The prior on the substitution rate at the root was
set to the value obtained by dividing the median distance
between the root and the tips in the ESTBRANCHES phylogram by
84 Myr. This yielded a rate of 0.0009 substitutions per site and
million years [S/(S!Myr)]. The prior for the Brownian motion
parameter, which controls the magnitude of autocorrelation
along the descending branches of the tree, was set to 1.11, with
a standard deviation of the same size. Markov chain Monte
Carlo (MCMC) samples were drawn for every 100th generation up to one million generations, with a burn-in of 100 000
cycles. Confidence in node ages was assessed using the
95 per cent credibility intervals calculated by MULTIDIVTIME.
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
Biogeography of the Cucurbitaceae
H. Schaefer et al.
Sicyos hillebrandii
Sicyos pachycarpus
Sicyos baderoa
Sicyos angulatus
Microsechium helleri
Sechiopsis tetraptera
Sicyosperma gracile
100
Sechium edule
Parasicyos dieterleae
Echinopepon wrightii
Echinopepon racemosus
89
Echinopepon paniculatus
Frantzia tacaco
Apatzingania arachoidea
Vaseyanthus insularis
Brandegea bigelovii
Pseudocyclanthera australis
98
100
Rytidostylis ciliata
99
Cyclanthera brachystachya
100
Hanburia mexicana
Elateriopsis oerstedii
78
Marah fabaceus
Marah macrocarpus
Echinocystis lobata
Linnaeosicyos amara
Trichosanthes montana
Trichosanthes kinabaluensis
Trichosanthes bracteata
99
Trichosanthes schlechteri
Gymnopetalum integrifolium
Trichosanthes pentaphylla
Gymnopetalum chinense
74
Hodgsonia heteroclita
Trichosanthes ovigera
61 97
Trichosanthes pendula
70
Trichosanthes villosa
Trichosanthes kirilowii
Trichosanthes postari
68 64
Trichosanthes reticulinervis
Trichosanthes cucumerina
Luffa operculata Mexico
Luffa operculata Senegal
86
Luffa echinata
Luffa aegyptiaca
Luffa graveolens
Luffa acutangula
Nothoalsomitra suberosa
100
Cucumis sativus
97
Cucumis hystrix
100
Cucumis maderaspatanus
Cucumis ritchiei
99
99
Cucumis javanicus
Cucumis melo
83
100
Cucumis zeyheri
Cucumis messorius
Cucumis sagittatus
69
Cucumis metuliferus
98
Cucumis asper
95
100
Cucumis oreosyce
Cucumis bryoniifolius
94
98
Cucumis humifructus
Cucumis hirsutus
98
Muellerargia jeffreyana
Muellerargia timorensis
76
Coccinia grandis
90
Coccinia sessilifolia
Coccinia rehmannii
Diplocyclos palmatus
Urceodiscus belensis
Papuasicyos papuana
Scopellaria marginata
Cucumeropsis mannii
79
Posadaea sphaerocarpa
72
Mslancium campestre
88
Melothria dulcis
Melothria pendula
80
Indomelothria chlorocarpa
Ruthalicia eglandulosa
Ruthalicia longipes
88
Trochomeria polymorpha
90
Trochomeria macrocarpa
99
Dactyliandra welwitschii
Ctenolepis cerasiformis
69
Zombitsia lucorum
88
Benincasa hispida
Praecitrullus fistulosus
88
Borneosicyos simplex
Solena heterophylla
77
Lemurosicyos variegatus
Cephalopentandra ecirrhosa
Raphidiocystis chrysocoma
64
Raphidiocystis phyllocalyx
Acanthosicyos horridus
98
Lagenaria breviflora
91
Lagenaria siceraria
Peponium caledonicum
Peponium vogelii
60
Citrullus colocynthis
Citrullus lanatus
100
Acanthosicyos naudinianus
Neoachmandra sphaerosperma
Neoachmandra japonica
Zehneria neocaledonica
Neoachmandra deltoidea
Neoachmandra samoensis
Zehneria baueriana
Zehneria grayana
84
Zehneria erythrobacca
Neoachmandra cunninghamii
Anangia macrosepala
Zehneria pallidinervia
64
Zehneria scabra
Zehneria bodinieri
82
Zehneria pisifera
94
Zehneria keayana
Zehneria minutiflora
Zehneria anomala
Cayaponia africana
Cayaponia
americana
100
Cayaponia podantha
96
Selysia prunifera
Abobra tenuifolia
Cionosicys macranthus
Schizocarpum filiforme
Schizocarpum palmeri
Schizocarpum reflexum
91
Tecunumania quetzalteca
Sicana odorifera
Calycophysum pedunculatum
Anacaona sphaerica
94
100 80
Penelopeia suburceolata
Cucurbita ficifolia
Cucurbita pepo
100
Cucurbita okeechobeensis subsp. martinezii
100
Cucurbita digitata subsp. palmata
Peponopsis adhaerens
98
Polyclathra cucumerina
100
Dieterlea aff. maxima
Dieterlea maxima
80
94
Dieterlea fusiformis
70
Ibervillea lindheimeri
Ibervillea millspaughii
100
Tumamoca macdougalii
Guraniopsis longipedicellata
Apodanthera mandonii
Ceratosanthes palmata
Halosicyos ragonesei
Kedrostis nana
Kedrostis africana
Melothrianthus smilacifolius
71
Gurania makoyana
Gurania spinulosa
Gurania tubulosa
Psiguria racemosa
Helmontia leptantha
90
Psiguria umbrosa
Wilbrandia verticillata
Doyerea emetocathartica
Apodanthera sagittifolia
Cucurbitella asperata
Corallocarpus
triangularis
99
96
100
Corallocarpus boehmii
Corallocarpus bainesii
96
Seyrigia humbertii
100
Trochomeriopsis diversifolia
83
Dendrosicyos socotranus
Eureiandra formosa
Bambekea racemosa
88
Edgaria darjeelingensis
96
Herpetospermum pedunculosum
99
Biswarea tonglensis
Schizopepon bryoniifolius
Austrobryonia pilbarensis
Austrobryonia argillicola
100
Austrobryonia micrantha
Austrobryonia centralis
100
83
Bryonia dioica
100
Bryonia alba
100
Bryonia verrucosa
Ecballium elaterium
Tricyclandra leandrii
Odosicyos bosseri
100
Ampelosicyos scandens
88
Ampelosicyos humblotii
100 100
Telfairia occidentalis
Telfairia pedata
88
Cogniauxia podolaena
Cogniauxia trilobata
92
Momordica cochinchinensis
Momordica cissoides
98
Momordica spinosa
85
Momordica cymbalaria
100
Momordica charantia
Momordica foetida
89
Momordica calantha
75
Thladiantha davidii
Thladiantha dubia
Thladiantha punctata
100
Thladiantha nudiflora
Thladiantha pustulata
97
Thladiantha hookeri
Sinobaijiania yunnanensis
100
Baijiania borneensis
Sinobaijiania smitinandii
72
Siraitia
siamensis
87
Siraitia grosvenorii
Microlagenaria africana
Indofevillea khasiana
100
Sicyeae
10 changes
100
100
Begoniaceae
Tetramelaceae
Datiscaceae
Coriariaceae
Corynocarpaceae
Anisophylleaceae
100
Bolbostemma paniculatum
Actinostemma tenerum
Xerosicyos perrieri
Xerosicyos danguyi
100
Zygosicyos tripartitus
88
Xerosicyos pubescens
100
Siolmatra brasiliensis
Zanonia indica
Gerrardanthus macrorhizus
76
100
Gerrardanthus paniculatus
Gerrardanthus grandiflorus
Sicydium diffusum
100
Sicydium tamnifolium
72
Chalema synanthera
92
Cyclantheropsis parviflora
90
Pseudosicydium
acariianthum
100
Pteropepon parodii
Fevillea trilobata
78
Fevillea pedatifolia
100
Fevillea anomalosperma
Fevillea pergamentacea
Neoalsomitra clavigera
100
99
Neoalsomitra integrifoliola
99
Neoalsomitra stephensiana
74
Neoalsomitra plena
100
Neoalsomitra trifoliolata
100
Neoalsomitra capricornica
99
Neoalsomitra angustipetala
100
99
Neoalsomitra podagrica
Neoalsomitra sarcophylla
99
Gynostemma pentaphyllum New Guinea
100
Gynostemma pentaphyllum
80
Hemsleya amabilis
100
Gomphogyne cirromitrata
Gomphogyne cissiformis
Bayabusua clarkei
68
Alsomitra macrocarpa
Begonia herbacea
100
Begonia oxyloba
Hillebrandia sandwichensis
100
Tetrameles nudiflora
Octomeles sumatrana
Datisca cannabina
Datisca glomerata
Coriaria myrtifolia
100
Coriaria nepalensis
Coriaria ruscifolia
Coriaria sarmentosa
Corynocarpus laevigatus
Anisophyllea fallax
Anisophyllea corneri
Combretoccarpus rotundatus
Cucurbiteae
Coniandreae
Actinostemmateae
Zanonieae
Fevilleeae
Gomphogyneae
outgroup
families
Schizopeponeae
Bryonieae
Telfairieae
Momodiceae
87
Thladiantheae
100
10 changes
Figure 1. (Caption overleaf.)
Proc. R. Soc. B (2009)
Siraitieae
Indofevilleeae
Cucurbitoideae
100
100
Cucurbitaceae
100
100
Benincaseae
845
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
846 H. Schaefer et al. Biogeography of the Cucurbitaceae
Figure 1. (Overleaf.). Best ML tree for Cucurbitaceae and relatives found with combined chloroplast gene, spacer and intron
sequences (5031 nucleotides) analysed under a GTRCGCI model with unlinked partitions for coding and non-coding regions.
Likelihood bootstrap values greater than 60% are given at the nodes. Rooting follows Zhang et al. (2006). The geographical
occurrence of genera is colour coded as follows: green, America (including Galapagos, Hawaii and the Caribbean); yellow,
mainland African; brown, Madagascar; red, Asia; blue, Australia/New Guinea/Polynesia; black, Europe.
To translate relative times into absolute times, the
Bayesian clock relied on the following simultaneous constraints. (i) The age of core Cucurbitales (the root node) was
constrained to maximally 84 Myr, based on the earliest fossils
of the sister group of Cucurbitales, the Fagales (Herendeen
et al. 1995). However, we also performed two runs in which
the root node was unconstrained or constrained to minimally
84 Myr. This yielded age estimates for the Cucurbitaceae
crown group that were older than the oldest angiosperm
fossils (132 Myr). In general, relaxed molecular clocks will
yield reliable ages only with at least one minimal and one
maximal constraint, the latter preferentially at or near the root
( Thorne et al. 1998). (ii) The age of the split between Datisca
and Octomeles/Tetrameles was set to minimally 68 Myr old or,
in an alternative run, to minimally 65.5 Myr old, based on the
fossil wood of Tetrameleoxylon prenudiflora from the Deccan
intertrappean beds at Mohgaonkalan in India (Lakhanpal &
Verma 1965; Lakhanpal 1970). These beds have been dated
to the Maastrichtian or Late Maastrichtian (Khajuria et al.
1994; Kar et al. 2003), and we therefore used either the
midpoint of the Maastrichtian (68 Myr) or the Maastrichtian/
Palaeocene border (65.5 Myr). (iii) The crown group of
Cucurbitaceae was set to minimally 65 Myr or, alternatively,
55.8 Myr, based on the seeds from the Palaeocene Felpham
flora (Collinson 1986; Collinson et al. 1993). These dates
span the upper and lower boundary of the Palaeocene.
(iv) The split between Linnaeosicyos, with tetracolpate–
reticulate pollen, and the remaining New World Sicyeae
(Schaefer et al. 2008a), which usually have polycolpate
pollen, was set to minimally 33.9 Myr or, alternatively,
23 Myr, based on Hexacolpites echinatus pollen from the
Oligocene of Cameroon, which is the oldest hexacolpate
Sicyeae-type pollen (Salard-Cheboldaeff 1978; Muller 1985).
Polycolpate pollen is not found in other Cucurbitaceae except
the African Neoachmandra peneyana ( Van der Ham &
Pruesapan 2006). The Oligocene epoch ranges from 33.9 to
23 Myr, and the stratum containing Hexacolpites has not been
precisely dated; therefore, in alternative analyses, we used the
upper or lower boundary as minimal constraints. (v) The split
between the Hispaniola endemics Anacaona and Penelopeia
was set to maximally 20 or 15 Myr, based on the age of
Dominican amber, which was produced by tropical trees and
provides a proxy for the presence of tropical forest on that
island; the amber age is estimated as 15–20 Myr (IturraldeVinent & MacPhee 1996).
To explore the sensitivity of the Bayesian clock to the
various priors, we performed alternative MCMC runs in
which we tested (a) the effect of a more clock-like Brownian
motion parameter of 0.4 instead of 1.11, (b) the effect of
using up to four data partitions, thereby allowing the genes
and spacer regions to have different rates, and (c) the effects of
varying the age constraints. For the latter exploration, we ran
an analysis in which all constraints were set to the lowest age
boundaries, another in which all constraints were set to the
highest age boundaries, and four analyses that used the
minimum age for one of the constraints and the maximum
age for the remaining constraints. All other parameters, such
Proc. R. Soc. B (2009)
as root rate and MCMC chain length, were constant between
these six runs. Finally, we performed a run (d) with mean ages
for four constraints, namely 66.8 Myr (constraint ii),
60.4 Myr (iii), 28.5 Myr (iv) and 17.5 Myr (v).
For the strict clock approach, rbcL branch lengths were
calculated under a GTRCGCIC clock model on the
preferred ML topology. The tree was imported into PAUP,
rooted on Coriaria (Zhang et al. 2006), and branch lengths
were then calculated under the ‘enforce clock’ option. The
distance between a calibration node and the present was
divided by the age of the calibration node to obtain a
substitution rate, and this rate was then used to calculate
the age of divergence events of interest. As calibration
nodes, we used either the age of the earliest Cucurbitaceae
seeds (constraint iii) or the oldest Sicyeae-type pollen
(constraint iv).
(e) Biogeographic analysis
For a dispersal-vicariance analysis, we used the 147-taxon
dataset also used for the clock runs and coded the distribution
ranges of all species in a binary matrix in MACCLADE. The
species were recorded as present in one of five regions: Asia;
Europe; Africa (including Madagascar); America (including
Caribbean, Galapagos, and Hawaii); and Australia (including
New Guinea and Polynesia). We then used the parsimonybased approach implemented in DIVA v. 1.1 (Ronquist 1996,
1997) to infer vicariance and dispersal events. The maximum
number of areas simultaneously occupied by hypothetical
ancestral lineages was experimentally constrained to 4, 3 or 2
because it is unlikely that an ancestral species would have
ranged over several continents.
3. RESULTS
Herbarium material yielded suitable DNA in more
than 95 per cent of the cases, even for 50–100-yearold collections.
The highest-scoring ML tree obtained for the 260-taxon
dataset (figure 1) shows Cucurbitaceae highly supported
as monophyletic and family relationships similar to those
found by Zhang et al. (2006) with a much larger amount
of sequence data. Within Cucurbitaceae, there are five
main clades (figure 1), namely: (i) a group of approximately 100 genera traditionally treated as subfamily
Cucurbitoideae (Kosteletzky 1833) and usually subdivided into several tribes (below); (ii) a clade of Asian
genera, including Alsomitra, Bayabusua and Neoalsomitra
that corresponds to the tribe Gomphogyneae of Bentham &
Hooker (1867); (iii) a clade of one African and five
Neotropical genera, including Fevillea and Sicydium, that
corresponds to the tribe Fevilleeae of Bentham & Hooker
(1867); (iv) a clade of a few genera from Madagascar,
continental Africa, Asia and South America corresponding
to the tribe Zanonieae of Blume (1826); and (v) a
clade consisting of the two Asian genera Actinostemma
and Bolbostemma.
Clades (ii–v) have been treated as subfamily Nhandiroboideae (an illegitimate name) or Zanonioideae
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
Biogeography of the Cucurbitaceae
(a taxonomic synonym of Fevilleoideae), but this subfamily is not supported as monophyletic by our data.
Clade (i), Cucurbitoideae, can be divided into geographically or morphologically more homogeneous groups that
correspond to the traditional tribes Herpetospermeae,
Bryonieae, Sicyeae, Coniandreae, Benincaseae and Cucurbiteae, plus a few clades of similar phylogenetic depth that
have not traditionally been ranked as tribes, such as the
Asian Thladiantha and Baijiania, the Asian/African Siraitia
and Microlagenaria, the African/Asian Momordica, the
African Telfairia, Cogniauxia, a group of Madagascan
genera and the Himalayan Indofevillea, which is sister to all
remaining Cucurbitoideae (figure 1). Well-known genera
found to be poly- or paraphyletic include Citrullus
(must include Acanthosicyos), Ampelosicyos (must include
Tricyclandra and Odosicyos), Gomphogyne (must include
Hemsleya), Xerosicyos (must include Zygosicyos), Apodanthera,
Psiguria and Trichosanthes.
MULTIDIVTIME dating runs with a Brownian motion
parameter of 0.4 instead of 1.11, yielding barely different
estimates for the ingroup nodes of interest. Runs that
allowed uncoupled rates for the two genes and the
spacers also yielded essentially identical estimates, and
final runs therefore modelled the data under a single
model. The estimates from the run in which all constraints
were set to their lowest boundaries differed significantly
from those obtained when all constraints were set to
their highest boundaries (Wilcoxon signed-rank test,
pZ0.0085; see fig. 1c,d in the electronic supplementary
material). Among the test runs in which one constraint
was set to the minimum age and the others to the
maximum age, only two yielded significantly different
results: the Tetrameleoxylon fossil set to the minimum age
and the Felpham flora seed set to the minimum age
(see fig. 1e in the electronic supplementary material).
However, all results were within the 95 per cent confidence intervals of the estimates obtained when the
constraints were set to mean ages (see fig. 1a,b in the
electronic supplementary material).
The substitution rates obtained under a strict clock
model calibrated with either the seed or the pollen fossil
(§2) were 0.00018 S/(S!Myr) (oldest Sycieae-type
pollen, constraint (iv)) or 0.00030 S/(S!Myr) (earliest
Cucurbitaceae seeds, constraint (iii)). An average rate of
0.00024 S/(S!Myr) yielded absolute times that for the
most part were older than those obtained with the relaxed
clock model (see table 2 in the electronic supplementary
material that lists the ages obtained with the relaxed
clock and with the strict clock model). The following
discussion focuses on the relaxed clock estimates
because they provide 95 per cent confidence intervals as
a measure of uncertainty.
The split between the two genera of Begoniaceae,
Begonia and Hillebrandia, is estimated as 29 (41–18) Myr
old, roughly the age of the Hawaiian archipelago
(ca 30 Myr), where Hillebrandia is endemic (Clement
et al. 2004). The deepest split in the Cucurbitaceae is ca
63 (69–61) Myr old, while crown group Cucurbitoideae
are 53 (60–48) Myr old, Gomphogyneae 56 (63–51)
Myr, Fevilleeae 46 (55–37) Myr and the Actinostemma/
Bolbostemma clade 52 (59–44) Myr. The Madagascan
Xerosicyos clade (split Zanonia–Xerosicyos) appears to be
49 (57–40) Myr old; the likewise Madagascan Ampelosicyos clade is 29 (39–19) Myr old. The Madagascar/
Proc. R. Soc. B (2009)
H. Schaefer et al.
847
Southeast Asia disjunction between the two species of
Muellerargia is only ca 12 (18–7) Myr old. The South
America/Asia disjunction in the Zanonia clade (Siolmatra–
Zanonia split), finally, is 24 (38–11) Myr old, and the
South America/Africa disjunction in the Fevillea clade
(Chalema–Cyclantheropsis split) 41 (51–31) Myr. Other
estimates of specific interest (Introduction; see table 2 in the
electronic supplementary material) concern Cayaponia
(stem age 10 (17–5) Myr), Luffa (stem age 35 (41–31)
Myr), Sicyos (stem age 8 (13–4) Myr) and the Lagenaria
crown group (8 (14–3) Myr; figure 2).
The DIVA analysis yielded Asia as the most likely
region of origin of the Cucurbitaceae (figure 3a). From
there, at least five lineages reached Africa, and 12 lineages
reached Australia. No fewer than seven lineages independently reached the American continent, and no fewer
than 13 lineages dispersed from Africa to Madagascar.
Some disjunctions are best explained as secondary
dispersals from Africa back to Asia (Coccinia and
Momordica). A few lineages reached Africa via LDD
from America (Cayaponia africana, Cucumeropsis, perhaps
Cyclantheropsis and Kedrostis).
4. DISCUSSION
(a) The geographical origins of the world’s
regional Cucurbitaceae floras
Our results suggest that Cucurbitaceae initially diversified
in Asia (specifically, the region north of the Tethys)
sometime in the Late Cretaceous. This fits with the
observation that India contains more deeply divergent
lineages of Cucurbitaceae than any other similar-sized
geographical area (Chakravarty 1946, 1959; this study). Of
the family’s Late Cretaceous radiations, two clades (the
Gomphogyneae and the Actinostemmateae) are now
almost restricted to subtropical Asia (figure 1). A third
clade, Fevilleeae, is mainly Neotropical except for a small
African ‘extension’, Cyclantheropsis. The ancestors of
Fevilleeae were probably more widely distributed in the
Laurasian tropics and reached the American continent by
dispersing across a still-narrow Atlantic (figure 3b; seeds of
Fevilleeae are wind- and water-dispersed). Cyclantheropsis
must result from a back dispersal from South America to
Africa in the Middle Eocene. The ancestors of the fourth
ancient clade, Zanonieae, apparently reached the African
continent early and from there dispersed to Madagascar (the
Early Eocene Xerosicyos lineage; figure 3b). Later, in the
Oligocene, at least two LDD events brought the Siolmatra
lineage to America and the Zanonia lineage back to tropical
Asia (figure 3c). The fifth and last ancient clade, the
Cucurbitoideae, diversified partly in Asia (e.g. Thladiantha,
Siraitia, Trichosanthes), and partly in Africa (e.g. Momordica,
Cucumis, Coccinia, Kedrostis). The cucumber (Cucumis
sativus) and its closest relatives (not all included in the
present study) evolved from a common ancestor ca 3 (6–1)
Myr ago. The wax gourd Benincasa and its sister group
Praecitrullus, an important vegetable in parts of India,
apparently split only 5 (10–1) Myr ago. Further dispersals
from Africa back to Asia are present within Momordica,
Coccinia, Kedrostis and Corallocarpus (figure 3c).
The native European cucurbit flora consists only of
Bryonia, with 10 species ( Volz & Renner 2008), and its
monotypic sister Ecballium, which probably represent a
lineage that spread along the Tethys border from Asia to
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
848 H. Schaefer et al. Biogeography of the Cucurbitaceae
73 ± 3
Tetrameles nudiflora
Octomeles sumatrana
Begonia herbacea
Hillebrandia sandwichensis
Datisca cannabina
Hemsleya amabilis
23 ± 5
Gomphogyne cissiformis
44 ± 5
Gynostemma pentaphyllum
39 ± 5
Neoalsomitra podagrica
51 ± 4
27 ± 5
Neoalsomitra capricornica
56 ± 3
Bayabusua clarkei
Alsomitra macrocarpa
Chalema synanthera
18 ± 5
Sicydium tamnifolium
39 ± 5
Pteropepon parodii
59 ± 3
34
±
6
41 ± 5
Pseudosicydium acariianthum
Cyclantheropsis parviflora
46 ± 5
Fevillea pergamentacea
23 ± 7
Fevillea trilobata
Xerosicyos pubescens
58 ± 3
17 ± 6
Xerosicyos danguyi
49 ± 4
Zanonia indica
24 ± 7
Siolmatra brasiliensis
54 ± 4
Gerrardanthus grandiflorus
28 ± 6
Gerrardanthus macrorhizus
52± 4
Bolbostemma paniculatum
32 ± 6
Actinostemma tenerum
Indofevillea khasiana
Gymnopetalum chinense
28 ± 3
Trichosanthes cucumerina
31 ± 3
Trichosanthes schlechteri
28 ± 3
Hodgsonia heteroclita
Hanburia mexicana
13 ± 3
Pseudocyclanthera australis
9±3
Cyclanthera brachystachya
15 ± 3
Echinopepon wrightii
12 ± 3
Vaseyanthus insularis
4±2
33 ± 2
Brandegea bigelovii
Sicyos baderoa
16 ± 3
2±2
Sicyos hillebrandii
8±2
Sechiopsis tetraptera
9±2
Sicyosperma gracile
10
±
3
17 ± 3
11
±
3
Sechium edule
35 ± 3
Parasicyos dieterleae
31 ± 2
Echinocystis lobata
15 ± 3
Marah macrocarpus
Linnaeosicyos amara
37 ± 3
Luffa operculata – Senegal
3±3
Luffa operculata – Mexico
11 ± 5
Luffa acutangula
Nothoalsomitra suberosa
Peponopsis adhaerens
16 ± 4
Cucurbita ficifolia
7
±
3
19 ± 4
Cucurbita pepo
Polyclathra cucumerina
Cionosicys macranthus
16 ± 3
Selysia prunifera
23 ± 4
10 ± 3
Cayaponia africana
17 ± 3
3±2
Cayaponia americana
19 ± 3
Schizocarpum reflexum
Tecunumania quetzalteca
21 ± 3
Penelopeia suburceolata
13 ± 3
Anacaona sphaerica
17 ± 3
Calycophysum pedunculatum
mannii
30 ± 4
1 ± 1 Cucumeropsis
Posadaea sphaerocarpa
8±3
Melancium campestre
10 ± 3
Melothria dulcis
7±3
15 ± 4
Melothria pendula
Indomelothria chlorocarpa GB
Ruthalicia eglandulosa
Muellerargia timorensis
12 ± 3
24 ± 4
Muellerargia jeffreyana
Cucumis sativus
16 ± 3
3±2
Cucumis hystrix
9±3
39 ± 3
14 ± 3
Cucumis melo
19
±
4
21 ± 4
Cucumis hirsutus
Coccinia grandis
15 ± 3
Diplocyclos palmatus
cunninghamii
1 ± 1 Neoachmandra
Neoachmandra japonica
9±3
11
±
3
Zehneria scabra
53 ± 3
20 ± 4
14 ± 4
Zehneria bodinieri
Zehneria anomala
Acanthosicyos horridus
Solena heterophylla
10 ± 3
Lemurosicyos variegatus
12
±
3
16
±
3
34 ± 3
Cephalopentandra ecirrhosa
Ctenolepis cerasiformis
19 ± 4
7±3
14 ± 3
Zombitsia lucorum
10 ± 3
Dactyliandra welwitschii
17 ± 3
7±3
13 ± 3
Trochomeria macrocarpa
Praecitrullus fistulosus
5±2
Benincasa hispida
Papuasicyos papuana
18 ± 3
Peponium vogelii
4±2
Peponium caledonicum
12 ± 3
Lagenaria siceraria
8±3
Lagenaria breviflora
16 ± 3
Citrullus lanatus
2±2
Citrullus colocynthis
41 ± 3
12 ± 4
Raphidiocystis chrysocoma
Bambekea racemosa
Corallocarpus boehmii
11 ± 3
38 ± 3
Cucurbitella asperata
13 ± 3
Halosicyos ragonesei
15 ± 4
Kedrostis africana
Wilbrandia verticillata
11
±
3
Gurania makoyana
17 ± 4
6±3
25 ± 4
Psiguria racemosa
13 ± 3
Doyerea emetocathartica
12 ± 3
Ceratosanthes palmata
15 ± 4
Melothrianthus smilacifolius
19 ± 4
Guraniopsis longipedicellata
14 ± 4
4±2
Apodanthera mandonii
13 ± 3
Dieterlea maxima
9
±
3
44 ± 3
22 ± 4
Ibervillea millspaughii
Trochomeriopsis diversifolia
15 ± 4
Seyrigia humbertii
Dendrosicyos socotranus
Biswarea tonglensis
14 ± 5
Herpetospermum pedunculosum
29 ± 4
Schizopepon bryoniifolius
Ecballium elaterium
23 ± 4
Bryonia dioica
7±3
Bryonia verrucosa
32 ± 4
46 ± 3
Austrobryonia micrantha
13 ± 4
Austrobryonia centralis
Ampelosicyos scandens
10 ± 5
Tricyclandra leandrii
5±3
Odosicyos bosseri
29 ± 5
Telfairia pedata
15 ± 5
49 ± 3
35 ± 5
Telfairia occidentalis
Cogniauxia podolaena
Momordica cochinchinensis
22 ± 6
Momordica cissoides
31 ± 5
Momordica charantia
51 ± 3
Siraitia grosvenorii
33 ± 8
Microlagenaria africana
Thladiantha hookeri
14 ± 5
Thladiantha dubia
45 ± 4
Baijiania borneensis
27 ± 6
Sinobaijiania yunnanensis
32 ± 6
Sinobaijiania smitinandii
26 ± 7
63 ± 4
29 ± 6
ii
i
81 ± 2
iv
63 ± 3
iii
v
80
60
40
20
0 age in myr
Figure 2. Chronogram obtained for Cucurbitaceae under a Bayesian autocorrelated rates relaxed clock model applied to the
combined data (5031 nucleotides) and calibrated with three minimal (yellow) and two maximal (orange) constraints as in run
(d ) of §2(d). Age estimates with their 95% confidence ranges shown in purple. Rooting follows Zhang et al. (2006). Green,
America (including Galapagos, Hawaii and the Caribbean); yellow, mainland Africa; brown, Madagascar; red, Asia; blue,
Australia/New Guinea/Polynesia; black, Europe.
Proc. R. Soc. B (2009)
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
Biogeography of the Cucurbitaceae
(a)
1
Cucurbitaceae
(b)
Fevilleeae
Zanonieae
2
3
s
Xerosicyo
(c)
Momordica 5
Siolmatra
Zanonia
4
6
Ampelosicyos
(d )
Cayaponia
9 Sicyos
7
88
Luffa
12x
11x
Figure 3. The biogeographic history of Cucurbitaceae as
inferred from the statistical approach described in the text
(coastlines drawn after Smith et al. 1994). (a) Late Cretaceous,
ca 70 Myr ago; (1) origin of the Cucurbitaceae in Asia.
(b) Palaeocene/Eocene, 60–40 Myr ago; (2) the ancestor of the
Zanonieae reaches Africa, the ancestor of the Xerosicyos
lineage reaches Madagascar, and (3) the ancestor of Fevilleae
reaches South America. (c) Oligocene, ca 30 Myr ago; (4) the
ancestor of the Siolmatra lineage disperses over the Atlantic
into South America, (5) the ancestors of Momordica cochinchinensis and Zanonia indica independently disperse from Africa
to Southeast Asia, and (6) the ancestor of the Ampelosicyos
lineage reaches Madagascar. (d ) Middle Miocene, ca 10 Myr
ago; (7) the ancestors of Neotropical Luffa disperse from
Africa to the Americas, (8) the ancestor of Cayaponia africana
disperses from South America to West Africa, and (9)
ancestors of several Sicyos species groups spread from South
America to Hawaii, Galapagos, New Zealand and Australia.
the Mediterranean 32 (41–24) Myr ago. The remaining
cucurbit species that occur in Europe are the result of
recent introductions (Echinocystis lobata, Sicyos angulatus
Proc. R. Soc. B (2009)
H. Schaefer et al.
849
and Thladiantha dubia) or casual escapes from cultivation
(Citrullus lanatus, Cucumis melo, C. sativus and Cucurbita
pepo). The closest extant relative of Bryonia and Ecballium
is the Australian genus Austrobryonia (four species), which
may have reached Australia from Asia ca 13 (21–6) Myr
ago (Schaefer et al. 2008b).
African Cucurbitoideae (25 genera) are the result of
five dispersals from Asia to Africa and two from America
to Africa (in the genera Cucumeropsis and Cayaponia).
The watermelon (Citrullus lanatus) and its sister species
(C. colocynthis) apparently evolved from a common
ancestor as recently as 2 (6–0.1) Myr ago. The lineage
leading to the cucumber tree, Dendrosicyos socotranus,
endemic on Socotra, some 350 km off the Arabian
peninsula, is estimated as 22 (30–14) Myr old, while the
Socotra archipelago is only some 10 Myr old (Ghebreab
1998). Dendrosicyos thus seems to be an island relict of a
progenitor lineage that went extinct on the mainland. This
example of a species that is twice as old as the island on
which it occurs cautions against using geological calibrations in molecular clock dating. Another supposed
example of a species being older than the island on
which it lives, that of the Hawaiian Hillebrandia sandwicensis (Clement et al. 2004), is not supported by our data
(see table 2 in the electronic supplementary material).
Madagascar has 16 native Cucurbitaceae genera with
50 species in total. From our data it appears that
Cucurbitaceae reached Madagascar at least 13 times,
apparently always from the African mainland, and that
these 13 ancestors then underwent local radiations,
giving rise to today’s 50 species. Using Madagascar as a
stepping stone, one of these clades, Peponium, later
reached the Seychelles (the endemic species there has
not yet been sequenced).
South America has approximately 350 species of
Cucurbitoideae in 47 genera that all descend from five
LDD events, mostly from Africa to South America. These
involved the ancestors of Cucurbiteae, Sicyinae, a clade of
Coniandreae, the Melothria clade and a subclade of Luffa.
Based on the tree topology (figure 1), Luffa originated in
the Old World or Australia, and one species then reached
the New World by LDD from Africa across the Atlantic as
suggested by Heiser & Schilling (1988). The fruit is dry
with fibrous tissue and probably well adapted to floating
(Ridley 1930). The Neotropical Melothria clade (figure 1)
appears to have crossed the Pacific because the sister
group of Melothria, Indomelothria, is endemic in Southeast
Asia. Today’s pumpkin and squash species (Cucurbita spp.
in the Cucurbiteae) apparently originated in Central or
South America, and the genus Cucurbita split from its
sister clade, Peponopsis, only some 16 (23–9) Myr ago.
North American Cucurbitaceae, finally, all descend from
seven expansions of Central and South American lineages
that occurred at widely different times (figure 2; see table 2
in the electronic supplementary material).
The indigenous Australian Cucurbitaceae flora consists
of 30 species in 12 genera of which two are endemic:
Nothoalsomitra, a liana species of Queensland’s rainforests;
and Austrobryonia, four species of trailers or creepers in
the dry regions of (mostly) Central Australia. This low
Australian species diversity is in marked contrast with the
minimally 12 independent dispersal events into Australia
(figures 1 and 3d ). The largest Australian ‘radiation’
comprises only four species (Schaefer et al. 2008b), even
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
850 H. Schaefer et al. Biogeography of the Cucurbitaceae
though the ecological conditions in Australian rainforests
and bushland are similar to those in Southeast Asian
forests and African bushland, where cucurbit diversity
is much higher.
Overall, Cucurbitaceae underwent at least 43 successful LDD events over the past 60 Myr, which would
translate into an average of seven LDDs every 10 Myr.
These events seem to have occurred throughout the
evolutionary past of the family, rather than being clustered
at a particular geological time (see table 2 in the electronic
supplementary material). The most striking case of a rapid
radiation following LDD is Sicyos, which reached Hawaii
only ca 2 (6–0.1) Myr ago (figure 3d ) and now has 15
Hawaiian species. Sicyos species often grow in seabird
colonies (www.botany.hawaii.edu/gradstud/eijzenga/
OIRC/lehua.htm# Vegetation), and the hooked barbs on
the fruits may help external dispersal on bird feathers.
Another extreme case is Muellerargia, the closest relative of
Cucumis (Renner & Schaefer 2008), which comprises just
two species, one endemic in Madagascar, and the other in
northeast Australia and Timor.
(b) Using museomics to achieve less biased taxon
sampling
Results of this study may have implications for the
conservation of wild Cucurbitaceae gene pools for future
crop improvement; this is almost certainly true for Luffa
and Cucumis. They also show that the gourd family, which
appears to be of Asian origin, has undergone numerous
natural LDD events, prior to anthropogenic LDD, and
that these events have played a large role in the build-up of
local cucurbit floras. More generally, this study illustrates
the great potential of herbarium material for generating
broad phylogenetic frameworks at relatively low costs and
in a time-efficient manner (no lengthy permit procedures).
With rapidly improving molecular techniques for DNA
isolation from tiny fragments of ancient collections, it is
now possible to study the phylogeny of worldwide lineages
without time-consuming and difficult field trips. And
broad geographical sampling in turn is the precondition
for assessing the full extent of LDD in the angiosperms at
different geological times and across different ocean
basins. Given that habitat destruction in the centres of
cucurbit diversity (Madagascar, Southeast Asia, West
Africa and Central America) is extremely high, the use
of herbarium material also may soon be the only option for
future research on the phylogenetics of this plant group.
We thank B. Duyfjes and W. J. de Wilde (L) for silica samples
and advice; M. Nee (NY ), M. Pignal (P), J. Wieringa (WAG)
and the curators of the herbaria AAU, ASU, B, BKF, BR,
BSC, CHR, CMU, E, EA, F, FTG, G, GIFU, K, KUN, L,
LE, LISC, M, MEXU, MO, NE, RSA, ULM, US and Z for
permission to sample DNA from specimens in their care;
E. Vosyka for help in the laboratory; and the German science
foundation for financial support (DFG RE603/3-1).
REFERENCES
Alsos, I. G., Eidesen, P. B., Ehrich, D., Skrede, I.,
Westergaard, K., Jacobsen, G. H., Landvik, J. Y.,
Taberlet, P. & Brochmann, C. 2007 Frequent longdistance plant colonization in the changing arctic. Science
316, 1606–1608. (doi:10.1126/science.1139178)
Bentham, G. & Hooker, J. D. 1867. Genera plantarum, vol.1.
London, UK: Reeve & Co.
Proc. R. Soc. B (2009)
Blume, C. L. von. 1826 Bijdragen tot de Flora van Nederlandsch
Indie 15. Batavia, IL: Ter Lands Drukkerij.
Chakravarty, H. L. 1946 Studies on Indian Cucurbitaceae
with special remarks on distribution and uses of economic
species. Indian J. Agr. Sci. 16, 1–89.
Chakravarty, H. L. 1959 Monograph on Indian Cucurbitaceae. Rec. Bot. Surv. India 17, 1–234.
Chen, J. F., Zhuang, F. Y., Liu, X. A. & Qian, C. T. 2004
Reciprocal differences of morphological and DNA characters in interspecific hybridization in Cucumis. Can.
J. Bot. 82, 16–21. (doi:10.1139/b03-107)
Clement, W. L., Tebbitt, M. C., Forrest, L. L., Blair, J. E.,
Brouillet, L., Eriksson, T. & Swensen, S. M. 2004
Phylogenetic position and biogeography of Hillebrandia
sandwicensis (Begoniaceae): a rare Hawaiian relict. Am.
J. Bot. 91, 905–917. (doi:10.3732/ajb.91.6.905)
Collinson, M. E. 1986 The Felpham flora—a preliminary
report. Tertiary Res. 8, 29–32.
Collinson, M. E., Boulter, M. C. & Holmes, P. R. 1993
Magnoliophyta (‘Angiospermae’). In The fossil record (ed.
M. J. Benton), pp. 809–841. London, UK: Chapman
and Hall.
De Wilde, W. J. J. O. & Duyfjes, B. E. E. 2007 The wild
species of Cucumis L. (Cucurbitaceae) in South-East Asia.
Adansonia, sér. 3 29, 239–248.
Dillehay, T. D., Rossen, J., Andres, T. C. & Williams, D. E.
2007 Preceramic adoption of peanut, squash, and cotton
in Northern Peru. Science 316, 1890–1893. (doi:10.1126/
science.1141395)
Galley, C., Bytebier, B., Bellstedt, D. U. & Linder, H. P. 2007
The Cape element in the Afrotemperate flora: from Cape
to Cairo? Proc. R. Soc. B 274, 535–543. (doi:10.1098/rspb.
2006.0046)
Gentry, A. H. 1991 The distribution and evolution of
climbing plants. In The biology of vines (eds F. E. Putz &
H. A. Mooney), pp. 3–42. Cambridge, UK: Cambridge
University Press.
Ghebreab, W. 1998 Tectonics of the Red Sea region
reassessed. Earth Sci. Rev. 45, 1–44. (doi:10.1016/
S0012-8252(98)00036-1)
Ghebretinsae, A. & Barber, J. C. 2006. Phylogenetic
relationships among species of Cucumis and Cucumella
(Cucurbitaceae): evidence from ITS, rpl16 and trnS-G
sequences. See www.2006.botanyconference.org/engine/
search/index.php?funcZdetail&aidZ498.
Heibl, C. & Cusimano, N. 2008. LAGOPUS—divergence
time estimation in R. See www.christophheibl.de/mdt/
mdtinr.html.
Heiser, C. B. & Schilling, E. E. 1988 Phylogeny and
distribution of Luffa (Cucurbitaceae). Biotropica 20,
185–191. (doi:10.2307/2388233)
Herendeen, P. S., Crane, P. R. & Drinnan, A. N. 1995
Fagaceous flowers, fruits, and capsules from the Campanian (Late Cretaceous) of central Georgia, USA. Int.
J. Plant Sci. 156, 93–116. (doi:10.1086/297231)
Iturralde-Vinent, M. E. & MacPhee, R. D. E. 1996 Age and
paleogeography of Dominican amber. Science 273,
1850–1852. (doi:10.1126/science.273.5283.1850)
Kar, R. K., Ambwani, K., Sahni, A. & Sharma, P. 2003
Unisexual flowers from the Deccan Intertrappean Bed of
Madhya Pradesh, India. Palaeobotanist 52, 73–79.
Khajuria, C. K., Prasad, G. V. R. & Manhas, B. K. 1994
Palaeontological constraints on the age of Deccan Traps,
Peninsular India. Newsletter on Stratigraphy 31, 21–32.
Kocyan, A., Zhang, L., Schaefer, H. & Renner, S. S. 2007
A multi-locus chloroplast phylogeny for the Cucurbitaceae
and its implications for character evolution and classification. Mol. Phyl. Evol. 44, 553–577. (doi:10.1016/
j.ympev.2006.12.022)
Downloaded from rspb.royalsocietypublishing.org on 8 March 2009
Biogeography of the Cucurbitaceae
Kosteletzky, V. F. 1833. Allgemeine medizinisch-pharmazeutische Flora. vol. 6. Prague, Czech Republic: Borrosch &
Andre.
Lakhanpal, R. N. 1970 Tertiary floras of India and their
bearing on the historical geology of the region. Taxon 19,
675–694. (doi:10.2307/1219280)
Lakhanpal, R. N. & Verma, J. K. 1965 Fossil wood of
Tetrameles from the Deccan Intertrappean beds of
Mohgaonkalan, Madhya Pradesh. Palaeobotanist 14,
209–213.
Lavin, M. & Beyra Matos, A. 2008 The impact of ecology and
biogeography on legume diversity, endemism, and phylogeny in the Caribbean region: a new direction in historical
biogeography. Bot. Rev. 74, 178–196. (doi:10.1007/
s12229-008-9006-8)
Maddison, W. P. & Maddison, D. K. 2003 MACCLADE:
analysis of phylogeny and character evolution. Sunderland,
MA: Sinauer Associates.
Muller, J. 1985 Significance of fossil pollen for angiosperm
history. Ann. Miss. Bot. Garden 71, 419–443. (doi:10.
2307/2399033)
Nee, M., Schaefer, H. & Renner, S. S. Submitted. A new,
anomalous species of Fevillea (Cucurbitaceae) from
Bolivia. Syst. Bot.
Renner, S. S. & Schaefer, H. 2008. Phylogenetics of Cucumis
(Cucurbitaceae) as understood in 2008. In Proc. IXth
EUCARPIA meeting on genetics and breeding of Cucurbitaceae (ed. M. Pitrat), pp. 53–58. See https://w3.avignon.
inra.fr/dspace/handle/2174/236.
Ridley, H. N. 1930 The dispersal of plants throughout the world.
Ashford, UK: L. Reeve.
Ronquist, F. 1996. DIVA v. 1.1. Computer program and
manual available by FTP from Uppsala University. See
ftp.uu.se or ftp.systbot.uu.se.
Ronquist, F. 1997 Dispersal-vicariance analysis: a new
approach to the quantification of historical biogeography.
Syst. Biol. 46, 195–203. (doi:10.2307/2413643)
Salard-Cheboldaeff, M. 1978 Sur la palynoflore Maestrichtienne et Tertiaire du bassin sédimentaire littoral du
Cameroun. Pollen et Spores 20, 215–260.
Sanjur, O. I., Piperno, D. R., Andres, T. C. & Wessel-Beaver,
L. 2002 Phylogenetic relationships among domesticated
and wild species of Cucurbita (Cucurbitaceae) inferred
from a mitochondrial gene: implications for crop plant
evolution and areas of origin. Proc. Natl Acad. Sci. USA 99,
535–540. (doi:10.1073/pnas.012577299)
Schaefer, H. & Renner, S. S. In press. Cucurbitaceae.
In Families and genera of vascular plants. (ed. K. Kubitzki).
Berlin, Germany: Springer.
Schaefer, H., Kocyan, A. & Renner, S. S. 2008a Linnaeosicyos
(Cucurbitaceae): a new genus for Trichosanthes amara, the
Caribbean sister species of all Sicyeae. Syst. Bot. 33,
349–355. (doi:10.1600/036364408784571707)
Proc. R. Soc. B (2009)
H. Schaefer et al.
851
Schaefer, H., Telford, I. R. H. & Renner, S. S. 2008b
Austrobryonia (Cucurbitaceae), a new Australian endemic
genus, is the closest living relative to the Eurasian and
Mediterranean Bryonia and Ecballium. Syst. Bot. 33,
125–132. (doi:10.1600/036364408783887357)
Smith, B. D. 1997 The initial domestication of Cucurbita pepo
in the Americas 10,000 years ago. Science 276, 932–934.
(doi:10.1126/science.276.5314.932)
Smith, A. G., Smith, D. G. & Funell, M. 1994 Atlas of
mesozoic and cenozoic coastlines. Cambridge, UK:
Cambridge University Press.
Stamatakis, A., Hoover, P. & Rougemont, J. 2008
A rapid bootstrap algorithm for the RAxML web-servers.
Syst. Biol. 75, 758–771. (doi:10.1080/106351508024
29642)
Thorne, J. L. & Kishino, H. 2002 Divergence time and
evolutionary rate estimation with multilocus data. Syst.
Biol. 51, 689–702. (doi:10.1080/10635150290102456)
Thorne, J. L., Kishino, H. & Painter, I. S. 1998 Estimating
the rate of evolution of the rate of molecular evolution.
Mol. Biol. Evol. 15, 1647–1657.
Van der Ham, R. & Pruesapan, K. 2006 Pollen morphology of
Zehneria s.l. (Cucurbitaceae). Grana 45, 241–248.
(doi:10.1080/00173130600874834)
Volz, S. M. & Renner, S. S. 2008 Hybridization, polyploidy,
and evolutionary transitions between monoecy and dioecy
in Bryonia (Cucurbitaceae). Am. J. Bot. 95, 1297–1306.
(doi:10.3732/ajb.0800187)
Wagner, W. L., Herbst, D. R. & Sohmer, S. H. 1990 Manual
of the flowering plants of Hawai’i. Honolulu, HI: University
of Hawaii Press.
Whitaker, T. W. & Carter, G. F. 1954 Oceanic drift of
gourds—experimental observations. Am. J. Bot. 41,
697–700. (doi:10.2307/2438952)
Winkworth, R. C., Wagstaff, S. J., Glenny, D. & Lockhart,
P. J. 2005 Evolution of the New Zealand mountain flora:
origins, diversification and dispersal. Org. Div. Evol. 5,
237–247. (doi:10.1016/j.ode.2004.12.001)
Yang, Z. 2007 PAML 4: a program package for phylogenetic
analysis by maximum likelihood. Mol. Biol. Evol. 24,
1586–1591. (doi:10.1093/molbev/msm088)
Zhang, L., Simmons, M. P., Kocyan, A. & Renner, S. S.
2006 Phylogeny of the Cucurbitales based on DNA
sequences of nine loci from three genomes: implications
for morphological and sexual system evolution.
Mol. Phyl. Evol. 39, 305–322. (doi:10.1016/j.ympev.
2005.10.002)
Zwickl, D. J. 2006. Genetic algorithm approaches for the
phylogenetic analysis of large biological sequence data sets
under the maximum likelihood criterion. PhD dissertation, The University of Texas, Austin, TX.