Divergence with gene flow in the rock-dwelling cichlids of Lake Malawi

A population genetic assessment of taxonomic species: The case of Lake Malawi cichlid fishes

Organisms sampled for population‐level research are typically assigned to species by morphological criteria. However, if those criteria are limited to one sex or life stage, or the organisms come from a complex of closely related forms, the species assignments may misdirect analyses. The impact of such sampling can be assessed from the correspondence of genetic clusters, identified only from patterns of genetic variation, to the species identified using only phenotypic criteria. We undertook this protocol with the rock‐dwelling mbuna cichlids of Lake Malawi, for which species within genera are usually identified using adult male coloration patterns. Given high local endemism of male colour patterns, and considerable allele sharing among species, there persists considerable taxonomic uncertainty in these fishes. Over 700 individuals from a single transect were photographed, genotyped and separately assigned: (a) to morphospecies using photographs; and (b) to genetic clusters using five widely used methods. Overall, the correspondence between clustering methods was strong for larger clusters, but methods varied widely in estimated number of clusters. The correspondence between morphospecies and genetic clusters was also strong for larger clusters, as well as some smaller clusters for some methods. These analyses generally affirm (a) adult male‐limited sampling and (b) the taxonomic status of Lake Malawi mbuna, as the species in our study largely appear to be well‐demarcated genetic entities. More generally, our analyses highlight the challenges for clustering methods when the number of populations is unknown, especially in cases of highly uneven sample sizes.

An allelic series at pax7a is associated with colour polymorphism diversity in Lake Malawi cichlid fish

Despite long-standing interest in the evolution and maintenance of discrete phenotypic polymorphisms, the molecular genetic basis of such polymorphism in the wild is largely unknown. Female sex-associated blotched colour polymorphisms found in cichlids of Lake Malawi, East Africa, represent a highly successful polymorphic phenotype, found and maintained in four genera across the geographic expanse of the lake. Previously, we identified an association with an allelic variant of the paired-box transcription factor gene pax7a and blotched colour morphs in Lake Malawi cichlid fishes. Although a diverse range of blotched phenotypes are present in Lake Malawi cichlid species, they all appeared to result from an allele of pax7a that produces increased levels of transcript. Here, we examine the developmental and genetic basis of variation among blotched morphs. First, we confirm that pax7a-associated blotch morphs result primarily from modulation of melanophore development and survival. From laboratory crosses and natural population studies, we identify at least three alleles of pax7a associated with discrete subtypes of blotched morphs, in addition to the ancestral pax7a allele. Genotypes at pax7a support initial evolution of a novel pax7a allele to produce the blotched class of morphs, followed by subsequent evolution of that pax7a blotched allele to produce additional alleles associated with discrete colour morphs. Variant alleles of pax7a produce different levels of pax7a transcript, correlating with pigmentation phenotype at the cellular level. This naturally selected allelic series should serve as a case study for understanding the molecular genetic control of pax7a expression and the evolution of sex-associated alleles.

Multivariate heritability of shape in June sucker (Chasmistes liorus) and Utah sucker (Catostomus ardens): shape as a functional trait for discriminating closely related species

Molecular genetic data suggest that June sucker (Chasmistes liorus) is only shallowly diverged from the co-occurring but phenotypically distinct Utah sucker (Catostomus ardens) in Utah Lake. Phenotypes representing both extreme morphologies (June sucker and Utah sucker) are observed in the small wild population, but relatively large numbers of intermediate phenotypes are also present. This relatively continuous variation between the two putative species could result from extensive hybridization (including reticulate evolutionary patterns) of genetically based phenotypes or incomplete divergence among lineages and extensive phenotypic plasticity with little genetic basis. To help inform the evolutionary history of June sucker and to provide critical information for management and restoration of June sucker populations, we evaluated the distribution of shape phenotypes among purebreds of each species and their hybrids and determined the heritability of shape and ecological performance between June sucker (C. liorus) and its sister species Utah sucker (C. ardens). Mouth shape of purebred June sucker and Utah sucker are located at the extremes, and hybrids are located midway between the purebreds. Multivariate heritability was relatively high for mouth shape at 0.27. Heritability for growth rate was high at 0.32-0.42, but variation was not associated with cross type. Genetically based variation in mouth shape has evolved fairly rapidly or has been maintained in the face of ongoing hybridization between the two species. Currently, there seems to be little evidence for differential selection between species that would maintain shape variation.

Identification of cichlid fishes from Lake Malawi using computer vision

BACKGROUND

The explosively radiating evolution of cichlid fishes of Lake Malawi has yielded an amazing number of haplochromine species estimated as many as 500 to 800 with a surprising degree of diversity not only in color and stripe pattern but also in the shape of jaw and body among them. As these morphological diversities have been a central subject of adaptive speciation and taxonomic classification, such high diversity could serve as a foundation for automation of species identification of cichlids.

METHODOLOGY/PRINCIPAL FINDING

