Persistence of neutral polymorphisms in Lake Victoria cichlid fish

Strong species structure but weak geographical structure in demersal Lake Victoria cichlids

Studying phenotypic and genetic differentiation between very young species can be very informative with regard to learning about processes of speciation. Identifying and characterizing genetic species structure and distinguishing it from spatial genetic structure within a species is a prerequisite for this and is often not given sufficient attention. Young radiations of cichlid fish are classical speciation study systems. However, it is only during the past decade that population genomics based on next-generation sequencing has begun to provide the power to resolve species and distinguish speciation from spatial population structure for the youngest of these radiations. The Lake Victoria haplochromine cichlids constitute the youngest large cichlid fish radiation, probably <20,000 years old. Earlier work showed that communities of rocky reef cichlids are composed of many reciprocally monophyletic species despite their very recent origins. Here, we build on this work by studying assemblages of offshore demersal cichlids, adding analyses of within-species spatial structure to the sympatric species structure. We sampled seven multispecies communities along a 6-km-long transect from one side of the Mwanza Gulf to the other side. We investigated whether phenotypically diagnosed putative species are reciprocally monophyletic and whether such monophyly is stable across species geographic ranges. We show that all species are genetically strongly differentiated in sympatry, that they are reciprocally monophyletic, and that monophyly is stable across distribution ranges. We found significant differentiation between geographically distinct populations in two species, but no or weak isolation by distance. We further found subtle but significant morphological differences between all species and a linear relationship between genomic and morphological distance which suggests that differences in morphology begin to accumulate after speciation has already affected genome-wide restrictions of gene flow.

Morphometric and Genetic Description of Trophic Adaptations in Cichlid Fishes

Since Darwin, biologists have sought to understand the evolution and origins of phenotypic adaptations. The skull is particularly diverse due to intense natural selection on feeding biomechanics. We investigated the genetic and molecular origins of trophic adaptation using Lake Malawi cichlids, which have undergone an exemplary evolutionary radiation. We analyzed morphological differences in the lateral and ventral head shape among an insectivore that eats by suction feeding, an obligate biting herbivore, and their F2 hybrids. We identified variation in a series of morphological traits—including mandible width, mandible length, and buccal length—that directly affect feeding kinematics and function. Using quantitative trait loci (QTL) mapping, we found that many genes of small effects influence these craniofacial adaptations. Intervals for some traits were enriched in genes related to potassium transport and sensory systems, the latter suggesting co-evolution of feeding structures and sensory adaptations for foraging. Despite these indications of co-evolution of structures, morphological traits did not show covariation. Furthermore, phenotypes largely mapped to distinct genetic intervals, suggesting that a common genetic basis does not generate coordinated changes in shape. Together, these suggest that craniofacial traits are mostly inherited as separate modules, which confers a high potential for the evolution of morphological diversity. Though these traits are not restricted by genetic pleiotropy, functional demands of feeding and sensory structures likely introduce constraints on variation. In all, we provide insights into the quantitative genetic basis of trophic adaptation, identify mechanisms that influence the direction of morphological evolution, and provide molecular inroads to craniofacial variation.

SINEs as Credible Signs to Prove Common Ancestry in the Tree of Life: A Brief Review of Pioneering Case Studies in Retroposon Systematics

Currently, the insertions of SINEs (and other retrotransposed elements) are regarded as one of the most reliable synapomorphies in molecular systematics. The methodological mainstream of molecular systematics is the calculation of nucleotide (or amino acid) sequence divergences under a suitable substitution model. In contrast, SINE insertion analysis does not require any complex model because SINE insertions are unidirectional and irreversible. This straightforward methodology was named the “SINE method,” which resolved various taxonomic issues that could not be settled by sequence comparison alone. The SINE method has challenged several traditional hypotheses proposed based on the fossil record and anatomy, prompting constructive discussions in the Evo/Devo era. Here, we review our pioneering SINE studies on salmon, cichlids, cetaceans, Afrotherian mammals, and birds. We emphasize the power of the SINE method in detecting incomplete lineage sorting by tracing the genealogy of specific genomic loci with minimal noise. Finally, in the context of the whole-genome era, we discuss how the SINE method can be applied to further our understanding of the tree of life.

Evolutionary divergence of novel open reading frames in cichlids speciation

Novel open reading frames (nORFs) with coding potential may arise from noncoding DNA. Not much is known about their emergence, functional role, fixation in a population or contribution to adaptive radiation. Cichlids fishes exhibit extensive phenotypic diversification and speciation. Encounters with new environments alone are not sufficient to explain this striking diversity of cichlid radiation because other taxa coexistent with the Cichlidae demonstrate lower species richness. Wagner et al. analyzed cichlid diversification in 46 African lakes and reported that both extrinsic environmental factors and intrinsic lineage-specific traits related to sexual selection have strongly influenced the cichlid radiation, which indicates the existence of unknown molecular mechanisms responsible for rapid phenotypic diversification, such as emergence of novel open reading frames (nORFs). In this study, we integrated transcriptomic and proteomic signatures from two tissues of two cichlids species, identified nORFs and performed evolutionary analysis on these nORF regions. Our results suggest that the time scale of speciation of the two species and evolutionary divergence of these nORF genomic regions are similar and indicate a potential role for these nORFs in speciation of the cichlid fishes.

Patterns of genomic differentiation between two Lake Victoria cichlid species, Haplochromis pyrrhocephalus and H. sp. 'macula'

BACKGROUND

The molecular basis of the incipient stage of speciation is still poorly understood. Cichlid fish species in Lake Victoria are a prime example of recent speciation events and a suitable system to study the adaptation and reproductive isolation of species.

RESULTS

Here, we report the pattern of genomic differentiation between two Lake Victoria cichlid species collected in sympatry, Haplochromis pyrrhocephalus and H. sp. 'macula,' based on the pooled genome sequences of 20 individuals of each species. Despite their ecological differences, population genomics analyses demonstrate that the two species are very close to a single panmictic population due to extensive gene flow. However, we identified 21 highly differentiated short genomic regions with fixed nucleotide differences. At least 15 of these regions contained genes with predicted roles in adaptation and reproductive isolation, such as visual adaptation, circadian clock, developmental processes, adaptation to hypoxia, and sexual selection. The nonsynonymous fixed differences in one of these genes, LWS, were reported as substitutions causing shift in absorption spectra of LWS pigments. Fixed differences were found in the promoter regions of four other differentially expressed genes, indicating that these substitutions may alter gene expression levels.

