Habitat use and its implications to functional morphology: niche partitioning and the evolution of locomotory morphology in Lake Tanganyikan cichlids (Perciformes: Cichlidae)

Animal locomotory morphology, i.e. morphological features involved in locomotion, is under the influence of a diverse set of ecological and behavioral factors. In teleost fish, habitat choice and foraging strategy are major determinants of locomotory morphology. In this study, we assess the influence of habitat use and foraging strategy on important locomotory traits, namely the size of the pectoral and caudal fins and the weight of the pectoral fin muscles, as applied to one of the most astonishing cases of adaptive radiation: the species flock of cichlid fishes in East African Lake Tanganyika. We also examine the course of niche partitioning along two main habitat axes, the benthic vs. limnetic and the sandy vs. rocky substrate axis. The results are then compared with available data on the cichlid adaptive radiation of neighbouring Lake Malawi. We find that pectoral fin size and muscle weight correlate with habitat use within the water column, as well as with substrate composition and foraging strategies. Niche partitioning along the benthic–limnetic axis in Lake Tanganyikan cichlids seems to follow a similar course as in Lake Malawi, while the course of habitat use with respect to substrate composition appears to differ between the cichlid assemblages of these two lakes.     

Unusually fine-scale genetic structuring found in rapidly speciating Malawi cichlid fishes

Mechanisms behind the explosive radiation of over 500 cichlid fish species from a single founding population in Lake Malawi during the last 700 000 years are poorly understood. Recent studies have suggested that the degree of population subdivision among the habitat patches within the lake may be responsible, but the evidence has been circumstantial: lack of a dispersal stage in haplochromine cichlids; genetic and colour variation among populations separated by large-scale geographical barriers; and fluctuating lake levels. One reason for the rapidity of speciation in these fishes may be that population subdivision is on a much finer scale than previously thought. Here we quantify the level of population subdivision and estimate migration at a scale of 700-1400 m, in order to investigate whether cichlid populations are sufficiently isolated from each other for allopatric divergence and perhaps speciation to take place. Using six microsatellite loci, we demonstrate the existence of highly significant genetic differentiation between subpopulations on adjacent headlands in each of four rock-dwelling haplochromine cichlid species. Our results suggest that these fish populations are divided into thousands of subunits among which genetic divergence is currently occurring, and that this may provide unprecedented opportunities for allopatric speciation. <br>     

Crater lake cichlids individually specialize along the benthic–limnetic axis

A common pattern of adaptive diversification in freshwater fishes is the repeated evolution of elongated open water (limnetic) species and high‐bodied shore (benthic) species from generalist ancestors. Studies on phenotype‐diet correlations have suggested that population‐wide individual specialization occurs at an early evolutionary and ecological stage of divergence and niche partitioning. This variable restricted niche use across individuals can provide the raw material for earliest stages of sympatric divergence. We investigated variation in morphology and diet as well as their correlations along the benthic‐limnetic axis in an extremely young Midas cichlid species, Amphilophus tolteca, endemic to the Nicaraguan crater lake Asososca Managua. We found that A. tolteca varied continuously in ecologically relevant traits such as body shape and lower pharyngeal jaw morphology. The correlation of these phenotypes with niche suggested that individuals are specialized along the benthic‐limnetic axis. No genetic differentiation within the crater lake was detected based on genotypes from 13 microsatellite loci. Overall, we found that individual specialization in this young crater lake species encompasses the limnetic‐ as well as the benthic macro‐habitat. Yet there is no evidence for any diversification within the species, making this a candidate system for studying what might be the early stages preceding sympatric divergence.     

Geographical ancestry of Lake Malawi's cichlid fish diversity

The Lake Malawi haplochromine cichlid flock is one of the largest vertebrate adaptive radiations. The geographical source of the radiation has been assumed to be rivers to the south and east of Lake Malawi, where extant representatives of the flock are now present. Here, we provide mitochondrial DNA evidence suggesting the sister taxon to the Lake Malawi radiation is within the Great Ruaha river in Tanzania, north of Lake Malawi. Estimates of the time of divergence between the Lake Malawi flock and this riverine sister taxon range from 2.13 to 6.76 Ma, prior to origins of the current radiation 1.20–4.06 Ma. These results are congruent with evaluations of 2–3.75 Ma fossil material that suggest past faunal connections between Lake Malawi and the Ruaha. We propose that ancestors of the Malawi radiation became isolated within the catchment during Pliocene rifting that formed both Lake Malawi and the Kipengere/Livingstone mountain range, before colonizing rivers to the south and east of the lake region and radiating within the lake basin. Identification of this sister taxon allows tests of whether standing genetic diversity has predisposed Lake Malawi cichlids to rapid speciation and adaptive radiation.     