Here we demonstrate a method for automatic classification of the Lake Malawi cichlids based on computer vision and geometric morphometrics. For this end we developed a pipeline that integrates multiple image processing tools to automatically extract informative features of color and stripe patterns from a large set of photographic images of wild cichlids. The extracted information was evaluated by statistical classifiers Support Vector Machine and Random Forests. Both classifiers performed better when body shape information was added to the feature of color and stripe. Besides the coloration and stripe pattern, body shape variables boosted the accuracy of classification by about 10%. The programs were able to classify 594 live cichlid individuals belonging to 12 different classes (species and sexes) with an average accuracy of 78%, contrasting to a mere 42% success rate by human eyes. The variables that contributed most to the accuracy were body height and the hue of the most frequent color.

CONCLUSIONS

Computer vision showed a notable performance in extracting information from the color and stripe patterns of Lake Malawi cichlids although the information was not enough for errorless species identification. Our results indicate that there appears an unavoidable difficulty in automatic species identification of cichlid fishes, which may arise from short divergence times and gene flow between closely related species.

Evolutionary History of Lake Tanganyika's Predatory Deepwater Cichlids

Hybridization among littoral cichlid species in Lake Tanganyika was inferred in several molecular phylogenetic studies. The phenomenon is generally attributed to the lake level-induced shoreline and habitat changes. These allow for allopatric divergence of geographically fragmented populations alternating with locally restricted secondary contact and introgression between incompletely isolated taxa. In contrast, the deepwater habitat is characterized by weak geographic structure and a high potential for gene flow, which may explain the lower species richness of deepwater than littoral lineages. For the same reason, divergent deepwater lineages should have evolved strong intrinsic reproductive isolation already in the incipient stages of diversification, and, consequently, hybridization among established lineages should have been less frequent than in littoral lineages. We test this hypothesis in the endemic Lake Tanganyika deepwater cichlid tribe Bathybatini by comparing phylogenetic trees of Hemibates and Bathybates species obtained with nuclear multilocus AFLP data with a phylogeny based on mitochondrial sequences. Consistent with our hypothesis, largely congruent tree topologies and negative tests for introgression provided no evidence for introgressive hybridization between the deepwater taxa. Together, the nuclear and mitochondrial data established a well-supported phylogeny and suggested ecological segregation during speciation.

Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background

Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales.

Results

Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R_ST = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R_ST = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R_ST = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m ~ 0.094 - 0.097) in the Volta populations.

Conclusions

This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

A METHODOLOGY FOR THE STOCHASTIC MODELING AND SIMULATION OF SYMPATRIC SPECIATION BY SEXUAL SELECTION

In the evolution literature, sympatric speciation is the origin of two, or more, species from a single local population. Many models have been developed to study the role of ecological competition and sexual selection in sympatric speciation.

In this paper we propose a methodology for systematically deriving efficient computational models to study speciation in populations evolving with overlapping generations. As a particular case, we consider sympatric speciation by sexual selection and we follow an individual based approach: a population is represented as a set of individuals that can mate and survive according to given probabilities.

We use our methodology to construct four different models for sympatric speciation, based on male traits and female preferences. These models differ in the genotypical representation of the individuals. Results of simulations in the different models are shown and discussed.

The study of the models show that sympatric speciation by sexual selection is unlikely, also with a favorable distribution of genotypes in the initial population.

Speciation within genomic networks: a case study based on Steatocranus cichlids of the lower Congo rapids

Hybridization in animals is a much more common phenomenon as previously thought and may have profound implications for speciation research. The cichlid genus Steatocranus (Teleostei: Cichlidae), a close relative to members of the East African cichlid radiations, radiated under riverine conditions in the lower Congo rapids and produced a small species flock. Previous phylogenetic analyses suggested that hybridization occurred and contributed to speciation in this genus. A re-analysis of an already published 2000 loci-AFLP data set explicitly testing for patterns of ancient gene flow provided strong evidence for a highly reticulate phylogenetic history of the genus. We provide, to our knowledge, the first example of a complex reticulate network in vertebrates, including multiple closely related species connected through ancient as well as recent gene flow. In this context, the limited validity of strictly bifurcating tree hypotheses as a phylogenetic basis for hypothesis testing in evolutionary biology is discussed.

The macroecology of rapid evolutionary radiation

A long-standing debate in ecology addresses whether community composition is the result of stochastic factors or assembly rules. Non-random, over-dispersed patterns of species co-occurrence have commonly been attributed to competition--a particularly important force in adaptive radiation. We thus examined the macroecology of the recently radiated cichlid rock-fish assemblage in Lake Malawi, Africa at a spectrum of increasingly fine spatial scales (entire lake to depth within rock-reef sites). Along this range of spatial scales, we observed a signal of community structure (decreased co-occurrence of species) at the largest and smallest scales, but not in between. Evidence suggests that the lakewide signature of structure is driven by extreme endemism and micro-allopatric speciation, while patterns of reduced co-occurrence with depth are indicative of species interactions. We identified a 'core' set of rock-reef species, found in combination throughout the lake, whose depth profiles exhibited replicated positive and negative correlation. Our results provide insight into how ecological communities may be structured differently at distinct spatial scales, re-emphasize the importance of local species interactions in community assembly, and further elucidate the processes shaping speciation in this model adaptive radiation.