CONCLUSIONS

These diverged short genomic regions may have contributed to the differentiation of two ecologically different species. Moreover, the origins of adaptive variants within the differentiated regions predate the geological formation of Lake Victoria; thus Lake Victoria cichlid species diversified via selection on standing genetic variation.

Phylogeographic analysis suggests two origins for the riparian azalea Rhododendron indicum (L.) Sweet

Ecological speciation is an important factor in the diversification of plants. The distribution of the woody species Rhododendron indicum, which grows along rivers and is able to withstand water flow when rivers flood (i.e. it is a rheophyte), is disjunct, in contrast to the widespread distribution of its relative, Rhododendron kaempferi. This study aimed to elucidate the phylogenetic relationships between R. indicum and R. kaempferi and the evolutionary processes that gave rise to them. The sequences of three non-coding chloroplast DNA regions (total length 1977 bp) were obtained from 21 populations covering the ranges of the two species. In addition, genome-wide SNPs were genotyped from 20 populations using a genotyping by sequencing method. Leaf morphologies were measured for eight representative populations. Two chloroplast DNA haplotypes, which were detected in R. indicum, were shared between the two species. Genome-wide SNPs identified two lineages in R. indicum and these lineages did not constitute a monophyletic group. Each of these two lineages was related to geographically close populations of R. kaempferi. Leaf morphology, which is a characteristic feature in rheophytes, was not differentiated between the two lineages in R. indicum. The morphological similarity between the two heterogeneous lineages may be a result of parallel evolution from R. kaempferi or of introgressive hybridization between the species due to strong selective pressure imposed by flooding.

Visual adaptation in Lake Victoria cichlid fishes: depth-related variation of color and scotopic opsins in species from sand/mud bottoms

Background

For Lake Victoria cichlid species inhabiting rocky substrates with differing light regimes, it has been proposed that adaptation of the long-wavelength-sensitive (LWS) opsin gene triggered speciation by sensory drive through color signal divergence. The extensive and continuous sand/mud substrates are also species-rich, and a correlation between male nuptial coloration and the absorption of LWS pigments has been reported. However, the factors driving genetic and functional diversity of LWS pigments in sand/mud habitats are still unresolved.

Results

To address this issue, nucleotide sequences of eight opsin genes were compared in ten Lake Victoria cichlid species collected from sand/mud bottoms. Among eight opsins, the LWS and rod-opsin (RH1) alleles were diversified and one particular allele was dominant or fixed in each species. Natural selection has acted on and fixed LWS alleles in each species. The functions of LWS and RH1 alleles were measured by absorption of reconstituted A1- and A2-derived visual pigments. The absorption of pigments from RH1 alleles most common in deep water were largely shifted toward red, whereas those of LWS alleles were largely shifted toward blue in both A1 and A2 pigments. In both RH1 and LWS pigments, A2-derived pigments were closer to the dominant light in deep water, suggesting the possibility of the adaptation of A2-derived pigments to depth-dependent light regimes.

Conclusions

The RH1 and LWS sequences may be diversified for adaptation of A2-derived pigments to different light environments in sand/mud substrates. Diversification of the LWS alleles may have originally taken place in riverine environments, with a new mutation occurring subsequently in Lake Victoria.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1040-x) contains supplementary material, which is available to authorized users.

Local Climate Heterogeneity Shapes Population Genetic Structure of Two Undifferentiated Insular Scutellaria Species

Spatial climate heterogeneity may not only affect adaptive gene frequencies but could also indirectly shape the genetic structure of neutral loci by impacting demographic dynamics. In this study, the effect of local climate on population genetic variation was tested in two phylogenetically close Scutellaria species in Taiwan. Scutellaria taipeiensis, which was originally assumed to be an endemic species of Taiwan Island, is shown to be part of the widespread species S. barbata based on the overlapping ranges of genetic variation and climatic niches as well as their morphological similarity. Rejection of the scenario of "early divergence with secondary contact" and the support for multiple origins of populations of S. taipeiensis from S. barbata provide strong evolutionary evidence for a taxonomic revision of the species combination. Further tests of a climatic effect on genetic variation were conducted. Regression analyses show nonlinear correlations among any pair of geographic, climatic, and genetic distances. However, significantly, the bioclimatic variables that represent the precipitation from late summer to early autumn explain roughly 13% of the genetic variation of our sampled populations. These results indicate that spatial differences of precipitation in the typhoon season may influence the regeneration rate and colonization rate of local populations. The periodic typhoon episodes explain the significant but nonlinear influence of climatic variables on population genetic differentiation. Although, the climatic difference does not lead to species divergence, the local climate variability indeed impacts the spatial genetic distribution at the population level.

Simple Biophysical Model Predicts Faster Accumulation of Hybrid Incompatibilities in Small Populations Under Stabilizing Selection

Speciation is fundamental to the process of generating the huge diversity of life on Earth. However, we are yet to have a clear understanding of its molecular-genetic basis. Here, we examine a computational model of reproductive isolation that explicitly incorporates a map from genotype to phenotype based on the biophysics of protein–DNA binding. In particular, we model the binding of a protein transcription factor to a DNA binding site and how their independent coevolution, in a stabilizing fitness landscape, of two allopatric lineages leads to incompatibilities. Complementing our previous coarse-grained theoretical results, our simulations give a new prediction for the monomorphic regime of evolution that smaller populations should develop incompatibilities more quickly. This arises as (1) smaller populations have a greater initial drift load, as there are more sequences that bind poorly than well, so fewer substitutions are needed to reach incompatible regions of phenotype space, and (2) slower divergence when the population size is larger than the inverse of discrete differences in fitness. Further, we find longer sequences develop incompatibilities more quickly at small population sizes, but more slowly at large population sizes. The biophysical model thus represents a robust mechanism of rapid reproductive isolation for small populations and large sequences that does not require peak shifts or positive selection. Finally, we show that the growth of DMIs with time is quadratic for small populations, agreeing with Orr’s model, but nonpower law for large populations, with a form consistent with our previous theoretical results.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