Evolution of a cichlid fish in a Lake Malawi satellite lake

Allopatric divergence in peripheral habitats may lead to rapid evolution of populations with novel phenotypes. In this study we provide the first evidence that isolation in peripheral habitats may have played a critical role in generation of Lake Malawi's cichlid fish diversity. We show that Lake Chilingali, a satellite lake 11.5 km from the shore of Lake Malawi, contains a breeding population of Rhamphochromis, a predatory genus previously thought to be restricted to Lake Malawi and permanently connected water bodies. The Lake Chilingali population is the smallest known Rhamphochromis, has a unique male nuptial colour pattern and has significant differentiation in mitochondrial DNA from Lake Malawi species. In laboratory mate choice trials with a candidate sister population from Lake Malawi, females showed a strong tendency to mate assortatively indicating that they are incipient biological species. These data support the hypothesis that isolation and reconnection of peripheral habitats due to lake level changes have contributed to the generation of cichlid diversity within African lakes. Such cycles of habitat isolation and reconnection may also have been important in evolutionary diversification of numerous other abundant and wide-ranging aquatic organisms, such as marine fishes and invertebrates.     

Mitochondrial genome primers for Lake Malawi cichlids

Resolving the evolutionary history of rapidly diversifying lineages like the Lake Malawi Cichlid Flock demands powerful phylogenetic tools. Although this clade of over 500 species of fish likely diversified in less than two million years, the availability of extensive sequence data sets, such as complete mitochondrial genomes, could help resolve evolutionary patterns in this group. Using a large number of newly developed primers, we generated whole mitochondrial genome sequences for 14 Lake Malawi cichlids. We compared sequence divergence across protein‐coding regions of the mitochondrial genome and also compared divergence in the mitochondrial loci to divergence at two nuclear protein‐coding loci, Mitfb and Dlx2. Despite the widespread sharing of haplotypes of identical sequences at individual loci, the combined use of all protein‐coding mitochondrial loci provided a bifurcating phylogenetic hypothesis for the exemplars of major lineages within the Lake Malawi cichlid radiation. The primers presented here could have substantial utility for evolutionary analyses of mitochondrial evolution and hybridization within this diverse clade.     

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits — including body shape, colour, lips and jaws — is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high‐bodied benthic species (Amphilophus astorquii). A total of 453 genome‐wide informative SNPs were identified in 240 F₂ hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy‐two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL‐linked genomic regions, genes that might contribute to divergence in body shape along the benthic–limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.     

Genetic population structure indicates sympatric speciation of Lake Malawi pelagic cichlids.

Allopatric processes of speciation have routinely been presented to explain the extraordinary radiation of the East African Great Lakes cichlid fish species flocks. The 21 or more species of pelagic cichlids within the Lake Malawi flock appear to have lake-wide distributions that challenge such a concept. Data from six microsatellite DNA loci indicate single, panmictic populations across the lake of three Diplotaxodon species. Levels of variability at these loci suggest that populations have been large and stable. Mitochondrial DNA sequence data (872 bp of control region + 981 bp of the NADH-2) from 90 species, representing all major clades within the Lake Malawi flock, indicate reciprocal monophyly of the pelagic clade. We suggest that these data support a hypothesis that speciation in sympatry is more plausible (and widespread) within the cichlid species flocks than previously thought.     

Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.)

Parallel phenotypic divergence in replicated adaptive radiations could either result from parallel genetic divergence in response to similar divergent selection regimes or from equivalent phenotypically plastic response to the repeated occurrence of contrasting environments. In post‐glacial fish, replicated divergence in phenotypes along the benthic‐limnetic habitat axis is commonly observed. Here, we use two benthic‐limnetic species pairs of whitefish from two Swiss lakes, raised in a common garden design, with reciprocal food treatments in one species pair, to experimentally measure whether feeding efficiency on benthic prey has a genetic basis or whether it underlies phenotypic plasticity (or both). To do so, we offered experimental fish mosquito larvae, partially burried in sand, and measured multiple feeding efficiency variables. Our results reveal both, genetic divergence as well as phenotypically plastic divergence in feeding efficiency, with the phenotypically benthic species raised on benthic food being the most efficient forager on benthic prey. This indicates that both, divergent natural selection on genetically heritable traits and adaptive phenotypic plasticity, are likely important mechanisms driving phenotypic divergence in adaptive radiation.     