Genetic isolation and morphological divergence mediated by high-energy rapids in two cichlid genera from the lower Congo rapids

Background

It is hypothesized that one of the mechanisms promoting diversification in cichlid fishes in the African Great Lakes has been the well-documented pattern of philopatry along shoreline habitats leading to high levels of genetic isolation among populations. However lake habitats are not the only centers of cichlid biodiversity - certain African rivers also contain large numbers of narrowly endemic species. Patterns of isolation and divergence in these systems have tended to be overlooked and are not well understood.

Results

We examined genetic and morphological divergence among populations of two narrowly endemic cichlid species, Teleogramma depressum and Lamprologus tigripictilis, from a 100 km stretch of the lower Congo River using both nDNA microsatellites and mtDNA markers along with coordinate-based morphological techniques. In L. tigripictilis, the strongest genetic break was concordant with measurable phenotypic divergence but no morphological disjunction was detected for T. depressum despite significant differentiation at mtDNA and nDNA microsatellite markers.

Conclusions

The genetic markers revealed patterns of philopatry and estimates of genetic isolation that are among the highest reported for any African cichlid species over a comparable geographic scale. We hypothesize that the high levels of philopatry observed are generated and maintained by the extreme hydrology of the lower Congo River.

Population differentiation in perch Perca fluviatilis: environmental effects on gene flow?

Environmental parameters were used to investigate barriers to gene flow and genetic differentiation in the Eurasian perch (Perca fluviatilis L.) at a small geographical scale in an archipelago system. Significant genetic differentiation was found among locations. Distance per se did not play a major role in the reduction of gene flow. Instead, the largest genetic differences between populations correlated with major changes in environmental conditions, such as temperature at time of spawning. The results show that genetic divergence can arise between populations in habitats thought to be highly connected, and that environmental variables can influence the level of gene flow between populations, including those that are at small spatial scales (tens of kilometres). The importance of a landscape approach when investigating genetic differentiation and defining barriers to gene flow is highlighted.

Multiple interacting loci control sex determination in lake Malawi cichlid fish

Several models have been proposed to suggest how the evolution of sex-determining mechanisms might contribute to speciation. Here, we describe the inheritance of sex in 19 fish species from the rapidly evolving flock of cichlids in Lake Malawi, Africa. We found that many of these species have a male heterogametic (XY) system on linkage group 7. Some species also segregate for a female heterogametic (ZW) system on linkage group 5 that is coincident with a dominant orange-blotch (OB) color pattern in females. The ZW system is epistatically dominant to the XY system when both are segregating within a family. Several lines of evidence suggest that additional sex-determining loci are segregating in some species. These results are consistent with the idea that genetic conflicts play an important role in the evolution of these species flocks and suggest that evolution of sex-determining mechanisms has contributed to the radiation of cichlid fish in East Africa.

Phylogeography and sympatric differentiation of the Arctic charr Salvelinus alpinus (L.) complex in Siberia as revealed by mtDNA sequence analysis

Sequence variation in the mtDNA control region of Arctic charr Salvelinus alpinus and Dolly Varden Salvelinus malma from 56 Siberian and North American populations was analysed to assess their phylogeographic relationships and the origins of sympatric forms. Phylogenetic trees confirm the integrity of phylogroups reported in previous mtDNA studies except that the Siberian group does not separate as a single cluster. Haplotype network analysis indicates the proximity of Siberian and Atlantic haplotypes. These are considered as one Eurasian group represented by the Atlantic, east Siberian (interior Siberia including Transbaikalia, Taimyr) and Eurosiberian (Finland, Spitsbergen, Taimyr) sub-groups. Salvelinus alpinus with presumably introgressed Bering group (malma) haplotypes were found along eastern Siberian coasts up to the Olenek Bay and the Lena Delta region, where they overlap with the Eurasian group and in the easternmost interior region. It is proposed that Siberia was colonized by S. alpinus in two stages: from the west by the Eurasian group and later from the east by the Bering group. The high diversity of Eurasian group haplotypes in Siberia indicates its earlier colonization by S. alpinus as compared with the European Alps. This colonization was rapid, proceeded from a diverse gene pool, and was followed by differential survival of ancestral mtDNA lineages in different basins and regions, and local mutational events in isolated populations. The results presented here support a northern origin of Transbaikalian S. alpinus, the dispersion of S. alpinus to the Lake Baikal Basin from the Lena Basin, segregation of S. alpinus between Lena tributaries and their restricted migration over the divides between sub-basins. These results also support sympatric origin of intralacustrine forms of S. alpinus.

Contrasting patterns of spatial genetic structure in sympatric rock-dwelling cichlid fishes