High prevalence of non-synonymous substitutions in mtDNA of cichlid fishes from Lake Victoria

When a population size is reduced, genetic drift may fix slightly deleterious mutations, and an increase in nonsynonymous substitution is expected. It has been suggested that past aridity has seriously affected and decreased the populations of cichlid fishes in Lake Victoria, while geographical studies have shown that the water levels in Lake Tanganyika and Lake Malawi have remained fairly constant. The comparably stable environments in the latter two lakes might have kept the populations of cichlid fishes large enough to remove slightly deleterious mutations. The difference in the stability of cichlid fish population sizes between Lake Victoria and the Lakes Tanganyika and Malawi is expected to have caused differences in the nonsynonymous/synonymous ratio, ω (=dN/dS), of the evolutionary rate. Here, we estimated ω and compared it between the cichlids of the three lakes for 13 mitochondrial protein-coding genes using maximum likelihood methods. We found that the lineages of the cichlids in Lake Victoria had a significantly higher ω for several mitochondrial loci. Moreover, positive selection was indicated for several codons in the mtDNA of the Lake Victoria cichlid lineage. Our results indicate that both adaptive and slightly deleterious molecular evolution has taken place in the Lake Victoria cichlids' mtDNA genes, whose nonsynonymous sites are generally conserved.

Do estimated and actual species phylogenies match? Evaluation of East African cichlid radiations

A large number of published phylogenetic estimates are based on a single locus or the concatenation of multiple loci, even though genealogies of single or concatenated loci may not accurately reflect the true history of species diversification (i.e., the species tree). The increased availability of genomic data, coupled with new computational methods, improves resolution of species relationships beyond what was possible in the past. Such developments will no doubt benefit future phylogenetic studies. It remains unclear how robust phylogenies that predate these developments (i.e., the bulk of phylogenetic studies) are to departures from the assumption of strict gene tree-species tree concordance. Here, we present a parametric bootstrap (PBST) approach that assesses the reliability of past phylogenetic estimates in which gene tree-species tree discord was ignored. We focus on a universal cause of discord-the random loss of gene lineages from genetic drift-and apply the method in a meta-analysis of East African cichlids, a group encompassing historical scenarios that are particularly challenging for phylogenetic estimation. Although we identify some evolutionary relationships that are robust to gene tree discord, many past phylogenetic estimates of cichlids are not. We discuss the utility of the PBST method for evaluating the robustness of gene tree-based phylogenetic estimations in general as well as for testing the clade-specific performance of species tree estimation methods and designing sampling strategies that increase the accuracy of estimated species relationships.

Multiple episodic evolution events in V1R receptor genes of East-African cichlids

Fish use olfaction to detect a variety of nonvolatile chemical signals, and thus, this sense is key to survival and communication. However, the contribution of the olfactory sense to social—especially reproductive—interactions in cichlids is still controversial. To obtain insights into this issue, we investigated the genes encoding V1Rs—possible candidates for reproductive pheromone receptors—among East-African cichlids. Interestingly, we found an excess of nonsynonymous over synonymous substitutions in four of six V1R genes in multiple cichlid lineages. First, we found that highly dimorphic V1R2 allele groups were shared among the cichlids inhabiting all East-African Great Lakes emerged through the episodic accumulation of the nonsynonymous substitutions prior to the radiation of the Lake Tanganyika species flock. We further detected such episodic events in V1R1 of the tribe Tropheini, and in V1R3 and V1R6 of the tribe Trematocarini. The excess of nonsynonymous substitutions in these examples were indicated as dN/dS > 1, which were all statistically significant by Fisher’s exact test. Furthermore, we speculate that the amino acid changes in these episodic events are likely functional switch because they occurred in the putative ligand-binding pocket. Our finding of the occurrence of multiple episodic events and the unexpected gene diversity in one unique gene family is suggestive of the contribution of the V1R to the species diversification and the social interaction in cichlids.

The diversity of male nuptial coloration leads to species diversity in Lake Victoria cichlids

The amazing coloration shown by diverse cichlid fish not only fascinates aquarium keepers, but also receives great attention from biologists interested in speciation because of its recently-revealed role in their adaptive radiation in an African lake. We review the important role of coloration in the speciation and adaptive evolution of Lake Victoria cichlids, which have experienced adaptive radiation during a very short evolutionary period. Mature male cichlids display their colors during mate choice. The color of their skin reflects light, and the reflected light forms a color signal that is received by the visual system of females. The adaptive divergence of visual perceptions shapes and diverges colorations, to match the adapted visual perceptions. The divergence of visual perception and coloration indicates that the divergence of color signals causes reproductive isolation between species, and this process leads to speciation. Differences in color signals among coexisting species act to maintain reproductive isolation by preventing hybridization. Thus, the diversity of coloration has caused speciation and has maintained species diversity in Lake Victoria cichlids.

Crater Lake Apoyo revisited--population genetics of an emerging species flock

The polytypic Nicaraguan Midas cichlids (Amphilophus cf. citrinellus) have been established as a model system for studying the mechanisms of speciation and patterns of diversification in allopatry and sympatry. The species assemblage in Crater Lake Apoyo has been accepted as a textbook example for sympatric speciation. Here, we present a first comprehensive data set of population genetic (mtDNA & AFLPs) proxies of species level differentiation for a representative set of individuals of all six endemic Amphilophus species occurring in Crater Lake Apoyo. AFLP genetic differentiation was partitioned into a neutral and non-neutral component based on outlier-loci detection approaches, and patterns of species divergence were explored with Bayesian clustering methods. Substantial levels of admixture between species were detected, indicating different levels of reproductive isolation between the six species. Analysis of neutral genetic variation revealed several A. zaliosus as being introgressed by an unknown contributor, hereby rendering the sympatrically evolving L. Apoyo flock polyphyletic. This is contrasted by the mtDNA analysis delivering a clear monophyly signal with Crater Lake Apoyo private haplotypes characterising all six described species, but also demonstrating different demographic histories as inferred from pairwise mismatch distributions.

Genetic structure of pelagic and littoral cichlid fishes from Lake Victoria

