Introduction

With their unique species assemblages, the Great Lakes of East Africa represent ideal model systems to study evolutionary processes. Lake Tanganyika is the oldest of the three East African Great Lakes with an estimated age of 9–12 myr (Cohen et al., 1993) and considered to be one of the most diverse freshwater ecosystems in the world. It is well known for its outstandingly diverse cichlid fish fauna that is characterized by an extraordinary high degree of endemism (Fryer & Iles, 1972). The species flock is older and the degree of ecological and morphological divergence among species is more extreme than in the other two Great Lakes (Greenwood, 1984). Many species show a high degree of population differentiation in terms of color and genetics along the shoreline (Sturmbauer & Meyer, 1992; Duftner et al., 2006; Sefc et al., 2007; Koblmüller et al., 2009). Thus, Lake Tanganyika is an ideal ecosystem to study genetic and phenotypic divergence and the formation of locally differentiated populations along its shoreline which in some parts is highly structured into rocky and sandy stretches. This is especially true for species that are restricted to a special type of habitat which is distributed discontinuously along the lake shores (Sturmbauer & Meyer, 1992; Kornfield & Smith, 2000; Duftner et al., 2006; Sefc et al., 2007).

The tribe Lamprologini is the most species-rich tribe of Lake Tanganyika comprising about 90 species in seven genera (Poll, 1956, 1986; Colombe & Allgayer, 1985; Stiassny, 1991; Schelly et al., 2006; Koblmüller et al., 2007; Day et al., 2007). This substrate-spawning lineage is monophyletic in Poll’s classification (Poll, 1986) and contains both endemic and non-endemic species (Sturmbauer et al., 1994). Non-endemic species have, besides inhabiting Lake Tanganyika itself, also colonized the Zaire River system (Sturmbauer et al., 1994) and the Malagarazi River (De Vos & Seegers, 1998). Outstanding diversity in morphology, behavior, and ecology makes members of this tribe particularly interesting when studying pathways of adaptive radiation (Sturmbauer et al., 1994). The genus Altolamprologus comprises two described species, Altolamprologus compressiceps and Altolamprologus calvus (Boulenger 1898), and possibly at least one gastropod shell breeding dwarf species (Konings, 1998). A. compressiceps is distributed all over the entire lake and lives in sympatry with A. calvus at the south-western part of the lake (Konings, 1998).

A mitochondrial phylogeny of the genus Altolamprologus showed that both described species are of quite recent origin and started to spread and diverge into color morphs very recently (Ernst, 2008). Furthermore, A. calvus is considered not to be completely genetically separated from A. compressiceps in the mtDNA phylogeny, in that lineage sorting toward populations of A. calvus is incomplete (Ernst, 2008). As the two entities are stable despite of sympatry, one can assume complete reproductive isolation. However, current data do not allow us to be certain that no genes are exchanged at very low frequency. A. compressiceps and A. calvus have a highly specialized morphology (Brichard, 1999), characterized by a very slender and high body and a tapered mouth. Jaws can be protruded for suction feeding and its teeth are forceps-like. In addition, eyes can be turned 90° out of their usual position. The genus Altolamprologus feeds specifically on shrimp of the family Atyidae (Yuma, 1994) and it thus seems highly likely that the peculiarly depressed body is a particular adaptation to this mode of foraging, since a laterally depressed body allows easier access to narrow rock crevices. Both species are restricted to rock- or cobble habitats, and in some cases gastropod shell beds.

It has been argued that distinct trophic adaptation and the strong bond to a specific type of habitat can facilitate population differentiation across ecological barriers and eventually lead to speciation (Sturmbauer, 1998). A. compressiceps is such a stenotopic and strongly habitat-bound species and thus is expected to exhibit a relatively high degree of differentiation among populations. Indeed, several color variants were found across the entire lake (Konings, 1998). Given its relatively young evolutionary age, a rather rapid local divergence is likely to have happened in this species.

As population structure can result in restricted gene flow, color variation and local adaptation in particular morphological features that indicate a rapid population divergence, population studies need to address both genetic and phenotypic divergence. Geometric morphometrics, a landmark-based approach for investigating body shape changes, has been shown to be able to detect fine-scale morphological differences in overall body shape. Furthermore, this technique has the potential to discriminate among assemblages with recent common ancestry (Adams et al., 2004; Maderbacher et al., 2008). Even populations living in only marginally different environmental conditions can exhibit morphological differences which can be detected by geometric morphometrics. In combination with neutral genetic markers, changes in the phenotype of an organism can be used to assess the intensity and direction of natural selection (Raeymaekers et al., 2007).