The cichlid fish of Lake Tanganyika in Eastern Africa are a celebrated example of both ecological and species diversification. Because population subdivision is likely to play an important role in the speciation process, understanding how habitat features interact with species' demographic, behavioral and ecological attributes to influence gene flow and population divergence may help explain the causes of high species richness in this and other systems. Here, we test the roles of isolation-by-habitat and isolation-by-distance in generating fine-scale population genetic structure in three sympatric species of habitat-restricted cichlids in Lake Tanganyika. Using multilocus microsatellite genotypes, we contrast patterns of population differentiation in these habitat specialists along a mosaic coastline of both favorable and unfavorable habitat. Despite their close phylogenetic relationship and shared habitat affinity, these species show striking differences in their pattern of genetic subdivision within the same geographical region, suggesting substantially different patterns of gene flow. In particular, two trophically specialized species exhibit much more restricted gene flow over sandy habitat than a trophically opportunistic species. This result suggests that ecological and behavioral traits have a strong influence on the scale and degree of population subdivision, a finding that has potentially important implications for understanding differential propensities for diversification among lineages and phylogenetic patterns of diversity.

Morphological and genetic structuring in the Utah Lake sucker complex

Population decline in the federally endangered June sucker (Chasmistes liorus), a lakesucker unique to Utah Lake, Utah, has been attributed in part to hybridization with the more widespread Utah sucker (Catostomus ardens). As a group, suckers in Utah Lake exhibit considerable external morphological variation. Meristic and morphological ambiguities, presumably the result of hybridization, create a continuum of intermediate forms between Chasmistes and Catostomus extremes and prevent definitive identification to species. Here we describe and evaluate the morphological and genetic variation in suckers in Utah Lake by comparing a morphological analysis with amplified fragment length polymorphism and microsatellite analyses. Suckers were morphologically differentiated using mouth characters associated with different feeding strategies: planktivory (June sucker) and benthivory (Utah sucker). Although we found no genetic evidence for a deep divergence between June and Utah morphs, significant, but slight population structuring accompanied the substantial morphological variation. Bayesian model-based genetic clustering analyses detected two sucker populations in Utah Lake; however, these clusters were not strongly concordant with morphological groupings or between marker systems. The suckers in Utah Lake present an interesting dilemma regarding conservation: should one conserve (breed and stock) a subset of the morphotypic variation in the Utah Lake sucker complex, focusing on the endangered June sucker morphotype, or should one conserve both June sucker and Utah sucker morphotypes in this complex, possibly maximizing evolutionary potential? We explore this question in the context of current genetic and morphological variation in the Utah Lake sucker complex as well as historical information on this complex and other lakesuckers.

Comparative analysis reveals signatures of differentiation amid genomic polymorphism in Lake Malawi cichlids

Low coverage survey sequencing shows that although Lake Malawi cichlids are phenotypically and behaviorally diverse, they appear genetically like a subdivided population.

Background

Cichlid fish from East Africa are remarkable for phenotypic and behavioral diversity on a backdrop of genomic similarity. In 2006, the Joint Genome Institute completed low coverage survey sequencing of the genomes of five phenotypically and ecologically diverse Lake Malawi species. We report a computational and comparative analysis of these data that provides insight into the mechanisms that make closely related species different from one another.

Results

We produced assemblies for the five species ranging in aggregate length from 68 to 79 megabase pairs, identified putative orthologs for more than 12,000 human genes, and predicted more than 32,000 cross-species single nucleotide polymorphisms (SNPs). Nucleotide diversity was lower than that found among laboratory strains of the zebrafish. We collected around 36,000 genotypes to validate a subset of SNPs within and among populations and across multiple individuals of about 75 Lake Malawi species. Notably, there were no fixed differences observed between focal species nor between major lineages. Roughly 3% to 5% of loci surveyed are statistical outliers for genetic differentiation (F_ST) within species, between species, and between major lineages. Outliers for F_ST are candidate genes that may have experienced a history of natural selection in the Malawi lineage.

Conclusion

We present a novel genome sequencing strategy, which is useful when evolutionary diversity is the question of interest. Lake Malawi cichlids are phenotypically and behaviorally diverse, but they appear genetically like a subdivided population. The unique structure of Lake Malawl cichlid genomes should facilitate conceptually new experiments, employing SNPs to identity genotype-phenotype association, using the entire species flock as a mapping panel.

Ecological separation in a polymorphic terrestrial salamander