The approximately 700 species of cichlids found in Lake Victoria in East Africa are thought to have evolved over a short period of time, and they represent one of the largest known examples of adaptive radiation. To understand the processes that are driving this spectacular radiation, we must determine the present genetic structure of these species and elucidate how this structure relates to the ecological conditions that caused their adaptation. We analyzed the genetic structure of two pelagic and seven littoral species sampled from the southeast area of Lake Victoria using sequences from the mtDNA control region and 12 microsatellite loci as markers. Using a Bayesian model-based clustering method to analyze the microsatellite data, we separated these nine species into four groups: one group composed of pelagic species and another three groups composed mainly of rocky-shore species. Furthermore, we found significant levels of genetic variation between species within each group at both marker loci using analysis of molecular variance (AMOVA), although the nine species often shared mtDNA haplotypes. We also found significant levels of genetic variation between populations within species. These results suggest that initial groupings, some of which appear to have been related to habitat differences, as well as divergence between species within groups took place among the cichlid species of Lake Victoria.

Origins of shared genetic variation in African cichlids

Cichlid fishes have evolved tremendous morphological and behavioral diversity in the waters of East Africa. Within each of the Great Lakes Tanganyika, Malawi, and Victoria, the phenomena of hybridization and retention of ancestral polymorphism explain allele sharing across species. Here, we explore the sharing of single nucleotide polymorphisms (SNPs) between the major East African cichlid assemblages. A set of approximately 200 genic and nongenic SNPs was ascertained in five Lake Malawi species and genotyped in a diverse collection of ∼160 species from across Africa. We observed segregating polymorphism outside of the Malawi lineage for more than 50% of these loci; this holds similarly for genic versus nongenic SNPs, as well as for SNPs at putative CpG versus non-CpG sites. Bayesian and principal component analyses of genetic structure in the data demonstrate that the Lake Malawi endemic flock is not monophyletic and that river species have likely contributed significantly to Malawi genomes. Coalescent simulations support the hypothesis that river cichlids have transported polymorphism between lake assemblages. We observed strong genetic differentiation between Malawi lineages for approximately 8% of loci, with contributions from both genic and nongenic SNPs. Notably, more than half of these outlier loci between Malawi groups are polymorphic outside of the lake. Cichlid fishes have evolved diversity in Lake Malawi as new mutations combined with standing genetic variation shared across East Africa.

Genome-wide RAD sequence data provide unprecedented resolution of species boundaries and relationships in the Lake Victoria cichlid adaptive radiation

Although population genomic studies using next generation sequencing (NGS) data are becoming increasingly common, studies focusing on phylogenetic inference using these data are in their infancy. Here, we use NGS data generated from reduced representation genomic libraries of restriction-site-associated DNA (RAD) markers to infer phylogenetic relationships among 16 species of cichlid fishes from a single rocky island community within Lake Victoria's cichlid adaptive radiation. Previous attempts at sequence-based phylogenetic analyses in Victoria cichlids have shown extensive sharing of genetic variation among species and no resolution of species or higher-level relationships. These patterns have generally been attributed to the very recent origin (<15,000 years) of the radiation, and ongoing hybridization between species. We show that as we increase the amount of sequence data used in phylogenetic analyses, we produce phylogenetic trees with unprecedented resolution for this group. In trees derived from our largest data supermatrices (3 to >5.8 million base pairs in width), species are reciprocally monophyletic with high bootstrap support, and the majority of internal branches on the tree have high support. Given the difficulty of the phylogenetic problem that the Lake Victoria cichlid adaptive radiation represents, these results are striking. The strict interpretation of the topologies we present here warrants caution because many questions remain about phylogenetic inference with very large genomic data set and because we can with the current analysis not distinguish between effects of shared ancestry and post-speciation gene flow. However, these results provide the first conclusive evidence for the monophyly of species in the Lake Victoria cichlid radiation and demonstrate the power that NGS data sets hold to resolve even the most difficult of phylogenetic challenges.

The research of W.E. Mayer (1953-2012): a spectrum of immune systems

Over a period of some 20 years, Werner Eugen Mayer played a significant role in establishing a framework for molecular studies of Mhc genes in multiple vertebrates. His work largely concerned gene isolation, sequencing, and related bioinformatic analyses both for the Mhc and for immune system genes of about 200 species, ranging from apes, monkeys, rodents, and marsupials, through to birds, bony fishes, and lampreys. In addition to his exploration of diverse Mhc genes, Werner is remembered for playing a critical role in the development of two important insights into the evolution of immune systems. His was among the first published DNA sequence-based descriptions of trans-species evolution of Mhc alleles, including the first description of the long-lived polymorphisms shared by humans and chimpanzees. This research opened the way for using Mhc polymorphisms in demographic analyses. The second important insight in which he played a prominent role involved the characterization of immune cells and their expressed genes in the lamprey, a jawless vertebrate. His findings helped to indicate the considerable degree to which extant immune mechanisms were co-opted in the creation of the adaptive immune system of jawed vertebrates.

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.

Correlation between nuptial colors and visual sensitivities tuned by opsins leads to species richness in sympatric Lake Victoria cichlid fishes

Reproductive isolation that prevents interspecific hybridization between closely related coexisting species maintains sympatric species diversity. One of the reproductive isolations is mate choice based on color signals (breeding color perceived by color vision). This is well known in several animal taxa, yet little is known about its genetic and molecular mechanism. Lake Victoria cichlid fishes are thought to be an example of sympatric species diversity. In the species inhabiting different light environments in rocky shore, speciation by sensory drive through color signals has been proposed by analyses of the long wavelength-sensitive (LWS) opsin gene and the male nuptial coloration. However, the genetic and molecular mechanism of how diversity of sympatric species occurring in the same habitat is maintained remains unknown. To address this issue, we determined nucleotide sequences of eight opsins of six sympatric species collected from a sandy-muddy shore--an ideal model system for studying sympatric species. Among eight opsins, the LWS and RH1 alleles were diversified and one particular allele is dominant or fixed in each species, and we propose that this is due to natural selection. The functions of their LWS alleles were also diversified as shown by absorption measurements of reconstituted visual pigments. To analyze the relationship between nuptial coloration and the absorption of LWS pigments, we systematically evaluated and defined nuptial coloration. We showed that the coloration was species specific with respect to hue and significantly differentiated by the index values of hue (dominant wavelength: λ(d)). The λ(d) value of the male nuptial coloration correlated with the absorption of LWS pigments from all the species, suggesting that reproductive isolation through mate choice using color signals may prevent sympatric interspecific hybridization, thereby maintaining the species diversity in sympatric species in Lake Victoria.