This study investigates to which extent morphological and genetic differentiation occurred among three geographically separated populations of the highly specialized shrimp-feeding rock-dweller A. compressiceps in the South of Lake Tanganyika.

Materials and methods

Sampling

119 individuals of A. compressiceps (54 males, 46 females, 20 unknown) were sampled from three different sites at the southern tip of Lake Tanganyika in Zambia (Mwese n = 39, [08°34′S, 30°47′E], Kasakalawe n = 40 [08°47′S, 31°04′E], Wonzye n = 40 [08°43′S, 31°08′E], see Fig. 1). Fin clips were taken from all individuals and preserved in 99% ethanol for DNA extraction. Whole specimens were stored in 10% formalin for geometric morphometric analysis.

Fig. 1
figure 1

Geographical location of Lake Tanganyika in East Africa, map of Lake Tanganyika with details of southern shore including sampling locations (Mwese, Kasakalawe and Wonzye)

Full size image

DNA extraction and microsatellite amplification

Whole genomic DNA was extracted from 119 specimens using a high-salt, ammonium acetate extraction protocol after a proteinase K digestion (modified after Bruford et al., 1998). Eleven polymorphic microsatellite loci developed for other cichlid species were screened for use in A. compressiceps. Five loci gave consistent results and exhibited enough variation to be used for population genetic analyses. Microsatellites were PCR-amplified in a 10 or 20 μl reaction mix containing 5–10 ng extracted DNA, 0.25 pM of each primer, 0.05 mM dNTP mix, 10 μg/μl BSA, 1.5 mM MgCl2, 10× Taq polymerase buffer and 0.5 U Taq DNA Polymerase. Forward primers were labeled with three different fluorescent dyes FAM, HEX, and NED. The PCR protocol started with an initial denaturation step at 94°C for 3 min followed by 40 cycles at 94°C for 30 s, annealing at 50–54°C for 30 s (specific annealing temperatures see Table 1) and extension at 72°C for 1 min followed by a final extension at 72°C for 7 min. PCR products were analyzed on an ABI 3130xl sequencer together with an ABI ROX500 size standard and allele calling was performed with GeneMapper software version 3.7.

Table 1 Characterization of five microsatellite loci: size range, number of alleles (N A), observed (H O) and expected (H E) heterozygosity and species in which the microsatellite locus was developed
Full size table

Data analysis

To detect scoring errors and genotyping artifacts such as null alleles, large allele dropout and stutter bands we used MICRO-CHECKER (van Oosterhout et al., 2004). Parameters of genetic diversity as number of alleles (N A), observed (H O) and expected (H E) heterozygosity and departures from Hardy–Weinberg equilibrium were calculated in Arlequin ver.3.11 (Excoffier et al., 2005). Arlequin was also used to quantify pairwise population differentiation calculating FST based on allele frequencies (assuming the infinite allele model IAM, Weir & Cockerham, 1984) and RST based on allele sizes assuming the stepwise mutation model (SMM, Slatkin, 1995). We also tested for isolation-by-distance performing a Mantel test.

Factorial correspondence analysis (FCA), a multivariable analysis to visualize genetic distances between genotypes, was performed in GENETIX ver.4.05.2 using the AFC3D procedure (Belkhir et al., 1996–2004). Factorial axes are ordered according to their eigenvalues, larger values explaining more of the general structure in the dataset. Using the software STRUCTURE version 2.2 (Pritchard et al., 2000), we conducted a population assignment test using Bayesian clustering based on allele frequencies to calculate the most probable number of populations (k) and to assign individuals to their population of origin. For each k (ranging form k = 1 to k = 3) we computed three iterations with a burn-in of 100,000 and 100,000 MCMC replicates using the admixture model (ancestry model) which assumes some degree of gene flow so that it may be possible that specimens stem from more than a single population. Furthermore, we chose the model of correlated allele frequencies which implies that all populations once originated from a single ancestral population and that observed genetic divergence is the result of genetic drift accumulating since the divergence (Martien et al., 2007). In addition, a recently implemented model for populations with low genetic divergence was used (Hubisz et al., 2009). For each value of k, we selected the MCMC run with the lowest estimated posterior probability (Ln P(D)).