Ancient hybridization and phenotypic novelty within Lake Malawi's cichlid fish radiation

Does hybridization play a broad innovative role in evolution? Many studies have shown hybrid origins of individual species, particularly in major adaptive radiations, but this may be a consequence, rather than a cause, of the existence of many closely related species. Cases of hybridization in the early stages of major adaptive radiations are comparatively rare. Here, we report phylogenetic evidence for ancient introgression between distinct lineages of the species-rich Lake Malawi haplochromine cichlid fishes. Mitochondrial DNA (mtDNA) sequences indicated surprisingly close relationships between the shallow-water rocky habitat "Mbuna" species and a group of dark-adapted "Deep-Benthic" species specialized for feeding in low-light conditions (dawn/dusk, under overhangs, and deep water). By contrast, analyses of nuclear amplified fragment length polymorphism data demonstrated that these Deep-Benthic cichlids were more closely related to shallow water "Shallow-Benthic" soft-sediment feeders, a group that shares similar head and body morphology. A coalescent-based computer simulation indicated that the mtDNA similarity of rocky habitat Mbuna species and dark-adapted Deep-Benthic species was due to hybridization rather than incomplete lineage sorting. Comparisons of morphology indicated that some Deep-Benthic species possessed novel morphology not present in other Lake Malawi species groups. Thus, these analyses support the hypothesis that ancient hybridization occurred within the Lake Malawi cichlid radiation, that the event occurred before the radiation of a species group adapted to low-light benthic habitats, and that this group went on to dominate the deep-water regions of Lake Malawi. The results of this study contribute to a growing literature consistent with a creative role of hybridization in the evolution of species diversity and adaptive radiations.     

Incipient sympatric speciation in Midas cichlid fish from the youngest and one of the smallest crater lakes in Nicaragua due to differential use of the benthic and limnetic habitats?

Understanding how speciation can occur without geographic isolation remains a central objective in evolutionary biology. Generally, some form of disruptive selection and assortative mating are necessary for sympatric speciation to occur. Disruptive selection can arise from intraspecific competition for resources. If this competition leads to the differential use of habitats and variation in relevant traits is genetically determined, then assortative mating can be an automatic consequence (i.e., habitat isolation). In this study, we caught Midas cichlid fish from the limnetic (middle of the lake) and benthic (shore) habitats of Crater Lake Asososca Managua to test whether some of the necessary conditions for sympatric speciation due to intraspecific competition and habitat isolation are given. Lake As. Managua is very small (<900 m in diameter), extremely young (maximally 1245 years of age), and completely isolated. It is inhabited by, probably, only a single endemic species of Midas cichlids, Amphilophus tolteca. We found that fish from the limnetic habitat were more elongated than fish collected from the benthic habitat, as would be predicted from ecomorphological considerations. Stable isotope analyses confirmed that the former also exhibit a more limnetic lifestyle than the latter. Furthermore, split‐brood design experiments in the laboratory suggest that phenotypic plasticity is unlikely to explain much of the observed differences in body elongation that we observed in the field. Yet, neutral markers (microsatellites) did not reveal any genetic clustering in the population. Interestingly, demographic inferences based on RAD‐seq data suggest that the apparent lack of genetic differentiation at neutral markers could simply be due to a lack of time, as intraspecific competition may only have begun a few hundred generations ago.     

Recent speciation between sympatric Tanganyikan cichlid colour morphs

Lake Tanganyika, Africa's oldest lake, harbours an impressive diversity of cichlid fishes. Although diversification in its radiating groups is thought to have been initially rapid, cichlids from Lake Tanganyika show little evidence for ongoing speciation. In contrast, examples of recent divergence among sympatric colour morphs are well known in haplochromine cichlids from Lakes Malawi and Victoria. Here, we report genetic evidence for recent divergence between two sympatric Tanganyikan cichlid colour morphs. These Petrochromis morphs share mitochondrial haplotypes, yet microsatellite loci reveal that their sympatric populations form distinct genetic groups. Nuclear divergence between the two morphs is equivalent to that which arises geographically within one of the morphs over short distances and is substantially smaller than that among other sympatric species in this genus. These patterns suggest that these morphs diverged only recently, yet that barriers to gene flow exist which prevent extensive admixture despite their sympatric distribution. The morphs studied here provide an unusual example of active diversification in Lake Tanganyika's generally ancient cichlid fauna and enable comparisons of speciation processes between Lake Tanganyika and other African lakes.     