1. When studying speciation, researchers commonly examine reproductive isolation in recently diverged populations. Polymorphic species provide an opportunity to examine the role of reproductive isolation in populations that may be in the process of divergence. 2. We examined a polymorphic population of Plethodon cinereus (red-backed salamanders) for evidence of sympatric ecological separation by colour morphology. Recent studies have correlated temperature and climate with colour morphology in this species, but no studies have looked at differences in diet or mate choice between colour morphs. We used artificial cover objects to assess salamander diet, mating preference and surface activity over a 2-year period at a field site in north-eastern Ohio. 3. We detected differences in diet between two colour morphs, striped and unstriped. The diets of striped individuals were significantly more diverse and were made up of more profitable prey than the diets of unstriped salamanders. 4. Opposite sex pairs were made up of individuals of the same colour morph and striped males were found more often with larger females than were unstriped males. 5. We corroborate findings of earlier studies suggesting that the unstriped form is adapted to warmer conditions. Unstriped individuals were the first to withdraw from the forest floor as temperatures fell in the late fall. We found no evidence that the colour morphs responded differently to abiotic factors such as soil moisture and relative humidity, and responses to surface temperatures were also equivocal. 6. We conclude that the two colour morphs exhibit some degree of ecological separation and tend to mate assortatively, but are unlikely to be undergoing divergence given the observed frequency of intermorph pairings.

Postglacial intra-lacustrine divergence of Icelandic threespine stickleback morphs in three neovolcanic lakes

The geographical context of divergence and local adaptation of lacustrine fish is controversial. Despite recent theoretical support for sympatric and parapatric divergence, empirical studies providing unequivocal support for this remain scant. An important component of such a case would be where multiple lakes have different morphs and a range of markers, both mitochondrial and nuclear, show monophyly within lakes. Here we describe such a situation for threespine sticklebacks in three lakes in Iceland. By analysing the variation at nuclear and mitochondrial markers in several freshwater and marine populations as well as three pairs of intra-lacustrine morphs we infer their phylogenetic relationships and colonization pattern. There were high levels of microsatellite variation in all populations and no evidence was found for either repeated colonization of marine fish or colonization from distinct glacial refugia. Intra-lacustrine threespine stickleback morphs in all three lakes show significant genetic divergence probably indicating restricted gene flow.

Species-Specific Population Structure in Rock-Specialized Sympatric Cichlid Species in Lake Tanganyika, East Africa

Species richness and geographical phenotypic variation in East African lacustrine cichlids are often correlated with ecological specializations and limited dispersal. This study compares mitochondrial and microsatellite genetic diversity and structure among three sympatric rock-dwelling cichlids of Lake Tanganyika, Eretmodus cyanostictus, Tropheus moorii, and Ophthalmotilapia ventralis. The species represent three endemic, phylogenetically distinct tribes (Eretmodini, Tropheini, and Ectodini), and display divergent ecomorphological and behavioral specialization. Sample locations span both continuous, rocky shoreline and a potential dispersal barrier in the form of a muddy bay. High genetic diversity and population differentiation were detected in T. moorii and E. cyanostictus, whereas much lower variation and structure were found in O. ventralis. In particular, while a 7-km-wide muddy bay curtails dispersal in all three species to a similar extent, gene flow along mostly continuous habitat appeared to be controlled by distance in E. cyanostictus, further restricted by site philopatry and/or minor habitat discontinuities in T. moorii, and unrestrained in O. ventralis. In contrast to the general pattern of high gene flow along continuous shorelines in rock-dwelling cichlids of Lake Malawi, our study identifies differences in population structure among stenotopic Lake Tanganyika species. The amount of genetic differentiation among populations was not related to the degree of geographical variation of body color, especially since more phenotypic variation is observed in O. ventralis than in the genetically highly structured E. cyanostictus.

Genetic population structure as indirect measure of dispersal ability in a Lake Tanganyika cichlid

Diversification and speciation processes are influenced by intrinsic (ecological specialization, dispersal) and extrinsic (habitat structure and instability) factors, but the effect of ecological characteristics on dispersal is difficult to assess. This study uses mitochondrial control region sequences to investigate the population structure and demographic history of the endemic Lake Tanganyika cichlid Neolamprologus caudopunctatus with a preference for the rock-sand interface along two stretches of continuous, rocky shoreline, and across a sandy bay representing a potential dispersal barrier. Populations along uninterrupted habitat were not differentiated; whereas, the sandy bay separated two reciprocally monophyletic clades. The split between the two clades between 170,000 and 260,000 years BP coincides with a period of rising water level following a major lowstand, and indicates that clades remained isolated throughout subsequent lake level fluctuations. Low long-term effective population sizes were inferred from modest genetic diversity estimates, and may be due to recent population expansions starting from small population sizes 45,000-60,000 years BP. Comparisons with available data from specialized rock-dwelling species of the same area suggest that habitat structure and lake level fluctuations determine phylogeographic patterns on large scales, while fine-scale population structure and demography are modulated by species-specific ecologies.

Studies of threespine stickleback developmental evolution: progress and promise