To investigate the historical population demography (population expansion), meta-population analyses were conducted. This was performed using the software BOTTLENECK ver.1.2.02 (Cornuet & Luikard, 1996) for detecting recent bottleneck events on the basis of allele frequency data. Populations that recently experienced a dramatic decline in effective population size are characterized by a significant excess in heterozygotes. We performed a Wilcoxon sign rank test to calculate probabilities for heterozygote excess in all three populations for three mutation models (IAM, TPM and SMM). To estimate population sizes and past migration rates among the three studied populations, we used MIGRATE ver.3.0.3 (Beerli & Felsenstein, 1999, 2001) applying a maximum likelihood inference. We used the stepwise mutation model for microsatellite markers and applied a MCMC approach, where we calculated 10 short chains (500 used trees out of sampled 10,000) and 3 long chains (5,000 used trees out of sampled 100,000). A migration matrix model with unequal population sizes and different migration rates was assumed (Beerli & Felsenstein, 1999, 2001).

Geometric morphometrics

113 specimens of A. compressiceps, which were also used for genetic analysis, were subjected to geometric morphometric analysis. The left body side of each preserved fish was scanned with a flatbed-scanner (Herler et al., 2007). On the basis of the landmark system introduced by Maderbacher et al. (2008) we developed a set of 17 homologous landmarks for A. compressiceps. Those landmarks were digitized for each specimen (Fig. 2) using tpsDIG (SUNY at Stony Brook, Rohlf 2003) to obtain coordinates for all 17 landmarks in all 113 specimens. Analysis of landmark data was carried out with IMP software (Sheets, 2003). Procrustes analysis was implemented to obtain shape information that is independent from position (orientation and rotation), size and location. For this study, partial procrustes superimposition was performed in CoordGen6f (IMP package, Sheets, 2003).

Fig. 2
figure 2

Landmark positions for geometric morphometric analysis. (1) Anterior tip of snout (2) groove at basis of forehead (3, 4) anterior and posterior insertion of dorsal fin (5, 7) upper and lower insertion of caudal fin (6) midpoint of origin of caudal fin (8, 9) posterior and anterior insertion of anal fin (10) insertion of ventral fin (11, 12) lower and upper insertion of pelvic fin (13) most posterior point of operculum (14, 15) extreme of orbit along the anteroposterior body axis (16) most posterior point of chin (17) most posterior point of mouth cleavage

Full size image

To evaluate the degree of variation among the three studied populations, a principal component analysis (PCA) was carried out in PCAGen6n (IMP package, Sheets, 2003) which creates a set of new independent variables that best describe the variance among studied populations. To study the overlap among populations a Canonical Variate Analysis (CVA) was carried out in CVAGen6j (IMP package, Sheets, 2003), that calculates partial warp scores to perform a CVA after a MANOVA (multivariate analysis of variance). Principal components maximize the variance of individual scores whereas canonical variates maximize the variance between the group means relative to the variance within the groups, thus leading to an optimal linear discrimination between the groups. For pairwise population comparison, deformation grids as well as Goodall’s F-test were calculated in TwoGroup6h (IMP package, Sheets, 2003).

Results

Genetic population differentiation

A total of 94 alleles were detected in five microsatellite loci. All five loci were highly polymorphic and all three populations displayed a high degree of genetic diversity (Mwese 0.780; Kasakalawe 0.783; Wonzye 0.820). Numbers of alleles per locus ranged from 11 (Pzeb1) to 28 (Pmv17) and populations were characterized by a total numbers of alleles as follows, 64 (Mwese), 66 (Kasakalawe), and 66 (Wonzye). A detailed characterization of all five loci (size range, number of alleles, observed and expected heterozygosity per locus) is given in Table 1. Parameters of genetic diversity in the three studied populations (number of alleles, observed and expected heterozygosity per locus and population) can be found in Table 2. One locus, Pzeb1, showed significant deviation from Hardy–Weinberg equilibrium in all three populations exhibiting an excess of heterozygotes. No evidence of null alleles or scoring errors was detected in Pzeb1 or any other locus. After excluding Pzeb1 from further analyses, F ST values still remained highly significant (P < 0.001).

Table 2 Parameters of genetic diversity in three populations
Full size table

