Mitochondrial phylogeny and phylogeography of East African squeaker catfishes (Siluriformes: Synodontis)

  • Stephan Koblmüller1,

    Affiliated with

    • Christian Sturmbauer1,

      Affiliated with

      • Erik Verheyen2,

        Affiliated with

        • Axel Meyer3 and

          Affiliated with

          • Walter Salzburger3Email author

            Affiliated with

            BMC Evolutionary Biology20066:49

            DOI: 10.1186/1471-2148-6-49

            Received: 10 April 2006

            Accepted: 19 June 2006

            Published: 19 June 2006

            Abstract

            Background

            Squeaker catfishes (Pisces, Mochokidae, Synodontis) are widely distributed throughout Africa and inhabit a biogeographic range similar to that of the exceptionally diverse cichlid fishes, including the three East African Great Lakes and their surrounding rivers. Since squeaker catfishes also prefer the same types of habitats as many of the cichlid species, we hypothesized that the East African Synodontis species provide an excellent model group for comparative evolutionary and phylogeographic analyses.

            Results

            Our analyses reveal the existence of six major lineages of Synodontis in East Africa that diversified about 20 MYA from a Central and/or West African ancestor. The six lineages show a clear geographic patterning. Two lineages are endemic to Lake Tanganyika (plus one non-endemic representative), and these are the only two Synodontis lineages that diversified further into a small array of species. One of these species is the cuckoo catfish (S. multipunctatus), a unique brood parasite of mouthbrooding haplochromine cichlids, which seems to have evolved in parallel with the radiation of its cichlid host lineage, the Tropheini. We also detect an accelerated rate of molecular evolution in S. multipunctatus, which might be the consequence of co-evolutionary dynamics.

            Conclusion

            We conclude that the ancestral lineage of today's East African squeaker catfish fauna has colonized the area before the Great Lakes have formed. This ancestor diversified rapidly into at least six lineages that inhabit lakes and rivers in East Africa. Lake Tanganyika is the only lake harboring a small species flock of squeaker catfishes.

            Background

            The Great Lakes in the East African Rift Valley (Fig. 1a) are home to an exceptionally diverse ichthyofauna [14]. The most famous elements of the lakes' faunas are the cichlid fishes that have formed species flocks of an unparalleled species-richness and degree of eco-morphological and behavioral complexity [1, 5, 6]. It has been estimated that almost 1,800 cichlid species inhabit Lakes Tanganyika, Malawi and Victoria [7, 8] and that these lake endemic species are likely to have evolved in the last few millions or as recently as the last thousands of years only [912]. It is, thus, not surprising that the cichlid species flocks from the East African Great Lakes have received considerable attention as model systems for the study of adaptive radiation and explosive speciation [6, 13, 14].
            http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-6-49/MediaObjects/12862_2006_Article_234_Fig1_HTML.jpg
            Figure 1

            The squeaker catfishes (Synodontis) in East Africa. (a) Map of southern and eastern Africa showingthe Great Lakes in East Africa and the main river systems in thearea. (b) The upside-down catfish (S. nigriventris) is characterized by its inverse swimming posture. Photo courtesy of E. Schraml. (c) The cuckoo catfish(S. multipunctatus) from Lake Tanganyika is a brood parasiteof mouthbrooding haplochromine cichlids. Here, the mouth-content ofa breeding female of Simochromis diagramma is shown: Wefound five larvae of S. diagramma plus one larger cuckoocatfish larvae (Photo: W. Salzburger).

            The family Cichlidae is, however, just one of several families of freshwater fishes that inhabit the East African Great Lakes [1, 15]. In Lake Tanganyika, the oldest lake in East Africa, more than 100 non-cichlid fish species in at least twenty families are recognized, about half of the species being endemic [1, 16, 17] (Table 1). Interestingly, all Tanganyikan fish families are known to also occur in the Congo River basin that is definitely older and holds a much larger number of fish families than any of the East African water-bodies [15, 18]. An initially West African to Central African origin has therefore been proposed for the East African representatives of many of these families [18]. So far, such a scenario has only been confirmed with phylogenetic analyses for the Cichlidae [9, 12, 1922], and Clupeidae (Wilson & Meyer, unpublished), while the Clariidae seem to have colonized Africa from the Arabian plate [23]. Eleven non-cichlid families with about 50 species have been described for Lake Malawi. Only one of these families is not represented in lakes Tanganyika and Victoria: A single anguillid species, the African mottled eel (Anguilla bengalensis labiata), migrates to the area of Lake Malawi via the Zambezi and Shire River. In Lake Victoria, representatives of thirteen non-cichlid fish families are known, two of which were introduced only recently. All other families are also found in Lake Tanganyika (Table 1). Aside of the Cichlidae, nine other families are known to live in all three Great Lakes: Bagridae, Characidae, Clariidae, Cyprinidae, Cyprinodontidae, Mastacembalidae, Mochokidae, Mormyridae, and Protopteridae.
            Table 1

            Freshwater fish families in the East African Great Lakes Victoria, Malawi and Tanganyika (according to refs. [1,16,17] and http://​www.​fishbase.​org)

            fish family

            Lake Victoria

            Lake Malawi

            Lake Tanganyika

            Amphiliidae

             

            Anabantidae

             

            Anguillidae

             

             

            Bagridae

            Centrarchidae

            √*

              

            Centropomidae

            √*

             

            Characidae

            Cichlidae

            Citharinidae

              

            Clariidae

            Clupeidae

              

            Cyprinidae

            Cyprinodontidae

            Kneriidae

              

            Malapteruridae

              

            Mastacembalidae

            Mochokidae

            Mormyridae

            Polypteridae

              

            Protopteridae

            Schilbeidae

             

            Tetraodontidae

              

            *introduced

            We were interested in the phylogeny and phylogeography of one genus of the catfish family Mochokidae, the squeaker catfishes (Synodontis), in East Africa. With about 190 species assigned to ten genera, the siluriform family Mochokidae is endemic to the African continent where these fishes inhabit freshwaters from small creeks to large streams and from small ponds to large lakes. The genus Synodontis Cuvier, 1816 (including the monotypic genera Brachysynodontis and Hemisynodontis) comprises about 120 species, and is the most species-rich and widespread genus of the Mochokidae. Synodontis species are found throughout Africa, except in the southernmost parts and the Maghreb, although most species occur in Central and West Africa [24]. They are called "squeakers" (also by many local tribes), because of the noise they produce when taken out of the water. Some larger Synodontis species are important food resources. Other species are traded as ornamental fish due to their coloration and because of the interesting behavior of some of the species. The upside-down catfishes (S. contractus, S. nigriventris), for example, show an inverse swimming posture (Fig. 1b). Even more bizarre is probably the brood parasitic cuckoo catfish from Lake Tanganyika (S. multipunctatus) that cedes its eggs among those of spawning mouthbrooding cichlids; the female cichlid inadvertently takes up the catfish-eggs into her buccal cavity, together with her own eggs, where the young catfishes nourish on the cichlid larvae (Fig. 1c) [25].

            Our interest in the East African Synodontis was founded on several reasons: Firstly, the Mochokidae are abundant in all three East African Great Lakes and are, as a consequence, an ideal system for comparative phylogenetic and phylogeographic analyses with the existing data on cichlid fishes. Secondly, of the nine non-cichlid families common to the three East African Great Lakes, the Mochokidae show the closest overlap in habitat in comparison to cichlids, in that they occur primarily in the littoral and sub-littoral zone [1, 16]. Thirdly, among the two East African mochokid genera, Synodontis is more species rich as compared to Chiloglanis, particularly in Lake Tanganyika, where seven species of Synodontis have been described. Finally, most Synodontis species can be considered as euryphagous, which enables them to cope with seasonal changes in food abundance and habitat change [26, 27]. Their feeding regime broadens their food niche and gives them a better ability to colonize different habitats as compared to more specialized fish species such as many lake-adapted cichlids. Thus, their colonization routes in East Africa should be good indicators of fish-accessible waterways in that region over evolutionary time spans.

            In order to obtain a better understanding of the evolution, diversity and phylogeography of East African squeaker catfish, we aimed to reconstruct the phylogeny of East African Synodontis with special emphasis on the sequence and timing of colonization events of the East African Great Lakes Tanganyika, Malawi and Victoria, and the evolution of "cuckooing" behavior in Lake Tanganyika. Therefore, we sequenced a portion of about 900 bp of the mitochondrial genome of 21 species of the genus Synodontis and the outgroup taxa Microsynodontis batesii and Chiloglanis sp. We included all but one described, plus one, until now, undescribed, species from Lake Tanganyika, the single described Lake Malawi Synodontis, one of the two species occurring in Lake Victoria, two species from East and South African river systems, as well as eleven representatives of Central and West African species (Table 2).
            Table 2

            List of specimens used for this study

            Taxon

            origin*

            sample locality

            collector**

            Taxon ID

            GenBank acc. nr. CR/ND6

            Synodontis njassae

            EA (LM)

            M'bamba Bay

            EV

            MbBay (1995)

            DQ663012

            DQ663077

               

            EV

            Mst45 (1995)

            DQ663023

            -

               

            EV

            M56M5 (1995)

            DQ663007

            DQ663070

               

            EV

            M5336 (1995)

            DQ663005

            DQ663068

               

            EV

            M5350 (1995)

            DQ663006

            DQ663069

               

            EV

            M7043 (1995)

            DQ663008

            DQ663071

               

            EV

            M9008 (1995)

            DQ663009

            DQ663072

               

            EV

            M9016 (1995)

            DQ663022

            -

               

            EV

            M9043 (1995)

            DQ663010

            DQ663073

               

            EV

            M9076 (1995)

            -

            DQ663074

               

            EV

            M9092 (1995)

            DQ663011

            DQ663075

               

            EV

            M9096 (1995)

            DQ663024

            -

               

            EV

            T995 (1995)

            DQ663025

            -

            Synodontis victoriae

            EA (LV, LK, M)

            Kisumu, LV

            WS

            5i8a (2004)

            DQ663013

            DQ663059

              

            Kisumu, LV

            WS

            5i8b (2004)

            DQ663014

            DQ663060

            Synodontis dhonti

            EA (LT)

            Tanganyika Lodge

            SK, CS

            3698

            DQ662994

            -

              

            Kasakalawe Lodge

            SK, CS

            3708

            DQ662995

            -

              

            Kasakalawe Lodge

            SK, CS

            3709

            DQ662996

            -

              

            Chituta Bay

            SK, CS

            3736

            DQ662997

            DQ663056

              

            Mtondwe Island

            SK, CS

            3738

            DQ662998

            -

              

            Katoto

            SK, CS

            3771

            DQ662999

            -

              

            Kasakalawe

            SK, CS

            3777

            DQ663000

            DQ663057

              

            Kasakalawe

            SK, CS

            3778

            DQ663001

            -

              

            Kalambo Lodge

            SK, CS

            4535

            DQ663002

            -

              

            Chituta Bay

            TV

            4536

            DQ663003

            -

              

            Chituta Bay

            TV

            4537

            DQ663004

            -

              

            Ulwile Island

            EV

            TB 127

            DQ663019

            DQ663065

              

            Ulwile Island

            EV

            TB 128

            DQ663020

            DQ663066

            Synodontis granulosus

            EA (LT)

            Chituta Bay

            SK, CS

            3735

            DQ662961

            -

              

            Mtondwe Island

            SK, CS

            3773

            DQ662962

            DQ663053

            Synodontis multipunctatus

            EA (LT)

            Muzumwa

            SK, CS

            3635

            DQ662963

            DQ663046

              

            Kalambo Lodge

            SK, CS

            3696

            DQ662964

            -

              

            south of Kalambo

            SK, CS

            3697

            DQ662965

            -

              

            Mpulungu

            SK, CS

            3712

            DQ662968

            -

              

            Mpulungu

            SK, CS

            3714

            DQ662969

            -

              

            Mpulungu

            SK, CS

            3715

            DQ662966

            -

              

            Mpulungu

            SK, CS

            3716

            DQ662967

            -

              

            Kasakalawe

            SK, CS

            3737

            DQ662970

            -

              

            ?***

            SK, CS

            3739

            DQ662971

            DQ664047

              

            Muzumwa

            SK, CS

            4112

            DQ662972

            -

              

            Muzumwa

            SK, CS

            4113

            DQ662973

            -

              

            Muzumwa

            SK, CS

            4114

            DQ662974

            -

              

            Chituta Bay

            TV

            4538

            DQ662975

            -

              

            Ulwile Island

            EV

            TB 073

            -

            DQ663062

              

            Ulwile Island

            EV

            TB 125

            DQ663015

            DQ663063

              

            Ulwile Island

            EV

            TB 126

            DQ663016

            DQ663064

              

            south of Mkangansi

            EV

            TB 872

            DQ663017

            DQ663067

            Synodontis petricola

            EA (LT)

            Kapembwa

            SK, CS

            3636

            DQ662976

            DQ663048

              

            Chituta Bay

            SK, CS

            3699

            DQ662977

            DQ663050

              

            Mpulungu

            SK, CS

            3731

            DQ662978

            -

              

            Mpulungu

            SK, CS

            3732

            DQ662979

            -

              

            Katoto

            SK, CS

            3769

            DQ662980

            -

              

            Katoto

            SK, CS

            3770

            DQ662981

            -

            Synodontis polli

            EA (LT)

            Kalambo Lodge

            SK, CS

            3695

            DQ662983

            DQ663049

              

            Kasakalawe Lodge

            SK, CS

            3710

            DQ662984

            -

              

            Kasakalawe Lodge

            SK, CS

            3711

            DQ662985

            -

              

            Mpulungu

            SK, CS

            3717

            DQ662986

            -

              

            Mpulungu

            SK, CS

            3733

            DQ662987

            DQ663052

              

            ?***

            SK, CS

            4030

            DQ662988

            -

              

            Mtondwe Island

            SK

            4539

            DQ662982

            -

              

            Ulwile Island

            EV

            TB 072

            DQ663018

            DQ663061

            Synodontis sp. nov.

            EA (LT)

            Mpulungu

            SK, CS

            3718

            DQ662989

            DQ663051

              

            Mpulungu

            SK, CS

            3775

            DQ662990

            DQ663054

              

            Gombe

            SK, CS

            3776

            DQ662991

            -

              

            ?***

            SK, CS

            4031

            DQ662992

            DQ663055

              

            Tanganyika Lodge

            SK, CS

            4532

            DQ662993

            -

            Synodontis nigromaculatus

            EA (LT, R)

            Lake Mweru

            RS

            4342

            DQ662954

            DQ663043

              

            Lake Mweru

            RS

            4339

            DQ662955

            DQ663045

            Synodontis zambezensis

            EA (R)

            Cinzombo Lodge; Luagwa R.

            CK

            3977

            DQ662952

            -

              

            Nsefu Lagoon; Luangwa R.

            CK

            3978

            DQ662953

            DQ663039

            Synodontis angelicus

            CA (Congo R.)

            ?***

            SK

            3744

            DQ662944

            DQ663033

              

            ?***

            SK

            3745

            DQ662945

            -

            Synodontis brichardi

            CA (Congo R.)

            ?***

            SK

            3993

            DQ662943

            DQ663035

            Synodontis eupterus

            WA, Nile

            ?***

            SK

            3861

            DQ662956

            DQ663036

              

            ?***

            SK

            4088

            DQ662957

            -

              

            ?***

            SK

            4108

            DQ662958

            DQ663037

            Synodontis flavitaeniatus

            CA (Congo R.)

            ?***

            JS

            75

            -

            DQ663077

            Synodontis nigrita

            WA, Nile

            ?***

            SK

            3868

            DQ662947

            DQ663038

              

            ?***

            SK

            4337

            DQ662948

            DQ663041

              

            ?***

            SK

            4338

            DQ662949

            DQ663042

            Synodontis nigriventris

            CA (Congo R.)

            ?***

            SK

            3700

            DQ662941

            DQ663030

              

            ?***

            SK

            3701

            DQ662942

            -

              

            ?***

            JS

            56

            DQ663021

            DQ663058

            Synodontis ocellifer

            WA

            ?***

            SK

            3994

            DQ662960

            DQ663031

            Synodontis polystigma

            CA (Congo R.)

            Lake Mweru

            RS

            4341

            DQ662951

            DQ663044

            Synodontis robbianus

            WA (Niger R.)

            ?***

            SK

            3702

            DQ662938

            DQ663029

              

            ?***

            SK

            3703

            DQ662939

            -

              

            ?***

            SK

            4089

            DQ662940

            -

            Synodontis robertsi

            CA (Congo R.)

            ?***

            SK

            4336

            DQ662946

            DQ663040

            Synodontis schoutedeni

            CA (Congo R.)

            ?***

            SK

            3860

            DQ662950

            DQ663034

            Synodontis velifer

            WA

            ?***

            WK

            3995

            DQ662959

            DQ663032

            Chiloglanis sp.

             

            Lunzua R.

            SK, CS

            3865

            DQ662935

            DQ663026

            Microsynodontis batesii

             

            ?***

            SK

            4024

            DQ662936

            DQ663027

              

            ?***

            SK

            4025

            DQ662937

            DQ663028

            *CA, Central Africa; LK, Lake Kivu; LM, Lake Malawi; LT, Lake Tanganyika; LV, Lake Victoria; M, Malagarazi River; R, East African rivers; WA, West Africa;

            **CS, Christian Sturmbauer; CK, Cyprian Katongo; EV, Erik Verheyen; JS, Jos Snoeks; RS, Robert Sinyinza; SK, Stephan Koblmüller; TV, Toby Veall; WK, Walter Köhldorfer; WS, Walter Salzburger

            ***samples obtained through ornamental fish trade.

            Coordinates of sampling localities (when available):

            LT – Chituta Bay, S 08°44' E 31°09'; Gombe, S 08°47' E 31°01'; Kalambo Lodge, S 08°39' E 31°13'; Kapembwa, S 08°37' E 30°51'; Kasakalawe, S 08°47' E 31°04'; Kasakalawe Lodge, S 08°47' E 31°06'; Katoto, S 08°48' E 31°01'; Mpulungu, S 08°46' E 31°06'; Mtondwe Island, S 08°43' E 31°07'; Muzumwa, S 08°42' E 31°12'; south of Kalambo, S 08°37' E 31°12'; south of Mkangansi, S 06°37' E 30°18'; Tanganyika Lodge, S 08°47' E 31°05'; Ulwile Island, S 07°72' E 30°34'Lunzua R. (near Mwanangwa), S 08°57' E 31°10'

            Results

            For this study, we analyzed 407 bp of the mitochondrial control region of 94 specimens and 426 bp of the NADH dehydrogenase subunit 6 gene (ND6) of 52 specimens (including the outgroup taxa). The likelihood mapping analysis yielded a fraction of 89.4% fully resolved quartets for the total combined data set (Fig. 2a), and 88.7% fully resolved quartets for the subset including the East African taxa only (Fig. 2b), indicating a fairly strong phylogenetic signal in both data sets. Within the genus Synodontis, pairwise distances (uncorrected p-distances) ranged up to 33.3% (S. shoutedeni versus S. multipunctatus) in the control region and up to 11.5 % (S. granulosus versus S. robertsi) in the ND6 gene segment. The largest uncorrected p-distances within the East African clade were 23.7% in the control region and 5.9% in ND6 (in both cases involving the species pair S. njassae and S. multipunctatus). The two gene segments analyzed, thus, corroborate the observation that in the mitochondrial genome the non-coding control region generally evolves about two to five times faster than protein coding genes [28].
            http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-6-49/MediaObjects/12862_2006_Article_234_Fig2_HTML.jpg
            Figure 2

            Results from the likelihood mapping analysis. (a) In theanalysis of the entire dataset, a percentage fraction of 89.4 % of all quartets was fully resolved (and 7.8 % were not resolved). (b) In the reduced dataset, 88.7 % of the quartets were resolved (7.9 % were not resolved). This points to a strong phylogenetic signal in the data.

            In our phylogenetic analyses, the internal branches interrelating major Central and West African lineages were short and, hence, the tree topologies remained somewhat unresolved and the branching topology varied depending on the tree-building algorithm used. This seems to suggest that several Central and West African lineages originated almost contemporaneously in the course of a major cladogenetic event. However, a much more thorough sampling and the analysis of a more slowly evolving gene segment (e.g., RAG1 [29]) would be necessary to shed light on this part of the evolutionary history of Synodontis in Africa. Our analyses consistently revealed that the Central and West African members of the genus Synodontis are ancestral to the East African representatives. The monophyly of the East African clade is clearly seen in the composite consensus tree (Fig. 3) based upon six most parsimonious trees [tree length, 49,875 steps; consistency index (CI) excluding uninformative characters, 0.47; retention index (RI), 0.74; rescaled consistency index (RC), 0.38], the neighbor joining (NJ) tree, the maximum likelihood (ML; Fig. 4) tree and the Bayesian inference phylogeny. The monophyly of the East African representatives was further confirmed by the individual analysis of the mitochondrial gene segments, the control region and the ND6 gene (trees not shown).
            http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-6-49/MediaObjects/12862_2006_Article_234_Fig3_HTML.jpg
            Figure 3

            Composite consensus tree of the phylogenetic analyses. The strictconsensus of the neighbor-joining tree, the most parsimonious trees, the optimal maximum likelihood topology (see Fig.4) and the Bayesian inference tree is shown. Numbers above the branches are neighbor-joining and maximum parsimony bootstrap values, numbers below the branches represent maximum likelihood bootstraps and Bayesian posterior probabilities. The grey box indicates the East African clade of Synodontis.

            http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-6-49/MediaObjects/12862_2006_Article_234_Fig4_HTML.jpg
            Figure 4

            Maximum likelihood tree. Maximum likelihood topology based on the K81uf+I+Γ model of molecular evolution [70] with nucleotide frequencies A, 0.3581, C, 0.2676, G, 0.1368, T, 0.2375, proportion of invariable sites (I), 0.2461, gamma shape parameter (α), 0.7306, and R-matrix A↔G, A↔T, C↔G and G↔T, 1.0000; A↔G, 7.7875 and C↔T, 1.2463. The blue box indicates the East African clade of Synodontis.

            The phylogenetic analyses of the East African Synodontis species revealed a similar outcome compared to that of the ancestral Central and West African lineages, with relatively short branches interrelating the main lineages (and low bootstrap support and posterior probabilities for these branches), indicating an almost contemporaneous origin of the main East African Synodontis lineages. Six distinct mitochondrial lineages can be distinguished within the East African clade of Synodontis. The species S. nigromaculatus, S. zambezensis, S. njassae and S. victoriae constitute separate clades, while the remaining Lake Tanganyika species can be divided into two distinct genetic lineages. The first Tanganyikan clade contains S. granulosus and the cuckoo catfish (S. multipunctatus); the second Tanganyikan clade includes S. dhonti, S. petricola, S. polli and S. sp. nov. In the maximum parsimony (MP) and NJ tree of the reduced dataset, the two clades of Lake Tanganyika Synodontis were resolved as a monophyletic group, albeit with low bootstrap support. In ML and BI, on the other hand, the two Tanganyikan clades were not resolved as sister groups, and S. victoriae was placed as sister group to the clade comprising S. dhonti, S. petricola, S. polli and S. sp. nov. – with similarly low support values. The four-cluster-likelihood-mapping (Fig. 5) favors a topology with a monophyletic origin of Lake Tanganyika Synodontis, and S. victoriae being more closely related to S. nigromaculatus, S. zambezensis and S. njassae. A Shimodaira-Hasegawa test [30] based on the tree topologies obtained from the reduced data set revealed that there are no significant differences between the topologies obtained by the different algorithms (P < 0.05).
            http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-6-49/MediaObjects/12862_2006_Article_234_Fig5_HTML.jpg
            Figure 5

            Four-cluster likelihood mapping analysis. The topology with a monophyly of the two endemic Tanganyikan clades receives the highest support. LM, Lake Malawi (S. njassae); LT I, first Lake Tanganyika clade; LT II, second Lake Tanganyika clade; LV, Lake Victoria (S. victoriae); R, riverine representatives.

            Within the second Tanganyikan clade the branching order among the four species S. dhonti, S. petricola, S. polli, S. sp. nov. differed with respect to the tree-building algorithm used. Together with the low bootstrap values and posterior probabilities this again points to an event of rapid cladogenesis [31]. In the other Tanganyikan clade, we observed relatively long branches (see Fig. 4), and also the uncorrected pairwise distances indicated a faster evolutionary rate in the cuckoo catfish (see above). For example, the maximum uncorrected pairwise distance in the control region in S. multipunctatus was 6.5%, whereas in S. njassae a maximum uncorrected pairwise distance of only 1.5% was found. Relative rate tests with RRTree indeed recovered a significantly higher (P < 0.01) substitution rate for S. multipunctatus compared to all other East African species. This was found for both gene segments.

            In order to tentatively date the major cladogenetic events in Synodontis in East Africa, we used r8s [32] and applied different calibration points for our analyses. In our first round of analysis, in which we applied all available dates, we obtained an age of 12.24 million years (MY) for the most recent common ancestor (MRCA) of S. granulosus and S. multipunctatus (minimum estimate: 11.29 MY; maximum estimate: 13.17 MY), and an age of 2.15 MY (min.: 0.97; max.: 2.37) for the MRCA of the cuckoo catfish (S. multipunctatus) (see Table 3). When applying two calibration points only, we obtained highly similar values plus an age of 7.5 MY for the second Tanganyikan clade (min.: 5.48; max.: 12.11). The MRCA estimates did not change substantially when only the fossil calibration was used (Table 3). The age for the entire genus was estimated with ca. 35 MY (min.: 22.87; max.: 40.12). This result should, however, be interpreted with caution since these dates lie outside the range of our calibration points, and since we might not have represented the entire diversity of Synodontis in Central, South and West Africa.
            Table 3

            Dating of the major cladogenetic events in Synodontis with r8s. We ran three independent analyses in r8s [32] with different calibration points using the maximum estimated age for the lacustrine habitat in lakes Malawi (1 MY [49]) and Tanganyika (6 MY [46–48]), as well as the minimum age of the East African clade of Synodontis as suggested by the oldest known Synodontis fossil in that region (>20 MY [45]). In the first analysis, all three calibrations were applied (1/6/20 calibration); in the second cycle, we used the Lake Malawi and the fossil calibration (1/20 calibration); in the third round, we only used the fossil based calibration (20 calibration). The numbers indicate the average value (in MY) obtained from a bootstrap approach with 30 replicates, the minimum and maximum values are depicted in round brackets (in italics). Square brackets indicate the time constraints used for the different r8s analyses and the range of the actual numbers used in the bootstrap replicates (in round brackets). The estimates for the age of the entire genus Synodontis should be interpreted with caution, as the values lie outside our range of calibration points

            Clade

            1/6/20 MY calibration

            1/20 MY calibration

            20 MY calibration

            Synodontis

            34.54 (23.42–38.14)

            35.39 (24.01–38.53)

            37.26 (22.87–40.12)

            East African clade

            [>20 (20–20.49)]

            [>20 (20–20.45)]

            [>20 (20–22.53)]

            L. Tanganyika I

            12.24 (11.29–13.17)

            12.36 (11.36–13.29)

            12.89 (11.24–14.81)

            L. Tanganyika II

            [<6 (4.77–6)]

            7.5 (5.48–12.11)

            7.88 (5.76–12.36)

            S. multipunctatus

            2.15 (0.97–2.37)

            2.21 (1.61–2.43)

            2.31 (1.64–2.58)

            L. Malawi (S. njassae)

            [<1 (0.99–1)]

            [<1 (0.99–1)]

            1.75 (1.34–2.18)

            Discussion

            Mitochondrial phylogeny of the East African Synodontis species

            Although the genus Synodontis is the most species-rich and widespread genus of the Mochokidae accounting for about a quarter of all African siluriform species, there is no molecular phylogeny available for this group. Earlier attempts based on a small number of sequences of the mitochondrial cytochrome b gene revealed a sister-group relationship between an East African clade (represented by S. zambezensis, S. nigromaculatus, S. njassae, and S. petricola) and some South-African representatives [33]. However, no Central and/or West African taxa have been included in this analysis that would allow placing the investigated species into a larger phylogenetic and phylogeographic context.

            We were interested in the phylogeny and phylogeography of the East African Synodontis species and, therefore, combined the mtDNA sequences of East African representatives with those of eleven Central and/or West African species. Our analyses of the combined dataset as well as the separate analyses of each mitochondrial locus, consistently revealed the monophyly of the East African clade and its derived position compared to the Central and/or West African taxa (Figs. 3, 4). While the phylogenetic relationships between most of the more ancestral Central and/or West African representatives could not be resolved, which is most likely due to the large genetic distances among these taxa, we could identify a sister-group relationship between a clade comprised by S. velifer and S. eupterus and S. ocellifer plus the East African clade (Figs. 3, 4).

            Within the East African clade, we could clearly identify six distinct lineages: S. zambezensis; S. nigromaculatus; S. victoriae (Lake Victoria); S. njassae (Lake Malawi); S. granulosus and S. multipunctatus (Lake Tanganyika I); S. dhonti, S. polli, S. petricola, S. sp. nov. (Lake Tanganyika II). However, it was not possible to resolve the phylogenetic relationships between these lineages, as indicated by the short branches interconnecting them, the lacking statistical support for these branches, and the fact that different tree building algorithms led to different phylogenetic hypotheses with regard to the relative position of these lineages. Such a lack of phylogenetic signal in ancestral branches is frequently observed in cases of relatively rapid lineage splitting events [12, 22, 29, 31, 3436], e.g., as consequence of the synchronized fragmentation of a founder population into allopatric units [31]. We, thus, suggest a scenario for the evolution of Synodontis in East Africa, in which one ancestral (colonizing) species was split into at least six sub-lineages within a relatively short period of time, and that – at least in some cases – these sub-lineages diversified further, albeit in isolation from each other. In this context, it is interesting to note that hybrids can artificially be produced between two of the lineages, the riverine taxa S. zambezensis and S. nigromaculatus [33, 37]. It is, however, unlikely that hybridization distorted our phylogenetic signal, as most of the lineages exclude each other geographically so that no natural hybridization was and is possible. The specimens of the only two species that partially overlap in range (S. zambezensis and S. nigromaculatus) showed reciprocal monophyly and appeared as distinct evolutionary units with an average sequence divergence of ca. 12% in the control region and ca. 4.5% in ND6.

            Lake Tanganyika is the only East African lake that harbors a small assemblage of endemic Synodontis species [16, 17]. Our phylogenetic analyses indicate the existence of two separate lineages of Synodontis in Lake Tanganyika (see Fig. 3) (note that also the non-endemic S. nigromaculatus is found in that lake). The first lineage comprises the relatively large Tanganyikan representative S. granulosus and the cuckoo catfish (S. multipunctatus), a brood parasite of mouth-brooding cichlids [25] (Fig. 1c). The second Tanganyikan clade includes the species S. dhonti, S. polli, S. petricola, and a yet undescribed species. There are some indications, e.g., the four-cluster likelihood mapping analysis (Fig. 5), for a monophyly of the two Lake Tanganyikan clades, although statistical support in the individual bootstrap analysis and in the Shimodaira-Hasegawa test is lacking. While in both Tanganyikan clades the monophyly of the taxa was supported with high bootstrap values and Bayesian posterior probabilities, the phylogenetic relationships between the four lineages of the second Tanganyikan clade could not be resolved with high confidence (Fig. 3). We interpret these short branches as the result of the relatively rapid lineage splitting at the beginning of the formation of Lake Tanganyika. It also seems possible that some hybridization took place in the early phases of species formation (see e.g., [38]), which could be tested with nuclear DNA sequences. However, it is unlikely that the observed lack of resolution in some branches interrelating main lineages (or taxa) is due to the markers used, since, overall, we observe a sufficiently strong phylogenetic signal in the dataset (Fig. 2), and since the control region is the fastest evolving section of the mitochondrial genome.

            Phylo-chronology of Synodontis in East Africa

            The reconstruction of divergence times from molecular data or phylogenies is not without potential pitfalls (see e.g., [39, 40]). Yet, such time estimates are useful for providing an approximate framework that puts diversification events into a temporal framework. Our phylo-chronological reconstruction of the major cladogenetic events in the East African clade of Synodontis (Table 3) is based on three independent calibration points, which were used in different combinations. Because these calibration points are derived from estimates, it is important to note that our approximations should be interpreted cautiously. For example, we used a minimum age of 20 MY for the MRCA of the East African clade, based on the oldest known fossils in that area [4145]. These Early Miocene fossils are, however, not particularly well dated (although 20 MY seems like a good approximation [43]), they have not been assigned to extant species, and the possibility always remains that still older fossils exist. However, together with the internal calibration pointes – the ages of lake basins – the estimated scenario for the evolution and phylogeography of Synodontis in East Africa is congruent with geological and paleolimnological events.

            With regard to the entire East African clade of Synodontis, the oldest known fossils from the Early Miocene are relatively close to the estimated age of the MRCA of that clade. In our bootstrap replicates, the age of the MRCA of the East African clade was estimated between 20 (the constrained minimum age in the analyses) and 22.53 MY. For the MRCA of the first Tanganyikan clade, with S. granulosus and S. multipunctatus, we obtained an age of 11.24–14.81. This is older than the age of the onset of a truly lacustrine deepwater habitat in Lake Tanganyika [4648], suggesting that the split between these two taxa took place before the ecosystem of a deep tropical lake has formed (e.g., in the riverine environment of the Congo and/or Proto-Malagarazi River) or, possibly more likely, a colonization in the initial period of lake formation when extensive swampy areas with shallow basins are likely to have existed in the area of present Lake Tanganyika. Thus, although both S. granulosus and, particularly, the cuckoo catfish (S. multipunctatus) are highly specialized lake species, their ancestors were riverine fish that seem to have independently achieved their lacustrine lifestyle. In contrast, the second Tanganyikan clade appears to have undergone its primary radiation at a later stage, most likely in a lacustrine (or semi-lacustrine) environment (ca. 5–8 MYA; although the maximum values from individual bootstrap replicates indicate that even older estimates are possible). The only independent estimation for the MRCA of S. njassae from Lake Malawi revealed an age of 1.34–2.18 MY. We find that the cuckoo catfish (S. multipunctatus) has diverged after about 2 MYA (0.97–2.58 MYA). Interestingly, this would be immediately after the host lineage of the cuckoo catfish, the mouthbrooding haplochromine cichlids, has initiated its radiation in Lake Tanganyika, according to a recent calibration [9].

            Our strategy of carrying out the analysis in three cycles also aimed to test the validity of some calibration points, such as the age of the MRCA within Lake Malawi and the second clade in Lake Tanganyika. While our analyses revealed highly consistent results for the MRCA of the first Tanganyikan clade and that of S. multipunctatus (Table 3), the estimates for the MRCA of S. njassae (Lake Malawi) and that of the second Tanganyikan clade are slightly older than expected based on the presumed formation of the truly lacustrine habitat in both lakes [4649]. This could, of course, be due to analytical limitations in estimating divergence times based on molecular phylogenies. Alternatively, it may indicate the actual coalescence age of the riverine MRCA of both lineages. However, it could also mean that the lake basins are somewhat older or that the divergence of Synodontis in the basins of lakes Malawi and Tanganyika began before the establishment of the present lacustrine habitat with its characteristic deep-water conditions, e.g., in more shallow pools in the transition phase between a swampy habitat of a proto-lake and a real lake.

            Phylogeography of Synodontis in East Africa

            Based on our phylogenetic and phylo-chronological analysis of Synodontis in East Africa, we can outline the following overall phylogeographic scenario: It seems that after an initial colonization of East Africa through a potentially widespread Central and/or West African representative (probably about 20 MYA), at least six lineages of Synodontis evolved that are each confined to particular geographic regions. The rapid spread of East African Synodontis can be explained by a colonization of the area via the ancient Congo-Malagarazi River connection, and incidents of drainage catchments towards the Nile systems (in the North) and between the Upper Congo and Zambezi drainage (in the South) (see ref. [50] and citations therein). For example, the two riverine species S. nigromaculatus and S. zambezensis are widely distributed in Southern and Eastern Africa (but are phylogenetically distinct from the South African clade [33]). The black-spotted squeaker catfish (S. nigromaculatus) is found in the Cunene River, the Okavango River, the upper reaches of the Congo River, the upper Zambezi basin, the Limpopo as well as in Lake Tanganyika, demonstrating the existence of fish accessible water-ways through-out eastern and southern Africa until recently (see e.g., [9, 35]). The plain squeaker (S. zambezensis) is found in the Zambezi and Limpopo system, as well as in Lake Rukwa. The maximum-likelihood analysis (Fig. 4) suggested a sister-group relationship between the easternmost species, S. zambezensis, and S. njassae from Lake Malawi, and between S. nigromaculatus and the first Tanganyikan clade. Both hypotheses would, phylogeographically, make sense. However, a more detailed sampling and, particularly, much longer sequences would be necessary to confirm this hypothesis.

            With our sampling strategy and phylogenetic analyses we are able to present a robust hypothesis for the evolution of Synodontis in East Africa (see above). However, we are not able to address questions related to the phylogeography of single species as has been undertaken for the clariid Clarias gariepinus [51]. It would, for example, be interesting to compare the phylogeography of S. nigromaculatus to one branch of haplochromine cichlids that have a comparable distribution range [9, 35, 50]. Similarly, a detailed phylogeographic scenario for S. victoriae and S. zambezensis could be used to test the existence of past fish-accessible connections between large East African lakes (see e.g., [9, 10, 18]). A much more thorough sampling of selected species and phylogeographic hypotheses for single species would be necessary in order to resolve these open questions.

            Comparative phylogeography between East African cichlids and squeaker catfishes

            One of our goals for this study was to compare the phylogeny and phylogeography of the East African squeaker catfishes to the cichlid fishes, for which an impressive amount of data exists (see e.g., [912, 35, 52, 53]). Our analyses reveal that, first of all, the initial diversification of Synodontis in East Africa is much older compared to the radiation of the extant cichlid species flocks in that area. (Note that there are representatives of more ancient cichlid lineages in East Africa; some of these have colonized East Africa secondarily [12, 21].) In fact, the main lineages of Synodontis are likely to have existed before the oldest lake, Lake Tanganyika, has formed. A Central to West African origin is, however, likely for both groups, and it seems plausible that the ancestors of both the East African cichlid radiations and the East African Synodontis clade colonized the area at about the same time – with the Congo River system playing a zoogeographic key-role for the spread of both groups.

            A main difference between the East African cichlids and Synodontis is also the evolutionary role that Lake Tanganyika has played. Lake Tanganyika is not only a reservoir for the most ancestral cichlid lineages in East Africa, but it has been identified as the cradle from which the most species-rich and widespread group of East African cichlids, the haplochromines, has evolved (see e.g., [9, 12, 54]). This "out of Tanganyika scenario" [9] implies that the species flocks of lakes Malawi and Victoria are derived from the older Tanganyikan radiation [9, 12, 54]. In Synodontis, we do not find evidence that a formerly endemic Tanganyikan lineage left the lake secondarily to colonize rivers and, subsequently, other lakes. Instead, Lakes Malawi and Victoria were seeded independently by ancestral lineages of Synodontis. Consequently, different and – in most cases – more ancient connections between water-systems were used by Synodontis to colonize East Africa.

            There are, however, also similarities in the distribution of cichlids and squeaker catfishes. The range of the black-spotted squeaker catfish (S. nigromaculatus), for example, is similar to that of one lineage of haplochromine cichlids, the Congolese-South African clade (according to ref. [9]), and S. zambezensis overlaps with Pseudocrenilabrus [9, 50]. A much more detailed sampling would be necessary to test for a putative congruence of time-scale and routes of colonization between these riverine cichlids and squeaker catfishes. Similarly, the occurrence of S. victoriae in lakes Kivu and Victoria seems to support the close relatedness of the fish faunas of both lakes [10]. A detailed phylogeographic scenario for S. victoriae would thus be of interest to further test for the commonality of geological and palaeolimnological events affecting several lineages of fishes.

            The evolution of brood parasitism in the cuckoo catfish (S. multipunctatus)

            The most fascinating behavioral adaptation in Synodontis is certainly the evolution of brood parasitism in S. multipunctatus from Lake Tanganyika [25] (Fig. 1c). In the course of our sampling we found larvae of S. multipunctatus in the mouth of Simochromis babaulti and S. diagramma females (with up to three catfishes per female).

            Our analyses resulted in two interesting findings regarding the evolution of cuckoo behavior in S. multipunctatus. First, our phylo-chronological analyses consistently revealed an age of about 2 MY for the MRCA of S. multipunctatus. This coincides with the MRCA of the mouthbrooding haplochromines of Lake Tanganyika [9], to which the current host-species of S. multipunctatus belong. It thus seems that brood parasitism evolved in S. multipunctatus at an early stage of the radiation of the Tropheini in Lake Tanganyika. Secondly, we find an accelerated rate of molecular evolution in both gene segments in the representatives of S. multipunctatus, as shown by longer branches in the phylogenies, larger genetic distances, and relative rate tests. A faster rate of DNA sequence evolution in brood parasites has been described in African finches, and it has been argued that this might be the consequence of a dynamic co-evolutionary history [55]. Iterative co-evolutionary processes between hosts and parasites might account for the observed significantly faster evolutionary rate in S. multipunctatus. Alternatively, there might be selective pressure to equalize the evolutionary rate in host-parasite systems. This could possibly explain why the mutation rate in the non-coding control region in S. multipunctatus approaches that observed in haplochromine cichlids [56] (6.5%; MRCA: 2 MYA), while the remaining representatives of the East African clade of Synodontis show considerably lower evolutionary rates (e.g., S. njassae : 1.5%; MRCA: 1–1.5 MYA).

            Conclusion

            Our phylogenetic and phylogeographic analyses of East African Synodontis species suggest the existence of at least six main lineages of Synodontis in East Africa that began to diversify about 20 MYA from a Central and/or West African ancestor. The six lineages show a clear-cut geographic pattern, and only the distributions of two riverine species, S. nigromaculatus and S. zambezensis, overlap. Representatives of three lineages are found in the oldest lake in the area, Lake Tanganyika, and two of these lineages are endemic and have undergone further diversification within the lake. The initial diversification of Synodontis in East Africa appears to be much older than that of the cichlid species flocks, yet, some taxa show similar geographic patterns compared to haplochromine cichlids, with which they often share the same habitat (in rivers and lakes). The endemic brood-parasitic cuckoo catfish (S. multipunctatus) from Lake Tanganyika seems to have diversified in parallel with its host species of the mouthbrooding modern haplochromines (the Tropheini). The significantly faster evolutionary rate in S. multipunctatus might be a consequence of co-evolutionary dynamics. Although Synodontis has colonized East Africa before the formation of the Great Lakes, which are home to the exceptionally diverse cichlid species flocks, and although the squeaker catfishes inhabit the same habitats as many of the cichlids, only a few species have evolved in one of the lakes. This, once more, illustrates the unparalleled, and still enigmatic, propensity of cichlid fishes to undergo explosive speciation.

            Methods

            Specimen information and DNA methods

            This study is based upon a total of 91 individuals of 21 described and one so far undescribed species of Synodontis, mainly from Lakes Tanganyika, Malawi and Victoria, and several African river systems. As outgroup taxa we used one individual of Chiloglanis sp. and two individuals of Microsynodontis batesii [57, 58]. Most of the specimens were sampled during several field expeditions from 1992 to 2004, while some additional samples were obtained from the aquarium trade (Table 2). Voucher specimens are available from the authors. Of all specimens, fin clips were taken and preserved in 96% ethanol.

            We analyzed a 407 bp segment of the most variable part of the mitochondrial control region (D-loop) of 94 specimens (including the outgroup taxa) and 426 bp of the more slowly evolving NADH dehydrogenase subunit 6 gene (ND6) of 52 specimens (including the outgroup taxa). We chose the mitochondrial control region to be able to directly compare our present results to the East African cichlid species flocks, for which the control region is the most common mitochondrial DNA marker [9, 10, 12, 35, 50, 59, 60]. The ND6 gene was chosen because it has previously been used in clariid catfishes [51]. We decided not to use the mitochondrial cytochrome b gene, because previous analyses indicated a reduced phylogenetic signal in that gene with respect to the East African Synodontis species [33].

            For DNA extraction we applied a proteinase K digestion followed by protein precipitation with ammonium acetate. As primers for amplification and sequencing of the most variable part of the mitochondrial control region we used L-Pro-F [61] and TDK-D [62]. For ND6 we used the primers ND5G and ND6L [63]. The PCR reactions were prepared for a total volume of 17 μl containing 0.085 μl of Taq DNA polymerase (BioThermTM), 1.7 μl of each primer (10 μM), 1.7 μl 10x dNTP mix, 1.7 μl 10x buffer, 7.62 μl high performance liquid chromatography (HPLC) water, and 2.5 μl of the extracted DNA. Amplification was performed on a GeneAmp 9700 PCR system (Applied Biosystems) under the following conditions: an initial denaturation phase at 94°C for 3 min followed 45 cycles with denaturation at 94° for 30 sec, primer annealing at 50°C for 30 sec, and extension at 72°C for 1 min 30 sec, with a final extension phase at 72°C for 10 min. The PCR-products were purified with ExoSAP-IT (Exonuclease I and Shrimp Alkaline Phosphatase in buffer; Amersham Biosciences) prior to being added as template for chain termination sequencing following the protocol described in ref. [64]. DNA fragments were purified with SephadexTM G-50 (Amersham Biosciences) following the manufacturer's instruction and subsequently visualized on an ABI 3100 capillary sequencer (Applied Biosystems). All sequences are available from GenBank under the accession numbers listed in Table 2.

            Phylogenetic analysis

            Alignment of the DNA sequences was performed using Clustal X [65] and improved by eye for the control region. Since all species were resolved as respectively monophyletic clades in a preliminary NJ tree (tree not shown) using all available sequences of the mitochondrial control region (calculated in PAUP* 4.0b10; [66]), we chose 52 taxa (including outgroup taxa) for amplification of the ND6 fragment. Several individuals had long insertions (> 100 bp) in the control region that were excluded for phylogenetic analyses. The combined dataset of control region and ND6, including 46 taxa, was used to reconstruct the phylogenetic relationships among the investigated Synodontis species. We applied neighbor-joining (NJ), maximum parsimony (MP) and maximum likelihood (ML) analyses in parallel, using the computer program PAUP*. Phylogenetic relationships were also reconstructed by Bayesian inference (BI) using MrBayes 3.1 [67].

            To assess the degree of saturation of transition (ti) and transversion (tv) mutations at each codon position of the ND6 gene, we plotted the number of mutations against the percentage of sequence divergence of all pairwise comparisons (not shown). Based on the estimated ti/tv ratio inferred from these pairwise comparisons we derived a proper weighting scheme for a weighted MP analysis. Due to the estimated ti/tv ratio of 6.81 for third codon positions of two-and threefold degenerate amino acids, 2.63 for third codon positions of fourfold degenerate amino acids, 2.00 for second codon positions, and 1.75 for first codon positions we applied the following weighting scheme (tv/ti): 70:10 for third codon positions of two- and threefold degenerate amino acids, 26:10 for third codon positions of fourfold degenerate amino acids, 70:35 for second codon positions, and 70:28 at first codon positions. C/T substitutions at the first codon position of leucine were treated as a fifth base and were down-weighted to the same weight as transitions at the third codon positions. To achieve a proper weighting scheme for the control region, we first performed a sliding window analysis to define regions exhibiting different degrees of variation (<10%, low variable; 10–20%, high variable; >20%, hyper variable) [68]. Despite considerable saturation, transition mutations contain substantial information necessary to resolve the relationships of evolutionary young splits, without obscuring the deep splits when they are down-weighted according to their frequency in relation to transversion mutations [22]. We thus evaluated the degree of saturation of ti and tv for each partition of the control region separately, and based upon the estimated ti/tv ratios of 1.40, 2.43 and 2.28 for the low, high and hyper variable regions we obtained following weighting schemes (tv/ti): low variable, 70:56; high and hyper variable, 70:28. To evaluate an appropriate substitution model for the ML analysis, we calculated hierarchical likelihood ratio test statistics using the program Modeltest 3.06 [69]. The model recovered was K81uf+I+Γ [70] with nucleotide frequencies A, 0.3581, C, 0.2676, G, 0.1368, T, 0.2375, proportion of invariable sites (I), 0.2461, gamma shape parameter (α), 0.7306, and R-matrix A↔C=G↔T+1; A↔T+C↔G=1.2463; and A↔G=C↔T=7.7875. For obtaining MP and ML topologies we applied heuristic search procedures with random addition of taxa and TBR branch swapping (1,000 replicates for MP; 100 replicates for ML). As standard measures of confidence we applied bootstrapping (1,000 pseudo-replicates for NJ and MP, 100 for ML). For the reconstruction of phylogenetic relationships by Bayesian inference analyses our data were partitioned by gene and codon position (only in ND6). Rate heterogeneity was set according to a gamma distribution with six rate categories (GTR model) for each data partition. Posterior probabilities were obtained from a 2,000,000-generation Metropolis-coupled Markov chain Monte Carlo simulation (four chains; chain temperature, 0.2; trees sampled every 100 generations), with parameters estimated from each data partition. A 50% majority-rule consensus tree was constructed after a burn-in of 1%, which is when likelihood values reached stationarity.

            To obtain a potentially better insight into the phylogenetic relationships among the East African member of the genus Synodontis we repeated the NJ, MP, ML and BI analyses with a reduced dataset including the species S. dhonti, S. granulosus, S. multipunctatus, S. nigromaculatus, S. njassae, S. petricola, S. polli, S. sp. nov., S. victoriae and S. zambezensis, using S. eupterus, S. velifer and S. ocellatus as outgroup taxa, based on the first round of phylogenetic analyses. Proper weighting schemes for MP analysis were derived from estimated ti/tv-ratios as described above. Since no tv were observed for third codon positions of two-and threefold degenerate amino acids, we applied the same weighting scheme as in the overall data set for these codon positions. Based on the estimated ti/tv ratio of 2.49 for third codon positions of fourfold degenerate amino acids, 2.28 for second codon positions, and 1.19 for first codon positions we applied following weighting scheme (tv/ti): 25:10 for third codon positions of fourfold degenerate amino acids, 70:31 for second codon positions, and 70:59 at first codon positions. C/T substitutions at the first codon position of leucine were again treated as a fifth base and down-weighted to the same weight as transitions at the third codon positions. For the control region we obtained ti/tv ratios of 2.28, 3.20 and 2.18 for low, high and hyper variable regions resulting in a weighting (tv/ti) of 70:31, 70:22 and 70:32, respectively. For NJ, ML and BI, the best fit substitution model, recovered by the hierarchical likelihood ratio tests implemented in Modeltest 3.06 [69], was HKY+I+Γ [71] with base frequencies A, 0.3565; C, 0.2814; G, 0.1308; T, 0.2313; ti/tv, 4.7844; I, 0.3046; and α, 0.5463. A four-cluster likelihood mapping analysis [72] using the program TREE-PUZZLE 5.1 [73] was applied to test the hypothesis of the monophyly of the two endemic Tanganyikan clades.

            To determine differences in substitution rates between East African representatives of the genus Synodontis and to test whether the cuckoo catfish (S. multipunctatus) indeed shows an accelerated rate of molecular evolution as indicated by its long branch, we conducted a relative rate test as implemented in RRTree [74]. We therefore tested S. multipunctatus against the remaining East African taxa, applying the K2P substitution model [75].

            We used r8s [32] to tentatively date the major cladogenetic events in Synodontis in East Africa. We used the maximum likelihood trees and applied the local molecular clock method with an optimization via the truncated newton method [76]. The analyses were performed in three cycles with different calibration points, and 30 bootstrap replicates for each cycle (the files for this bootstrapping approach were generated with Mesquite [77]; the maximum likelihood trees were generated with PAUP* [66] as described above). In the first cycle, we used all three available calibration points: (i) the maximum age of about one million years for the lacustrine radiation of the Lake Malawi species S. njassae, based on the age of the truly lacustrine habitat [49] (see also refs. [9, 56]); (ii) the maximum age of about six million years for the lacustrine radiation of the second Tanganyikan clade including S. dhonti, S. petricola, S. polli and S. sp. nov., based on the age of the truly lacustrine habitat in Lake Tanganyika [4648] (see also refs. [9, 22, 64]) (Note that we did not use the first Tanganyikan clade because of the apparently deviating evolutionary rate in that clade and the possibility that S. granulosus might represent a separate lineage that existed before the formation of the lacustrine habitat in the Lake Tanganyika basin); and (iii) the minimum age of about 20 million years for the most recent common ancestor (MRCA) of the East African clade of Synodontis, based on the oldest Synodontis fossils in that area [4145]. We note that these East African Synodontis fossils from the Early Miocene have not been assigned to extant species (see e.g., ref. [45]), which is difficult due to the diversity of the family and the fact that the criteria used for the identification of extant species are not applicable to the fossil specimens [78]. However, because these Synodontis fossils resemble extant East African species, and because of the rich fossil record since the Miocene, it is highly likely that these specimens belong to the East African clade. The second and third cycle of analyses aimed to corroborate internal calibrations. We, therefore, ran an analysis with a single calibration point at 20 million years for the minimum age of the East African clade, and an analysis, in which we constrained the maximum age for S. njassae to one million year and the minimum age for the East African clade to 20 million years.

            Abbreviations

            BI: 

            Bayesian inference

            ML: 

            maximum likelihood

            MP: 

            maximum parsimony

            MRCA: 

            most recent common ancestor

            mtDNA: 

            mitochondrial DNA

            MY: 

            million years

            MYA: 

            million years ago

            ND6: 

            NADH Dehydrogenase Subunit 6

            NJ: 

            neighbor joining

            ti: 

            transition mutation

            tv: 

            transversion mutation.

            Declarations

            Acknowledgements

            We thank C. Kapasa, P. Ngalande, H. Phiri, L. Makasa, R. Shapola, D. Sinyinza, and the team at the Mpulungu Station of the Ministry of Agriculture and Cooperatives, Republic of Zambia as well as L. Mumba, G. Mutenda, and C. Katongo, from the University of Zambia in Lusaka for their cooperation during fieldwork. We are also grateful to C. Danady, C. Katongo, W. Köhldorfer, R. Sinyinza, J. Snoeks, and T. Veall for providing samples, to E. Hespeler and J. Hauck for technical assistance in the laboratory of AM, and to K. Sefc for valuable comments on the manuscript. This study was financed by the Austrian Science Foundation (grant P17968) to CS, the Deutsche Forschungsgemeinschaft (DFG) to AM and WS. SK was further supported by the University of Graz and the Austrian Academy of Sciences (DOC-fellowship).

            Authors’ Affiliations

            (1)
            Department of Zoology, Karl-Franzens-University Graz
            (2)
            Vertebrate Department, Royal Belgian Institute of Natural Sciences
            (3)
            Lehrstuhl für Zoologie und Evolutionsbiologie, Department of Biology, University of Konstanz

            References

            1. Fryer G, Iles TD: The cichlid fishes of the Great Lakes of Africa: Their biology and Evolution. Edinburgh: Oliver & Boyd 1972.
            2. Coulter GW, ed: Lake Tanganyika and its Life. London: Oxford University Press 1991.
            3. Greenwood PH: Morphology, endemism and speciation in African cichlid fishes. Verhandlungen der Deutschen Zoologischen Gesellschaft 1973, 66: 115–124.
            4. Boulenger G: Catalogue of the fresh-water fishes of Africa in the Brithish Museum (Natural History). London: British Museum for Natural History 1898.
            5. Stiassny MLJ, Meyer A: Cichlids of the Rift Lakes. Scientific American 1999, 280: 64–69.View Article
            6. Kocher TD: Adaptive evolution and explosive speciation: the cichlid fish model. Nature Rev Genet 2004, 5: 288–298.View ArticlePubMed
            7. Snoeks J, Rüber L, Verheyen E: The Tanganyika problem: comments on the taxonomy and distribution patterns of its cichlid fauna. Speciation in Ancient Lakes (Edited by: Martens K, Goddeeris G, Coulter GW). Stuttgart: Schweizerbart'sche Verlagsdruckerei 1994, 355–372.
            8. Turner GF, Seehausen O, Knight ME, Allender CJ, Robinson RL: How many species of cichlid fishes are there in African lakes? Mol Ecol 2001, 10: 793–806.View ArticlePubMed
            9. Salzburger W, Mack T, Verheyen E, Meyer A: Out of Tanganyika: Genesis, explosive speciation, key-innovations and phylogeography of the haplochromine cichlid fishes. BMC Evolutionary Biology 2005, 5: 17.View ArticlePubMed
            10. 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.View ArticlePubMed
            11. 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.View ArticlePubMed
            12. 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.View ArticlePubMed
            13. 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.View ArticlePubMed
            14. Kornfield I, Smith PF: African Cichlid Fishes: Model systems for evolutionary biology. Annu Rev Ecol Syst 2000, 31: 163–196.View Article
            15. Poll M: Revision de la faune ichthyologique du Lac Tanganika. Annales du Musée r. du Congo belge, Bruxelles, Zoologie, Series 1 1946, 4: 141–364.
            16. Coulter GW: Composition of the flora and fauna. Lake Tanganyika and its life (Edited by: Coulter GW). New York: Oxford University Press 1991, 200–274.
            17. Brichard P: Cichlids and all other fishes of Lake Tanganyika. Neptune City, NJ: T H F Publications 1989.
            18. Coulter GW: Zoogeography, affinities and evolution, with special regard to the fish. Lake Tanganyika and its Life (Edited by: Coulter GW). New York: Oxford University Press 1991, 275–305.
            19. Zardoya R, Vollmer DM, Craddock C, Streelman JT, Karl S, Meyer A: Evolutionary conservation of microsatellite flanking regions and their use in resolving the phylogeny of cichlid fishes (Pisces: Perciformes). Proc R Soc Lond B Biol Sci 1996, 263: 1589–98.View Article
            20. Mayer WE, Tichy H, Klein J: Phylogeny of African cichlid fishes as revealed by molecular markers. Heredity 1998, 80: 702–714.View ArticlePubMed
            21. 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.PubMed
            22. Koblmuller S, Duftner N, Katongo C, Phiri H, Sturmbauer C: Ancient divergence in bathypelagic Lake Tanganyika deepwater cichlids: mitochondrial phylogeny of the tribe Bathybatini. J Mol Evol 2005, 60: 297–314.View ArticlePubMed
            23. Agnese JF, Teugels GG: Insight into the phylogeny of African Clariidae (Teleostei, Siluriformes): implications for their body shape evolution, biogeography, and taxonomy. Mol Phylogenet Evol 2005, 36: 546–53.View ArticlePubMed
            24. Poll M: Révision des Synodontis africains (famille Mochokidae). Mus R Afr Centr Tervuren (Belg) Zool 1971, 191: 1–497.
            25. Sato T: A brood parasitic catfish of mouthbrooding cichlid fishes in Lake Tanganyika. Nature 1986, 323: 58–59.View ArticlePubMed
            26. Winemiller KO, Kelso LC-Winemiller: Comparative ecology of catfishes of the Upper Zambezi River floodplain. J Fish Biol 1996, 49: 1043–1061.View Article
            27. Sanyanga RA: Food composition and selectivity of Synodontis zambezensis (Pisces: Mochokidae) in Lake Kariba, and the ecological implications. Hydrobiologia 1998, 361: 89–99.View Article
            28. Meyer A: Evolution of mitochondrial DNA in fishes. Biochemistry and molecular biology of fishes (Edited by: Hochachka PW, Mommsen TP). Amsterdam: Elsevier Science Publishers 1993, 2: 1–33.
            29. Clabaut C, Salzburger W, Meyer A: Comparative phylogenetic analyses of the adaptive radiation of Lake Tanganyika cichlid fish: nuclear sequences are less homoplasious but also less informative than mitochondrial DNA. J Mol Evol 2005, 61: 666–81.View ArticlePubMed
            30. Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol 1999, 16: 1114–1116.
            31. 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.View Article
            32. Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics 2003, 19: 301–2.View ArticlePubMed
            33. Bruwer EE, van der Bank FH, Steyn GJ, Wink M: Phylogenetic relationships between Synodontis species: some preliminary results. Afr J of Aquat Sci 2000, 25: 89–92.View Article
            34. Sturmbauer C, Koblmuller S, Sefc KM, Duftner N: Phylogeographic history of the genus Tropheus , a lineage of rock-dwelling cichlid fishes endemic to Lake Tanganyika. Hydrobiologia 2005, 542: 335–366.View Article
            35. Joyce DA, Lunt DH, Bills R, Turner GF, Katongo C, Duftner N, Sturmbauer C, Seehausen O: An extant cichlid fish radiation emerged in an extinct Pleistocene lake. Nature 2005, 435: 90–5.View ArticlePubMed
            36. 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.View ArticlePubMed
            37. Bruwer EE, van der Bank FH: Biochemical genetic markers (hybrizymes) to identify Synodontis zambezensis x S. nigromaculatus hybrids (Teleostei: Mochokidae). Biochem Syst Ecol 2002, 30: 641–649.View Article
            38. Seehausen O: Hybridization and adaptive radiation. Trends in Ecology and Evolution 2004, 19: 198–207.View ArticlePubMed
            39. Page RDM, Holmes EC: Molecular evolution: a phylogenetic approach. Oxford: Blackwell Sciences 1998.
            40. Gillespie JH: The causes of Molecular Evolution. Oxford: Oxford University Press 1991.
            41. Greenwood PH: Fish remains from Miocene deposits of Rusinga Island and Kavirondo Province, Kenya. Annals and Magazine of Natural History 1951, 12: 1192–1201.
            42. Greenwood PH: Review of cenozoic freshwater fish faunas in Africa. Annals of The Geological Survey of Egypt 1974, 4: 221–232.
            43. Van Couvering JAH: Early records of freshwater fishes in Africa. Copeia 1977, 1977: 163–166.View Article
            44. Van Neer W: New Late Tertiary fish fossils from the Sinda region, Eastern Zaire. African Study Monographs 1992, (Suppl 17) : 27–47.
            45. Stewart KM: The freshwater fish of Neogene Africa (Miocene-Pleistocene): systematics and biogeography. Fish and Fisheries 2001, 2: 177–230.View Article
            46. Tiercelin JJ, Mondeguer A: The geology of the Tanganyika Trough. Lake Tanganyika and its Life (Edited by: Coulter GW). New York: Oxford University Press 1991, 7–48.
            47. 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 African Rift system, deduced from multichannel and high-resolution reflection seismic data and piston core evidence. Basin Research 1996, 8: 1–28.View Article
            48. 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.View Article
            49. 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.
            50. Katongo C, Koblmüller S, Duftner N, Makasa L, Sturmbauer C: Phylogeography and speciation in the Pseudocrenilabrus philander species complex in Zambian Rivers. Hydrobiologia 2005, 542: 221–233.View Article
            51. Giddelo CS, Arndt AD, Volckaert FAM: Impact of rifting and hydrography on the genetic structure of Clarias gariepinus in eastern Africa. J Fish Biol 2002, 60: 1252–1266.View Article
            52. 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.View ArticlePubMed
            53. 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.View ArticlePubMed
            54. 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
            55. Sorenson MD, Payne RB: A single ancient origin of brood parasitism in African finches: implications for host-parasite coevolution. Evolution Int J Org Evolution 2001, 55: 2550–67.
            56. 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.PubMed
            57. De Pinna MCC: Higher-level phylogeny of Siluriformes, with a new classification of the order (Teleostei, Ostariophysi). PhD Thesis New York: The City University of New York 1993.
            58. M Hardman: The phylogenetic relationships among non-diplomystid catfishes as inferred from mitochondrial cytochrome b sequences; the search for the ictalurid sister taxon (Otophysi: Siluriformes). Mol Phylogenet Evol 2005, 37: 700–20.View ArticlePubMed
            59. C Sturmbauer, A Meyer: Mitochondrial phylogeny of the endemic mouthbrooding lineages of cichlid fishes of Lake Tanganyika, East Africa. Mol Biol Evol 1993, 10: 751–768.PubMed
            60. 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
            61. Meyer A, Morrissey JM, Schartl M: Recurrent origin of a sexually selected trait in Xiphophorus fishes inferred from a molecular phylogeny. Nature 1994, 368: 539–42.View ArticlePubMed
            62. Lee WJ, Conroy J, Howell WH, Kocher TD: Structure and evolution of teleost mitochondrial control regions. J Mol Evol 1995, 41: 54–66.View ArticlePubMed
            63. Park LK, Brainard MA, Dightman DA, Winans GA: Low levels of intraspecific variation in the mitochondrial DNA of chum salmon (Oncorhynchus keta). Mol Mar Biol Biotechnol 1993, 2: 362–70.PubMed
            64. Duftner N, Koblmuller S, Sturmbauer C: Evolutionary relationships of the Limnochromini, a tribe of benthic deepwater cichlid fish endemic to Lake Tanganyika, East Africa. J Mol Evol 2005, 60: 277–89.View ArticlePubMed
            65. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 1997, 25: 4876–82.View ArticlePubMed
            66. 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, 4.0 ed City: Sinauer 2003.
            67. Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Bioinformatics 2001, 17: 754–755.View ArticlePubMed
            68. Sturmbauer C, Meyer A: Genetic divergence, speciation and morphological stasis in a lineage of African cichlid fishes. Nature 1992, 358: 578–81.View ArticlePubMed
            69. Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics 1998, 14: 817–8.View ArticlePubMed
            70. Kimura M: Estimation of evolutionary distances between homologous nucleotide sequences. Proc Natl Acad Sci U S A 1981, 78: 454–8.View ArticlePubMed
            71. Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol 1985, 22: 160–174.View ArticlePubMed
            72. Strimmer K, von Haeseler A: Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proc Natl Acad Sci USA 1997, 94: 6815–9.View ArticlePubMed
            73. Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: a maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics 2002, 18: 502–504.View ArticlePubMed
            74. Robinson-Rechavi M, Huchon D: RRTree: relative-rate tests between groups of sequences on a phylogenetic tree. Bioinformatics 2000, 16: 296–7.View ArticlePubMed
            75. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 1980, 16: 111–20.View ArticlePubMed
            76. Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol 2002, 19: 101–109.PubMed
            77. Maddison WP, Maddison DR: Mesquite: A modular system for evolutionary analysis. [http://​mesquiteproject.​org] Book Mesquite: A modular system for evolutionary analysis, 1.06 ed. City 2005.
            78. Pinton A, Fara E, Otero O: Spine anatomy reveals the diversity of catfish through time: a case study of Synodontis (Siluriformes). Naturwissenschaften 2006, 93: 22–6.View ArticlePubMed

            Copyright

            © Koblmüller et al. 2006

            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.

            Advertisement