A promising route for understanding the origin and diversification of organismal form is through studies at the intersection of evolution and development (evo-devo). While much has been learned over the last two decades concerning macroevolutionary patterns of developmental change, a fundamental gap in the evo-devo synthesis is the integration of mathematical population and quantitative genetics with studies of how genetic variation in natural populations affects developmental processes. This micro-evo-devo synthesis requires model organisms with which to ask empirical questions. Threespine stickleback fish (Gasterosteus aculeatus), long a model for studying behavior, ecology and evolution, is emerging as a prominent model micro-evo-devo system. Research on stickleback over the last decade has begun to address the genetic basis of morphological variation and sex determination, and much of this work has important implications for understanding the genetics of speciation. In this paper we review recent threespine stickleback micro-evo-devo results, and outline the resources that have been developed to make this synthesis possible. The prospects for stickleback research to speed the micro-(and macro-) evo-devo syntheses are great, and this workhorse model system is well situated to continue contributing to our understanding of the generation of diversity in organismal form for many more decades.

Genetic differentiation and reproductive isolation in a Cyprinodon fish species flock from Laguna Chichancanab, Mexico

The Cyprinodon species flock from Laguna Chichancanab, aged 8000 years, provides another potential case of sympatric speciation. The flock consists of seven morphologically distinct species, each within partially different trophic niches, and a group of specimens which cannot unequivocally be assigned to one of these species. Genetic analyses, based on mtDNA and five microsatellite loci, revealed significant genetic differentiation of one species, C. maya, from other members of the species flock, providing strong evidence for reproductive isolation. For the remaining members of the flock significant genetic structuring was detected, with some evidence of gene flow with the most abundant species C. beltrani. These analyses suggest that speciation proceeds with ongoing hybridisation, and further suggest that the morphologically unidentifiable specimens found in the lake are probably hybrids. I propose that in the Cyprinodon species flock besides disruptive selection sexual selection plays an important part in achieving and maintaining reproductive isolation.

Nuclear gene variation and molecular dating of the cichlid species flock of Lake Malawi

The cichlid fishes of Lake Malawi are famously diverse. However, evolutionary studies have been difficult because of their recent and uncertain phylogenetic history. Portions of 12 nuclear loci were sequenced in nine rock-dwelling species (mbuna) and three representatives of pelagic nonmbuna species. In contrast to the pattern of variation at mitochondrial genes, which do provide phylogenetic resolution at the level of mbuna versus nonmbuna, and among some genera, the nuclear loci were virtually devoid of phylogenetic signal. Only a small minority of variable positions were phylogenetically informative, and no phylogenetic branches are supported by more than one site. From the nuclear gene perspective the Malawian radiation appears to be a star phylogeny, as if the founding of the lake was accompanied by a partial bottleneck. The pattern is different from that found in Lake Victoria, in which nuclear loci share large amounts of ancestral variation. In the case of nuclear genes of Lake Malawi, the absence of phylogenetically informative variation suggests a relative absence of ancestral variation. Nuclear genes also differed from the mitochondria in having nearly twice the amount of divergence from Oreochromis (tilapia). An approximate splitting time between mbuna and nonmbuna lineages was estimated as 0.7 Myr. Oreochromis is estimated to have diverged from the cichlids in Lake Malawi and Lake Tanganyika about 18 MYA.

Behavioral genomics and the study of speciation at a porous species boundary

Porous species boundaries are characterized by differential gene flow, where some regions of the genome experience divergent evolution while others experience the homogenizing effects of gene flow. If species can arise or remain distinct despite gene flow between them, speciation can only be understood on a gene by gene level. To understand the genetics of speciation, we therefore must identify the targets of selection that cause divergent evolution and identify the genetic architecture underlying such "speciation phenotypes". This will enable characterization of genomic regions that are "free to flow" between species, and those that diverge in the face of gene flow. We discuss this problem in the genus Laupala, a morphologically cryptic, flightless group of crickets that has radiated in Hawaii. Because songs are used in courtship and always distinguish close relatives of Laupala as well as species in sympatry, we argue that songs in Laupala are speciation phenotypes. Here, we present our approaches to identify the underlying genomic regions and song genes that differentiate closely related species. We discuss what is known about the genetic basis of this species difference derived from classic quantitative genetics and quantitative trait locus mapping experiments. We also present a model of the molecular expression of cricket song to assist in our goal to identify the genes involved in song variation. As most species are sympatric and exchange genes with congeners, we discuss the importance of understanding the genetic and genomic architecture of song as a speciation phenotype that must be characterized to identify differential patterns of gene flow at porous species boundaries.

On the origin of Lake Malawi cichlid species: a population genetic analysis of divergence

The cichlid fishes of Lake Malawi are famously diverse. However, phylogenetic and population genetic studies of their history have been difficult because of the great amount of genetic variation that is shared between species. We apply a recently developed method for fitting the "isolation with migration" divergence model to a data set of specially designed compound loci to develop portraits of cichlid species divergence. Outgroup sequences from a cichlid from Lake Tanganyika permit model parameter estimates in units of years and effective population sizes. Estimated speciation times range from 1,000 to 17,000 years for species in the genus Tropheops. These exceptionally recent dates suggest that Malawi cichlids as a group experience a very active and dynamic diversification process. Current effective population size estimates range form 2,000 to near 40,000, and to >120,000 for estimates of ancestral population sizes. It appears that very recent speciation and gene flow are among the reasons why it has been difficult to discern the phylogenetic history of Malawi cichlids.