Factorial correspondence analysis of microsatellite data showed a tendency of population separation with some individuals overlapping (Fig. 3). The first three factorial axes explain 11% of the total genetic variance. Eigenvalues of the axes are 0.34, 0.31, and 0.30. Wonzye was separated from Mwese through factorial component 1 (3.8%) whereas Kasakalawe was separated from Wonzye and Mwese through factorial component 2 (3.5%).

Fig. 3
figure 3

Results of the factorial correspondence analysis (FCA) of microsatellite data

Full size image

Pairwise population differentiation between all three populations based on allele frequencies (F ST) showed relatively low values (0.017–0.027) but all were highly significant (Table 3). Population differentiation based on allele sizes (R ST) also exhibited rather low values (0.006–0.028) with two population comparisons showing significant differentiation (Mwese–Wonzye and Kasakalawe–Wonzye). Population comparison between Mwese and Kasakalawe could not address significant differentiation (Table 3). Genetic differentiation was highest between Kasakalawe and Wonzye despite the smallest geographical distance. Therefore, a Mantel test did not show any significant correlations between geographic and genetic distances among populations that indicates no isolation by distance. Population assignment tests using a Bayesian clustering method based on allele frequencies conducted with STRUCTURE could not infer the right number of populations and were not able to assign individuals to their population of origin. Even using a new, recently implemented model in STRUCTURE which compensates low levels of genetic divergence (Hubisz et al., 2009), no population structure was detected most likely due to very low F ST values obtained from pairwise population comparisons.

Table 3 Pairwise population differentiation among three populations: F ST/R ST values below diagonal, Goodall’s F values above diagonal
Full size table

No recent genetic bottleneck was inferred by BOTTLENECK as no significant heterozygosity excess was detected except for Mwese under the infinite allele model. Estimations of population sizes and past migration rates among the three studied populations performed in MIGRATE (Beerli & Felsenstein, 1999, 2001) showed that Theta, describing population sizes, ranged from 7.11 (Mwese) to 10.65 (Wonzye). Migration rates (M) were almost the same among all populations but Mwese showed slightly higher immigration rates than Kasakalawe and Wonzye (Table 4). Immigration rates from Kasakalawe and Wonzye to Mwese were 2.0 and 1.9 individuals per generation, respectively.

Table 4 Results of the migration rate analysis: Pairwise population comparisons of population sizes (Theta) and migration rates (M) among the three studied populations
Full size table

Geometric morphometrics

Principal component analysis (PCA) could not clearly separate distinct clusters of populations but nevertheless showed a tendency of population differentiation (Fig. 4). Specimens from Mwese and Kasakalawe were partly overlapping on Principal Component 2 (describing 9.8% of total variance). Specimens from Wonzye completely overlapped with the other two populations. Principal Component 1 described 65.5% of total variance but along this axis there is no population separation observable.

Fig. 4
figure 4

Results of the principal component analysis (PCA)

Full size image

Canonical variate analysis (CVA) of 17 landmarks showed a stronger tendency of population separation because of a better linear discrimination of groups (Fig. 5). Some individuals were overlapping which was observed more frequently between specimens from Mwese and Wonzye than between other populations. Mwese was separated from Wonzye through Canonical Variate 1 (Eigenvalue 2.1196) whereas Kasakalawe showed separation from Mwese and Wonzye through Canonical Variate 2 (Eigenvalue 1.0361). Pairwise population differentiation based on shape change (Goodall’s F test) showed highly significant results. Highest morphological differentiation was found between Mwese and Kasakalawe, whereas F ST scores obtained from comparisons between Mwese–Wonzye and Kasakalawe–Wonzye were very similar (Table 3). Deformation grids showing pairwise population comparisons indicate that shape changes mainly occurred in the head region of A. compressiceps in all three population comparisons (Fig. 6a–c).

Fig. 5
figure 5

Results of the canonical variate analysis (CVA) of 17 landmarks

Full size image
Fig. 6
figure 6

Deformation grids: a Kasakalawe and Mwese, b Mwese and Wonzye, and c Kasakalawe and Wonzye. Arrow size and direction refer to the amount of divergence and the trend of morphological differentiation between the two compared populations

Full size image

Discussion

Genetic population differentiation