Ancient origin and recent divergence of a haplochromine cichlid lineage from isolated water bodies in the East African Rift system

Phylogenetic analysis identified haplochromine cichlids from isolated water bodies in the eastern branch of the East African Rift system as an ancient lineage separated from their western sister group in the course of the South Kenyan-North Tanzanian rift system formation. Within this lineage, the close phylogenetic relatedness among taxa indicates a recent common ancestry and historical connections between now separated water bodies. In connection with a total lack of local genetic diversity attributable to population bottlenecks, the data suggest cycles of extinction and colonization in the unstable habitat provided by the small lakes and rivers in this geologically highly active area.

Multilocus phylogeny of Crenicichla (Teleostei: Cichlidae), with biogeography of the C. lacustris group: species flocks as a model for sympatric speciation in rivers

First multilocus analysis of the largest Neotropical cichlid genus Crenicichla combining mitochondrial (cytb, ND2, 16S) and nuclear (S7 intron 1) genes and comprising 602 sequences of 169 specimens yields a robust phylogenetic hypothesis. The best marker in the combined analysis is the ND2 gene which contributes throughout the whole range of hierarchical levels in the tree and shows weak effects of saturation at the 3rd codon position. The 16S locus exerts almost no influence on the inferred phylogeny. The nuclear S7 intron 1 resolves mainly deeper nodes. Crenicichla is split into two main clades: (1) Teleocichla, the Crenicichla wallacii group, and the Crenicichla lugubris-Crenicichla saxatilis groups ("the TWLuS clade"); (2) the Crenicichla reticulata group and the Crenicichla lacustris group-Crenicichla macrophthalma ("the RMLa clade"). Our study confirms the monophyly of the C. lacustris species group with very high support. The biogeographic reconstruction of the C. lacustris group using dispersal-vicariance analysis underlines the importance of ancient barriers between the middle and upper Paraná River (the Guaíra Falls) and between the middle and upper Uruguay River (the Moconá Falls). Our phylogeny recovers two endemic species flocks within the C. lacustris group, the Crenicichla missioneira species flock and the herein discovered Crenicichla mandelburgeri species flock from the Uruguay and Paraná/Iguazú Rivers, respectively. We discuss putative sympatric diversification of trophic traits (morphology of jaws and lips, dentition) and propose these species flocks as models for studying sympatric speciation in complex riverine systems. The possible role of hybridization as a mechanism of speciation is mentioned with a recorded example (Crenicichla scottii).

Reconstructing the evolutionary history of the radiation of the land snail genus Xerocrassa on Crete based on mitochondrial sequences and AFLP markers

BACKGROUND

A non-adaptive radiation triggered by sexual selection resulted in ten endemic land snail species of the genus Xerocrassa on Crete. Only five of these species and a more widespread species are monophyletic in a mitochondrial gene tree. The reconstruction of the evolutionary history of such closely related species can be complicated by incomplete lineage sorting, introgression or inadequate taxonomy. To distinguish between the reasons for the nonmonophyly of several species in the mitochondrial gene tree we analysed nuclear AFLP markers.

RESULTS

Whereas six of the eleven morphologically delimited Xerocrassa species from Crete are monophyletic in the mitochondrial gene tree, nine of these species are monophyletic in the tree based on AFLP markers. Only two morphologically delimited species could not be distinguished with the multilocus data and might have diverged very recently or might represent extreme forms of a single species. The nonmonophyly of X. rhithymna with respect to X. kydonia is probably the result of incomplete lineage sorting, because there is no evidence for admixture in the AFLP data and the mitochondrial haplotype groups of these species coalesce deeply. The same is true for the main haplotype groups of X. mesostena. The nonmonophyly of X. franciscoi might be the result of mitochondrial introgression, because the coalescences of the haplotypes of this species with some X. mesostena haplotypes are shallow and there is admixture with neighbouring X. mesostena.

CONCLUSION

The most likely causes for the nonmonophyly of species in the mitochondrial gene tree of the Xerocrassa radiation on Crete could be inferred using AFLP data by a combination of several criteria, namely the depth of the coalescences in the gene tree, the geographical distribution of shared genetic markers, and concordance with results of admixture analyses of nuclear multilocus markers. The strongly subdivided population structure increases the effective population size of land snail species and, thus, the likelihood of a long persistence of ancestral polymorphisms. Our study suggests that ancestral polymorphisms are a frequent cause for nonmonophyly of species with a strongly subdivided population structure in gene trees.

Genetic variation and demographic history of the Haplochromis laparogramma group of Lake Victoria-An analysis based on SINEs and mitochondrial DNA

More than 500 endemic haplochromine cichlid species inhabit Lake Victoria. This striking species diversity is a classical example of recent explosive adaptive radiation thought to have happened within the last approximately 15,000 years. In this study, we examined the population structure and historical demography of 3 pelagic haplochromine cichlid species that resemble in morphology and have similar niche, Haplochromis (Yssichromis) laparogramma, Haplochromis (Y.) pyrrhocephalus, and Haplochromis (Y.) sp. "glaucocephalus". We investigated the sequences of the mitochondrial DNA control region and the insertion patterns of short interspersed elements (SINEs) of 759 individuals. We show that sympatric forms are genetically differentiated in 4 of 6 cases, but we also found apparent weakening of the genetic differentiation in areas with turbid water. We estimated the timings of population expansion and species divergence to coincide with the refilling of the lake at the Pleistocene/Holocene boundary. We also found that estimates can be altered significantly by the choice of the shape of the molecular clock. If we employ the nonlinear clock model of evolutionary rates in which the rates are higher towards the recent, the population expansion was dated at around the event of desiccation of the lake ca. 17,000 YBP. Thus, we succeeded in clarifying the species and population structure of closely related Lake Victoria cichlids and in showing the importance of applying appropriate clock calibrations in elucidating recent evolutionary events.