Age and spread of the haplochromine cichlid fishes in Africa

The Haplochromini are by far the most species-rich cichlid fish tribe that originated along with the so-called primary radiation of the Lake Tanganyika cichlid species flock, i.e. at the same time during which the majority of the endemic Lake Tanganyika cichlid tribes emerged. Unlike the other tribes, the haplochromines are not restricted to Lake Tanganyika but distributed throughout Africa, except for the northwestern part of the continent. Haplochromine cichlids seeded the adaptive radiation of cichlid fishes in Lakes Malawi, Kivu, Victoria, Turkana, as well as in the now extinct paleo-Lake Makgadikgadi. Here we present a comprehensive phylogenetic and phylogeographic analysis of haplochromine cichlids that is based upon DNA sequences of two mitochondrial gene segments of riverine taxa covering all major African biogeographic regions where haplochromines are found. Our analysis revealed that six lineages of haplochromines originated within a short period of time, about 5.3-4.4 MYA. These haplochromine lineages show a highly complex phylogeographic pattern, probably severely influenced by climate- and/or geology-induced changes of the environment, with river capture events most likely playing an important role for species dispersal.     

Niche divergence facilitated by fine‐scale ecological partitioning in a recent cichlid fish adaptive radiation

Ecomorphological differentiation is a key feature of adaptive radiations, with a general trend for specialization and niche expansion following divergence. Ecological opportunity afforded by invasion of a new habitat is thought to act as an ecological release, facilitating divergence, and speciation. Here, we investigate trophic adaptive morphology and ecology of an endemic clade of oreochromine cichlid fishes (Alcolapia) that radiated along a herbivorous trophic axis following colonization of an isolated lacustrine environment, and demonstrate phenotype‐environment correlation. Ecological and morphological divergence of the Alcolapia species flock are examined in a phylogenomic context, to infer ecological niche occupation within the radiation. Species divergence is observed in both ecology and morphology, supporting the importance of ecological speciation within the radiation. Comparison with an outgroup taxon reveals large‐scale ecomorphological divergence but shallow genomic differentiation within the Alcolapia adaptive radiation. Ancestral morphological reconstruction suggests lake colonization by a generalist oreochromine phenotype that diverged in Lake Natron to varied herbivorous morphologies akin to specialist herbivores in Lakes Tanganyika and Malawi.     

Signal and preference divergence among populations of the non‐endemic basal Lake Malawi cichlid fish Astatotilapia calliptera (Perciformes: Cichlidae)

In Lake Malawi and Lake Victoria, cichlid fishes have diversified into hundreds of species, many reproductively isolated by mate choice. Territorial males tend to be more aggressive to similar‐coloured males, facilitating coexistence of divergent colour morphs or species. Behavioural mate choice and aggression biases of species and allopatric populations of specialized rocky shore cichlids are influenced by divergent signals such as male colour. Believed to be basal to the Lake Malawi haplochromine radiation, and inhabiting shallow weedy areas of the lake and neighbouring water bodies, Astatotilapia calliptera also shows allopatric variation in colour. Here, it is demonstrated that such signal divergence is associated with tendencies of females to mate with males of their own population and also for males to preferentially attack males of their own population, indicating that preference divergence related to signal divergence in allopatry may have operated throughout the adaptive radiation of the Malawian cichlids. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 180–188.     

Peripheral Isolate Speciation of a Lake Malawi Cichlid Fish from Shallow Muddy Habitats

