Out of Tanganyika: Genesis, explosive speciation, key-innovations and phylogeography of the haplochromine cichlid fishes
© Salzburger et al; licensee BioMed Central Ltd. 2005
Received: 19 November 2004
Accepted: 21 February 2005
Published: 21 February 2005
The adaptive radiations of cichlid fishes in East Africa are well known for their spectacular diversity and their astonishingly fast rates of speciation. About 80% of all 2,500 cichlid species in East Africa, and virtually all cichlid species from Lakes Victoria (~500 species) and Malawi (~1,000 species) are haplochromines. Here, we present the most extensive phylogenetic and phylogeographic analysis so far that includes about 100 species and is based on about 2,000 bp of the mitochondrial DNA.
Our analyses revealed that all haplochromine lineages are ultimately derived from Lake Tanganyika endemics. We find that the three most ancestral lineages of the haplochromines sensu lato are relatively species poor, albeit widely distributed in Africa, whereas a fourth newly defined lineage – the 'modern haplochromines' – contains an unparalleled diversity that makes up more than 7% of the worlds' ~25,000 teleost species. The modern haplochromines' ancestor, most likely a riverine generalist, repeatedly gave rise to similar ecomorphs now found in several of the species flocks. Also, the Tanganyikan Tropheini are derived from that riverine ancestor suggesting that they successfully re-colonized Lake Tanganyika and speciated in parallel to an already established cichlid adaptive radiation. In contrast to most other known examples of adaptive radiations, these generalist ancestors were derived from highly diverse and specialized endemics from Lake Tanganyika. A reconstruction of life-history traits revealed that in an ancestral lineage leading to the modern haplochromines the characteristic egg-spots on anal fins of male individuals evolved.
We conclude that Lake Tanganyika is the geographic and genetic cradle of all haplochromine lineages. In the ancestors of the replicate adaptive radiations of the 'modern haplochromines', behavioral (maternal mouthbrooding), morphological (egg-spots) and sexually selected (color polymorphism) key-innovations arose. These might be – together with the ecological opportunity that the habitat diversity of the large lakes provides – responsible for their evolutionary success and their propensity for explosive speciation.
"At some stage in the past the waterways of Africa were, from the fishes' point of view, accessibly interconnected." P. H. Greenwood (1983)
With estimated numbers of about 1,000 and 500 species respectively the assemblages of cichlid fishes of East African lakes Malawi (LM) and Victoria (LV) are, by far, the most species-rich species flocks [1–6]. In Lake Tanganyika (LT), Africa's oldest lake, "only" about 200 to 250 cichlid species occur, but they are phenotypically, behaviorally and genetically more diverse [7–9]. Remarkably, almost 100% of the species of all of these flocks are endemics and the three East African Great Lakes do not have a single cichlid species in common . While in LT at least twelve eco-morphologically distinct tribes can be clearly distinguished phylogenetically [7, 9, 11], the species flocks of LM and LV are entirely comprised of cichlids assigned to only one of these tribes, the Haplochromini [1, 4, 10, 12]. Additionally, several less species-rich flocks of haplochromines are found in smaller lakes in East Africa , some of which have been combined with the LV radiation into a "superflock" of closely related species [4, 12, 13].
That the haplochromines have a particularly strong propensity for explosive speciation among cichlids is best illustrated by their unparalleled species-richness and diversity. With the exception of two LT lineages (Lamprologini: ~70–100 species; Ectodini: ~25–30 species), discrete adaptive radiations of non-haplochromines are comprised by about a dozen species at most. Thus, which evolutionary novelties might be causally related to the explosive speciation that distinguish the haplochromines from all other cichlid lineages remains a crucial question for the understanding of the explosive patterns of cichlid evolution in East Africa. Haplochromines furthermore represent prime examples for parallel evolution, and it is particularly the lacustrine haplochromine species flocks that independently evolved morphologies and color-patterns that are convergent between species of different lakes [2, 17, 18].
In order to gain a deeper understanding of the adaptive radiations of cichlids in general and the formation of East African cichlid species flocks in particular, several crucial questions still remain to be answered: (i) What evolutionary lineages from which geographic ranges make up the diversity of haplochromines? (ii) Which were the founding lineages of the lacustrine adaptive radiations and were they riverine generalists? (iii) Can particular behavioral and/or morphological key-innovations be identified that might be causally related with the diversification of haplochromines? With the aim of addressing these questions, we conducted the most extensive phylogenetic and phylogeographic study of haplochromine cichlids so far, analyzing a portion of up to 2,000 bp of the mitochondrial genome of about 100 species. We included representatives of relevant cichlid tribes from LT [7, 9, 15] as well as members of all major riverine and lacustrine haplochromine lineages and all but one haplochromine genus (following ).
In the maximum likelihood tree of the dataset including 100 taxa (Fig. 3) the LT Eretmodini were placed as sister group to the LT Lamprologini, followed by a clade comprised by the two representatives of the Limnochromini (Limnochromis auritus, Triglachromis otostigma) plus Cyphotilapia frontosa, and all remaining taxa. Among these, a clade that includes four LT tribes (Cyprichromini, Ectodini, Limnochromini, Perissodini) plus the Orthochromis species from the Malagarasi drainage, which had so far been considered to belong to the Haplochromini, was recovered as sister group to the remaining haplochromine representatives. The latter, the haplochromines sensu lato, clustered into four distinct groups, with a Congolese/South-African lineage (CSA) sister to the Pseudocrenilabrus-, the Astatoreochromis-lineage, and the 'modern haplochromines' (i.e. a clade comprised by the LT Tropheini sister to the LM representatives, several East-African riverine lineages and the members of the LV region superflock). Bayesian inference revealed the same branching order for the different lineages. In the strict consensus topology of the 78,617 most parsimonious trees (unweighted tree length: 4400; tree not shown), the Eretmodini were again resolved as sister group to the Lamprologini, the Limnochromini plus C. frontosa, four LT tribes plus the Malagarasi Orthochromis and the haplochromines sensu lato. Here, the Astatoreochromis-lineage was resolved as most ancestral lineage, and as sister group to a clade comprised by the Pseudocrenilabrus- and the CSA lineage, and the modern haplochromines. Also, the neighbor-joining analysis recovered a tree (not shown) with the Eretmodini as sister-group to the Lamprologini, the Limnochromini plus C. frontosa, the four LT tribes including the Malagarasi River Orthochromis species, and the haplochromines sensu lato. In the neighbor-joining tree, the Pseudocrenilabrus-lineage occupied the most ancestral branch in the haplochromines sensu lato, followed by the Congolese/South-African clade, the Astatoreochromis-lineage and the modern haplochromines. A Shimodaira-Hasegawa test  revealed that there are no significant differences between the topologies obtained with the different algorithms (P < 0.05). Similarly, in the four-cluster likelihood mapping analysis  none of the three possible alternative branching orders among the four main lineages of the haplochromines sensu lato received support greater than 50%.
Phylogeny and evolutionary origin of haplochromines
Mitochondrial DNA (mtDNA) sequence data have a long and successful history in the study of East African cichlid evolution (see e.g., [1, 4, 9, 17, 24, 25]). Limitations with mtDNA have only been encountered when focusing on the phylogeny among extremely closely related species due to the possibility of the persistence of ancestral polymorphism (see e.g., ), or because of hybridization events (see e.g., [27–29]). Nevertheless, mtDNA sequences proved to be particularly suitable for the reconstruction of the East African cichlid phylogeny at the tribal level and for tribal assignments [9, 15]; for phylogenetic reconstruction within older tribes [30, 31]; and for phylogeographic analyses [4, 32]. Also, because of the extremely fast rate of lineage formation in cichlids, nuclear and even some mitochondrial genes  are too slowly evolving to contain phylogenetic information (reviewed in ).
The different phylogenetic algorithms, with which we analyzed our data, revealed largely congruent results. In all analyses, and in agreement with previous results , we found that the Eretmodini are placed as sister group to the substrate spawning Lamprologini – with an estimated number of up to 100 species the most species-rich tribe in LT – plus several LT tribes and all haplochromine representatives (see Fig. 3 for the maximum likelihood tree). The molecular phylogenies thus corroborate that all haplochromines are ultimately derived from LT cichlids and that their ancestor(s) are likely to have left LT secondarily.
We consider all species that belong to the monophyletic group descending from this ancestor (asterisk in Fig. 3), to being haplochromines sensu lato. These are further divided into four distinct groups, a Congolese/Southern African lineage (CSA), the genera Pseudocrenilabrus (Ps.) and Astatoreochromis (As.), and the modern haplochromines (MH). While the respective monophyly of these four lineages was supported by high bootstrap values and Bayesian posterior probabilities, our analyses could not unambiguously resolve the exact relationships between these four lineages [maximum likelihood and Bayesian inference (see Fig. 3): (CSA, (Ps., (As., (MH)))); maximum parsimony: (As., ((Ps., CSA), (MH))); neighbor joining: (Ps., (CSA, (As., (MH))))]. An evaluation of these alternative hypotheses by means of a Shimodaira-Hasegawa test  and a four-cluster likelihood mapping analysis  indicated that none of these alternative branching orders receives significantly more support than the others. This suggests that the four lineages of the haplochromines sensu lato evolved almost contemporaneously from a common ancestor. This is further supported by the observation of relatively short branches interrelating these four lineages and the generally low bootstrap values and Bayesian posterior probabilities supporting the corresponding relationships proposed by the different algorithms.
The CSA lineage is composed of several widespread and moderately species-rich groups of the Congo drainage and Southern Africa and consists of two main clades: A clade with a species from the Congo drainage (Ctenochromis oligacanthus) ancestral to Southern African genera (Serranochromis, Sargochromis, Pharyngochromis) was resolved as sister group to a clade comprised by solely Congo drainage taxa (Orthochromis, Haplochromis sp. nov., Cyclopharynx and Thoracochromis). Within the Pseudocrenilabrus- and Astatoreochromis-clades, branch lengths were relatively short and a more detailed phylogeographic sampling would be necessary to resolve the relationships between the different geographic lineages.
The modern haplochromines consist of species flocks of an unparalleled diversity. They include the endemic LT tribe Tropheini (~25 species) sister-group to a clade comprised by the entire species flock of LM (~1,000 species) and several East African riverine and lacustrine lineages (~200 species) plus the LV region superflock (~600 species). With approximately 1,800 species this – here phylogenetically defined – monophyletic lineage makes up about 7% of all known teleost fish.
Astatotilapia, Thorachochromis, and Orthochromis are polyphyletic genera
Our phylogenies show that several genera are in fact polyphyletic, and major taxonomic revisions will be required in the future to take our phylogenetic results into consideration. For example, and in agreement with previous studies (see e.g., [4, 14, 33]), Astatotilapia emerges as polyphyletic genus, with representatives assigned to both the East African riverine clade in the modern haplochromines and the LV region superflock.
The genus Thoracochromis, represented in our analysis by T. brauschi, has also been shown to be polyphyletic before , with T. brauschi from the Congo drainage as a more ancestral lineage, and T. petronius and T. pharyngealis from the Nile drainage with affinities to the LV region superflock (note that the Nile River Thoracochromis of  are consequently listed as Haplochromis in ). The placement of T. brauschi as sistergroup to T. petronius and T. pharyngealis plus the remaining representatives of the LV region superflock in the AFLP based phylogeny of  seems to contradict our mtDNA based results in which T. brauschi was identified as member of the CSA lineage. However, the reported branching order did not receive considerable bootstrap support in . Also, the choice of Astatoreochromis alluaudi (mislabeled as Astatotilapia alluaudi in ) as single outgroup species seems problematic, as our present analyses (see above) and former results  indicate that A. alluaudi and not T. brauschi is more closely related to the modern haplochromines (and, thus, also to the LV region superflock). Further analyses including nuclear DNA sequence data and more taxa assigned to the genus will be necessary to address this problem.
Based on the phylogeny it is apparent that the Orthochromis lineage, which is confined to the Malagarasi River system and two isolated rivers, East of LT , is not part of the radiation of haplochromine cichlids. This is further supported by the breeding behavior of these fish: While the Malagarasi River Orthochromis are biparental caregivers , the haplochromines sensu lato are all maternal mouthbrooders. Thus, the Orthochromis-species from the Malagarasi area, a group of exclusively riverine fish, should be placed into its own tribe. The name Orthochromis is, however, also used for riverine species from the Congo drainage. In our analyses, Orthochromis polyacanthus and O. stormsi from the Congo River system fall – according to their distribution – into the CSA-clade leaving also the genus Orthochromis polyphyletic (see also ). We suggest to using Orthochromis for the CSA lineage representatives, since O. polyacanthus was the first species of the genus to be described , and Schwetzochromini (as tribe name) and Schwetzochromis (as genus name) for the Malagarasi area species, since this name was repeatedly used for some species of that complex (see e.g., [9, 34, 35]).
Phylogeography and phylochronology
We note that the application of a molecular clock for estimating divergence times has the potential of not being without problems for several reasons (see e.g., [37, 38]). However, a molecular-clock-based time estimate does surely provide an approximate framework for phylogeographic inferences. Our phylogeographic scenario (Fig. 4a), which is derived from the maximum likelihood phylogenetic and molecular clock analyses, suggests that several lineages independently left LT to colonize surrounding river systems and consequently other lakes in East Africa. The molecular clock calibration based on the chronogram generated with r8s  (Fig. 4b,c) yielded about 2.4 MYA (1.22 – 4.02 MYA) for the split of the CSA lineage from the common ancestor of the haplochromines sensu lato, and about 2 MYA (1.15 – 3.89 MYA) for the first branching events within the CSA clade. The spread of Congo drainage taxa into southern river systems occurred at a later stage, most likely at the relatively shallow watershed between upper branches of the Congo River and the Zambezi River – a scenario that is also supported by the placement of Serranochromis sp. (from Lake Mweru in the upper Congo) as sister group to the Zambezi/Southern African genera Sargochromis and Pharyngochromis in our phylogenies. However, further sampling in that area would be necessary to reconstruct the southward spread of the CSA lineage.
At essentially the same time as the CSA lineage, the ancestors of the Pseudocrenilabrus- and the Astatoreochromis- lineage diverged from their common ancestor. Despite their large distributional ranges – they also colonized the LV (both lineages) and LM (Pseudocrenilabrus) region – the genera Astatoreochromis and Pseudocrenilabrus never underwent considerable speciation. The three described Astatoreochromis species occur in the LV region including Lakes Edward and George (A. alluaudi), in the rivers Rusizi and Lukuga (A. straeleni), and in the Malagarasi River (A. vanderhorsti). The three species of Pseudocrenilabrus occur from the Nile system to the LV region (P. multicolor), in Eastern- and Southern Africa including LM (P. philander), and in the central Congo basin (P. nicholsi). All analyzed representatives are relatively closely related suggesting a recent spread of these lineages in East Africa. However, we did not include P. nicholsi, which is morphologically different from P. multicolor and P. philander and would – if it really belonged to the Pseudocrenilabrus-lineage – represent the only haplochromine in the Congo drainage that is not a member of the CSA lineage.
The most recent common ancestor of the modern haplochromines was dated to have existed about 1.8 MYA (0.66 – 3.78 MYA) (Fig. 4). This ancestral lineage forms the crucial phylogenetic and biogeographic link between the species flocks of all three East African Great Lakes, and its discovery documents the existence of much earlier hypothesized fish-accessible waterways between these waterbodies [40, 41]. Apparently, the Malagarasi River (and possibly the Rusizi) played a major role for the dispersal of these fishes, since many modern haplochromine lineages occur in these drainages and in lakes South-Eastern and North of LT exclusively, which argues against the view that LM haplochromines originated from Zambezi River stocks . Whether or not Lake Rukwa has ever acted as link between the faunas of LT and LM  cannot be answered by our data. Lake Rukwa seems to have overflowed at its maximum levels into LT several times. However, Lake Rukwa has also become very shallow in recent geological times and it might have dried up completely  eradicating its original fauna. At present, Lake Rukwa harbors haplochromines that belong to the East-African riverine clade in Figs. 2, 3[4, 25].
Our analyses also recovered another closely related lineage to the LV region superflock in the East-African riverine clade, in addition to Haplochromis gracilior form Lake Kivu . This lineage includes H. paludinosus that occurs in the Malagarasi (which was already suggested by ), as well as undescribed species from Tanzania and Lake Edward (Figs. 2, 3). It is, however, unclear by which waterway haplochromine cichlids once colonized Lake Kivu. The flow of the Rusizi, presently from Lake Kivu into LT with the Panzi falls as strong barrier for fish migration, might actually have been reversed before the uplift of the Virunga volcanoes north of Lake Kivu as suggested by deposits of fossil LT mollusks and fluviatile sands in the upper Rusizi valley [41, 42]. This connection could possibly explain how haplochromines of LT origin were able to colonize Lake Kivu about 1.5 million years ago (Fig. 4).
Evolutionary key-innovations of haplochromines
One of only few synapomorphies of the haplochromines sensu lato is the particular type of cranial apophysis for the upper pharyngeal bones . The distinctive organization of the pharyngeal apophysis, a second set of jaws that is functionally decoupled from the oral ones , is characteristic to all cichlids and has been interpreted as prominent feature that – because of its adaptability – contributes to the cichlids' evolutionary success [2, 3, 10, 43]. It is, however, not evident how the relatively minor morphological modification of part of that structure in the haplochromines  might function as an evolutionary key-innovation. Interestingly, however, all haplochromines sensu lato are maternal mouthbrooders with the females alone incubating the eggs in their buccal cavities [10, 12]. Mouthbrooding, which is regarded as an adaptation to predation pressure [44–46], has evolved several times independently and in diverse behavioral modes in cichlids [10, 22, 47, 48]. The characteristic maternal mouthbrooding behavior displayed by haplochromines is believed to being a derived character state [35, 46, 47]. Mouthbrooding strongly limits the number of eggs and fry that can be raised and might have led to generally much smaller population sizes, which has, for example, population genetic implications on fixation of alleles, and might result in smaller effective population sizes. Furthermore, mouthbrooding species may be considered to being promising colonizers of new habitats, since only a single mouthbrooding female is necessary for the founding of a new population.
An eminent feature of several female mouthbrooding cichlid genera is the occurrence of egg-spots on the anal fins of males. In some species also females show such ovoid markings, but these are smaller and much less conspicuous than in males. Also, some species of the modern haplochromines, e.g., some deep-water lineages of LM, have lost their egg-spots secondarily. In mimicking real eggs to attract females, these egg-spots function as natural releasers [22, 23], or intra-specific sexual advertisement , apparently serving to ensure a greater fertilization success of the eggs by bringing about greater proximity of the female's mouth to the male's genital opening. Based on the molecular phylogeny, we could trace the origin of the characteristic egg-spots (ocelli) [10, 23] to the common ancestor of the Astatoreochromis-lineage and the modern haplochromines. There are other cichlid species in which males display yellow or red marks on their pelvic, dorsal or anal fins, but only in these lineages true egg-spots on the males' anal fins with a yellow, orange or red center and a colorless/transparent outer ring [10, 22] are found. Interestingly, the branch leading to the Astatoreochromis-lineage and the modern haplochromines is the one with a pronounced potential for an increased rate of speciation (see Fig. 4c). Based on the character state reconstructions (Fig. 5) it seems likely that this ancestor was riverine. Thus, it may be concluded that the egg-spots first evolved in a haplochromine cichlid that inhabited a turbid riverine environment, where these conspicuous markings would seem to be particularly effective and necessary for intra-specific communication.
Another innovation that further distinguishes the exceptionally species-rich modern haplochromines from all other cichlids is the occurrence of numerous color morphs, often accompanied by sexual color dimorphism. Inter- and intra-specific polychromatism combined with maternal mouthbrooding involving egg-spots as releasers can be hypothesized to being permissive for sexual selection through female choice and, hence, the haplochromines' propensity for species formation, as sexual selection is probably a major causal factor in the origin of isolating mechanisms and the maintenance of reproductive isolation [18, 49–53]. These distinctive features of the modern haplochromines, that have arisen just in their ancestor, in combination with the numerous ecological niches that are provided by the large East African lakes might thus have induced a considerable increase of the haplochromines' evolutionary potential. The importance of large waterbodies for the evolution of the modern haplochromines is reflected by the fact that these cichlids only radiated in lakes (and species number rather correlates to the size, but not to the age, of a lake), whereas the riverine lineages are all species-poor albeit often widespread (Figs. 1, 3).
Replicate adaptive radiations of the 'modern haplochromines'
A common feature of many adaptive radiations is that their founders are believed to have had a more generalist's lifestyle, while adaptive radiations themselves are defined by being composed of highly specialized species with narrower niche widths [54, 55]. Theory predicts that generalists more likely have better dispersal abilities and are expected to be able to adapt readily to novel environmental settings . A generalist ancestor scenario fits well with the diversification of haplochromine cichlids. Generalist riverine species of the genera Astatoreochromis, Astatotilapia, Pseudocrenilabrus, and Haplochromis (e.g., bloyeti), are ancestral to the adaptive radiations of the Tropheini of LT, and/or the radiations of LM and the LV region superflock. These genera are widely distributed and not confined to Eastern Africa, and they are the only ones that could inhabit the waterways that – over geological time spans – connected the lakes of Eastern Africa.
The phylogeny presented here (Fig. 3) reveals that modern haplochromines gave rise to several major adaptive radiations; the most prominent ones are those of LM and LV. Interestingly, it uncovers that also the radiation of the Tropheini from LT  must now be considered as an additional radiation of the modern haplochromines, corroborating the much older perception that LT accommodates several independent species flocks . It is further suggested by mapping the fishes' lifestyle onto our molecular phylogeny that the highly specialized Tropheini are descendents of a river-living species. This implies that the ancestor of the Tropheini successfully re-entered the lake habitat and evolved into the presently dominant group in the rocky littoral zone of LT. Thus, this lineage of modern haplochromines managed to occupy "empty niches" in an apparently "full" ecosystem, as all remaining tribes, which now account for about 200 species, had already been established when the ancestor of the Tropheini secondarily entered LT (Figs. 2, 3). The observation that these fish underwent an independent adaptive radiation in LT underlines the haplochromines' propensity for speciation.
In an apparent contrast to most other known examples of adaptive radiations  is the finding that the generalist ancestors of the haplochromine species flocks were derived from already highly diverse and specialized LT endemics (Fig. 3). Therefore, specialization may not be an "evolutionary one-way street", but rather some lineages have reversed their level of specialization, i.e., generalists arose from highly specialized lineages, yet, apparently retained their high propensity for speciation and level of evolvability (see ). The faunal revolution of LT's radiation of cichlids was thus not confined to the lake habitat itself (see also [8, 9, 58]), but it effectively involved large parts of Africa via the intermediate step of repeatedly evolving generalist riverine lineages – in much the same way as the adaptive radiation in LV produced haplochromine species that secondarily colonized surrounding rivers .
Our phylogeny of haplochromines provides strong support for replicate adaptive radiations in East African cichlids. The concept of replicate radiations, in which the same sequence of adaptations to ecological niches evolved repeatedly in lineages that inhabit similar environments, has been developed based on sympatric species pairs of fishes in postglacial lakes and on the Anolis lizard ecomorphs on different islands of the Greater Antilles [59–61]. Our inclusive phylogenetic and phylogeographic study shows that similar ecological types of cichlids in the different East African lakes evolved independently (see also [2, 17]), yet it also shows that the convergent ecotypes in the species flocks of LM, LV, Lake Kivu as well as in the Tropheini [2, 10, 12, 13, 17] arose from the same ancestral phenotype in the ancestor of the modern haplochromines. We suggest that a combination of behavioral (maternal mouthbrooding) and morphological innovations (egg-spots, color polymorphisms, pronounced sexual dichromatism) as well as ecological opportunities (after the colonization of large lakes) might have predestined this particular lineage to give rise to these replicate adaptive radiations.
It has been noted before that lineages of LT origin have left the lake secondarily (there are, for example, about five lamprologine species that are found in the Congo and Malagarasi Rivers) [9, 58]. Here, we show that the entire haplochromine diversity has its origin in LT corroborating the view that ancient lakes not only preserve biodiversity but also act as biodiversity hotspots, genetic reservoirs and cradles from which new lineages evolve [4, 8, 9]. What remains to be answered is where the LT cichlids originated and to what extent a proposed and meanwhile desiccated Pliocene lake in the Congo plains [41, 42] was the source of the ancient LT lineages, pushing back even further the onset of replicate adaptive radiations in East African cichlids.
Our phylogenetic analyses that include representatives of all major haplochromine lineages show that all haplochromines are derived from Lake Tanganyika cichlids. While the Orthochromis species of the Malagarasi area apparently do not belong to the radiation of the haplochromine cichlids and should be placed into a new tribe, we defined four new lineages within the haplochromines sensu lato: A clade combining Congolese and South-African genera (CSA-lineage), the Pseudocrenilabrus-, the Astatoreochromis-lineage, and the exceptionally species-rich modern haplochromines. The ca. 1,800 species of modern haplochromines are comprised of the entire haplochromine species flocks of Lake Malawi and the Lake Victoria region, some 200 riverine and lacustrine species, as well as the Tanganyikan Tropheini, which are likely to have evolved from a riverine ancestor and secondarily colonized Lake Tanganyika. This proposed "out of Tanganyika" scenario of haplochromine evolution is in agreement with the geological and palaeo-geological history of East Africa. Based on a character-state reconstruction from this new phylogeny, we were able to discover the evolution of several key-innovations that arose in the lineage leading to the modern haplochromines. These character reconstructions suggest that a combination of behavioral (maternal mouthbrooding) and morphological characteristics (egg-spots, color polymorphisms, pronounced sexual dichromatism) as well as ecological opportunities (after the colonization of large lakes) might have predestined this particular lineage to give rise to replicate adaptive radiations and, therefore, be causally related to the extraordinary success of these particular cichlid fish radiations.
Specimen information and DNA methods
For this study, a total of 304 specimens were analyzed. We combined all available GenBank entries from previous studies [4, 9, 15, 25, 62] with 180 newly determined DNA sequences. Table 1 [see Additional file 1] lists specimen information, geographic origin of the specimens, names of collectors, and GenBank accession numbers. When available, voucher specimens have been deposited at the Royal Museum for Central Africa, Tervuren, Belgium.
Sample preparation, polymerase chain reaction (PCR) amplification and DNA sequencing have been performed as described elsewhere  for both mitochondrial DNA segments the complete non-coding control region and the entire NADH Dehydrogenase Subunit II (ND2) gene. Forward and reverse sequences have been assembled using the computer programs Sequence Navigator (Applied Biosystems, USA) and Sequencher (GenCodes, USA).
Phylogenetic reconstruction and hypotheses testing
The complete sequences of the mitochondrial control region (895 bp) and the ND2 gene (1,047 bp) were aligned using the computer program Clustal W ; alignments have been further adjusted by eye. Up to 34 gaps had to be included in the control region alignment, which were coded as indels. Due to missing data on the 5'-end of the control region in about one quarter of the taxa, a terminal section of 34 bp has been excluded from the phylogenetic analysis, leading to an alignment of 1,908 bp for the combined dataset.
For phylogenetic reconstruction we performed maximum likelihood, maximum parsimony and neighbor-joining methods in parallel using the computer program PAUP* 4.0b10 . Two taxa belonging to tilapiine cichlids were included as outgroup based on previous phylogenetic analyses using mitochondrial DNA [9, 15, 16], nuclear DNA  and SINE insertion patterns [66, 67] demonstrating that the Tilapia-lineages are ancestral to all remaining LT tribes (but excluding Tylochromis, which is ancestral to the Tilapiini). We did not include representatives of three LT tribes – the Bathybatini, Trematocarini and Tylochromini -, as these ancient lineages were shown to have existed before the primary radiation of mouthbrooders in LT [9, 15, 67]. After an initial neighbor-joining analysis including the control region sequences of all 304 specimens, we reduced the dataset to 100 taxa based on the obtained topology. This reduction was necessary to allow computational feasibility for maximum likelihood and maximum parsimony analyses. The optimal model of molecular evolution for the maximum likelihood analysis was determined in a likelihood-ratio test running the computer program Modeltest v3.06 . For the heuristic maximum likelihood search of the combined dataset we used the general-time-reversible model of molecular evolution, with a gamma shape parameter of 0.7937 and a proportion of invariable sites of 0.3426. Due to the many closely related taxa in the dataset, maximum parsimony analyses were completed for 1010 rearrangements. For neighbor joining, we used the HKY model and conducted a bootstrap analysis with 5,000 replicates. Bootstrap analyses for maximum likelihood were performed with 100 replicates and for maximum parsimony with 1,000 replicates. We also applied Bayesian inference of phylogeny with Mr. Bayes 3.0b4  running four Metropolis Coupled Monte-Carlo-Markov-Chains in parallel for 250.000 generations, using the general-time-reversible model with gamma correction (six types of substitutions), and excluding 5 % of the trees as burn-in. The obtained topologies from the different phylogenetic algorithms were evaluated by means of a nonparametric Shimodaira-Hasegawa test  under a resampling-estimated log-likelihood as implemented in PAUP* . To estimate the support for distinct internal branches critical for our interpretations, we performed a four-cluster likelihood mapping analysis with the program PUZZLE 5.0 [20, 70] in which we grouped the taxa of the haplochromines sensu lato into four groups according to the four lineages that were recovered from the phylogenetic analyses (CSA-, Pseudocrenilabrus-, Astatoreochromis-lineage, and modern haplochromines).
To test whether or not shifts in the probabilities of speciation occurred on certain branches of the obtained phylogeny, we quantitatively tested the fit of this tree to a Markovian null model in which the probability of speciation is equally distributed along branches. Therefore, we calculated the cumulative probability p c  for relevant branches based on the different trees. We then mapped lifestyles (riverine versus lacustrine), breeding behavior, and the occurrence of true egg-spots on the maximum likelihood tree using MacClade 4.0 (Sinauer, Sunderland, MA) for maximum parsimony character state reconstructions and Mesquite  for maximum likelihood reconstructions. To tentatively date the major cladogenetic events in the haplochromines, we constructed a chronogram based on a maximum likelihood tree with constraint molecular clock (see above for search parameters). Therefore, we used the computer program r8s  applying the local molecular clock method and an optimization via the truncated newton method . Confidence intervals were assessed by means of a bootstrap approach. We simulated 25 bootstrap matrices with Mesquite and, for each matrix, constructed a maximum likelihood tree (general-time-reversible model; model parameters, gamma shape correction and proportion of invariable sites being estimated from each matrix; rearrangements limited to 1,000). The resulting trees were then analyzed with r8s as described above, and in addition with the minimum and maximum values of geological datings separately, in order to define upper and lower bounds. As calibration points, we used an age estimate for the LM species flock of about one million years [74, 75], the maximum age of 200.000 years for the LV region superflock [1, 4, 25] as well as the time window for the Lukuga connection between LT and the Congo system (between the minima at 3.5 MYA and 1.1 MYA) [76, 77]. An escape of LT lineages into river systems draining into the lake was possible at any time, whereas the only connection to the Congo system, the Lukuga, was available at periods of high lake level stands only. The Lukuga channel, the only outlet of LT, was dry when first seen in 1874 by Cameron, but four years later the lake overflowed. Since then, the Lukuga was repeatedly flooded. Intermittent connections with the Congo River system in the late Pliocene/early Pleistocene were suggested on the basis of large lake level fluctuations in LT during periods with increased precipitation [41, 76–79]. The resulting molecular clock rate fits well with previously used rates in East African cichlids [4, 75].
List of abbreviations used
million years ago
NADH Dehydrogenase Subunit II
This paper is dedicated to the memory of the eminent African ichthyologist L. De Vos (1957–2003). We thank R. Abila, C. Katongo, I. Kornfield, R. Paul, T. Reuter, U. Schliewen, E. Schraml, L. Seegers, O. Seehausen, J. Snoeks, C. Sturmbauer, and L. De Vos for providing specimens and/or sharing unpublished sequences; K. Lezzar and J.J. Tiercelin for discussion; A. McArthur and M. J. Sanderson for technical assistance; M. Barluenga, P. Bunje, G. Fryer and an anonymous reviewer for valuable comments on the manuscript. This work was funded by the European Union (Marie Curie fellowship) and the Landesstiftung Baden-Württemberg to W. S.; a grant from The Belgian Federal Office for Scientific, Technical and Cultural Affairs and support through the Ecological Genetics Research Network of the Fund for Scientific Research-Flanders to E. V.; and grants of the Deutsche Forschungsgemeinschaft to A. M.
- Meyer A, Kocher TD, Basasibwaki P, Wilson AC: Monophyletic origin of Lake Victoria cichlid fishes suggested by mitochondrial DNA sequences. Nature. 1990, 347: 550-3. 10.1038/347550a0.View ArticlePubMed
- Stiassny MLJ, Meyer A: Cichlids of the Rift Lakes. Scientific American. 1999, 280: 64-69.View Article
- Kornfield I, Smith PF: African Cichlid Fishes: Model systems for evolutionary biology. Annu Rev Ecol Syst. 2000, 31: 163-196. 10.1146/annurev.ecolsys.31.1.163.View Article
- Verheyen E, Salzburger W, Snoeks J, Meyer A: Origin of the superflock of cichlid fishes from Lake Victoria, East Africa. Science. 2003, 300: 325-9. 10.1126/science.1080699.View ArticlePubMed
- Kocher TD: Adaptive evolution and explosive speciation: the cichlid fish model. Nature Reviews Genetics. 2004, 5: 288-298. 10.1038/nrg1316.View ArticlePubMed
- Salzburger W, Meyer A: The species flocks of East African cichlid fishes: recent advances in molecular phylogenetics and population genetics. Naturwissenschaften. 2004, 91: 277-90. 10.1007/s00114-004-0528-6.PubMed
- Poll M: Classification des Cichlidae du lac Tanganika, Tribus, genres et espéces,. T. XLV – Fascicule. 1986, Brussels: Académie Royale de Belgique, 2
- Nishida M: Lake Tanganyika as an evolutionary reservoir of old lineages of East African cichlid fishes: Inferences from allozyme data. Experentia. 1991, 47: 974-979.View Article
- Salzburger W, Meyer A, Baric S, Verheyen E, Sturmbauer C: Phylogeny of the Lake Tanganyika cichlid species flock and its relationship to the Central and East African haplochromine cichlid fish faunas. Syst Biol. 2002, 51: 113-35. 10.1080/106351502753475907.View ArticlePubMed
- Fryer G, Iles TD: The cichlid fishes of the Great Lakes of Africa: Their biology and Evolution. 1972, Edinburgh: Oliver & Boyd
- Takahashi T: Systematics of Tanganyikan cichlid fishes (Teleostei: Perciformes). Ichthyological Research. 2003, 50: 367-382. 10.1007/s10228-003-0181-7.View Article
- Greenwood PH: Towards a phyletic classification of the 'genus' Haplochromis (Pisces, Cichlidae) and related taxa. Bull Br Mus nat Hist (Zool). 1979, 35: 265-322.
- Greenwood PH: Towards a phyletic classification of the 'genus' Haplochromis (Pisces, Cichlidae) and related taxa. Part 2. Bull Br Mus nat Hist (Zool). 1980, 39: 1-99.
- Seehausen O, Koetsier E, Schneider MV, Chapman LJ, Chapman CA, Knight ME, Turner GF, van Alphen JJ, Bills R: Nuclear markers reveal unexpected genetic variation and a Congolese-Nilotic origin of the Lake Victoria cichlid species flock. Proc R Soc Lond B Biol Sci. 2003, 270: 129-37. 10.1098/rspb.2002.2153.View Article
- 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-32. 10.1006/mpev.1995.1039.View ArticlePubMed
- Klett V, Meyer A: What, if anything, is a Tilapia? Mitochondrial ND2 phylogeny of tilapiines and the evolution of parental care systems in the African cichlid fishes. Mol Biol Evol. 2002, 19: 865-83.View ArticlePubMed
- Kocher TD, Conroy JA, McKaye KR, Stauffer JR: Similar morphologies of cichlid fish in lakes Tanganyika and Malawi are due to convergence. Mol Phylogenet Evol. 1993, 2: 158-165. 10.1006/mpev.1993.1016.View ArticlePubMed
- Meyer A: Phylogenetic relationships and evolutionary processes in East African cichlids. Trends in Ecology and Evolution. 1993, 8: 279-284. 10.1016/0169-5347(93)90255-N.View ArticlePubMed
- Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol. 1999, 16: 1114-1116.View Article
- Strimmer K, von Haeseler A: Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proc Natl Acad Sci U S A. 1997, 94: 6815-9. 10.1073/pnas.94.13.6815.PubMed CentralView ArticlePubMed
- Slowinski JB, Guyer C: Testing the stochasticity of patterns of organismal diversity: An improved null model. American Naturalist. 1989, 134: 907-921. 10.1086/285021.View Article
- Wickler W: Zur Stammesgeschichte funktionell korrelierter Organ- und Verhaltensmerkmale: Ei-Attrappen und Maulbrüten bei afrikanischen Cichliden. Zeitschrift für Tierpsychologie. 1962, 19: 129-164.View Article
- Wickler W: 'Egg-dummies' as natural releasers in mouth-breeding cichlids. Nature. 1962, 194: 1092-1093.View Article
- Sturmbauer C, Meyer A: Genetic divergence, speciation and morphological stasis in a lineage of African cichlid fishes. Nature. 1992, 358: 578-81. 10.1038/358578a0.View ArticlePubMed
- Nagl S, Tichy H, Mayer WE, Takezaki N, Takahata N, Klein J: The origin and age of haplochromine fishes in Lake Victoria, east Africa. Proc R Soc Lond B Biol Sci. 2000, 267: 1049-61. 10.1098/rspb.2000.1109.View Article
- Moran P, Kornfield I: Retention of an ancestral polymorphism in the Mbuna Species Flock (Teleostei: Cichlidae) of Lake Malawi. Mol Biol Evol. 1993, 10: 1015-1029.
- Salzburger W, Baric S, Sturmbauer C: Speciation via introgressive hybridization in East African cichlids?. Mol Ecol. 2002, 11: 619-25. 10.1046/j.0962-1083.2001.01438.x.View ArticlePubMed
- Rüber L, Meyer A, Sturmbauer C, Verheyen E: Population structure in two sympatric species of the Lake Tanganyika cichlid tribe Eretmodini: evidence for introgression. Mol Ecol. 2001, 10: 1207-25. 10.1046/j.1365-294X.2001.01259.x.View ArticlePubMed
- Seehausen O: Hybridization and adaptive radiation. Trends in Ecology and Evolution. 2004, 19: 198-207. 10.1016/j.tree.2004.01.003.View ArticlePubMed
- Sturmbauer C, Hainz U, Baric S, Verheyen E, Salzburger W: Evolution of the tribe Tropheini from Lake Tanganyika: synchronized explosive speciation producing multiple evolutionary parallelism. Hydrobiologia. 2003, 500: 51-64. 10.1023/A:1024680201436.View Article
- Koblmüller S, Salzburger W, Sturmbauer C: Evolutionary relationships in the sand-dwelling cichlid lineage of lake Tanganyika suggest multiple colonization of rocky habitats and convergent origin of biparental mouthbrooding. J Mol Evol. 2004, 58: 79-96. 10.1007/s00239-003-2527-1.View ArticlePubMed
- Baric S, Salzburger W, Sturmbauer C: Phylogeography and evolution of the tanganyikan cichlid genus Tropheus based upon mitochondrial DNA sequences. J Mol Evol. 2003, 56: 54-68. 10.1007/s00239-002-2380-7.View ArticlePubMed
- Meyer A, Montero C, Spreinat A: Evolutionary history of the cichlid fish species flocks of the East African great lakes inferred from molecular phylogenetic data. Advances In Limnology. 1994, 44: 409-425.
- De Vos L, Seegers L: Seven new Orthochromis species (Teleostei: Cichlidae) from the Malagarasi, Luiche and Rugufu basins (Lake Tanganyika drainage), with notes on their reproductive biology. Ichthyol Explor Freshwaters. 1998, 9: 371-420.
- Boulenger G: Catalogue of the fresh-water fishes of Africa in the Brithish Museum (Natural History). 1898, London: Brithish Museum for Natural History
- Page RDM, Holmes EC: Molecular evolution: a phylogenetic approach. 1998, Oxford: Blackwell Sciences
- Gillespie JH: The causes of Molecular Evolution. 1991, Oxford: Oxford University Press
- Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19: 301-2. 10.1093/bioinformatics/19.2.301.View ArticlePubMed
- Greenwood PH: The zoogeography of African freshwater fishes: bioaccountancy or biogeography?. Evolution, Time and Space: The Emergence of the Biosphere. Edited by: Sims RW, Price JH, Whalley PES. 1983, London: Academic Press, 23: 179-199.
- Coulter GW: Zoogeography, affinities and evolution, with special regard to the fish. Lake Tanganyika and its Life. Edited by: Coulter GW. 1991, New York: Oxford University Press, 275-305.
- Beadle L: The Inland Waters of Tropical Africa. 1974, London: Longman
- Liem KF: Evolutionary strategies and morphological innovations: cichlid pharyngeal jaws. Syst Zool. 1973, 22: 425-441.View Article
- Hert E: The function of egg-spots in an African mouth-brooding cichlid fish. Anim Behav. 1989, 37: 726-732. 10.1016/0003-3472(89)90058-4.View Article
- Mrowka W: Egg stealing in a mouthbrooding cichlid fish. Anim Behav. 1987, 35: 923-925.View Article
- Kuwamura T: The evolution of parental care and mating systems among Tanganyikan cichlids. Fish Communities in Lake Tanganyika. Edited by: Kawanabe H, Hori M, Nagoshi M. 1997, Kyoto, Japan: Kyoto University Press, 59-86.
- Goodwin NB, Balshine-Earn S, Reynolds JD: Evolutionary transitions in parental care in cichlid fish. Proc R Soc Lond B Biol Sci. 1998, 265: 2265-2272. 10.1098/rspb.1998.0569.View Article
- Barlow GW: The Cichlid Fishes. Nature's Grand Experiment in Evolution. 2000, Cambridge, MA: Perseus Publishing
- Mayr E: Evolution of fish species flocks: a commentary. Evolution of Fish Species Flocks. Edited by: Echelle AA, Kornfield I. 1984, Orono, Maine: University of Maine at Orono Press, 3-12.
- Dominey WJ: Effects of sexual selection and life history on speciation: Species flocks in African cichlids and Hawaiian Drosophila. Evolution of Fish Species Flocks. Edited by: Echelle AA, Kornfield I. 1984, Orono, Maine: University of Maine at Orono Press, 231-250.
- Danley PD, Kocher TD: Speciation in rapidly diverging systems: lessons from Lake Malawi. Mol Ecol. 2001, 10: 1075-86. 10.1046/j.1365-294X.2001.01283.x.View ArticlePubMed
- Panhuis TM, Butlin R, Zuk M, Tregenza T: Sexual selection and speciation. Trends in Ecology and Evolution. 2001, 16: 364-371. 10.1016/S0169-5347(01)02160-7.View ArticlePubMed
- Turner GF, Burrows MT: A model of sympatric speciation by sexual selection. Proc R Soc Lond B Biol Sci. 1995, 260: 287-292.View Article
- Simpson GG: The major features of evolution. 1953, New York: Columbia University Press
- Schluter D: The Ecology of Adaptive Radiation. 2000, New York: Oxford University Press
- Kosswig C: Ways of speciation in cichlids. Copeia. 1963, 1963: 238-244.View Article
- Kirschner M, Gerhart J: Evolvability. Proc Natl Acad Sci U S A. 1998, 95: 8420-8427. 10.1073/pnas.95.15.8420.PubMed CentralView ArticlePubMed
- Sturmbauer C, Verheyen E, Meyer A: Mitochondrial phylogeny of the Lamprologini, the major substrate spawning lineage of cichlid fishes from Lake Tanganyika in Eastern Africa. Mol Biol Evol. 1994, 11: 691-703.PubMed
- Schluter D, McPhail JD: Character displacement and replicate adaptive radiation. Trends in Ecology and Evolution. 1993, 8: 197-200. 10.1016/0169-5347(93)90098-A.View ArticlePubMed
- Losos JB, Jackmann TR, Larson A, De Queiroz K, Rodrigues-Schettino L: Contingency and determinism in replicated adaptive radiations of island lizards. Science. 1998, 279: 2115-2118. 10.1126/science.279.5359.2115.View ArticlePubMed
- Bernatchez L, Chouinard A, Guoquing L: Integrating molecular genetics and ecology in studies of adaptive radiation: whitefish, Coregonus sp., as a case study. Biological Journal of the Linnean Society. 1999, 68: 173-194. 10.1006/bijl.1999.0337.View Article
- Shaw PW, Turner GF, Idid MR, Robinson RL, Carvalho GR: Genetic population structure indicates sympatric speciation of Lake Malawi pelagic cichlids. Proc R Soc Lond B Biol Sci. 2000, 267: 2273-80. 10.1098/rspb.2000.1279.View Article
- Swofford DL: PAUP* – Phylogenetic Analyses Using Parsinomy and other methods, version 4.0. Book PAUP* – Phylogenetic Analyses Using Parsinomy and other methods, version 4.0 (Editor ed.^eds.), 4.0 ed. 2003, City: Sinauer
- Mayer WE, Tichy H, Klein J: Phylogeny of African cichlid fishes as revealed by molecular markers. Heredity. 1998, 80: 702-714.View ArticlePubMed
- Takahashi K, Terai Y, Nishida M, Okada N: Phylogenetic relationships and ancient incomplete lineage sorting among cichlid fishes in Lake Tanganyika as revealed by analysis of the insertion of retroposons. Mol Biol Evol. 2001, 18: 2057-66.View ArticlePubMed
- Terai Y, Takahashi K, Nishida M, Sato T, Okada N: Using SINEs to probe ancient explosive speciation: "hidden" radiation of African cichlids?. Mol Biol Evol. 2003, 20: 924-30. 10.1093/molbev/msg104.View ArticlePubMed
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-8. 10.1093/bioinformatics/14.9.817.View ArticlePubMed
- Strimmer K, von Haeseler A: Quartet puzzling: A quartet maximum-likelihood method for reconstructing tree topologies. Mol Biol Evol. 1996, 13: 964-969.View Article
- Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol. 2002, 19: 101-109.View ArticlePubMed
- Delvaux D: Age of Lake Malawi (Nyasa) and water level fluctuations. Mus R Afr Centr Tervuren (Belg) Dept Geol Min Rapp Ann. 1995, 1995–1996: 99-108.
- Sturmbauer C, Baric S, Salzburger W, Rüber L, Verheyen E: Lake level fluctuations synchronize genetic divergences of cichlid fishes in African lakes. Mol Biol Evol. 2001, 18: 144-54.View ArticlePubMed
- Cohen AS, Lezzar KE, Tiercelin JJ, Soreghan M: New paleogeographic and lake-level reconstructions of Lake Tanganyika: implications for tectonic, climatic and biological evolution in a rift lake. Basin Research. 1997, 7: 107-132. 10.1046/j.1365-2117.1997.00038.x.View Article
- Lezzar KE, Tiercelin JJ, De Batist M, Cohen AS, Bandora T, Van Regensbergen P, Le Turdu C, Mifundu W, Klerkx J: New seismic stratigraphy and Late Tertiary history of the North Tanganyika Basin, East Africn Rift system, deduced from multichannel and high-resolution reflection seismic data and piston core evidence. Basin Research. 1996, 8: 1-28.View Article
- Tiercelin JJ, Mondeguer A: The geology of the Tanganyika Trough. Lake Tanganyika and its Life. Edited by: Coulter GW. 1991, New York: Oxford University Press, 7-48.
- Cohen AS, Soreghan M, Scholz CA: Estimating the age of formation of lakes: An example from Lake Tanganyika, East African Rift system. Geology. 1993, 21: 511-514. 10.1130/0091-7613(1993)021<0511:ETAOFO>2.3.CO;2.View Article
- Lee WJ, Conroy J, Howell WH, Kocher TD: Structure and evolution of teleost mitochondrial control regions. J Mol Evol. 1995, 41: 54-66. 10.1007/BF00174041.View ArticlePubMed
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.