Pleistocene desiccation in East Africa bottlenecked but did not extirpate the adaptive radiation of Lake Victoria haplochromine cichlid fishes

The Great Lakes region of East Africa, including Lake Victoria, is the center of diversity of the mega-diverse cichlid fishes (Perciformes: Teleostei). Paleolimnological evidence indicates dramatic desiccation of this lake ca. 18,000-15,000 years ago. Consequently, the hundreds of extant endemic haplochromine species in the lake must have either evolved since then or refugia must have existed, within that lake basin or elsewhere, from which Lake Victoria was recolonized. We studied the population history of the Lake Victoria region superflock (LVRS) of haplochromine cichlids based on nuclear genetic analysis (12 microsatellite loci from 400 haplochomines) of populations from Lake Kivu, Lake Victoria, and the connected and surrounding rivers and lakes. Population genetic analyses confirmed that Lake Kivu haplochromines colonized Lake Victoria. Coalescent analyses show a 30- to 50-fold decline in the haplochromine populations of Lake Victoria, Lake Kivu, and the region ca. 18,000-15,000 years ago. We suggest that this coincides with drastic climatic and geological changes in the late Pleistocene. The most recent common ancestor of the Lake Victoria region haplochromines was estimated to have existed about 4.5 million years ago, which corresponds to the first radiation of cichlids in Lake Tanganyika and the origin of the tribe Haplochrominii. This relatively old evolutionary origin may explain the high levels of polymorphism still found in modern haplochromines. This degree of polymorphism might have acted as a "genetic reservoir" that permitted the explosive radiation of hundreds of haplochromines and their array of contemporary adaptive morphologies.

Population structure of two closely related pelagic cichlids in Lake Victoria, Haplochromis pyrrhocephalus and H. laparogramma

Cichlid fishes in Lake Victoria show spectacular diversification that is thought to be recent. Therefore, by investigating those fishes, we may be able to elucidate recently completed or ongoing speciation processes. We studied the population structures of two closely related pelagic cichlid species, Haplochromis pyrrhocephalus and H. laparogramma, using a mitochondrial DNA locus and 12 nuclear microsatellite loci as putative neutral markers. Ten and two populations of H. pyrrhocephalus and H. laparogramma, respectively, were sampled from the southern part of Lake Victoria. We grouped those 12 populations into four mutually differentiated regional populations, one of which consisted of the two H. laparogramma populations. The levels of differentiation were substantial at the mitochondrial locus (F(ST) = 0.03-0.54), but very low at microsatellite loci (R(ST) = 0.008-0.116). The data from both types of loci indicated that the regional population of H. laparogramma was first separated from those of H. pyrrhocephalus if we set aside one erratic population of H. pyrrhocephalus. The data also suggested recent population expansions of the two species, the time scales for which were estimated to be on the order of 10(4)-10(5) years. These data suggested that dynamic speciation processes accompanied occasional spawning of new species and population size changes in this lake.

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.

Non-equilibrium estimates of gene flow inferred from nuclear genealogies suggest that Iberian and North African wall lizards (Podarcis spp.) are an assemblage of incipient species

Background

The study of recently-diverged species offers significant challenges both in the definition of evolutionary entities and in the estimation of gene flow among them. Iberian and North African wall lizards (Podarcis) constitute a cryptic species complex for which previous assessments of mitochondrial DNA (mtDNA) and allozyme variation are concordant in describing the existence of several highly differentiated evolutionary units. However, these studies report important differences suggesting the occurrence of gene flow among forms. Here we study sequence variation in two nuclear introns, β-fibint7 and 6-Pgdint7, to further investigate overall evolutionary dynamics and test hypotheses related to species delimitation within this complex.

Results

Both nuclear gene genealogies fail to define species as monophyletic. To discriminate between the effects of incomplete lineage sorting and gene flow in setting this pattern, we estimated migration rates among species using both F_ST-based estimators of gene flow, which assume migration-drift equilibrium, and a coalescent approach based on a model of divergence with gene flow. Equilibrium estimates of gene flow suggest widespread introgression between species, but coalescent estimates describe virtually zero admixture between most (but not all) species pairs. This suggests that although gene flow among forms may have occurred the main cause for species polyphyly is incomplete lineage sorting, implying that most forms have been isolated since their divergence. This observation is therefore in accordance with previous reports of strong differentiation based on mtDNA and allozyme data.

Conclusion

These results corroborate most forms of Iberian and North African Podarcis as differentiated, although incipient, species, supporting a gradual view of speciation, according to which species may persist as distinct despite some permeability to genetic exchange and without having clearly definable genetic boundaries. Additionally, this study constitutes a warning against the misuse of equilibrium estimates of migration among recently-diverged groups.

Evolutionary history and biogeographic affinities of the serranochromine cichlids in Zambian rivers

Zambian rivers belong to two major drainages: the Luapula-Congo and the Zambezi River system. Tectonic activity repeatedly altered drainage systems in Africa, so that current fish faunas can only be understood in the context of historic drainage capture events. We use phylogenetic relationships of one widespread lineage of haplochromine cichlids, the serranochromines, to trace their biogeographic expansion and diversification in Zambia. The mitochondrial DNA phylogeny suggests five ancient clades, and their common ancestor was likely to have invaded from the Lower Congo River. The branching intervals in the linearized tree analysis suggest three major cladogenesis events and two periods of faunal exchange. The five clades originating in the Congo River drainage diversified further; one stayed in the Congo drainage, one diversified in the Zambezi system only, while the three clades of mixed distribution underwent diversification in the Zambezi system, to re-enter the Congo drainage very recently, as indicated by sometimes zero mutation differences among taxa from different drainages. Our hypothesis is consistent with the suggested radiation in the extinct Lake palaeo-Makgadikgadi, so that we propose that the Zambian serranochromine fauna in part represents survivors of the extinct lacustrine flock plus several novel species that originated in situ.

Nuclear and mitochondrial data reveal different evolutionary processes in the Lake Tanganyika cichlid genus Tropheus

BACKGROUND