Much of the exceptional diversity of cichlid fishes in the African Great Lakes can be explained by geographic variation among isolated populations of species specialised to live on patchily distributed rocky habitat. However, there are also many endemic species that are not specialised for rocky shores. These appear to experience weaker geographic isolation. Major decreases in lake volume may have segregated such populations in isolated refugia in the distant past, but subsequent range changes have likely eliminated most of the phylogeographic signal of these events. Divergence in currently isolated peripheral water bodies may be more amenable to the study of recent processes of allopatric divergence. We investigate a haplochromine cichlid fish, here referred to as Lethrinops sp. ‘chilingali’, isolated in a small satellite lake near Lake Malawi, and the candidate sister taxon, Lethrinops lethrinus, which inhabits shallow muddy habitats in the main lake and associated water bodies. The satellite lake form from Lake Chilingali showed significant morphological differentiation, with a less ventrally-placed mouth and shorter snout, associated with a shift in diet from a diverse range of benthic invertebrates towards specialisation on mid-water chaoborid larvae and pupae. The Lake Chilingali population showed substantially reduced mitochondrial DNA diversity and no haplotype sharing was observed with populations from the main lake system. In laboratory experiments, putative species showed a high degree of assortative mating and territorial males were significantly more aggressive towards intruders of their own population. This study adds to the evidence that rapid evolution of novel phenotypes in peripheral habitats can add to the diversity of lacustrine cichlids through the evolution of at least partial reproductive isolation in allopatry.     

Experimental evidence for trait utility of gill raker number in adaptive radiation of a north temperate fish

North temperate fish in post‐glacial lakes are textbook examples for rapid parallel adaptive radiation into multiple trophic specialists within individual lakes. Speciation repeatedly proceeded along the benthic–limnetic habitat axis, and benthic–limnetic sister species diverge in the number of gill rakers. Yet, the utility of different numbers of gill rakers for consuming benthic vs. limnetic food has only very rarely been experimentally demonstrated. We bred and raised families of a benthic–limnetic species pair of whitefish under common garden conditions to test whether these species (i) show heritable differentiation in feeding efficiency on zooplankton, and (ii) whether variation in feeding efficiency is predicted by variation in gill raker numbers. We used zooplankton of three different size classes to investigate prey size dependency of divergence in feeding efficiency and to investigate the effect strength of variation in the number of gill rakers. Our results show strong interspecific differences in feeding efficiency. These differences are largest when fish were tested with the smallest zooplankton. Importantly, feeding efficiency is significantly positively correlated with the number of gill rakers when using small zooplankton, also when species identity is statistically controlled for. Our results support the hypothesis that a larger number of gill rakers are of adaptive significance for feeding on zooplankton and provide one of the first experimental demonstrations of trait utility of gill raker number when fish feed on zooplankton. These results are consistent with the suggested importance of divergent selection driven feeding adaptation during adaptive radiation of fish in post‐glacial lakes.     

Foraging environment determines the genetic architecture and evolutionary potential of trophic morphology in cichlid fishes

Phenotypic plasticity allows organisms to change their phenotype in response to shifts in the environment. While a central topic in current discussions of evolutionary potential, a comprehensive understanding of the genetic underpinnings of plasticity is lacking in systems undergoing adaptive diversification. Here, we investigate the genetic basis of phenotypic plasticity in a textbook adaptive radiation, Lake Malawi cichlid fishes. Specifically, we crossed two divergent species to generate an F₃ hybrid mapping population. At early juvenile stages, hybrid families were split and reared in alternate foraging environments that mimicked benthic/scraping or limnetic/sucking modes of feeding. These alternate treatments produced a variation in morphology that was broadly similar to the major axis of divergence among Malawi cichlids, providing support for the flexible stem theory of adaptive radiation. Next, we found that the genetic architecture of several morphological traits was highly sensitive to the environment. In particular, of 22 significant quantitative trait loci (QTL), only one was shared between the environments. In addition, we identified QTL acting across environments with alternate alleles being differentially sensitive to the environment. Thus, our data suggest that while plasticity is largely determined by loci specific to a given environment, it may also be influenced by loci operating across environments. Finally, our mapping data provide evidence for the evolution of plasticity via genetic assimilation at an important regulatory locus, ptch1. In all, our data address long‐standing discussions about the genetic basis and evolution of plasticity. They also underscore the importance of the environment in affecting developmental outcomes, genetic architectures, morphological diversity and evolutionary potential.     

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

Hybrid origin of a cichlid population in Lake Malawi: implications for genetic variation and species diversity

The importance of species recognition to taxonomic diversity among Lake Malawi cichlids has been frequently discussed. Hybridization - the apparent breakdown of species recognition - has been observed sporadically among cichlids and has been viewed as both a constructive and a destructive force with respect to species diversity. Here we provide genetic evidence of a natural hybrid cichlid population with a unique colour phenotype and elevated levels of genetic variation. We discuss the potential evolutionary consequences of interspecific hybridization in Lake Malawi cichlids and propose that the role of hybridization in generating both genetic variability and species diversity of Lake Malawi cichlids warrants further consideration.