Analyses of amplified fragment length polymorphisms (AFLP) indicate rapid radiation of Diospyros species (Ebenaceae) endemic to New Caledonia
© Turner et al.; licensee BioMed Central Ltd. 2013
Received: 14 August 2013
Accepted: 9 December 2013
Published: 12 December 2013
Radiation in some plant groups has occurred on islands and due to the characteristic rapid pace of phenotypic evolution, standard molecular markers often provide insufficient variation for phylogenetic reconstruction. To resolve relationships within a clade of 21 closely related New Caledonian Diospyros species and evaluate species boundaries we analysed genome-wide DNA variation via amplified fragment length polymorphisms (AFLP).
A neighbour-joining (NJ) dendrogram based on Dice distances shows all species except D. minimifolia, D. parviflora and D. vieillardii to form unique clusters of genetically similar accessions. However, there was little variation between these species clusters, resulting in unresolved species relationships and a star-like general NJ topology. Correspondingly, analyses of molecular variance showed more variation within species than between them. A Bayesian analysis with BEAST produced a similar result. Another Bayesian method, this time a clustering method, Structure, demonstrated the presence of two groups, highly congruent with those observed in a principal coordinate analysis (PCO). Molecular divergence between the two groups is low and does not correspond to any hypothesised taxonomic, ecological or geographical patterns.
We hypothesise that such a pattern could have been produced by rapid and complex evolution involving a widespread progenitor for which an initial split into two groups was followed by subsequent fragmentation into many diverging populations, which was followed by range expansion of then divergent entities. Overall, this process resulted in an opportunistic pattern of phenotypic diversification. The time since divergence was probably insufficient for some species to become genetically well-differentiated, resulting in progenitor/derivative relationships being exhibited in a few cases. In other cases, our analyses may have revealed evidence for the existence of cryptic species, for which more study of morphology and ecology are now required.
KeywordsCryptic species Island flora Morphological diversification Progenitor/derivative relationships Species radiation Woody plants
Island floras are often characterized by high levels of endemism and groups of closely related but morphological and ecological divergent species that are mostly the result of single colonisation events followed by radiation e.g. [1, 2]. New Caledonia was cited as one of the 34 biodiversity hotspots recognized by Conservation International [3, 4]. Nearly 75% of the native flora is endemic , which is the fourth highest for an island . While the continental part of New Caledonia (mainly Grande Terre) was entirely submerged during the Eocene (until 37 mya), a thick layer of heavy-metal-rich oceanic mantle accumulated . Today, around one-third of the main island, Grande Terre, is still overlaid with ultramafic substrates. Generally, Grande Terre is a substrate mosaic , which is cited as one reason for the high level of endemism found there e.g. . The climate in New Caledonia ranges from tropical to subtropical, and the main island is split by a mountain range into a humid eastern and a dry western part with prevailing winds and rain coming from the south-east. Taking climatic and geological factors together, Grande Terre has a wide range of environmental diversity. The main vegetation types in New Caledonia are evergreen humid forests, maquis, dry forests, littoral vegetation, and savannah .
Occurrence of Diospyros species in different habitats in New Caledonia
Humid mountain forest
D. parviflora, D. trisulca
D. flavocarpa, D. labillardierei
Humid low land forest
D. glans, D. pancheri, D. parviflora, D. umbrosa
D. minimifolia, D. pustuala, D. tridentata
D. cherrieri, D. erudita, D. minimifolia, D. perplexa, D. pustulata
D. calciphila, D. inexplorata
Diospyros colonised New Caledonia via long-distance dispersal at least four times . In previous studies based on low-copy nuclear and/or multiple plastid markers [12, 13], it was possible to resolve phylogenetic relationships for the majority of Diospyros species, except for one group of endemics from New Caledonia. Of the 31 New Caledonian Diospyros species, 24 belong to this clade of closely related endemics. In previous analyses, this strongly supported group is related to species found on islands throughout the Indian and Pacific Oceans as far east as Hawai’i [12, 13]. However, due to extremely low levels of sequence divergence, it was not possible to tease apart relationships between these species (they formed a hard polytomy in most individual trees, and there was little informative variation that permitted clustering of pairs or groups of species). Most of these closely related species are morphologically and ecologically clearly differentiated (for examples see ), and several species are narrow endemics restricted to small areas.
Amplified fragment length polymorphism (AFLP; ) is a fingerprinting technique that has proven to be useful for revealing phylogenetic relationship among closely related taxa (e.g. Hypochaeris, ; Lactuca, ; Phylica, ; Trollius, ; Ranunculus alpestris, ; Puya, [20, 21]; Araucaria, ). In contrast to standard phylogenetic markers, AFLP variation is spread across the whole genome, spanning both coding and non-coding DNA regions and may therefore be more representative of overall genetic patterns present as well as being highly informative for phylogenetic analyses at the low phylogenetic level [23, 24]. Compared to other fingerprinting techniques AFLP shows increased reproducibility and does not require any prior knowledge of the analysed genomes. However, there are some detrimental issues to consider when working with AFLP data; these include potential non-homology and non-independence of fragments, asymmetry in the probability of loss/gain of fragments, and problems in distinguishing heterozygote from homozygote bands e.g. [23, 25]. Despite these difficulties, several authors have used AFLPs to reveal phylogenetic relationships corroborated by analyses of other types of data, especially for species that have diverged recently or radiated within a short period of time e.g. [15, 17, 23, 26].
After excluding 186 replicates the final matrix used for analyses contained 192 individuals and 792 fragments. The AFLP profiles showed good reproducibility with a mean error-rate of 2.4% across all replicated samples. Because the focus of this study was on the phylogenetic relationships between species and species limits rather than intra-specific population genetics, we are presenting and discussing mostly the results of inter-specific relationships. We are presenting here only unrooted trees due to the low resolution of their backbone. We analysed the data using neighbour-joining (NJ) dendrograms and principal coordinate analysis (PCO) with different distance methods, and in both cases the Dice distance gave the highest resolution of relationships between species.
Results of different AMOVAs conducted
Percentage of variation
Among populations within groups
Structure no admixed
The average gene diversity over loci within populations ranged from 0.03 in D. erudita (population 4) to 0.12 in D. parviflora (population 22). Contrary to predictions, the highest number of polymorphic sites, pair-wise differences and average gene diversity were not found in the admixed populations (according to the Structure results) but in D. parviflora (for details see Additional file 2).
“Explosive” radiations featuring rapid opportunistic morphological and ecological diversification are phenomena previously reported for some islands (e.g.  and references therein). Extreme ancestral bottlenecks, together with on-going hybridization and incomplete lineage sorting, can prevent phylogenetic reconstruction in cases of island radiations if they have been recent and produced many species . However, a good understanding of phylogenetic relationships within radiating groups is key for further evolutionary studies into mechanisms and whether change is adaptive, due to drift in small populations or other phenomena .
For the endemic New Caledonian Diospyros species, previous studies, based on multiple plastid  and low-copy nuclear  markers, showed 21 species to be closely related (Figure 1) and were not able to clearly resolve phylogenetic relationships among them. In the combined data set (plastid and nuclear markers; ) only seven of the 21 species included were found to form highly supported groups of accessions from single species. Individuals belonging to each of the remaining 14 species failed to cluster according to their taxonomic circumscription. Dating analysis based on plastid and low-copy nuclear markers showed that the common ancestor of this clade of endemic New Caledonian Diospyros species has arrived in New Caledonia around nine million years ago . Diospyros vieillardii has been shown to be sister to the rest of this endemic clade and separated from the rest of the species around 7.2 million years ago.
Results of the current study using genome-wide AFLP markers reveal that most species form unique groups paralleling recognised species. Around one-third (eight species, NJ dendrogram, Figure 2A) and one-half (11 species, Bayesian tree, Figure 2B) of the species, are genetically distinct with high support (Figure 2). However, the overall AFLP results prove unable to clearly resolve the backbone of trees, similar to previous results obtained from analyses of DNA sequence data . Intra-specific variation was greater (~80%) than that found at inter-specific level (~20%). This low ratio of among- versus within-species divergence in the context of considerable morphological and ecological divergence is indicative of a recent diversification . Such a process can explain why we were able to get clear species boundaries for most species but were unable to clearly resolve phylogenetic relationships between them.
Two species that did not form well-defined clades (D. minimifolia and D. parviflora) were previously considered by White  to show variability in leaf morphology that may indicate that they are in fact a collection of several species. For D. minimifolia White  mentioned that the type population (close to population 15 of this study) has smaller leaves compared to other populations of this species. In our results this population clusters together with the majority of the D. minimifolia accessions; the population that is separated from the rest (population 16) is from Gaji. According to White D. parviflora is a wide-spread species, showing considerable variability of leaf morphology even within populations, making it impossible to differentiate these into different species. Our results show all accessions of D. parviflora, except those from Plateau de Tango (population 24), to form a group. All included accessions from D. parviflora are from ultramafic localities.
To our surprise, the AFLP results do not show any significant grouping according to ecological (edaphic, climatic, elevational), geographical or morphological factors (Additional file 3). The two weakly differentiated groups revealed by Structure and PCO also do not correspond to any conspicuous phenotypic characteristics. The allele-frequency divergence between the two groups found by Structure is low, which explains why we did not observe the two groups in the Bayesian and NJ tree-building results. Taken together, these results indicate that positive selection has perhaps acted on few genomic regions  and has resulted in phenotypic diversification of New Caledonian Diospyros. Variation in copy number of specific genomic regions may be an additional aspect of molecular variation that, although invisible to AFLP markers, could form the basis of adaptation to different environmental conditions .
The individuals of D. vieillardii, D. umbrosa and D. flavocarpa form a minimally isolated group (squares in the grey group) in the PCO (Figure 3). Previous phylogenetic analyses (Figure 1) showed these three species to be sister to the rest of the taxa. Due to its morphological and ecological features D. sp. Pic N’ga from Île des Pins could be a hybrid between D. calciphila and D. vieillardii, but D. vieillardii is now not known from this island. In PCO, individuals of this putative species are located between individuals of D. calciphila and D. vieillardii (Figure 3). The split between the two groups observed (Figures 3 and 4) could be relatively old, separating two lineages that developed in isolated regions. For instance, dry periods of the Pleistocene caused aridification in many areas, and some vegetation types persisted only in local refugia e.g. [34–36]. After climatic conditions became more favourable, the two groups probably expanded rapidly into newly suitable habitats where they overlapped; the time scale of these fluctuations (ca. 0.02 – 0.1 myr; ) was probably not enough to allow woody species with long generation time such as Diospyros to diverge and become permanently reproductively isolated . There are a few admixed individuals in the Structure analysis (Figure 4), which implies that hybridization might have played a role in evolution of this group.
Accelerated rates of evolution at few genes as a result of positive selection could have resulted in the morphological and ecological diversification apparent today in this group of New Caledonian Diospyros species. Furthermore, in addition to retention of ancestral polymorphisms, frequent gene flow could have acted against genome-wide genetic differentiation between the species. Barriers to gene flow between these species may be highly porous, with only few genes responsible for ecological and morphological adaptations evolving on distinct trajectories under strong selection, which leaves the rest of their genomes open to gene flow . Finding these few genes with AFLP is realistically improbable because they are a miniscule component in comparison the rest of these genomes. In the case of a recent and rapid radiation in plants, it could be argued that the bulk of regions sampled by AFLP have not evolved quickly enough to accumulate substitutions that could indicate species relationships. Our results are similar to those found in various other island genera (e.g. Araucaria in New Caledonia, ; Ourisia in New Zealand, ).
Diospyros vieillardii, which is sister to the rest of the taxa belonging to this group of New Caledonian endemics [12, 13], is confined to ultramafic soils, which supports the hypothesis of this being an exaptation of the progenitor of this New Caledonian Diospyros clade to ultramafic soils when the whole island was still covered by heavy-metal-rich substrates; similar findings have been made in other plant groups in New Caledonia e.g. . Later, erosion reduced the extent of this geological layer to one third of the island , and existing species began to move onto other substrates where they subsequently diverged, forming distinct species. Such observations have been made in various other New Caledonian groups (e.g. Araucaria, ; Spiraeanthemum, ; Codia, ). A few studies have examined the adaptive basis and processes involved in rapid radiations in New Caledonia e.g.  and Hawai’i (e.g. lobeliads, ; silverswords, ). Linking ecological parameters and/or phenotypic traits associated with speciation has to be done with caution because range alterations, subsequent evolution, and species extinctions might have erased initial signals found in only a few genes. Therefore, the associations observed today may be misleading, and the specific conditions/traits that were indeed linked to speciation, if any, may no longer be present .
Further work involving common garden experiments would provide insights into the effect of environmental conditions on morphological traits and therefore plasticity of genomes of the New Caledonian Diospyros species. Unfortunately, such experiments are time and cost intensive. It is difficult to obtain ripe fruits of all Diospyros species, and in addition it is difficult to germinate and grow them, which is a crucial aspect of conducting such experiments. Reciprocal transplantation of seedlings across environments are of course more easily conducted than common garden experiments, but they are still time consuming and costly; in addition species adapted to one soil type often will not survive when transplanted to other soil types.
Although New Caledonian Diospyros are morphologically and ecologically diverse, they show little genetic divergence (based on DNA sequences and AFLP data). In this case of the endemic clade of New Caledonian Diospyros, AFLP data did not provide enough information to resolve phylogenetic relationships between the species, but it was sensitive enough for testing for the presence of genetic species boundaries. However, the AFLP results exhibit a good correlation with morphology-based species concepts. Further studies of this New Caledonian Diospyros group with deeper sampling of the genome using next generation sequencing methods are needed to get a clearer picture of the processes that formed this group.
Table of accessions; showing all individuals used in this study
D. calciphila F.White
26, littoral forest
(MPU, NOU, P)
D. cherrieri F.White
21, dry forest
24, dry forest
D. erudita F.White
21, dry forest
22, dry forest
WU062858, Chambrey & Turner 20 (NOU)
D. flavocarpa (Vieill. ex P.Parm.) F.White
10, humid mountain forest
12, dense humid mountain forest
D. glans F.White
1, forest near river
NOU053705, NOU030755, WU062846
5, dense forest near road
BT082, BT084, BT087, BT093-BT094
6, forest near river
D. impolita F.White
7, mesophyll forest near beach
D. inexplorata F.White
25, littoral forest
D. labillardierei F.White
10, river edge in mountain forest
13, river edge
D. minimifolia F.White
11, dry forest
18, mesophyll forest near beach
BT263-BT264, BT266-267, BT269-BT270
21, dry forest
D. pancheri Kosterm.
2, forest near road
JM6619, JM6620 (NOU)
5, dense forest near road
D. parviflora (Schltr.) Bakh.
3, wet forest
BT059, BT062-BT063, BT068
4, wet dense forest
BT080, BT085, BT089-BT090
6, forest near river
20, humid forest at low elevation
tree no. 23109
23, mountain forest
D. perplexa F.White
11, forest near river
D. pustulata F.White
8, dry forest
11, dry forest
BT257-BT258, BT265, BT268, BT271-BT272
21, dry forest
D. revolutissima F.White
D. tridentata F.White
15, dry forest at low elevation
D. trisulca F.White
BT185, BT192, BT197, BT199-BT201
14, mountain forest
NOU031344, JM6637 (NOU)
D. umbrosa F.White
BT061, BT065-BT066, BT071, BT073
4, wet dense forest
13, dense humid forest
BT246-BT247, BT251, BT254, BT256
20, humid forest at low elevation
D. veillonii F.White
18, mesophyll forest near beach
D. vieillardii (Hiern) Kosterm.
1, forest near river
4, dry open forest
BT088, BT091-BT092, BT098, BT100
6, forest near river
28, forest near river
D. yahouensis (Schltr.) Kosterm.
19, mesophyll forest
D. sp. Pic N'ga
Preparation and amplification of fragments followed the protocol of Vos et al.  with some modifications. Restriction of genomic DNA with two restriction enzymes and ligation of double-stranded adaptors to the resulting restricted fragments were performed in one step in a thermal cycler (Veriti, AB, Life Technologies; 37°C for 2 h followed by a 30 min hold at 17°C). Reactions comprised 1.1 μL 10x T4 DNA ligase buffer (Promega), 1.1 μL 0.5 M NaCl, 0.55 μL BSA (1 mg/ mL; New England BioLabs), 50 μM MseI adaptors (genXpress), 5 μM EcoRI adaptors (genXpress), 1 U MseI restriction endonuclease (New England BioLabs), 5 U EcoRI restriction endonuclease (New England BioLabs), 1 U T4 DNA ligase (Promega), and 0.5 μg DNA and were made up to a total volume of 11 μL with water. Ligated DNA fragments were diluted 10-fold with sterile water. Preselective amplification reactions contained 1.14 μL 10x RedTaq PCR reaction buffer (Sigma), 0.2 U RedTaq DNA polymerase (Sigma), 0.22 μL dNTPs (10 mM; AB, Life Technologies), 0.58 μL preselective primer pairs (EcoRI-A and MseI-CT, each 5 μM; Sigma), and 2 μL diluted restriction-ligation product, and were brought with water to a total volume of 10 μL. Amplification was carried out in the same machine used for restriction-ligation with the following profile: 2 min at 72°C, 20 cycles of 10 sec denaturing at 94°C, 30 sec annealing at 56°C, 2 min extension at 72°C, and a final extension step for 30 min at 60°C. The preselective PCR products were diluted 10-fold with sterile water. Reactions for selective amplification contained 0.5 μL 10x RedTaq PCR reaction buffer (Sigma), 0.1 U RedTaq DNA polymerase (Sigma), 0.11 μL dNTPs (10 mM; AB, Life Technologies), 0.27 μL MseI-primer (5 μM; Sigma), 0.27 μL EcoRI-primer (1 μM; Sigma), and 1 μL diluted preselective amplification product and were brought to a total volume of 5 μL with water. They were carried out in a GeneAmp PCR System 9700 (AB, Life Technologies) with the following profile: 1 min at 94°C, 9 cycles of 1 sec at 94°C, 30 sec at 65-57°C (reducing the temperature at 1°C per cycle), 2 min at 72°C, 25 cycles of 1 sec at 94°C, 30 sec at 56°C, 2 min at 72°C and a final extension for 30 min at 60°C. The selective PCR products were purified using Sephadex G-50 Superfine (GE Healthcare Bio-Sciences) applied to a MultiScreen-HV 96-Well Plate (Millipore) in three steps of 200 μL each and settled at 750 g (1, 1 and 5 min, respectively). The same speed was used for centrifugation of the samples (5 μL of each selective PCR product), again for 5 min. Two microliters of the eluate were combined with 10 μL HiDi and 0.1 μL GeneScan 500 ROX (AB, Life Technologies) and denatured for 3 min at 95°C before running them on a capillary sequencer (3130xl Genetic Analyzer, AB, Life Technologies).
The selective primer pairs (6Fam-EcoRI-AGC/MseI-CTGA, Vic-EcoRI-ATG/MseI-CTCG and Ned-EcoRI-ATC/MseI-CTGA) were chosen because they generated clear and not too many bands (thus decreasing the risk of fragments co-migrating by chance), with sufficient variability in preliminary tests. Although the genome size of the New Caledonian Diospyros species (1C-value: 1.5 – 2.3 pg; ) is smaller than the mean 1C-value of eudicots (2.7 pg, ), we found the AFLP profiles generated with Msel primers with four selective bases much clearer than those obtained from primers with just three selective bases.
Reproducibility was checked by repeating ca. 80% of the samples. This high number of repetitions was necessary because of initial difficulties with fragment sizing.
Scoring and phylogenetic analysis
Sizing and scoring of the data was performed with GeneMarker v2.2.0 (SoftGenetics). After pre-analysis using default settings, sizing profiles of all samples were checked and where necessary manually corrected. Most of these corrections concerned one of the following peaks of the size standard: 35 bp, 50 bp and 139 bp. These peaks were often not correctly recognized by the GeneMarker program. High-quality sizing-profiles (score > 90) were obtained for all samples. A panel of scorable fragments was established for each primer combination, and fragments between 65 – 510 bp were scored. The relative fluorescent unit (RFU) threshold was set at 40. Automatic scoring was conducted using Local Southern peak call, peak saturation, base line subtraction, spike removal, pull up correction, and a stutter peak filter of 5% (as described in ). The results were exported as presence/absence matrix. The outcome of the automatic scoring was manually checked and corrected for errors. These errors mostly concerned peaks for which shape was atypical. In total 486 samples corresponding to 231 individuals were scored. From 186 individuals replicate samples were performed (between two and five replicates per individual). The differences between the different samples (replicates) were counted and divided by the total number of phenotypic comparisons to get the error rate (calculated according to Bonin et al. ). After initial analysis (neighbour-joining, NJ) of the complete data set, replicates of samples and obviously failed samples were excluded from further analyses. As replicated samples of the corresponding individuals mostly clustered together, selection of samples from each individual for further analyses was random and not according to any pattern or protocol. For the final analyses we ended up with 192 individuals.
All three primer-combinations were combined in a single matrix and analysed together. Different distance measures were tested for their power to resolve relationships with our data set. Distance matrixes were calculated in PAUP* v4b10 (; Nei-Li distance) and SplitsTree v4.12.6 (; uncorrected P, Dice, corrected and uncorrected Hamming). Phylogenetic relationships based on previously mentioned distance matrices were reconstructed using SplitsTree v4.12.6  to create unrooted NJ dendrograms. To assess robustness of branches NJ-bootstrap (NJ-BS) analyses were performed using SplitsTree v4.12.6  and PAUP* v4b10 . Bayesian inference (BI) was conducted with BEAST v1.7.5 , with two runs each 20 million generations, sampling every 1,000th generation and removal of the first 30% of trees as burn in.
To visualise the pattern of genetic clustering of individuals and populations, we plotted principal coordinate analysis (PCO) using the R-package scatterplot3d  based on an individual Dice distance matrix, and respectively, on AMOVA-derived pair-wise F ST distances calculated with Arlequin v126.96.36.199 . To investigate further significant groupings of the included individuals we used the program Structure v2.3.3 [55, 56] on the Bioportal computing cluster of the University Oslo . We ran Structure for K = 1–23 with 10 replicates each and a model based on admixture and independent allelic frequencies, without taking into account information regarding sampling localities. Each run had 3 million iterations with 10% additional burn in. The calculation of deltaK (∆K; ) and preparation of the input file for Clumpp was done with Harvester . Production of a combined file from the ten replicates of the best K was perfomed using Clumpp v1.1.2  with the full search algorithm. The graphical representation of Structure results was prepared with Distruct v1.1 .
Both non-hierarchical and hierarchical analyses of molecular variance (AMOVA) and calculations of population statistics were conducted using Arlequin v188.8.131.52 . For hierarchical AMOVAs groups have been defined based on different possible clusterings (Additional file 4) according to Structure results, taxonomy, distribution patterns and ecological traits.
Availability of supporting data
AFLP presence/absence matrix and phylogenetic analyses are deposited in treeBASE under study 14798 (http://purl.org/phylo/treebase/phylows/study/TB2:S14798).
This work was funded by a grant from the Austrian Science Fund (FWF, Project-Number: P 22159-B16) awarded to R. Samuel. The authors thank the team of the Department of Systematic and Evolutionary Botany of the University of Vienna as well as the team of the IRD AMAP Noumea for support with this study. Special thank goes to V Klenja for the lab work. Thanks to the following persons for their help with lab work, field work and ideas to improve this work: F Balao Robles, J-P Butin, C Chambrey, G Dagostini, E Grasserbauer, V Hequet, G Kohl, D & I Létocart, W Nigote and H Vandrot. Specimens are deposited in the herbaria of Noumea (NOU), University of Montpellier II (MPU) and the University of Vienna (WU).
- Givnish TJ, Millam KC, Mast AR, Paterson TB, Theim TJ, Hipp AL, Henss JM, Smith JF, Wood KR, Sytsma KJ: Origin, adaptive radiation and diversification of the Hawaiian lobeliads (Asterales: Campanulaceae). Proc R Soc B. 2009, 276: 407-146. 10.1098/rspb.2008.1204.PubMed CentralPubMedView ArticleGoogle Scholar
- Knope ML, Morden CW, Funk VA, Fukami T: Area and the rapid radiation of Hawaiian Bidens (Asteraceae). J Biogeogr. 2012, 39: 1206-1216. 10.1111/j.1365-2699.2012.02687.x.View ArticleGoogle Scholar
- Mittermeier RA, Gil PR, Hoffmann M, Pilgrim J, Brooks T, Mittermeier CG, Lamoreux J, da Fonseca GAB: Hotspots revisited: Earth’s biologically richest and most endangered terrestrial ecoregions. 2004, Mexico City: CEMEXGoogle Scholar
- Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J: Biodiversity hotspots for conservation priorities. Nature. 2000, 403: 853-858. 10.1038/35002501.PubMedView ArticleGoogle Scholar
- Morat P, Jaffré T, Tronchet F, Munzinger J, Pillon Y, Veillon J-M, Chalopin M: Le référentiel taxonomique Florical et les caractéristiques de la flore vasculaire indigène de la Nouvelle-Calédonie. Adansonia. 2012, 34: 177-219.View ArticleGoogle Scholar
- Lowry PP: Diversity, endemism and extinction in the flora of New Caledonia: a review. Rare, threatened, and endangered floras of Asia and the Pacific rim. Edited by: Peng CF, Lowry PPII. 1998, Taiwan: Institute of Botany, Taipei, 181-206.Google Scholar
- Pelletier B, Payri C, Richer De Forges B: Geology of the New Caledonia region and its implications for the study of the New Caledonian biodiversity. Compendium of marine species from New Caledonia, Documents Scientifiques et Techniques II4. 2006, New Caledonia: Institut de Recherche pour le Développement Nouméa, 17-30.Google Scholar
- Maurizot P, Vendé-Leclerc M: New Caledonia geological map, scale 1/500000. Direction de l’Industrie, des Mines et de l’Energie - Service de la Géologie de Nouvelle-Calédonie, Bureau de Recherches Géologiques et Minières. 2009Google Scholar
- Pillon Y, Munzinger J, Amir H, Lebrun M: Ultramafic soils and species sorting in the flora of New Caledonia. J Ecol. 2010, 98: 1108-1116. 10.1111/j.1365-2745.2010.01689.x.View ArticleGoogle Scholar
- Jaffré T, Rigault F, Munzinger J: La végétation. Atlas de la Nouvelle-Calédonie. Edited by: Bonvallot J, Gay J-C, Habert E. 2012, Nouméa: IRD Editions, 77-80.Google Scholar
- Duangjai S, Wallnöfer B, Samuel R, Munzinger J, Chase MW: Generic delimitation and relationships in Ebenaceae sensu lato: evidence from six plastid DNA regions. Am J Bot. 2006, 93: 1808-1827. 10.3732/ajb.93.12.1808.PubMedView ArticleGoogle Scholar
- Duangjai S, Samuel R, Munzinger J, Forest F, Wallnöfer B, Barfuss MHJ, Fischer G, Chase MW: A multi-locus plastid phylogenetic analysis of the pantropical genus Diospyros (Ebenaceae), with an emphasis on the radiation and biogeographic origins of the New Caledonian endemic species. Mol Phylogenet Evol. 2009, 52: 602-620. 10.1016/j.ympev.2009.04.021.PubMedView ArticleGoogle Scholar
- Turner B, Munzinger J, Duangjai S, Temsch EM, Stockenhuber R, Barfuss MHJ, Chase MW, Samuel R: Molecular phylogenetic of New Caledonian Diospyros (Ebenaceae) using plastid and nuclear markers. Mol Phylogenet Evol. 2013, 69: 740-763. 10.1016/j.ympev.2013.07.002.PubMed CentralPubMedView ArticleGoogle Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, Van de Lee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M, Zabeau M: AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res. 1995, 23: 4407-4414. 10.1093/nar/23.21.4407.PubMed CentralPubMedView ArticleGoogle Scholar
- Tremetsberger K, Stuessy TF, Kadlec G, Urtubey E, Baeza CM, Beck SG, Valdebenito HA, Ruas CF, Matzenbacher NI: AFLP phylogeny of South American species of Hypochaeris (Asteraceae, Lactuceae). Syst Bot. 2006, 31: 610-626. 10.1600/036364406778388520.View ArticleGoogle Scholar
- Koopman WJM, Zevenbergen MJ, van den Berg RG: Species relationships in Lactuca s.l. (Lactuceae, Asteraceae) inferred from AFLP fingerprints. Am J Bot. 2001, 88: 1881-1887. 10.2307/3558364.PubMedView ArticleGoogle Scholar
- Richardson JE, Fay MF, Cronk QCB, Chase MW: Species delimitation and the origin of populations in island representatives of Phylica (Rhamnaceae). Evolution. 2003, 57: 816-827.PubMedView ArticleGoogle Scholar
- Despré L, Giells L, Redoutet B, Taberlet P: Using AFLP to resolve phylogenetic relationships in a morphologically diversified plant species complex when nuclear and chloroplast sequences fail to reveal variability. Mol Phylogenet Evol. 2003, 27: 185-196. 10.1016/S1055-7903(02)00445-1.View ArticleGoogle Scholar
- Paun O, Schönswetter P, Winkler M, Tribsch A, IntraBioDiv Consortium: Historical divergence vs. contemporary gene flow: evolutionary history of the calcicole Ranunculus alpestris group (Ranunculaceae) in the European Alps and the Carpathians. Mol Ecol. 2008, 17: 4263-4275. 10.1111/j.1365-294X.2008.03908.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Schulte K, Silvestro D, Kiehlmann E, Vesely S, Novoa P, Zizka G: Detection of recent hybridization between sympatric Chilean Puya species (Bromeliaceae) using AFLP markers and reconstruction of complex relationships. Mol Phylogenet Evol. 2010, 57: 1105-1119. 10.1016/j.ympev.2010.09.001.PubMedView ArticleGoogle Scholar
- Jabaily RS, Sytsma KJ: Historical biogeography and life-history evolution of Andean Puya (Bromeliaceae). Bot J Linn Soc. 2012, 171: 201-224.View ArticleGoogle Scholar
- Gaudeul M, Rouhan G, Gardner MF, Hollingsworth PM: AFLP markers provide insights into the evolutionary relationships and diversification of New Caledonian Araucaria species (Araucariaceae). Am J Bot. 2012, 99: 68-81. 10.3732/ajb.1100321.PubMedView ArticleGoogle Scholar
- Koopman WJM: Phylogenetic signal in AFLP data sets. Syst Biol. 2005, 54: 197-217. 10.1080/10635150590924181.PubMedView ArticleGoogle Scholar
- Degnan JH, Rosenberg NA: Discordance of species trees with their most likely gene trees. PLoS Genet. 2006, 2: 762-768.View ArticleGoogle Scholar
- Meudt HM, Clarke AC: Almost forgotten or latest practice? AFLP applications, analyses and advances. Trends Plant Sci. 2007, 12: 106-117. 10.1016/j.tplants.2007.02.001.PubMedView ArticleGoogle Scholar
- Bussell JD, Waycott M, Chappill JA: Arbitrarily amplified DNA markers as characters for phylogenetic inference. Perspect Plant Ecol, Evol Systematics. 2005, 7: 3-26. 10.1016/j.ppees.2004.07.001.View ArticleGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software Structure: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedView ArticleGoogle Scholar
- Vigouroux Y, Glaubitz JC, Matsouka Y, Goddman MM, Sánchez GJ, Doebley J: Population structure and genetic diversity of New World maize races assessed by DNA microsatellites. Am J Bot. 2008, 94: 1240-1253.View ArticleGoogle Scholar
- Glor RE: Phylogenetic insights on adaptive radiation. Ann Rev Ecol, Evol Systematics. 2010, 41: 251-270. 10.1146/annurev.ecolsys.39.110707.173447.View ArticleGoogle Scholar
- Lerner HRL, Meyer M, James HF, Hofreiter M, Fleischer RC: Multilocus resolution of phylogeny and timescale in the extant adaptive radiation of Hawaiian honeycreepers. Curr Biol. 2011, 21: 1838-1844. 10.1016/j.cub.2011.09.039.PubMedView ArticleGoogle Scholar
- White F: Flore de la Nouvelle-Calédonie et Dépendances. 19. Ébénacées. 1993, Paris: Muséum National d’Histoire NaturelleGoogle Scholar
- Kapralov MV, Votintseva AA, Filatov DA: Molecular adaptation during a rapid adaptive radiation. Mol Biol Evol. 2013, 30: 1051-1059. 10.1093/molbev/mst013.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmidt JM, Good RT, Appleton B, Sherrard J, Raymant GC: Copy number variation and transposable elements feature in recent ongoing adaptation at the Cyp6g1 locus. PLoS Genet. 2010, 6: e1000998-10.1371/journal.pgen.1000998.PubMed CentralPubMedView ArticleGoogle Scholar
- Pintaud J-C, Tanguy J, Puig H: Chorology of New Caledonian palms and possible evidence of Pleistocene rain forest refugia. C R Acad Sci. 2011, 324: 453-463.View ArticleGoogle Scholar
- Pillon Y, Hopkins HC, Munzinger J, Amir H, Chase MW: Cryptic species, gene recombination and hybridization in the genus Spiraeanthemum (Cunoniaceae) from New Caledonia. Bot J Linn Soc. 2009, 161: 137-152. 10.1111/j.1095-8339.2009.00997.x.View ArticleGoogle Scholar
- Poncet V, Munoz F, Munzinger J, Pillon Y, Gomez C, Couderc M, Tranchant-Dubreuil C, Hamon S, de Kochko A: Phylogeography and niche modelling of the relict plant Amborella trichopoda (Amborellaceae) reveal multiple Pleistocene refugia in New Caledonia. Mol Ecol. 2013, doi:10.1111/mec.12554Google Scholar
- Bennett KD: Milankovitch cycles and their effects on species in ecological and evolutionary time. Paleobiology. 1990, 16: 11-21.Google Scholar
- Kane NC, King MG, Barker MS, Raduski A, Karrenberg S, Yatabe Y, Knapp SJ, Rieseberg LH: Comparative genomic and population genetic analyses indicate highly porous genomes and high levels of gene flow between divergent Helianthus species. Evolution. 2009, 63: 2061-2075. 10.1111/j.1558-5646.2009.00703.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Meudt HM, Lockhart PJ, Bryant D: Species delimitation and phylogeny of a New Zealand plant species radiation. BMC Evol Biol. 2009, 9: 111-10.1186/1471-2148-9-111.PubMed CentralPubMedView ArticleGoogle Scholar
- Pillon Y, Munzinger J, Amir H, Hopkins HC, Chase MW: Reticulate evolution on a mosaic of soils: diversification of the New Caledonian endemic genus Codia (Cunoniaceae). Mol Ecol. 2009, 18: 2263-2275. 10.1111/j.1365-294X.2009.04178.x.PubMedView ArticleGoogle Scholar
- Murienne J, Guilbert E, Grandcolas P: Species diversity in the New Caledonian endemic genera Cephalidiosus and Nobarnus (Insecta: Heteroptera: Tingidae), an approach using phylogeny and species distribution modelling. Biol J Linn Soc. 2009, 97: 177-184. 10.1111/j.1095-8312.2008.01184.x.View ArticleGoogle Scholar
- Givnish TJ, Montgomery RA, Goldstein G: Adaptive radiation of photosynthetic physiology in the Hawaiian lobeliads: light regimes, static light responses, and whole-plant compensation points. Am J Bot. 2004, 91: 228-246. 10.3732/ajb.91.2.228.PubMedView ArticleGoogle Scholar
- Baldwin BG, Sanderson MJ: Age and rate of diversification of the Hawaiian silversword alliance (Compositae). Proc Natl Acad Sci USA. 1998, 95: 9402-9406. 10.1073/pnas.95.16.9402.PubMed CentralPubMedView ArticleGoogle Scholar
- Barraclough TG: What can phylogenetics tell us about speciation in the Cape flora?. Divers Distrib. 2006, 12: 21-26. 10.1111/j.1366-9516.2006.00208.x.View ArticleGoogle Scholar
- Ibanez T, Munzinger J, Dagostini G, Hequet V, Rigault F, Jaffré T, Birnbaum P: Structural and floristic diversity of mixed tropical rainforest in New Caledonia: New data from the New Caledonian Plant Inventory and Permanent Plot Network (NC-PIPPN). Appl Veg Sci. doi:10.1111/avsc.1270Google Scholar
- Tel-Zur N, Abbo S, Myslabodski D, Mizrahi Y: Modified CTAB procedure for DNA isolation from epiphytic cacti of genera Hylocereus and Selenicereus (Cactaceae). Plant Mol Biol Report. 1999, 17: 249-254. 10.1023/A:1007656315275.View ArticleGoogle Scholar
- Bennett MD, Leitch IJ: Angiosperm DNA C-values database (release 8.0, Dec. 2012). http://www.kew.org/cvalues/,
- Safer S, Tremetsberger K, Guo Y-P, Kohl G, Samuel MR, Stuessy TF, Stuppner H: Phylogenetic relationships in the genus Leontopodium (Asteraceae: Gnaphalieae) based on AFLP data. Bot J Linn Soc. 2011, 165: 364-377. 10.1111/j.1095-8339.2011.01117.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Bonin A, Bellemain E, Bronken Eidesen P, Pompanon F, Brochmann C, Taberlet P: How to track and assess genotyping errors in population genetic studies. Mol Ecol. 2004, 13: 3261-3273. 10.1111/j.1365-294X.2004.02346.x.PubMedView ArticleGoogle Scholar
- Swofford DL: PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. 2003, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23: 254-267.PubMedView ArticleGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.PubMed CentralPubMedView ArticleGoogle Scholar
- Ligges U, Mächler M: Scatterplot3d - an R package for visualizing multivariate data. J Stat Softw. 2003, 8: 1-20.View ArticleGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005, 1: 47-50.Google Scholar
- Pritchard JK, Stephens M, Donnely P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Hubisz MJ, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009, 9: 1322-1332. 10.1111/j.1755-0998.2009.02591.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Lifeportal. http://www.uio.no/english/services/it/research/hpc/lifeportal/,
- Earl DA, Von Holdt BM: Structure Harvester: a website and program for visualizing Structure output and implementing the Evanno method. Conserv Genet Resour. 2012, 4: 359-361. 10.1007/s12686-011-9548-7.View ArticleGoogle Scholar
- Jakobsson M, Rosenberg NA: CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007, 23: 1801-1806. 10.1093/bioinformatics/btm233.PubMedView ArticleGoogle Scholar
- Rosenberg NA: Distruct: a program for the graphical display of population structure. Mol Ecol Notes. 2004, 4: 137-138.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.