Cichlid fishes are notorious for their wealth of intra- and interspecific colour pattern diversity. In Lake Tanganyika, the endemic genus Tropheus represents the most impressive example for geographic variation in the pattern and hue of integument colouration, but the taxonomy of the over 100 mostly allopatric colour morphs remains to a large degree unresolved. Previous studies of mitochondrial DNA sequence data revealed polyphyly of the six nominally described species and complex phylogeographic patterns influenced by lake level fluctuations and population admixture, and suggested the parallel evolution of similar colour patterns in divergent evolutionary lineages. A gene tree of a rapidly radiating group may be subject to incomplete and stochastic lineage sorting, and to overcome this problem we used multi-locus, nuclear AFLP data in comparison with mtDNA sequences to study diversification, migration and introgression in Tropheus colour morphs in Lake Tanganyika.

RESULTS

Significant incongruence between phylogenetic reconstructions from mitochondrial and AFLP data suggested incomplete sorting of mitochondrial haplotypes as well as frequent introgression between differentiated lineages. In contrast to the mitochondrial phylogeny, the AFLP phenogram was largely congruent with species classifications, colour pattern similarities, and in many cases also with the current geographic distribution of populations, and did not produce evidence of convergent colour pattern evolution. Homoplasy in the AFLP data was used to identify populations that were strongly affected by introgression.

CONCLUSION

Different evolutionary processes were distinguished by the combination of mitochondrial and AFLP data. Mitochondrial phylogeographic patterns retained signals of large-scale migration events triggered by historical, major lake level fluctuations, whereas AFLP data indicated genetic cohesion among local groups of populations resulting from secondary contact of adjacent populations in the course of the more frequently occurring, minor lake level fluctuations. There was no support for the parallel evolution of similar colour patterns in the AFLP data. Genetic signatures of introgression and hybridisation detected in several populations suggest that lake level fluctuations drove the stunning diversification of Tropheus morphs not only through population fragmentation, but also by promoting hybridisation between differentiated morphs in secondary contact.

Gene flow between species of Lake Victoria haplochromine fishes

The haplochromine cichlid fishes of Lake Victoria (LV), East Africa, are a textbook example of adaptive radiation-a rapid divergence of multiple morphologically distinguishable forms from a few founding lineages. The forms are generally believed to constitute a "flock" of several hundred reproductively isolated species in a dozen or so genera. This belief has, until now, not been subjected to a test, however. Here, we compare genetic variation at 11 loci in 10 haplochromine populations of 6 different species. Although the genetic diversity in the populations is quite high, using a variety of statistical tests, we find no evidence of genetic differentiation among the populations of LV haplochromines. On genetic distance trees, populations of the same species intermingle with those of different species. At the molecular level, the species are indistinguishable from one another. Genetic comparisons with closely related species in 2 crater lakes indicate that the species within LV continue exchanging genes. These observations have important implications for phylogenetic reconstruction. The approach used in this study is applicable to other instances of adaptive radiation.

SINEs of speciation: tracking lineages with retroposons

The value of short interspersed elements (SINEs) for diagnosing common ancestry is being expanded to examine the differential sorting of lineages through the course of speciation events. Because most SINEs are neutral markers of identical descent, are not precisely excised from the genome and have a known ancestral condition, they are advantageous for reconciling gene trees and species trees with minimal phylogenetic error. A population perspective on SINE evolution combined with coalescence theory provides a context for investigating the phenomenon of ancestral polymorphism and its role in producing incongruent SINE insertion patterns among multiple loci. Studies of human Alu repeats demonstrate the value of young polymorphic SINEs for assessing human genomic diversity and tracking ancient demographics of human populations, whereas incongruent insertion patterns revealed by older fixed SINE loci, such as those in African cichlid fishes, contain information that might help identify ancient radiations that are otherwise obscured by accumulated mutations in sequence data. Here, we review the utility of retroposons for inferring common ancestry, discuss limits to the method, and clarify confusion by providing examples from the literature that illustrate how discordant multi-locus insertion patterns of retroelements can indicate lineage-sorting events that should not be misinterpreted as phylogenetic noise.

Widespread Genealogical Nonmonophyly in Species of Pinus Subgenus Strobus

Phylogenetic relationships among Pinus species from subgenus Strobus remain unresolved despite combined efforts based on nrITS and cpDNA. To provide greater resolution among these taxa, a 900-bp intron from a late embryogenesis abundant (LEA)-like gene (IFG8612)was sequenced from 39 pine species, with two or more alleles representing 33 species. Nineteen of 33 species exhibited allelic nonmonphyly in the strict consensus tree, and 10 deviated significantly from allelic monophyly based on topology incongruence tests. Intraspecific nucleotide diversity ranged from 0.0 to 0.0211, and analysis of variance shows that nucleotide diversity was strongly associated (P < 0.0001)with the degree of species monophyly. Although species nonmonophyly complicates phylogenetic interpretations, this nuclear locus offers greater topological support than previously observed for cpDNA or nrITS. Lacking evidence for hybridization, recombination, or imperfect taxonomy, we feel that incomplete lineage sorting remains the best explanation for the polymorphisms shared among species. Depending on the species, coalescent expectations indicate that reciprocal monophyly will be more likely than paraphyly in 1.71 to 24.0 x 10(6) years, and that complete genome-wide coalescence in these species may require up to 76.3 x 10(6) years. The absence of allelic coalescence is a severe constraint in the application of phylogenetic methods in Pinus, and taxa sharing similar life history traits with Pinus are likely to show species nonmonophyly using nuclear markers.

Evolution of Pleopsidium (Lichenized Ascomycota) S943 Group I Introns and the Phylogeography of an Intron-Encoded Putative Homing Endonuclease