Male-male competition and nuptial-colour displacement as a diversifying force in Lake Victoria cichlid fishes

We propose a new mechanism for diversification of male nuptial-colour patterns in the rapidly speciating cichlid fishes of Lake Victoria. Sympatric closely related species often display nuptial colours at opposite ends of the spectrum with males either blue or yellow to red. Colour polymorphisms within single populations are common too. We propose that competition between males for breeding sites promotes such colour diversification, and thereby speciation. We hypothesize that male aggression is primarily directed towards males of the common colour, and that rare colour morphs enjoy a negatively frequency-dependent fitness advantage. We test our hypothesis with a large dataset on the distributions and nuptial colorations of 52 species on 47 habitat islands in Lake Victoria, and with a smaller dataset on the within-spawning-site distributions of males with different coloration. We report that territories of males of the same colour are negatively associated on the spawning site, and that the distribution of closely related species over habitat islands is determined by nuptial coloration in the fashion predicted by our hypothesis. Whereas among unrelated species those with similar nuptial colour are positively associated, among closely related species those with similar colour are negatively associated and those with different colour are positively associated. This implies that negatively frequency-dependent selection on nuptial coloration among closely related species is a sufficiently strong force to override other effects on species distributions. We suggest that male-male competition is an important and previously neglected agent of diversification among haplochromine cichlid fishes.

Variation in habitat preference and population structure among three species of the Lake Malawi cichlid genus Protomelas

Several studies have demonstrated strong population structuring over small distances in the rocky-shore mbuna cichlid fishes from Lake Malawi, suggesting the potential for allopatric speciation. However, many endemic Lake Malawi cichlids are neither mbuna, nor confined to rocky shores. Using microsatellites, we investigated the population structure in three species of the non-mbuna genus Protomelas. The rocky-shore P. taeniolatus showed high levels of population structure even over distances of less than 1 km, while the sandy-shore species P. similis showed no significant structure over distances up to 21 km. Protomelas fenestratus, which is generally found at the interface between rocks and sand, also showed low levels of population structure. Our results suggest that the model of allopatric speciation based on habitat fragmentation within the current lake basin may be equally applicable to rocky-shore non-mbuna as to mbuna, but that an alternative model is required to explain speciation among sandy-shore species as well as the deep-water and pelagic species.

Population-structure and genetic diversity in a haplochromine fish cichlid of a satellite lake of Lake Victoria

The approximately 500 species of the cichlid fish species flock of Lake Victoria, East Africa, have evolved in a record-setting 100,000 years and represent one of the largest adaptive radiations. We examined the population structure of the endangered cichlid species Xystichromis phytophagus from Lake Kanyaboli, a satellite lake to Lake Victoria in the Kenyan Yala wetlands. Two sets of molecular markers were analysed--sequences of the mitochondrial control region as well as six microsatellite loci--and revealed surprisingly high levels of genetic variability in this species. Mitochondrial DNA sequences failed to detect population structuring among the three sample populations. A model-based population assignment test based on microsatellite data revealed that the three populations most probably aggregate into a larger panmictic population. However, values of population pairwise FST indicated moderate levels of genetic differentiation for one population. Eleven distinct mitochondrial haplotypes were found among 205 specimens of X. phytophagus, a relatively high number compared to the total number of 54 haplotypes that were recovered from hundreds of specimens of the entire cichlid species flock of Lake Victoria. Most of the X. phytophagus mitochondrial DNA haplotypes were absent from the main Lake Victoria, corroborating the putative importance of satellite lakes as refugia for haplochromine cichlids that went extinct from the main lake in the last decades and possibly during the Late Pleistocene desiccation of Lake Victoria.

The species flocks of East African cichlid fishes: recent advances in molecular phylogenetics and population genetics

With more than 3,000 species, the fish family Cichlidae is one of the most species-rich families of vertebrates. Cichlids occur in southern and central America, Africa, Madagascar, and India. The hotspot of their biodiversity is East Africa, where they form adaptive radiations composed of hundreds of endemic species in several lakes of various sizes and ages. The unparalleled species richness of East African cichlids has been something of a conundrum for evolutionary biologists and ecologists, since it has been in doubt whether these hundreds of species arose by allopatric speciation or whether it is necessary to invoke somewhat less traditional models of speciation, such as micro-allopatric, peripatric, or even sympatric speciation or evolution through sexual selection mediated by female choice. Ernst Mayr's analyses of these evolutionary uniquely diverse species assemblages have contributed to a more direct approach to this problem and have led to a deeper understanding of the patterns and processes that caused the formation of these huge groups of species. We review here recent molecular data on population differentiation and phylogenetics, which have helped to unravel, to some extent, the patterns and processes that led to the formation and ecological maintenance of cichlid species flocks. It is becoming apparent that sexually selected traits do play an important role in speciation in micro-allopatric or even sympatric settings. Species richness seems to be roughly correlated with the surface area, but not the age, of the lakes. We observe that the oldest lineages of a species flock of cichlids are often less species-rich and live in the open water or deepwater habitats. While the species flocks of the Lake Malawai and the Lake Victoria areas were shown to be monophyletic, the cichlid assemblage of Lake Tanganyika seems to consist of several independent species flocks. Cichlids emerge as an evolutionary model system in which many fundamental questions in evolution and ecology can be tested successfully, yet for other fish species flocks the relative importance of alternative mechanisms of speciation is likely to differ from that in cichlid fish.

