Age estimates for an adaptive lake fish radiation, its mitochondrial introgression, and an unexpected sister group: Sailfin silversides of the Malili Lakes system in Sulawesi
© Stelbrink et al.; licensee BioMed Central Ltd. 2014
Received: 17 April 2014
Accepted: 22 April 2014
Published: 3 May 2014
The Malili Lakes system in central Sulawesi (Indonesia) is a hotspot of freshwater biodiversity in the Wallacea, characterized by endemic species flocks like the sailfin silversides (Teleostei: Atherinomorpha: Telmatherinidae) radiation. Phylogenetic reconstructions of these freshwater fishes have previously revealed two Lake Matano Telmatherina lineages (sharpfins and roundfins) forming an ancient monophyletic group, which is however masked by introgressive hybridization of sharpfins with riverine populations. The present study uses mitochondrial data, newly included taxa, and different external calibration points, to estimate the age of speciation and hybridization processes, and to test for phylogeographic relationships between Kalyptatherina from ancient islands off New Guinea, Marosatherina from SW Sulawesi, and the Malili Lakes flock.
Contrary to previous expectations, Kalyptatherina is the closest relative to the Malili Lakes Telmatherinidae, and Marosatherina is the sister to this clade. Palaeogeographic reconstructions of Sulawesi suggest that the closer relationship of the Malili Lakes radiation to Kalyptatherina might be explained by a 'terrane-rafting’ scenario, while proto-Marosatherina might have colonized Sulawesi by marine dispersal. The most plausible analysis conducted here implies an age of c. 1.9 My for the onset of divergence between the two major clades endemic to Lake Matano. Diversification within both lineages is apparently considerably more recent (c. 1.0 My); stream haplotypes present in the sharpfins are of even more recent origin (c. 0.4 My).
Sulawesi’s Telmatherinidae have most likely originated in the Sahul Shelf area, have possibly reached the island by both, marine dispersal and island/terrane-rafting, and have colonized the Malili Lakes system from rivers. Estimates for the split between the epibenthic sharpfins and the predominantly pelagic to benthopelagic roundfins in Lake Matano widely coincide with geological age estimates of this rift lake. Diversification within both clades clearly predates hybridization events with stream populations. For Lake Matano, these results support a scenario of initial benthic-pelagic divergence after colonization of the lake by riverine populations, followed by rapid radiation within both clades within the last 1 My. Secondary hybridization of stream populations with the sharpfins occurred more recently, and has thus most likely not contributed to the initial divergence of this benthic species flock.
KeywordsMolecular clock Mitochondrial DNA Southeast Asia Sulawesi Biogeography Adaptive radiation Introgressive hybridization
Adaptive radiations of organisms restricted to habitat islands are among the prime model systems for investigating speciation processes in nature [1–4]. Estimations of rates of radiation, as well as the reconstruction of past geographic scenarios of divergence, require the incorporation of a temporal axis into phylogenetic reconstructions [5–7]. Likewise, the dating of phylogenetic splits may contribute to the understanding of other processes potentially contributing to evolutionary divergence (and adaptive radiation), such as introgressive hybridization [8, 9].
The ancient lakes of Sulawesi’s central highlands are a hotspot of aquatic diversity, strongly dominated by endemic species flocks, including radiations of freshwater fishes . The Malili Lakes system (Matano, Mahalona, and Towuti) – a hydrological chain of three main lakes interconnected by rivers – is the main habitat of the sailfin silversides radiation (Atherinomorpha: Telmatherinidae). About 30 morphospecies of these sexually dimorphic fishes that typically possess bright male colourations are distinguished in the Malili Lakes drainage and a few surrounding rivers . Sailfin silversides show conspicuous colour polymorphisms [12, 13], and there are clear indications that ecological speciation has shaped their adaptive radiation [14–16]. Phylogenetic analyses suggest that the two lineages of Telmatherina radiating in the hydrological head of the lakes system, ancient graben-lake Matano, form an ancient monophyletic group that was later compromised by introgressive hybridization from stream populations . This introgression has affected only the “sharpfins”, a lineage of predominantly epibenthic sailfin silversides, whereas their rather pelagic sister group, the “roundfins”, show no indications of genetic exchange with stream populations ([17, 18]; see  for discussion).
Morphological data support a clade composed of Marosatherina ladigesi from south-west Sulawesi, and Kalyptatherina helodes from the islands Batanta and Misool off the Vogelkop peninsula (Birds’ Head, New Guinea), as sister group of the Telmatherinidae in and close by the Malili Lakes . However, only Marosatherina has been considered as an outgroup to the lakes radiation in genetic studies so far [11, 21]. Likewise, Telmatherinidae are represented only by Marosatherina, a species available worldwide by the aquarium trade, in most phylogenetic studies targeting relationships within the Atherinomorpha (e.g., [22–24]).
Here, we use mitochondrial data to (i) assess the relationship of Kalyptatherina relative to Marosatherina and the Malili Lakes Telmatherinidae, (ii) estimate the age of the sailfin silverside radiation of the Malili Lakes, and (iii) provide an estimation of the age of the mitochondrial introgression present in Lake Matano’s sharpfin Telmatherina. For this, we combine sequence data of the Telmatherinidae with data of representatives of the Melanotaeniidae, the closely related rainbowfishes from Australia and New Guinea [22–25], and estimate divergence times using both indirect and geological calibration points, and substitution rates suggested by  as telmatherinid fossils are missing.
DNA extraction, amplification and sequencing
Primers used in the present study (see Methods for PCR conditions)
5′-CAT ACC CCA AAC ATG TTG GT-3′
5′-GTS GST TTT CAC TCC CGC TTA-3′
5′-CTA CCT GAA GAG ATC AAA AC-3′
5′-CGC GTT TAG CTG TTA ACT AA-3′
5′-TGA AGG AGG ATT TAG CAG TAA G-3′
5′-TCT CTG TGG CAA AAG AGT-3′
5′-AAG TGA TTG CGC TAC CTT CGC AC-3′
Phylogenetic analyses and estimation of divergence times
Single gene sequences were aligned using MAFFT (; default settings; http://www.ebi.ac.uk/Tools/msa/mafft) and corrected by eye; ambiguous sites at the 5′ end of the 12S-16S fragment were removed manually. Both single gene alignments were concatenated using SequenceMatrix v. 1.7.8.  resulting in a total alignment of 2,102 bp. The GTR + G substitution model was used for maximum likelihood (ML) analyses using RAxML BlackBox (; partition model; 100 bootstrap replicates) and Bayesian inference (BI) using MrBayes v. 3.1.2 (; partition model; ngen = 10,000,000; nchains = 4; samplefreq = 500; burnin = 10,001). Iso rhotophilus was used as outgroup in all analyses as suggested in ; phylograms were visualized and re-rooted in FigTree v. 1.4 (available at: http://tree.bio.ed.ac.uk/software/figtree).
Divergence time estimates for selected nodes
Mean age (lower and upper 95% HPD) [My]
103.09 (82.70, 122.15)
70.36 (39.70, 106.48)
33.15 (23.43, 43.20)
31.31 (16.40, 48.90)
93.46 (69.89, 116.89)
64.19 (36.83, 95.98)
30.46 (22.53, 39.22)
28.65 (14.92, 43.70)
83.62 (59.92, 107.04)
57.69 (33.54, 85.19)
27.47 (20.40, 34.95)
26.03 (13.56, 39.08)
67.20 (40.38, 93.53)
46.10 (22.62, 70.57)
21.88 (13.85, 31.38)
20.40 (9.07, 32.21)
41.45 (18.04, 65.92)
22.64 (8.27, 39.89)
10.44 (3.89, 17.11)
9.95 (3.14, 18.18)
75.85 (53.25, 98.86)
52.00 (29.69, 76.73)
24.94 (18.31, 31.63)
23.56 (12.17, 35.75)
54.59 (35.08, 75.21)
39.49 (24.69, 56.03)
18.74 (13.30, 24.86)
17.41 (9.49, 26.50)
40.78 (24.86, 58.33)
30.34 (20.55, 40.64)
13.09 (9.22, 13.73) | 32.7 (28.4, 37.4)
12.74 (7.54, 18.43)
34.50 (20.65, 49.72)
25.85 (18.83, 32.88)
11.04 (7.69, 14.78) | 27.0 (23.8, 30.8)
10.87 (6.70, 14.98)
64.76 (43.47, 86.22)
43.85 (25.27, 65.90)
21.35 (15.20, 27.93)
20.19 (10.38, 30.97)
57.80 (37.80, 78.33)
38.95 (22.54, 59.33)
19.10 (13.38, 25.44)
18.11 (9.05, 27.98)
37.88 (18.74, 56.83)
25.00 (11.00, 39.86)
13.15 (6.95, 19.24)
12.30 (4.90, 20.65)
42.86 (25.60, 62.34)
27.17 (14.37, 42.94)
13.60 (8.89, 18.66)
12.93 (6.18, 20.36)
28.90 (15.01, 43.48)
19.15 (9.56, 30.43)
9.11 (5.50, 13.06)
8.45 (3.79, 13.45)
17.43 (8.21, 27.65)
11.99 (5.78, 19.06)
5.35 (3.22, 8.12)
5.18 (2.15, 8.53)
9.59 (3.52, 17.08)
7.13 (3.21, 12.20)
3.12 (1.47, 4.91)
3.00 (1.01, 5.28)
6.43 (2.40, 11.60)
4.45 (1.95, 7.62)
1.98 (0.90, 3.20)
1.86 (0.63, 3.27)
3.66 (1.18, 7.03)
2.45 (0.81, 4.31)
1.04 (0.39, 1.82)
1.01 (0.31, 1.92)
3.05 (0.90, 5.81)
2.24 (0.73, 4.15)
0.99 (0.36, 1.78)
0.93 (0.27, 1.80)
13.94 (6.35, 22.44)
10.15 (4.84, 16.54)
4.40 (2.32, 6.70)
4.32 (1.72, 7.17)
3.60 (0.73, 7.77)
2.92 (0.64, 6.29)
1.27 (0.33, 2.68)
1.22 (0.25, 2.68)
1.21 (0.22, 2.61)
0.92 (0.17, 1.98)
0.40 (0.09, 0.83)
0.37 (0.07, 0.83)
1.20 (0.08, 2.97)
0.96 (0.08, 2.44)
0.41 (0.04, 1.00)
0.38 (0.02, 0.95)
10.23 (4.58, 17.24)
7.24 (3.09, 11.98)
3.18 (1.54, 5.07)
3.03 (1.19, 5.33)
6.14 (2.08, 11.43)
4.39 (1.44, 8.03)
1.92 (0.73, 3.42)
1.82 (0.48, 3.45)
2.36 (0.64, 4.63)
1.68 (0.50, 3.26)
0.73 (0.23, 1.34)
0.70 (0.18, 1.37)
5.05 (1.54, 9.36)
3.84 (1.09, 7.02)
1.67 (0.61, 3.02)
1.61 (0.48, 3.07)
Resulting rate [% / My]
Phylogenetic relationships and age estimates
The Melanotaeniidae form a highly supported monophyletic group, including Rhadinocentrus ornatus from Queensland, and all species of the genera Melanotaenia, Chilatherina and Glossolepis analyzed. In this group, R. ornatus represents the sister taxon to three distinct and highly supported clades from New Guinea (and surrounding islands), and northern Australia (Figure 1). Interestingly, Glossolepis dorityi and G. pseudoincisus share one haplotype (Figure 1), which is most likely due to gene flow between the two species inhabiting the very same river system in the northern lowlands of New Guinea. Two Melanotaenia species from the Bird’s Head (M. angfa) and Batanta Island (M. batanta) form a monophyletic group, being sister to two separate clades comprising Melanotaenia (sub)species from southern New Guinea and northern Australia, and species of the genera Chilatherina, Glossolepis and Melanotaenia from northern New Guinea (cf. “western”, “southern”, and “northern” clade in ; Figure 1). This is largely congruent with findings of . The inferred mean age estimates for the western–northern/southern split (node 8 in Figure 1) and the northern–southern split (node 9 in Figure 1) range from 12.7-40.8 and 10.9-34.5 My, respectively, depending on if node 9 was enforced in the respective analysis (see Table 2, analyses [B and D]).
Phylogenetic reconstructions (ML and BI) show a well to highly supported clade comprising members of the genus Pseudomugil (P. signifer, P. reticulatus, P. furcatus, and P. pellucidus; note that the Pseudomugilidae are paraphyletic; see also below), and the Telmatherindae. Within the Telmatherinidae, Marosatherina ladigesi from the Maros karst (southwest Sulawesi) is basal to the clade including Kalyptatherina helodes (Batanta and Misool) plus the sailfin silversides from the Malili Lakes system in Central Sulawesi, Indonesia (node 14; mean age range: 8.5-28.9 My). Estimated divergence times range from 18.1-57.8 My for the split of Telmatherinidae and the most recent Pseudomugilidae clade (node 11), and are quite similar for the first diversification events within each family (12.3-37.9 My for node 12 for the most recent Pseudomugil clade, and 12.9-42.9 My for node 13, respectively).
Three different genera belong to the Malili Lakes species flock, namely Paratherina, Tominanga, and Telmatherina (node 15 = sailfin silversides split: 5.2-17.4 My). The morphologically and nuclear distinct Paratherina (Lakes Mahalona and Towuti; node 24 = TMRCA, time to most recent common ancestor of Paratherina and Telmatherina: 3.0-10.2 My, node 27 = 'speciation’ , the first intra-clade diversification within Paratherina: 1.6-5.1 My) and Tominanga (Lake Mahalona; node 21 = TMRCA of Tominanga and Telmatherina: 1.2-3.6 My, node 23 = 'speciation’: 0.4-1.2 My) are both clearly supported as monophyletic. In contrast, species of the genus Telmatherina fall in three distinct mitochondrial clades. Node 25 (1.8-6.1 My) represents the MRCA of a clade (from now on we refer to MRCAs of particular clades when discussing nodes) is composed of the lake-dwelling Telmatherina celebensis from Lakes Mahalona and Towuti, and several populations of stream-dwelling Telmatherina bonti. Node 22 (0.4-1.2 My) comprises sharpfin specimens from Lake Matano, together with stream-dwellers; node 16 contains the remaining Telmatherina from Lake Matano, namely sharpfins (node 19), roundfins (node 18), and two stream Telmatherina from River Nuha (north of Matano) and River Tominanga (draining Lake Mahalona to Towuti). Taken together, these results are consistent with those reported by . In the light of nuclear and morphological data (cf. ), the mitochondrial signatures provide evidence for substantial hybridization between lake- and stream-dwelling Telmatherina. This results in two haplotype clades, the “original” Matano haplotypes (node 19) being sister to the endemic roundfins (node 18), and the introduced haplotypes of (node 22) within Lake Matano’s sharpfins.
The endemic sailfin silversides radiation of the Malili Lakes serves as a model system in evolutionary ecology research (see [10, 19, 41] for reviews). The detailed and temporal reconstruction of the phylogenetic history of the Telmatherinidae is crucial for exploring the likely conditions underlying speciation processes, including the spatial origin and morphological traits of ancestral and introgressed populations. It remained, however, unclear if it is justified to consider Marosatherina from SW Sulawesi as the sister species to the lakes radiation, and if the age of the lakes radiation in fact falls into preliminary lake age estimations proposed by geologists and limnologists, which are however not yet fully reliable . These preliminary seismic data suggest an age of at least 600,000-700,000 years for Lake Towuti , while the geological fault formation, in which Lake Matano is embedded, might be comparatively older possibly starting around 4 Mya . An onset of the Malili Lakes system formation in the early Pleistocene c. 1–2 Mya is plausible (Robert Hall, pers. comm.).
The closest relative to the sailfin silversides radiation
The phylogenetic reconstructions presented here are based on mitochondrial markers only, and hence do not allow inferences about cyto-nuclear discordances within the sailfin silversides species flock (cf. ). However, this mtDNA dataset allows inferring a first hypothesis for reconstructing phylogenetic relationships among disjunctly distributed telmatherinid species, whose mtDNA haplotypes most likely have preserved the vicariant phylogenetic signal.
Our analyses clearly suggest that mtDNA haplotypes of the Malili Lakes radiation are more closely related to Kalyptatherina helodes, the only telmatherinid species occurring on islands east of Sulawesi (Misool and Batanta), rather than to Marosatherina ladigesi from SW Sulawesi, a species previously considered the sister taxon to the Malili flock. Its inclusion into the Telmatherinidae is highly plausible in the light of morphological data , but its placement as the sister taxon to the lacustrine flocks appears surprising. However, reconstructions of the complex geological history of Sulawesi and adjacent islands may provide explanations for these findings (see section below).
The phylogenetic relationships of the families within Atheriniformes, and the composition of these families, remain partially controversial, most likely due to substantial differences in taxon sampling and the methods applied (Figure 2). Two important taxa could not be included in the present dataset, namely Cairnsichthys rhombosomoides (Queensland, Australia), and the rather widespread brackish water family Phallostethidae (priapumfishes; recorded from Sundaland, Luzon, Palawan and Southwest Sulawesi; [44, 45]. Cairnsichthys is suggested to be basal to Pseudomugil and Telmatherina according to the molecular phylogeny by , while the phallostethids are placed as the sister to Pseudomugil according to . Morphological studies by  suggested Melanotaeniidae as sister to Madagascar’s Bedotiidae, which are in turn most closely related to a clade composed of Telmatherinidae, Pseudomugilidae (blue eyes), and three other atheriniform families. The combined analysis of molecular and morphological data by  supported the sister group relationship of monophyletic Melanotaeniidae and Bedotiidae; a clade composed of Marosatherina (Telmatherinidae) and two Pseudomugil species (Pseudomugilidae) represent the sister clade to Melanotaeniidae and Bedotiidae in that study. A comprehensive study by  reported contrasting results, based on the combined analysis of one mitochondrial and nuclear marker. The single representatives of Pseudomugil and Telmatherina analyzed by  were also supported as most closely related, but nested within a clade composed of Melanotaeniidae and Phallostethidae. Most recently,  confirmed the inclusion of Phallostethidae within the atherinomorphs as distant sister to a monophyletic group of Pseudomugil species, based on 10 nuclear markers. Unmack et al.  found in a comprehensive multilocus molecular study (seven mtDNA markers and one nuclear marker) that Melanotaeniidae are the monophyletic sister group to Madagascar’s Rheocles and Bedotia; Melanotaeniidae were again supported as sister to a clade composed of Marosatherina and Pseudomugil. The results of the present study clearly support the close relationship between Telmatherinidae and Pseudomugilidae reported by previous workers. However, the Pseudomugilidae are also clearly paraphyletic, with P. signifer being basal to all remaining Pseudomugilidae and Telmatherinidae (Figure 1). This finding is consistent with the results of  (see also Figure 2), and highlights the need for a detailed study addressing the systematic position of several Pseudomugil species, especially that of P. signifer.
Divergence time estimates of the split of Kalyptatherina and the central Sulawesi Malili flock of 8.5-28.9 Mya renders a scenario of terrane-rafting the most plausible explanation for the present-day distribution pattern, given the region’s geological history. In contrast, and considering its present local distribution and estimated node ages, the ancestor of Marosatherina most likely colonized western Sulawesi by dispersal from the Sahul Shelf, the current distribution of Pseudomugil and Melanotaeniidae.
A scenario of 'terrane-rafting’ provides a plausible explanation for the sister group relationship between Kalyptatherina, endemic to the small islands off the Vogelkop Peninsula of New Guinea, and the Malili Lakes sailfin silversides. Geological elements formerly belonging to the Sula Spur were in proximity to the Australian margin – including old offshore islands like Batanta and Misool – before this promontory was extended, moved westwards, and finally collided with West Sulawesi (e.g., [37, 47]). It appears most plausible that the population ancestral to the Malili sailfin silversides originates from the Sahul Shelf area, and was dispersed on such a 'terrane raft’ when the Sula Spur was extended and moved westwards until this fragment ('raft’) collided with West Sulawesi. However, given the temporal uncertainties in both the separation of the Sula Spur from the Sahul Shelf (c. 15 Mya; see also ) and divergence time estimates among the four analyses, it remains difficult to test this hypothesis. Mean ages and credibility intervals suggest that this scenario might be plausible for analyses [A] and [B] (15.0-43.5 My and 9.6-30.4 My; see Table 2), while the credibility intervals obtained from analyses [C] and [D] (5.5-13.1 and 3.8-13.5) would slightly post-date the estimated age of the Sula Spur separation. However, the lack of fossil remains requires denser sampling and the incorporation of multiple markers to explain the relationships of present-day geographically distant groups, which probably have been in vicinity in the past.
A marine dispersal explanation for the sister group relationship between Marosatherina, a monotypic genus endemic to the Maros karst in SW Sulawesi, and all the remaining sailfin silversides, appears most likely based on its current distribution and the divergence time estimates inferred (node 13: 12.9-42.9 My). West Sulawesi and the remaining geological parts of the island are of different origin (Sunda Shelf and Sahul Shelf), and amalgamated not until in the Miocene due to tectonic movements. It appears plausible to assume that the population ancestral to Marosatherina might have originated in the Sahul Shelf area and colonized present-day West Sulawesi by marine dispersal; this requires, however, the assumption that the dispersing proto-Marosatherina had a tolerance for marine conditions. Such a salt tolerance might indeed be a plesiomorphic character of sailfin silversides. The occurrence of Kalyptatherina, Neostethus (Phallostethidae, present with one species on Sulawesi, Neostethus djajaorum; ), and also some Pseudomugil species in brackish habitats, provides support for this assumption. Likewise, Marosatherina as well as Telmatherina bonti tolerate brackish waters under aquarium conditions (F.H., pers. obs.). However, none of the Sulawesi sailfin silversides has ever been reported from such habitats in nature, and it remains thus unclear, if the assumption that ancestral sailfin silversides were able to cross marine barriers, is in fact realistic. The endemism of all recent species of the family provides a substantial argument against profound abilities for marine dispersal. As an alternative, partly complementary scenario, the ancestral population of Marosatherina might have been widespread across both Sunda and Sahul Shelf areas, followed by extinction across large extents during periods of sea-level fluctuations. In that case, Marosatherina would represent a relict species as several areas of West Sulawesi remained above sea level or at least were covered only by shallow water during the island’s history according to palaeogeographic reconstructions [47, 52]. However, this again assumes that the ancestral population might have been, at least to some degree, saltwater-tolerant.
In line with earlier workers , the present results support the monophyly of the rainbow fishes (Melanotaeniidae). As expected, its genera Chilatherina and Glossolepis are nested within Melanotaenia, and the three major geographic clades recovered correspond to the expected freshwater ecoregions of “southern”, “northern” and “western” New Guinea [24, 53]. The island’s Central Highlands are the major barrier putatively separating the “northern” and “southern” clades, and provide an opportunity to estimate divergence rates also within the Telmatherinidae (see below).
Divergence time estimation: The Malili Lakes radiation
Depending on the method applied, molecular clock approaches estimate the onset of the Malili Lakes radiation to 5.2-17.4 My, but the youngest estimate appears by far the most plausible, given the estimated timeframe for Sulawesi’s final amalgamation (see above). This implies ages of about 1.9 My for the initial split of the benthic sharpfins and the predominantly pelagic to benthopelagic roundfins, estimates of c. 1.0 My for divergence within these two lineages inside Lake Matano, and substantially less (0.4 My) for the lineage of haplotypes introgressed by stream populations into Matano’s sharpfins.
 provided a first age estimation for the divergence among the three mitochondrial haplotype clades present in Lake Matano’s Telmatherina radiation. This approach was based on a constant rate of evolution, and applied a genetic distance-age ratio of 1-2%/My. For the sailfin silversides endemic to Lake Matano, the deepest and according to geological data (see  and Robert Hall, pers. comm.) oldest lake of the system, these analyses suggested an age of 0.95-1.9 My separating roundfins (“Clade I” in ; see  for the identity of these clades) and sharpfins (“Clade II”). Haplotypes originating from streams and rivers, present in Lake Matano’s sharpfins due to introgressive hybridization , diverged from the lacustrine haplotypes (sharpfins + roundfins) in Roy et al.’s  analyses 1.85-3.7 My ago. Their ingroup dataset did however not include the other sailfin silversides species occurring outside of Lake Matano, from the remaining lakes, rivers, and streams of the Malili Lakes; it appears accordingly unclear if this framework is suited for providing reliable estimates for the relevant splits. In the present study, we combine geological and indirect calibration points, as well as recently suggested substitution rates, to estimate and carefully discuss the timing of the most relevant splits within the sailfin silversides radiation in a relaxed molecular clock framework.
Analysis [A] – indirect calibration for the Iso–Melanotaeniasplit
We first used the split between Iso hawaiiensis and Melanotaenia lacustris as an indirect calibration point, obtained from a study of ricefishes (Adrianichthyidae) by , in analysis [A]. This approach provides an age of 17.4 My for the Malili Lakes radiation (node 15); Lake Matano’s roundfins are estimated to an age of 3.7 Mya (node 18) in that analysis (Table 2). Under analysis [A], the age of the Malili Lakes radiation (node 15) significantly predates the proposed age for the formation of present-day Sulawesi, namely the final amalgamation of the North and West arms of Sulawesi with the Sula Spur (c. 10–20 Mya; [47, 50], as well as the geological evidence for the age of the Lake Matano.
Some technical issues might account for the observed inconsistencies between the proposed ages of these nodes, and the geological reconstruction of Sulawesi’s history and its ancient lakes. Dating based on indirect calibration points presupposes that adequate calibration points were used in the source analysis. The 21 fossil and six biogeographic calibration points used by  are remarkable, but all fall outside the atherinomorphs. These priors, plus possible saturation effects, may have led to the bias of an overestimation of splits within the atherinomorphs. A recently published, and more comprehensive dated phylogeny by , might justify this inference. There, the atherinomorphs are estimated to have originated c. 80 Mya, while the genera Iso and Melanotaenia are more recent compared to the split in Setiamarga et al.’s study (: 93–113 My vs. : c. 50 My); however, some basal nodes remain unresolved in , and thus both genera do not form a sister group, as in . Thus, we conclude that the inferred node ages from analysis [A] very likely substantially overestimate the age of the Malili Lakes radiation, as well as that of the intralacustrine speciation and hybridization processes.
Analysis [B] and [C] – indirect calibration for Melanotaeniaclades and rate of 1%/My
Estimated divergence times for the New Guinean northern and southern Melanotaenia clades obtained from Unmack et al.’s  study have been used as indirect calibration points in analysis [B]. In addition, we used the substitution rate of 1%/My suggested by  in their rainbowfish dataset (analysis [C]). Interestingly, analysis [B] shows comparable ages for node 8 and node 9 compared to , however, the resulting ucld.mean (uncorrelated lognormal relaxed clock mean) is 0.47%/My for this analysis. Accordingly, age estimates inferred from analysis [B] clearly predate the hypothesized starting point of the New Guinean Highlands uplift (node 9: 25.9 My). However, the MRCA of the Melanotaenia clade is 39.5 My, and is therefore quite similar to the divergence time estimates of the single Melanotaenia lineage in .
Likewise, divergence times differ remarkably between the few comparable nodes of analysis [C] and Unmack et al.’s  study (see Table 2). The considerable mismatch between analysis [C] and the study by  might be related to the different prior conditions in the respective analyses. With respect to the target question, the taxon sampling available is more complete in the present study. Differences in the resulting estimates may also be due to the genetic markers applied: substitution rates vary across genetic markers and among taxa, and might therefore cause deviating node ages though using the same prior substitution rate.
Analysis [D] – geological calibration using the Central Highlands formation
In a final analysis, we used age estimations for the Central Highlands in New Guinea, and the clades of rainbow fish separated by this major barrier, for estimating the age of the relevant nodes in the sailfin silversides phylogeny (analysis [D]). The Central Highlands in New Guinea are a massive mountain range, spanning almost the entire island with a maximum altitude of c. 5,000 m above sea level. The continuous formation of the Central Highlands since the late Miocene (c. 10–14 Mya; ) clearly resulted in the interruption of former watersheds, and led to the present-day regions of freshwater biodiversity. This explains the presence of two well-resolved clades of rainbow fishes north and south of this mountain range ([24, 54], this study); the ancestral populations became most likely continuously separated. However, given the nature of geological events, which should in most cases be considered rather as continuous processes than as distinct, precise events, it appears problematic to reconstruct the exact timeframe in which the proposed ancestral melanotaeniid population in New Guinea was initially separated, giving raise to the two distinct clades observed to date.
Compared to the indirect calibration approach, the prior age of node 9 is substantially younger in analysis [D] (mean age: 10.9 My) than the posterior age inferred from analysis [A] (mean age: 34.5 My). This results in a likewise substantially more recent age of both, the onset of the Malili Lakes radiation, and that of its radiating clades (see Table 2 and Figure 3). Accordingly, node ages derived from analysis [D] appear more plausible than those from analyses [A] and [B]. However, the root height representing the age of Iso considerably underestimates the divergence time inferred from  by about 70 My, while in contrast, Iso is estimated to be c. 15 My younger in analysis [D] compared to the study by . This underestimation may not only be due to conflicting topologies, the position of Iso within the atherinomorphs (see Discussion above), and the comparatively recent calibration point used (node 9), but might be also correlated with saturation effects in basal nodes by solely using mitochondrial markers.
Telmatherina of Lake Matano’s endemic sharpfin radiation carry either mitochondrial haplotypes closely related to those of the lakes’ roundfins, or those introgressed by riverine populations [17–19]. The age of the haplotypes introgressed into the Matano flock (mainly node 22) is comparatively young, estimated in analysis [D] to less than 400,000 years. In contrast, age estimates for the “native” Matano sharpfin haplotypes (node 19: 0.9 My), its sister – the roundfins (node 18: 1.0 My) –, and the clade of T. celebensis from the lower lakes of the system (node 26: 0.7 My) are comparatively older, and appear largely congruent to the general age estimates for the Malili Lakes (see above). It appears therefore very likely that the mitochondrial introgression observed has occurred rather recently, in comparison to the age of the lakes’ native haplotypes (see Figure 3). Shared haplotypes in highly distinct lake- and stream-dwelling Telmatherina are also observed in Lakes Towuti’s and Mahalona’s T. celebensis, and several populations of T. bonti (see also Figure 1 and supplementary figure one in , incorporating more stream populations). Node 26, comprising these lacustrine and riverine populations, is estimated to c. 700,000 years ago in analysis [D], a time that coincides with the proposed age for Lake Towuti .
Similarly to the Lake Matano Telmatherina radiation, Paratherina represents a small, monophyletic radiation, occurring in off- and inshore waters of Lakes Mahalona and Towuti. Analysis [D] suggests that the first diversification event within Paratherina (node 27) has occurred c. 1.6 Mya. It appears reasonable to assume that the onset of diversification may have taken place in the older Lake Mahalona, from where the putatively younger Lake Towuti has been colonized. The Paratherina populations of both lakes are likely still connected via Tominanga River, as suggested by the recent microsatellite study of . Further support for possible riverine dispersal of Paratherina comes from the historical presence of species shared with Lake Towuti in the small hill-lake Lontoa (or Wawontoa; see ), connected to the large lake by rivers. This lake has however undergone substantial degradation, and the presence of Paratherina could not be confirmed during recent surveys (F.H., pers. obs.).
Divergence times inferred for the Malili Lakes radiation clearly predate both the final formation of Sulawesi and any suggested age estimates for the Malili Lakes for the majority of analyses performed. Hence, node ages derived from analysis [D], based on geological calibration by the New Guinean highland barrier, seem most plausible to us. This means that divergence times obtained from other sources, i.e., [24, 36], possibly overestimate telmatherinid and melanotaeniid clade ages, again highlighting the issues related with molecular clock analyses (see Discussion above).
Some concluding remarks can be made on speciation and hybridization processes in the sailfin silversides radiation, based on the – in our view – most plausible analysis [D]. Based on the present topology, riverine Telmatherina bonti populations not only cluster within the lacustrine clades, indicating hybridization events; some of these riverine lineages also appear basal to some lacustrine populations. In line with similar results from an earlier phylogenetic study , this clearly indicates that the Malili Lakes were colonized by riverine populations, which appears highly plausible, and meets patterns observed also in invertebrate radiations, like e.g., the pachychilid snails .
The present study indicates that the Sulawesi telmatherinids might have originated c. 3–5 Mya, a period when present-day Sulawesi was being formed through a series of tectonic events such as e.g., mountain uplifts in West Sulawesi and the Matano fault formation [42, 52]. The Matano fault gave rise to the rift lake Matano, and probably also initiated the formation of the remaining lakes of the Malili Lakes system. Preliminary geological reconstructions suggest that the Malili Lakes are generally 1–2 My old, while preliminary seismic data support with 600,000-700,000 years a younger age for Lake Towuti. This geological and seismic evidence is in line with the present phylogenetic reconstruction, which shows comparatively recent diversification and hybridization events within Telmatherina celebensis and their riverine relatives, compared to diversification and introgression in Lake Matano (node 26; see Figures 1 and 3). Age estimates suggest diversification along a benthic-pelagic axis, into sharpfins and roundfins, c. 1.9 Mya after Lake Matano was colonized by stream populations, followed by a rapid radiation in both of these clades in the last 1 My. Secondary hybridization did probably not affect initial divergence within Lake Matano’s sharpfin radiation, as the age of the introgressed haplotypes clearly postdates the initial diversification by about 600,000 years.
Availability of supporting data
We thank the Research Center for Biology, Indonesian Institute of Sciences (LIPI) and Kementerian Riset dan Teknologi Republik Indonesia (RISTEK) for the permit to conduct research in Indonesia. PT. INCO / PT. VALE Indonesia Tbk., provided outstanding logistic support in Sulawesi. We are grateful to Johannes Graf for extensive samples of rainbowfish species and Robert Shelly for providing several atheriniform samples. We thank the Axios Review editor and reviewers (Michael Matschiner and two anonymous reviewers), and Jobst Pfaender for valuable comments on an earlier version of the manuscript. Fieldwork was funded by research grants of the Deutsche Forschungsgemeinschaft, DFG (to UKS: 567/3-1-3; to FH: HE 5707/2-1).
- Schluter D: The Ecology of Adaptive Radiation. 2000, Oxford: Oxford University Press, 1-284.
- Losos JB, Mahler DL: Adaptive radiation: The interaction of ecological opportunity, adaptation, and speciation. Evolution since Darwin: the first 150 years. Edited by: Bell MA, Futuyma DJ, Eanes WF, Levinton JS. 2010, Sunderland, MA: Sinauer Associates, 381-420.
- Martin CH, Wainwright PC: Multiple fitness peaks on the adaptive landscape drive adaptive radiation in the wild. Science. 2013, 339: 208-211. 10.1126/science.1227710.PubMedView Article
- Santos ME, Salzburger W: How cichlids diversify. Science. 2012, 338: 619-621. 10.1126/science.1224818.PubMedView Article
- Bess EC, Catanach TA, Johnson KP: The importance of molecular dating analyses for inferring Hawaiian biogeographical history: A case study with bark lice (Psocidae: Ptycta). J Biogeogr. 2014, 1: 158-167.View Article
- 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 Article
- Friedman M, Keck BP, Dornburg A, Eytan RI, Martin CH, Darrin C, Wainwright PC, Near TJ, Hulsey CD: Molecular and fossil evidence place the origin of cichlid fishes long after Gondwanan rifting. Proc R Soc London B. 2013, 280: 1-8.
- Schwarzer J, Misof B, Ifuta SN, Schliewen UK: Time and origin of cichlid colonization of the lower Congo rapids. PLoS One. 2011, 6: e22380-10.1371/journal.pone.0022380.PubMedPubMed CentralView Article
- Willis SC, Farias IP, Ortí G: Testing mitochondrial capture and deep coalescence in Amazonian cichlid fishes (Cichlidae: Cichla). Evolution. 2014, 68: 256-268. 10.1111/evo.12230.PubMedView Article
- von Rintelen T, von Rintelen K, Glaubrecht M, Schubart CD, Herder F: Aquatic biodiversity hotspots in Wallacea: the species flocks in the ancient lakes of Sulawesi, Indonesia. Biotic Evolution and Environmental Change in Southeast Asia. Edited by: Gower DJ, Johnson KG, Richardson JE, Rosen BR, Rüber L, Williams ST. 2012, Cambridge: Cambridge University Press, 290-315.View Article
- Herder F, Schwarzer J, Pfaender J, Hadiaty RK, Schliewen UK: Preliminary checklist of sailfin silversides (Teleostei: Telmatherinidae) in the Malili Lakes of Sulawesi (Indonesia), with a synopsis of systematics and threats. Verhandlungen der Gesellschaft für Ichthyologie. 2006, 5: 139-163.
- Gray SM, Dill LM, Tantu FY, Loew ER, Herder F, McKinnon JS: Environment-contingent sexual selection in a colour polymorphic fish. Proc R Soc London B. 2008, 275: 1785-1791. 10.1098/rspb.2008.0283.View Article
- Pfaender J, Gray SM, Rick IP, Chapuis S, Hadiaty RK FH: Spectral data reveal unexpected cryptic colour polymorphism in female sailfin silverside fish from ancient Lake Matano. Hydrobiologia. 2013, doi:10.1007/s10750-013-1557-y
- Pfaender J, Schliewen UK, Herder F: Phenotypic traits meet patterns of resource use in the radiation of “sharpfin” sailfin silverside fish in Lake Matano. Evol Ecol. 2010, 24: 957-974. 10.1007/s10682-009-9332-2.View Article
- Pfaender J, Miesen FW, Hadiaty RK, Herder F: Adaptive speciation and sexual dimorphism contribute to diversity in form and function in the adaptive radiation of Lake Matano’s sympatric roundfin sailfin silversides. J Evol Biol. 2011, 24: 2329-2345. 10.1111/j.1420-9101.2011.02357.x.PubMedView Article
- Cerwenka AF, Wedekind JD, Hadiaty RK, Schliewen UK, Herder F: Alternative egg-feeding tactics in Telmatherina sarasinorum, a trophic specialist of Lake Matano’s evolving sailfin silversides fish radiation. Hydrobiologia. 2012, 693: 131-139. 10.1007/s10750-012-1099-8.View Article
- Herder F, Nolte AW, Pfaender J, Schwarzer J, Hadiaty RK, Schliewen UK: Adaptive radiation and hybridization in Wallace’s dreamponds: evidence from sailfin silversides in the Malili lakes of Sulawesi. Proc R Soc London B. 2006, 273: 2209-2217. 10.1098/rspb.2006.3558.View Article
- Schwarzer J, Herder F, Misof B, Hadiaty RK, Schliewen UK: Gene flow at the margin of Lake Matano’s adaptive sailfin silverside radiation: Telmatherinidae of River Petea in Sulawesi. Hydrobiologia. 2008, 615: 201-213. 10.1007/s10750-008-9561-3.View Article
- Herder F, Schliewen UK: Beyond sympatric speciation: radiation of sailfin silverside fishes in the Malili lakes (Sulawesi). Evolution in Action. Edited by: Glaubrecht M. 2010, Berlin: Springer, 465-483.View Article
- Aarn WI, Ivantsoff W, Kottelat M: Phylogenetic analysis of telmatherinidae (Teleostei: Atherinomorpha), with description of Marosatherina, a new genus from Sulawesi. Ichthyol Explor Freshw. 1998, 9: 311-323.
- Roy D, Paterson G, Hamilton PB, Heath DD, Haffner GD: Resource-based adaptive divergence in the freshwater fish Telmatherina from Lake Matano, Indonesia. Mol Ecol. 2007, 16: 35-48.PubMedView Article
- Sparks JS, Smith WL: Phylogeny and biogeography of the Malagasy and Australasian rainbowfishes (Teleostei: melanotaenioidei): gondwanan vicariance and evolution in freshwater. Mol Phylogenet Evol. 2004, 33: 719-734. 10.1016/j.ympev.2004.07.002.PubMedView Article
- Setiamarga DHE, Miya M, Yamanoue Y, Mabuchi K, Satoh TP, Inoue JG, Nishida M: Interrelationships of Atherinomorpha (medakas, flyingfishes, killifishes, silversides, and their relatives): the first evidence based on whole mitogenome sequences. Mol Phylogenet Evol. 2008, 49: 598-605. 10.1016/j.ympev.2008.08.008.PubMedView Article
- Unmack PJ, Allen GR, Johnson JB: Phylogeny and biogeography of rainbowfishes (Melanotaeniidae) from Australia and New Guinea. Mol Phylogenet Evol. 2013, 67: 15-27. 10.1016/j.ympev.2012.12.019.PubMedView Article
- Parenti LR: Relationships of atherinomorph fishes (Teleostei). Bull Mar Sci. 1993, 52: 170-196.
- Kocher TD, Conroy JA, McKaye KR, Stauffer JR, Lockwood SF: Evolution of NADH dehydrogenase subunit 2 in East African cichlid fish. Mol Phylogenet Evol. 1995, 4: 420-432. 10.1006/mpev.1995.1039.PubMedView Article
- Rüber L, Van Tassell JL, Zardoya R: Rapid speciation and ecological divergence in the American seven-spined gobies (Gobiidae, Gobiosomatini) inferred from a molecular phylogeny. Evolution. 2003, 57: 1584-1598. 10.1111/j.0014-3820.2003.tb00366.x.PubMedView Article
- Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008, 9: 286-298. 10.1093/bib/bbn013.PubMedView Article
- Vaidya G, Lohman DJ, Meier R: SequenceMatrix: concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics. 2011, 27: 171-180. 10.1111/j.1096-0031.2010.00329.x.View Article
- Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML web servers. Syst Biol. 2008, 57: 758-771. 10.1080/10635150802429642.PubMedView Article
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedView Article
- Bloom DD, Unmack PJ, Gosztonyi AE, Piller KR, Lovejoy NR: It’s a family matter: molecular phylogenetics of Atheriniformes and the polyphyly of the surf silversides (Family: Notocheiridae). Mol Phylogenet Evol. 2012, 62: 1025-1030. 10.1016/j.ympev.2011.12.006.PubMedView Article
- Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.PubMedView Article
- 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.PubMedPubMed CentralView Article
- Miller MA, Pfeiffer W, Schwartz T: Creating the CIPRES Science Gateway for inference of large phylogenetic trees. Proceedings of the Gateway Computing Environments Workshop (GCE). 2010, New Orleans, 1-8.View Article
- Setiamarga DHE, Miya M, Yamanoue Y, Azuma Y, Inoue JG, Ishiguro NB, Mabuchi K, Nishida M: Divergence time of the two regional medaka populations in Japan as a new time scale for comparative genomics of vertebrates. Biol Lett. 2009, 5: 812-816. 10.1098/rsbl.2009.0419.PubMedPubMed CentralView Article
- Hill KC, Hall R: Mesozoic-Cenozoic evolution of Australia’s New Guinea margin in a west Pacific context. Geol Soc Aust Spec Publ Soc Am Spec Pap. 2003, 22/372: 265-289.
- Rambaut A, Drummond AJ: Tracer v. 1.5. 2007, Available from http://tree.bio.ed.ac.uk/software/tracer,
- Saeed B, Ivantsoff W, Crowley LELM: Systematic relationships of atheriniform fishes within division 1 of the series Atherinomorpha (Actinopterygii) with relevant historical perspectives. Vopr Ichtiologii. 1994, 34: 1-32.
- Near TJ, Dornburg A, Eytan RI, Keck BP, Smith WL, Kuhn KL, Moore JA, Price SA, Burbrink FT, Friedman M, Wainwright PC: Phylogeny and tempo of diversification in the superradiation of spiny-rayed fishes. Proc Natl Acad Sci USA. 2013, 110: 12738-12743. 10.1073/pnas.1304661110.PubMedPubMed CentralView Article
- Vaillant JJ, Haffner GD, Cristecu ME: The ancient lakes of Indonesia: towards integrated research on speciation. Integr Comp Biol. 2011, 51: 634-643. 10.1093/icb/icr101.PubMedView Article
- Wilson MEJ, Moss SJ: Cenozoic palaeogeographic evolution of Sulawesi and Borneo. Palaeogeogr Palaeoclimatol Palaeoecol. 1999, 145: 303-337. 10.1016/S0031-0182(98)00127-8.View Article
- Saeed B, Ivantsoff W: Kalyptatherina, the first telmatherinid genus known outside of Sulawesi. Ichthyol Explor Freshw. 1991, 2: 227-238.
- Parenti LR: Phylogenetic systematics and biogeography of phallostethid fishes (Atherinomorpha, Phallostethidae) of Northwestern Borneo, with description of a new species. Copeia. 1996, 1996: 703-712. 10.2307/1447535.View Article
- Parenti LR, Louie KD: Neostethus djajaorum, new species, from Sulawesi, Indonesia, the first phallostethid fish (Teleostei: Atherinomorpha) known from east of Wallace’s line. Raffles Bull Zool. 1998, 43: 139-150.
- Moss SJ, Wilson MEJ: Biogeographic implications of the Tertiary palaeogeographic evolution of Sulawesi and Borneo. Biogeography and Geological Evolution of SE Asia. Edited by: Hall R, Holloway JD. 1998, Leiden: Backhuys Publishers, 133-163.
- Hall R: Southeast Asia’s changing palaeogeography. Blumea. 2009, 54: 148-161. 10.3767/000651909X475941.View Article
- Spakman W, Hall R: Surface deformation and slab–mantle interaction during Banda arc subduction rollback. Nat Geosci. 2010, 3: 562-566. 10.1038/ngeo917.View Article
- Lohman DJ, de Bruyn M, Page T, von Rintelen K, Hall R, Ng PKL, Shih HT, Carvalho GR, von Rintelen T: Biogeography of the Indo-Australian archipelago. Annu Rev Ecol Evol Syst. 2011, 42: 205-226. 10.1146/annurev-ecolsys-102710-145001.View Article
- Stelbrink B, Albrecht C, Hall R, von Rintelen T: The biogeography of Sulawesi revisited: is there evidence for a vicariant origin of taxa on Wallace’s “anomalous island”?. Evolution. 2012, 66: 2252-2271. 10.1111/j.1558-5646.2012.01588.x.PubMedView Article
- Hall R: Late Jurassic-Cenozoic reconstructions of the Indonesian region and the Indian Ocean. Tectonophysics. 2012, 570–571: 1-41.View Article
- Hall R: Sundaland and Wallacea: geology, plate tectonics and palaeogeography. Biotic Evolution and Environmental Change in Southeast Asia. Edited by: Gower DJ, Johnson KG, Richardson JE, Rosen BR, Rüber L, Williams ST. 2012, Cambridge: Cambridge University Press, 32-78.View Article
- Abell R, Thieme ML, Revenga C, Bryer M, Kottelat M, Bogutskaya N, Coad B, Mandrak N, Contreras Balderas S, Bussing W, Stiassny MLJ, Skelton P, Allen GR, Unmack P, Naseka A, Ng R, Sindorf N, Robertson J, Armijo E, Higgins JV, Heibel TJ, Wikramanayake E, Olson D, Lopéz HL, Reis RE, Lundberg JG, Sabaj Pérez MH, Petry P: Freshwater ecoregions of the world: a new map of biogeographic units for freshwater biodiversity conservation. Bioscience. 2008, 58: 403-414. 10.1641/B580507.View Article
- McGuigan K: Phylogenetic relationships and historical biogeography of melanotaeniid fishes in Australia and New Guinea. Mar Freshw Res. 2000, 51: 713-723. 10.1071/MF99159.View Article
- Walter RP, Hogan JD, Haffner GD, Heath DD: Genetic structure and connectivity among lake populations of threatened Paratherina sailfin silversides from Sulawesi, Indonesia. Conserv Genet. 2011, 12: 1387-1393. 10.1007/s10592-011-0238-5.View Article
- Kottelat M: Sailfin silversides (Pisces: Telmatherinidae) of Lakes Towuti, Mahalona and Wawontoa (Sulawesi, Indonesia) with descriptions of two new genera and two new species. Ichthyol Explor Freshw. 1990, 1: 35-54.
- von Rintelen T, Wilson AB, Meyer A, Glaubrecht M: Escalation and trophic specialization drive adaptive radiation of viviparous freshwater gastropods in the ancient lakes on Sulawesi, Indonesia. Proc R Soc London B. 2004, 271: 2541-2549. 10.1098/rspb.2004.2842.View Article
- Stelbrink B, Stöger I, Hadiaty RK, Schliewen UK, Herder F: Age estimates for an adaptive lake fish radiation, its mitochondrial introgression, and an unexpected sister group: Sailfin silversides of the Malili Lakes system in Sulawesi. Dryad Digital Repository. 2014, doi:10.5061/dryad.8dh7r at http://datadryad.org,
- Miya M, Takeshima H, Endo H, Ishiguro NB, Inoue JG, Mukai T, Satoh TP, Yamaguchi M, Kawaguchi A, Mabuchi K, Shirai SM, Nishida M: Major patterns of higher teleostean phylogenies: a new perspective based on 100 complete mitochondrial DNA sequences. Mol Phylogenet Evol. 2003, 26: 121-138. 10.1016/S1055-7903(02)00332-9.PubMedView Article
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.