Species exhibiting intraspecific variation provide ideal opportunities to study evolutionary processes responsible for diversification. Natural selection, sexual selection, genetic drift, or a combination of these factors, may act as diversifying forces. In this respect, coloration has received much attention as it responds to natural and sexual selection (Endler, 1983; Seehausen & van Alphen, 1999), but may also evolve through drift alone (Hoffman et al., 2006). Moreover, species-specific characteristics like site fidelity, stenotopy through specialization to particular trophic niches, and territorial organization strongly influence the degree of genetic differentiation among local populations.

Pairwise population differentiation between all three populations based on allele frequencies (F ST) showed relatively low values (1.75–2.74) but all were highly significant (Table 3). However, the observed values are well within the range among populations of other cichlid species in Lake Tanganyika. Variabilichromis moorii, a stenotopic herbivorous habitat specialist confined to rocky habitat without obvious morphological differentiation among sites showed at least twice as much genetic differentiation than A. compressiceps (genetic divergence ranging from 4.0–16.0%, Duftner et al., 2006). Tropheus moorii, a highly phylopatric species of which about 120 color morphs are described (Konings, 1998; Schupke, 2003), showed a higher degree of genetic population divergence than A. compressiceps with F ST values of 3.1 to 6.8% (Sefc et al., 2007). Higher values of population differentiation (F ST ranging from 4.7–18.8%) were found in Eretmodus cyanostictus, a species exhibiting a high degree of site philopatry (Sefc et al., 2007) but hardly any color differentiation. Moreover, the less highly philopatric and more mobile cichlid fish species Ophthalmotilapia ventralis (FST 0.9–1.4%; Sefc et al., 2007) exhibited less genetic differentiation among populations than A. compressiceps. With a level of 1.8 to 2.7% of population differentiation A. compressiceps lies well above the values for more mobile species but remains in the lower range of more philopatric species. However, ongoing gene flow among populations might also explain rather low F ST values in a stenotopic species like A. compressiceps (see also Martien et al., 2007).

R ST values were generally lower than F ST values, with two population comparisons showing significant results. As the stepwise mutation model represents a rather strict mutation model, such low values are not surprising given the overall low population differentiation among the three studied populations.

Despite its description as a habitat-bound species, analysis of past migration rates showed ongoing migration among all three studied populations even across larger geographic distances that could prevent more pronounced genetic differentiation. As little as one individual per generation can be sufficient to prevent substantial population differentiation (Wright, 1931). Taken together, the clear but low degree of population differentiation observed in A. compressiceps can best be explained by (1) its young evolutionary age and consequently recent spread along the southern shores of Lake Tanganyika and (2) by its rather stenotopic and habitat-bound life style. The question about the actual degree of ongoing gene flow among populations cannot be answered, but the existence of several color morphs which are either following the distribution of continuous rocky shore lines or rocky outcrops suggests very low levels of gene flow across habitat breakpoints. Such breakpoints are sandy stretches along the shoreline or river estuaries. They have been shown to promote population differentiation by preventing gene flow in species adapted to rocky habitats (Sefc et al., 2007).

As already mentioned, low levels of population differentiation over large geographic distances compared to other cichlids with similar life history traits indicates a relatively recent separation of A. compressiceps populations in the southern basin of Lake Tanganyika. This assumption is in accordance with results of a mitochondrial phylogeny suggesting a recent divergence of the genus Altolamprologus (Ernst, 2008) based on low levels of sequence divergence, frequently shared haplotypes and incomplete lineage sorting among A. compressiceps and A. calvus. According to a TCS haplotype network, Altolamprologus seems to have evolved from a local ancestor that underwent a second diversification. A. calvus seems to have split from a local A. compressiceps population, presumably somewhere in the southwest of the lake, and remained distributed quite locally. Based on mtDNA data, A. compressiceps and A. calvus are considered two closely related species. MtDNA data also showed that A. compressiceps first evolved in an unknown locality and subsequently dispersed across the entire lake (Ernst, 2008). Populations separated by large geographic distances (e.g., coast to coast) are expected to exhibit a higher level of genetic differentiation than adjacent populations. Even a large geographic distance between two populations in this study (about 40 km between Mwese and Wonzye) did not result in higher population differentiation than within smaller geographic ranges in A. compressiceps populations. This may possibly be due to the relative recent separation of the existing populations from a single ancestral population (Slatkin, 1993). Non-significant isolation by distance based on our microsatellite data might also be attributed to sampling at only three sites and using too few microsatellite loci. Even with an excess of heterozygotes at one locus, no bottleneck events were detected in any of the three populations.