Divergent selection during speciation of Lake Malawi cichlid fishes inferred from parallel radiations in nuptial coloration

Repeated evolution of the same phenotypic difference during independent episodes of speciation is strong evidence for selection during speciation. More than 1,000 species of cichlids, >10% of the world's freshwater fish species, have arisen within the past million years in Lakes Malawi and Victoria in eastern Africa. Many pairs of closely related sympatric species differ in their nuptial coloration in very similar ways. Nuptial coloration is important in their mate choice, and speciation by sexual selection on genetically or ecologically constrained variation in nuptial coloration had been proposed, which would repeatedly produce similar nuptial types in different populations, a prediction that was difficult to test in the absence of population-level phylogenies. We measured genetic similarity between individuals within and between populations, species, and lake regions by typing 59 individuals at >2,000 polymorphic genetic loci. From these data, we reconstructed, to our knowledge, the first larger species level phylogeny for the most diverse group of Lake Malawi cichlids. We used the genetic and phylogenetic data to test the divergent selection scenario against colonization, character displacement, and hybridization scenarios that could also explain diverse communities. Diversity has arisen by replicated radiations into the same color types, resulting in phenotypically very different, yet closely related, species within and phenotypically highly similar yet unrelated sets of species between regions, which is consistent with divergent selection during speciation and is inconsistent with colonization and character displacement models.

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.

Phylogeography of Lake Malawi cichlids of the genus Pseudotropheus: significance of allopatric colour variation

One of the most compelling features of the cichlid fishes of the African Great Lakes is the seemingly endless diversity of male coloration. Colour diversification has been implicated as an important factor driving cichlid speciation. Colour has also been central to cichlid taxonomy and, thus, to our concept of species diversity. We undertook a phylogeographical examination of several allopatric populations of the Lake Malawi cichlid Pseudotropheus zebra in order to reconstruct the evolutionary history of the populations, which exhibit one of two dorsal fin colours. We present evidence that populations with red dorsal fins (RT) are not monophyletic. The RT population defining the northern limit of the distribution has evidently originated independently of the southern RT populations, which share a common ancestry. This evidence of species-level colour convergence is an important discovery in our understanding of cichlid evolution. It implies that divergence in coloration may accompany speciation, and that allopatric populations with similar coloration cannot be assumed to be conspecific. In addition to this finding, we have observed evidence for introgression, contributing to current evidence that this phenomenon may be extremely widespread. Thus, in species-level phylogenetic reconstructions, including our own, consideration must be given to the potential effects of introgression.

The estimation of population differentiation with microsatellite markers

Microsatellite markers are routinely used to investigate the genetic structuring of natural populations. The knowledge of how genetic variation is partitioned among populations may have important implications not only in evolutionary biology and ecology, but also in conservation biology. Hence, reliable estimates of population differentiation are crucial to understand the connectivity among populations and represent important tools to develop conservation strategies. The estimation of differentiation is c from Wright's FST and/or Slatkin's RST, an FST -analogue assuming a stepwise mutation model. Both these statistics have their drawbacks. Furthermore, there is no clear consensus over their relative accuracy. In this review, we first discuss the consequences of different temporal and spatial sampling strategies on differentiation estimation. Then, we move to statistical problems directly associated with the estimation of population structuring itself, with particular emphasis on the effects of high mutation rates and mutation patterns of microsatellite loci. Finally, we discuss the biological interpretation of population structuring estimates.

Genetic evidence for adaptation-driven incipient speciation of Drosophila melanogaster along a microclimatic contrast in "Evolution Canyon," Israel

Substantial genetic differentiation, as great as among species, exists between populations of Drosophila melanogaster inhabiting opposite slopes of a small canyon. Previous work has shown that prezygotic sexual isolation and numerous differences in stress-related phenotypes have evolved between D. melanogaster populations in "Evolution Canyon," Israel, in which slopes 100-400 m apart differ dramatically in aridity, solar radiation, and associated vegetation. Because the canyon's width is well within flies' dispersal capabilities, we examined genetic changes associated with local adaptation and incipient speciation in the absence of geographical isolation. Here we report remarkable genetic differentiation of microsatellites and divergence in the regulatory region of hsp70Ba which encodes the major inducible heat shock protein of Drosophila, in the two populations. Additionally, an analysis of microsatellites suggests a limited exchange of migrants and lack of recent population bottlenecks. We hypothesize that adaptation to the contrasting microclimates overwhelms gene flow and is responsible for the genetic and phenotypic divergence between the populations.