The sporadic distribution of nuclear group I introns among different fungal lineages can be explained by vertical inheritance of the introns followed by successive losses, or horizontal transfers from one lineage to another through intron homing or reverse splicing. Homing is mediated by an intron-encoded homing endonuclease (HE) and recent studies suggest that the introns and their associated HE gene (HEG) follow a recurrent cyclical model of invasion, degeneration, loss, and reinvasion. The purpose of this study was to compare this model to the evolution of HEGs found in the group I intron at position S943 of the nuclear ribosomal DNA of the lichen-forming fungus Pleopsidium. Forty-eight S943 introns were found in the 64 Pleopsidium samples from a worldwide screen, 22 of which contained a full-length HEG that encodes a putative 256-amino acid HE, and 2 contained HE pseudogenes. The HEGs are divided into two closely related types (as are the introns that encode them) that differ by 22.6% in their nucleotide sequences. The evolution of the Pleopsidium intron-HEG element shows strong evidence for a cyclical model of evolution. The intron was likely acquired twice in the genus and then transmitted via two or three interspecific horizontal transfers. Close geographical proximity plays an important role in intron-HEG horizontal transfer because most of these mobile elements were found in Europe. Once acquired in a lineage, the intron-HEG element was also vertically transmitted, and occasionally degenerated or was lost.

Reticulate phylogeny of gastropod-shell-breeding cichlids from Lake Tanganyika--the result of repeated introgressive hybridization

Background

The tribe Lamprologini is the major substrate breeding lineage of Lake Tanganyika's cichlid species flock. Among several different life history strategies found in lamprologines, the adaptation to live and breed in empty gastropod shells is probably the most peculiar. Although shell-breeding arose several times in the evolutionary history of the lamprologines, all obligatory and most facultative shell-breeders belong to the so called "ossified group", a monophyletic lineage within the lamprologine cichlids. Since their distinctive life style enables these species to live and breed in closest vicinity, we hypothesized that these cichlids might be particularly prone to accidental hybridization, and that introgression might have affected the evolutionary history of this cichlid lineage.

Results

Our analyses revealed discrepancies between phylogenetic hypotheses based on mitochondrial and nuclear (AFLP) data. While the nuclear phylogeny was congruent with morphological, behavioral and ecological characteristics, several species – usually highly specialized shell-breeders – were placed at contradicting positions in the mitochondrial phylogeny. The discordant phylogenies strongly suggest repeated incidents of introgressive hybridization between several distantly related shell-breeding species, which reticulated the phylogeny of this group of cichlids. Long interior branches and high bootstrap support for many interior nodes in the mitochondrial phylogeny argue against a major effect of ancient incomplete lineage sorting on the phylogenetic reconstruction. Moreover, we provide morphological and genetic (mtDNA and microsatellites) evidence for ongoing hybridization among distantly related shell-breeders. In these cases, the territorial males of the inferred paternal species are too large to enter the shells of their mate, such that they have to release their sperm over the entrance of the shell to fertilize the eggs. With sperm dispersal by water currents and wave action, trans-specific fertilization of clutches in neighboring shells seem inevitable, when post-zygotic isolation is incomplete.

Conclusion

From the direct observation of hybrids we conclude that hybridization between distantly related gastropod-shell-breeding cichlids of Lake Tanganyika follows inevitably from their ecological specialization. Moreover, the observed incongruence between mtDNA and nuclear multilocus phylogeny suggests that repeated hybridization events among quite distantly related taxa affected the diversification of this group, and introduced reticulation into their phylogeny.

Evolution of the complementary sex-determination gene of honey bees: balancing selection and trans-species polymorphisms

The mechanism of sex determination varies substantively among evolutionary lineages. One important mode of genetic sex determination is haplodiploidy, which is used by approximately 20% of all animal species, including >200,000 species of the entire insect order Hymenoptera. In the honey bee Apis mellifera, a hymenopteran model organism, females are heterozygous at the csd (complementary sex determination) locus, whereas males are hemizygous (from unfertilized eggs). Fertilized homozygotes develop into sterile males that are eaten before maturity. Because homozygotes have zero fitness and because common alleles are more likely than rare ones to form homozygotes, csd should be subject to strong overdominant selection and negative frequency-dependent selection. Under these selective forces, together known as balancing selection, csd is expected to exhibit a high degree of intraspecific polymorphism, with long-lived alleles that may be even older than the species. Here we sequence the csd genes as well as randomly selected neutral genomic regions from individuals of three closely related species, A. mellifera, Apis cerana, and Apis dorsata. The polymorphic level is approximately seven times higher in csd than in the neutral regions. Gene genealogies reveal trans-species polymorphisms at csd but not at any neutral regions. Consistent with the prediction of rare-allele advantage, nonsynonymous mutations are found to be positively selected in csd only in early stages after their appearances. Surprisingly, three different hypervariable repetitive regions in csd are present in the three species, suggesting variable mechanisms underlying allelic specificities. Our results provide a definitive demonstration of balancing selection acting at the honey bee csd gene, offer insights into the molecular determinants of csd allelic specificities, and help avoid homozygosity in bee breeding.

African cichlid fish: a model system in adaptive radiation research

The African cichlid fish radiations are the most diverse extant animal radiations and provide a unique system to test predictions of speciation and adaptive radiation theory. The past few years have seen major advances in the phylogenetics, evolutionary biogeography and ecology of cichlid fish. Most of this work has concentrated on the most diverse radiations. Unfortunately, a large number of small radiations and 'non-radiations' have been overlooked, potentially limiting the contribution of the cichlid system to our understanding of speciation and adaptive radiation. I have reviewed the literature to identify 33 intralacustrine radiations and 76 failed radiations. For as many as possible I collected information on lake size, age and phylogenetic relationships. I use these data to address two questions: (i) whether the rate of speciation and the resulting species richness are related to temporal and spatial variation in ecological opportunity and (ii) whether the likelihood of undergoing adaptive radiation is similar for different African cichlid lineages. The former is a key prediction of the ecological theory of adaptive radiation that has been presumed true but remains untested for cichlid radiations. The second is based on the hypothesis that the propensity of cichlids to radiate is due to a key evolutionary innovation shared by all African cichlids. The evidence suggests that speciation rate declines through time as niches get filled up during adaptive radiation: young radiations and early stages of old radiations are characterized by high rates of speciation, whereas at least 0.5 Myr into a radiation speciation becomes a lot less frequent. The number of species in cichlid radiations increases with lake size, supporting the prediction that species diversity increases with habitat heterogeneity, but also with opportunity for isolation by distance. Finally, the data suggest that the propensity to radiate within lakes is a derived property that evolved during the evolutionary history of some African cichlids, and the appearance of which does not coincide with the appearance of proposed key innovations in morphology and life history.