This observation is interesting as significant demographic effects of the most recent water level fluctuations in Lake Tanganyika (32–14.000 years ago, Lezzar et al., 1996; Cohen et al., 1993, 1997; McGlue et al., 2008) have been suggested to be the major factor for observed demographic structure in other sympatric cichlid species (Verheyen et al., 1996; Sturmbauer et al., 2001). It may well be that the spread of Altolamprologus was too recent to show the same patterns as in species with pre-existing populations.

Morphological population differentiation

For the study of subtle morphological differences in overall body shape, landmark-based morphometric approaches are better able to visualize and quantify changes among populations than traditional methods (Adams et al., 2004; Maderbacher et al., 2008). In contrast to the nuclear genetic results (highest genetic differentiation between Kasakalawe and Wonzye), geometric morphometric analysis revealed highest phenotypic differentiation between Kasakalawe and Mwese (Goodall’s F of 33.24) which are separated by a geographical distance of about 40 km (Fig. 1). The deformation grid comparing these two populations (Fig. 6a) indicated that differences in morphology mostly occur in the head region. As head morphology was shown to be a strong indicator for trophic specialization, slight habitat differences may enforce slightly differential selection directions to promote local adaptation. If gene flow continues to be restricted, these adaptations concerning feeding efficiency could subsequently lead to greater differences in the trophic apparatus and might promote speciation. As A. calvus lives in sympatry in a relatively short stretch of coast line, its origin might have followed this scenario.

Maderbacher et al. (2008) studied the overall morphology of T. moorii applying traditional and geometric morphometrics, where the latter revealed significant differentiation among populations in body shape. Shape changes primarily occurred in the head region which is considered to be the most variable part in the cichlid body due to trophic specializations (Fryer & Iles, 1972). As previous studies showed (e.g., Albertson & Kocher, 2001; Albertson et al., 2003; Postl et al., 2008), the oral jaw apparatus plays a fundamental role in cichlid evolution and diversification. Postl et al. (2008) assessed the same populations of T. moorii as Maderbacher et al. (2008) using geometric morphometric methods investigating viscerocranial bone elements of the oral jaw apparatus. Results showed that most of these structures were highly conserved and did not show significant variation among populations, but the dentary, the tooth-bearing anterior part of the lower jaw, showed significant differentiation among studied populations. A. compressiceps, a feeding specialist for shrimp shows highly specialized morphological adaptations in the head region, a protrudable mouth for efficient suction feeding, as well as in overall body shape in the form of a highly compressed body (Yuma, 1994). More detailed studies investigating the oral jaw apparatus in A. compressiceps would be necessary to determine the degree of differentiation in these bony elements.

Despite the relatively young divergence in the genus Altolamprologus, significant population differentiation on a small geographic scale was found in nuclear DNA as well as overall morphology of A. compressiceps. Assuming an only recent dispersal across the southern lake basin based on mtDNA data, significant population differentiation may be due to special habitat and diet requirements, as well as restricted gene flow across ecological barriers. Although both genetic and morphological markers showed statistically significant population differentiation, levels of differentiation among populations are not consistent between markers. Overall body shape and in this case particularly head features might be affected by natural selection. In contrast, neutral genetic markers like microsatellites are considered to evolve independently from the highly complex mechanisms of natural selection so that differences in these loci can accumulate over time. Other evolutionary forces like mutation and genetic drift are more likely to be responsible for the genetic divergence among populations causing dissimilar pattern in genetic and morphological differentiation (Beebee & Rowe, 2004; Halliburton, 2004). In wild populations, P ST (Raeymaekers et al., 2007) represents an equivalent to Q ST (Spitze, 1993) quantifying the proportion of among-population phenotypic variance. Comparing neutral genetic population differentiation (F ST) to phenotypic differentiation (P ST) could help separate effects of natural selection from effects of genetic drift contributing to the process of population differentiation in A. compressiceps but a sensible analysis of F ST/P ST among wild populations would require more than the number of populations sampled for this study.

Clearly, further studies are needed to gain a better understanding of the processes driving or affecting speciation in the genus Altolamprologus. As described in Egger et al. (2007), nuclear and mitochondrial analyses may reveal different processes contributing to the formation of species. Additional nuclear genetic marker systems like microsatellites or AFLPs comparing sympatric and allopatric A. compressiceps and A. calvus over the entire distribution range should be the next step to shed further light on the evolutionary relationships among species and populations.