Annotation of expressed sequence tags for the East African cichlid fish Astatotilapia burtoni and evolutionary analyses of cichlid ORFs

The puzzling phylogeography of the haplochromine cichlid fish Astatotilapia burtoni

Astatotilapia burtoni is a member of the "modern haplochromines," the most species-rich lineage within the family of cichlid fishes. Although the species has been in use as research model in various fields of research since almost seven decades, including developmental biology, neurobiology, genetics and genomics, and behavioral biology, little is known about its spatial distribution and phylogeography. Here, we examine the population structure and phylogeographic history of A. burtoni throughout its entire distribution range in the Lake Tanganyika basin. In addition, we include several A. burtoni laboratory strains to trace back their origin from wild populations. To this end, we reconstruct phylogenetic relationships based on sequences of the mitochondrial DNA (mtDNA) control region (d-loop) as well as thousands of genomewide single nucleotide polymorphisms (SNPs) derived from restriction-associated DNA sequencing. Our analyses reveal high population structure and deep divergence among several lineages, however, with discordant nuclear and mtDNA phylogenetic inferences. Whereas the SNP-based phylogenetic hypothesis uncovers an unexpectedly deep split in A. burtoni, separating the populations in the southern part of the Lake Tanganyika basin from those in the northern part, analyses of the mtDNA control region suggest deep divergence between populations from the southwestern shoreline and populations from the northern and southeastern shorelines of Lake Tanganyika. This phylogeographic pattern and mitochondrial haplotype sharing between populations from the very North and the very South of Lake Tanganyika can only partly be explained by introgression linked to lake-level fluctuations leading to past contact zones between otherwise isolated populations and large-scale migration events.

Gene expression signatures of mating system evolution

The diversity of mating systems among animals is astounding. Importantly, similar mating systems have evolved even across distantly related taxa. However, our understanding of the mechanisms underlying these convergently evolved phenotypes is limited. Here, we examine on a genomic scale the neuromolecular basis of social organization in cichlids of the tribe Ectodini from Lake Tanganyika. Using field-collected males and females of four closely related species representing two independent evolutionary transitions from polygyny to monogamy, we take a comparative transcriptomic approach to test the hypothesis that these independent transitions have recruited similar gene sets. Our results demonstrate that while lineage and species exert a strong influence on neural gene expression profiles, social phenotype can also drive gene expression evolution. Specifically, 331 genes (∼6% of those assayed) were associated with monogamous mating systems independent of species or sex. Among these genes, we find a strong bias (4:1 ratio) toward genes with increased expression in monogamous individuals. A highly conserved nonapeptide system known to be involved in the regulation of social behavior across animals was not associated with mating system in our analysis. Overall, our findings suggest deep molecular homologies underlying the convergent or parallel evolution of monogamy in different cichlid lineages of Ectodini.

Parental investment matters for maternal and offspring immune defense in the mouthbrooding cichlid Astatotilapia burtoni

BACKGROUND

Parental care, while increasing parental fitness through offspring survival, also bears cost to the care-giving parent. Consequentially, trade offs between parental care and other vitally important traits, such as the immune system seem evident. In co-occurring phases of parental care and immunological challenges negative consequences through a resource allocation trade off on both the parental and the offspring conditions can be predicted. While the immune system reflects parental stress conditions, parental immunological investments also boost offspring survival via the transfer of immunological substances (trans-generational immune priming). We investigated this relationship in the mouthbrooding East African cichlid Astotatilapia burtoni. Prior to mating, females were exposed to an immunological activation, while others remained immunologically naïve. Correspondingly, the immunological status of females was either examined directly after reproduction or after mouthbrooding had ceased. Offspring from both groups were exposed to immunological challenges to assess the extent of trans-generational immune priming. As proxy for immune status, cellular immunological activity and gene expression were determined.

RESULTS

Both reproducing and mouthbrooding females allocate their resources towards reproduction. While upon reproduction the innate immune system was impeded, mouthbrooding females showed an attenuation of inflammatory components. Juveniles from immune challenged mouthbrooding females showed downregulation of immune and life history candidate genes, implying a limitation of trans-generational plasticity when parents experience stress during the costly reproductive phase.

CONCLUSION

Our results provide evidence that both parental investment via mouthbrooding and the rise of the immunological activity upon an immune challenge are costly traits. If applied simultaneously, not only mothers seem to be impacted in their performance, but also offspring are impeded in their ability to react upon a potentially virulent pathogen exposure.

Variations on a theme: Genomics of sex determination in the cichlid fish Astatotilapia burtoni

Background

Sex chromosomes change more frequently in fish than in mammals or birds. However, certain chromosomes or genes are repeatedly used as sex determinants in different members of the teleostean lineage. East African cichlids are an enigmatic model system in evolutionary biology representing some of the most diverse extant vertebrate adaptive radiations. How sex is determined and if different sex-determining mechanisms contribute to speciation is unknown for almost all of the over 1,500 cichlid species of the Great Lakes. Here, we investigated the genetic basis of sex determination in a cichlid from Lake Tanganyika, Astatotilapia burtoni, a member of the most species-rich cichlid lineage, the haplochromines.

Results

We used RAD-sequencing of crosses for two populations of A. burtoni, a lab strain and fish caught at the south of Lake Tanganyika. Using association mapping and comparative genomics, we confirmed male heterogamety in A. burtoni and identified different sex chromosomes (LG5 and LG18) in the two populations of the same species. LG5, the sex chromosome of the lab strain, is a fusion chromosome in A. burtoni. Wnt4 is located on this chromosome, representing the best candidate identified so far for the master sex-determining gene in our lab strain of A. burtoni.

Conclusions

Cichlids exemplify the high turnover rate of sex chromosomes in fish with two different chromosomes, LG5 and LG18, containing major sex-determining loci in the two populations of A. burtoni examined here. However, they also illustrate that particular chromosomes are more likely to be used as sex chromosomes. Chromosome 5 is such a chromosome, which has evolved several times as a sex chromosome, both in haplochromine cichlids from all Great Lakes and also in other teleost fishes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3178-0) contains supplementary material, which is available to authorized users.

Using an integrative approach to investigate the evolution of behaviour

Behaviour is a central focus of interest in biology because it has an impact on several aspects of an organism's life. Evolutionary biologists have realised the advantage of an integrative approach that jointly studies the molecular, cellular and physiological levels of an individual to link them with the organismal behavioural phenotype. First, this mechanistic information helps in understanding physiological and evolutionary constraints acting on the behavioural response to the environment and its evolution. Second, it furthers our understanding of the process of molecular convergent evolution. Finally, we learn about natural variation in molecular, cellular and physiological traits present in wild populations and their underlying genetic basis, which can be a substrate for selection to act on. I illustrate these points using our work on behaviour variation in fishes. The information on the mechanistic bases of behaviour variation in various species and behaviours will contribute to an ecological annotation of genes and to uncover new mechanisms implicated in how this astonishing behavioural diversity arose, is maintained and will evolve.

Social odors conveying dominance and reproductive information induce rapid physiological and neuromolecular changes in a cichlid fish

Background

Social plasticity is a pervasive feature of animal behavior. Animals adjust the expression of their social behavior to the daily changes in social life and to transitions between life-history stages, and this ability has an impact in their Darwinian fitness. This behavioral plasticity may be achieved either by rewiring or by biochemically switching nodes of the neural network underlying social behavior in response to perceived social information. Independent of the proximate mechanisms, at the neuromolecular level social plasticity relies on the regulation of gene expression, such that different neurogenomic states emerge in response to different social stimuli and the switches between states are orchestrated by signaling pathways that interface the social environment and the genotype. Here, we test this hypothesis by characterizing the changes in the brain profile of gene expression in response to social odors in the Mozambique Tilapia, Oreochromis mossambicus. This species has a rich repertoire of social behaviors during which both visual and chemical information are conveyed to conspecifics. Specifically, dominant males increase their urination frequency during agonist encounters and during courtship to convey chemical information reflecting their dominance status.

Results

We recorded electro-olfactograms to test the extent to which the olfactory epithelium can discriminate between olfactory information from dominant and subordinate males as well as from pre- and post-spawning females. We then performed a genome-scale gene expression analysis of the olfactory bulb and the olfactory cortex homolog in order to identify the neuromolecular systems involved in processing these social stimuli.

Conclusions

Our results show that different olfactory stimuli from conspecifics’ have a major impact in the brain transcriptome, with different chemical social cues eliciting specific patterns of gene expression in the brain. These results confirm the role of rapid changes in gene expression in the brain as a genomic mechanism underlying behavioral plasticity and reinforce the idea of an extensive transcriptional plasticity of cichlid genomes, especially in response to rapid changes in their social environment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1255-4) contains supplementary material, which is available to authorized users.

Characterizing the transcriptome of yellow-cheek carp (Elopichthys bambusa) enables evolutionary analyses within endemic East Asian Cyprinidae

The identification of genes that may be responsible for the divergence of closely related species is one of the central goals of evolutionary biology. The species of endemic East Asian Cyprinidae diverged less than 8millionyears ago, and the morphological differences among these species are great. However, the genetic basis of their divergence remains unknown. In this report, we investigated the transcriptome of one endemic East Asian cyprinid - the yellow-cheek carp Elopichthys bambusa. A comparison with the publicly available transcriptomes of other endemic East Asian cyprinids, including the silver carp (Hypophthalmichthys molitrix) and blunt-nose black bream (Megalobrama amblycephala), revealed a number of candidate adaptive genes in each species, such as zona pellucida glycoprotein 2 in E. bambusa and zebrafish vitelline envelope protein in M. amblycephala. An enrichment test showed the enrichment of some specific gene ontology (GO) terms for these putatively adaptive genes. Taken together, our work is the first step toward elucidating the genes that may be related to the divergence of endemic East Asian Cyprinidae, and these genes identified as being probably under positive selection should be good candidates for subsequent evolutionary and functional studies.

Evolution of genomic structural variation and genomic architecture in the adaptive radiations of African cichlid fishes

African cichlid fishes are an ideal system for studying explosive rates of speciation and the origin of diversity in adaptive radiation. Within the last few million years, more than 2000 species have evolved in the Great Lakes of East Africa, the largest adaptive radiation in vertebrates. These young species show spectacular diversity in their coloration, morphology and behavior. However, little is known about the genomic basis of this astonishing diversity. Recently, five African cichlid genomes were sequenced, including that of the Nile Tilapia (Oreochromis niloticus), a basal and only relatively moderately diversified lineage, and the genomes of four representative endemic species of the adaptive radiations, Neolamprologus brichardi, Astatotilapia burtoni, Metriaclima zebra, and Pundamila nyererei. Using the Tilapia genome as a reference genome, we generated a high-resolution genomic variation map, consisting of single nucleotide polymorphisms (SNPs), short insertions and deletions (indels), inversions and deletions. In total, around 18.8, 17.7, 17.0, and 17.0 million SNPs, 2.3, 2.2, 1.4, and 1.9 million indels, 262, 306, 162, and 154 inversions, and 3509, 2705, 2710, and 2634 deletions were inferred to have evolved in N. brichardi, A. burtoni, P. nyererei, and M. zebra, respectively. Many of these variations affected the annotated gene regions in the genome. Different patterns of genetic variation were detected during the adaptive radiation of African cichlid fishes. For SNPs, the highest rate of evolution was detected in the common ancestor of N. brichardi, A. burtoni, P. nyererei, and M. zebra. However, for the evolution of inversions and deletions, we found that the rates at the terminal taxa are substantially higher than the rates at the ancestral lineages. The high-resolution map provides an ideal opportunity to understand the genomic bases of the adaptive radiation of African cichlid fishes.

The evolutionary genomics of cichlid fishes: explosive speciation and adaptation in the postgenomic era

With more than 1,500 species, cichlid fishes provide textbook examples of recent and diverse adaptive radiations, rapid rates of speciation, and the parallel evolution of adaptive phenotypes among both recently and distantly related lineages. This extraordinary diversity has attracted considerable interest from researchers across several biological disciplines. Their broad phenotypic variation coupled with recent divergence makes cichlids an ideal model system for understanding speciation, adaptation, and phenotypic diversification. Genetic mapping, genome-wide analyses, and genome projects have flourished in the past decade and have added new insights on the question of why there are so many cichlids. These recent findings also show that the sharing of older DNA polymorphisms is extensive and suggest that linage sorting is incomplete and that adaptive introgression played a role in the African radiation. Here, we review the results of genetic and genomic research on cichlids in the past decade and suggest some potential avenues to further exploit the potential of the cichlid model system to provide a better understanding of the genomics of adaptation and speciation.

Phenotype-environment association of the oxygen transport system in trimorphic European whitefish (Coregonus lavaretus) populations

Replicated adaptive radiation events, typified by phenotypic divergence across resource axes, provide important insight into the eco-evolutionary dynamics that lead to the formation of new species. Here, we show that in trimorphic adaptive radiations of European whitefish (Coregonus lavaretus), divergence of the oxygen transport system has occurred across the pelagic/littoral (shallow)-profundal (deep) resource axis, and at multiple biological scales. Profundal whitefish exhibited significantly larger red blood cells (RBCs), a greater proportion of cathodic hemoglobin protein components, and higher hemoglobin transcript abundance in kidney compared to littoral and pelagic morphs. Hemoglobin transcript abundance in brain and gill, but not kidney, and anodic hemoglobin protein component diversity in blood were also linked to variation at an intronic single nucleotide polymorphism (SNP). As the whitefish morphs differ in population genetic structure at this SNP, hemoglobin transcript and protein divergence between profundal and pelagic/littoral morphs is likely being driven by genetic divergence. Our findings, along with our previous work on lake whitefish, highlight the importance of the oxygen transport system to the postglacial colonization of novel lacustrine environments by whitefish throughout the northern hemisphere.

Gene duplication in an African cichlid adaptive radiation

BACKGROUND

Gene duplication is a source of evolutionary innovation and can contribute to the divergence of lineages; however, the relative importance of this process remains to be determined. The explosive divergence of the African cichlid adaptive radiations provides both a model for studying the general role of gene duplication in the divergence of lineages and also an exciting foray into the identification of genomic features that underlie the dramatic phenotypic and ecological diversification in this particular lineage. We present the first genome-wide study of gene duplication in African cichlid fishes, identifying gene duplicates in three species belonging to the Lake Malawi adaptive radiation (Metriaclima estherae, Protomelas similis, Rhamphochromis "chilingali") and one closely related species from a non-radiated riverine lineage (Astatotilapia tweddlei).

RESULTS

Using Astatotilapia burtoni as reference, microarray comparative genomic hybridization analysis of 5689 genes reveals 134 duplicated genes among the four cichlid species tested. Between 51 and 55 genes were identified as duplicated in each of the three species from the Lake Malawi radiation, representing a 38%-49% increase in number of duplicated genes relative to the non-radiated lineage (37 genes). Duplicated genes include several that are involved in immune response, ATP metabolism and detoxification.

CONCLUSIONS

These results contribute to our understanding of the abundance and type of gene duplicates present in cichlid fish lineages. The duplicated genes identified in this study provide candidates for the analysis of functional relevance with regard to phenotype and divergence. Comparative sequence analysis of gene duplicates can address the role of positive selection and adaptive evolution by gene duplication, while further study across the phylogenetic range of cichlid radiations (and more generally in other adaptive radiations) will determine whether the patterns of gene duplication seen in this study consistently accompany rapid radiation.

Tol2-mediated generation of a transgenic haplochromine cichlid, Astatotilapia burtoni

Cichlid fishes represent one of the most species-rich and rapid radiations of a vertebrate family. These ~2200 species, predominantly found in the East African Great Lakes, exhibit dramatic differences in anatomy, physiology, and behavior. However, the genetic bases for this radiation, and for the control of their divergent traits, are unknown. A flood of genomic and transcriptomic data promises to suggest mechanisms underlying the diversity, but transgenic technology will be needed to rigorously test the hypotheses generated. Here we demonstrate the successful use of the Tol2 transposon system to generate transgenic Astatotilapia burtoni, a haplochromine cichlid from Lake Tanganyika, carrying the GFP transgene under the control of the ubiquitous EF1α promoter. The transgene integrates into the genome, is successfully passed through the germline, and the widespread GFP expression pattern is stable across siblings and multiple generations. The stable inheritance and expression patterns indicate that the Tol2 system can be applied to generate A. burtoni transgenic lines. Transgenesis has proven to be a powerful technology for manipulating genes and cells in other model organisms and we anticipate that transgenic A. burtoni and other cichlids will be used to test the mechanisms underlying behavior and speciation.

Female preference for males depends on reproductive physiology in the African cichlid fish Astatotilapia burtoni

Mate choice is fundamental to sexual selection, yet little is known about underlying physiological mechanisms that influence female mating decisions. We investigated the endocrine underpinnings of female mate choice in the African cichlid Astatotilapia burtoni, a non-seasonal breeder. In addition to profiling behavioral and hormonal changes across the female reproductive cycle, we tested two hypotheses regarding possible factors influencing female mate choice. We first asked whether female mate choice is influenced by male visual and/or chemical cues. A. burtoni females were housed for one full reproductive cycle in the center of a dichotomous choice apparatus with a large (attractive) or small (unattractive) conspecific male on either side. Females associated mostly with small, less attractive males, but on the day of spawning reversed their preference to large, attractive males, with whom they mated almost exclusively, although this choice depended on the relative amount of androgens released into the water by small males. We next asked whether male behavior or androgen levels change in relation to the stimulus females' reproductive state. We found that stimulus male aggression decreased and reproductive displays increased as the day of spawning approached. Moreover male testosterone levels changed throughout the females' reproductive cycle, with larger males releasing more testosterone into the water than small males. Our data suggest that female association in a dichotomous choice assay is only indicative of the actual mate choice on the day of spawning. Furthermore, we show that male behavior and hormone levels are dependent on the reproductive state of conspecific females.

The ecological and genetic basis of convergent thick-lipped phenotypes in cichlid fishes

The evolution of convergent phenotypes is one of the most interesting outcomes of replicate adaptive radiations. Remarkable cases of convergence involve the thick-lipped phenotype found across cichlid species flocks in the East African Great Lakes. Unlike most other convergent forms in cichlids, which are restricted to East Africa, the thick-lipped phenotype also occurs elsewhere, for example in the Central American Midas Cichlid assemblage. Here, we use an ecological genomic approach to study the function, the evolution and the genetic basis of this phenotype in two independent cichlid adaptive radiations on two continents. We applied phylogenetic, demographic, geometric morphometric and stomach content analyses to an African (Lobochilotes labiatus) and a Central American (Amphilophus labiatus) thick-lipped species. We found that similar morphological adaptations occur in both thick-lipped species and that the 'fleshy' lips are associated with hard-shelled prey in the form of molluscs and invertebrates. We then used comparative Illumina RNA sequencing of thick vs. normal lip tissue in East African cichlids and identified a set of 141 candidate genes that appear to be involved in the morphogenesis of this trait. A more detailed analysis of six of these genes led to three strong candidates: Actb, Cldn7 and Copb. The function of these genes can be linked to the loose connective tissue constituting the fleshy lips. Similar trends in gene expression between African and Central American thick-lipped species appear to indicate that an overlapping set of genes was independently recruited to build this particular phenotype in both lineages.

Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish

BACKGROUND

The availability of a large number of recently sequenced vertebrate genomes opens new avenues to integrate cytogenetics and genomics in comparative and evolutionary studies. Cytogenetic mapping can offer alternative means to identify conserved synteny shared by distinct genomes and also to define genome regions that are still not fine characterized even after wide-ranging nucleotide sequence efforts. An efficient way to perform comparative cytogenetic mapping is based on BAC clones mapping by fluorescence in situ hybridization. In this report, to address the knowledge gap on the genome evolution in cichlid fishes, BAC clones of an Oreochromis niloticus library covering the linkage groups (LG) 1, 3, 5, and 7 were mapped onto the chromosomes of 9 African cichlid species. The cytogenetic mapping data were also integrated with BAC-end sequences information of O. niloticus and comparatively analyzed against the genome of other fish species and vertebrates.

RESULTS

The location of BACs from LG1, 3, 5, and 7 revealed a strong chromosomal conservation among the analyzed cichlid species genomes, which evidenced a synteny of the markers of each LG. Comparative in silico analysis also identified large genomic blocks that were conserved in distantly related fish groups and also in other vertebrates.

CONCLUSIONS

Although it has been suggested that fishes contain plastic genomes with high rates of chromosomal rearrangements and probably low rates of synteny conservation, our results evidence that large syntenic chromosome segments have been maintained conserved during evolution, at least for the considered markers. Additionally, our current cytogenetic mapping efforts integrated with genomic approaches conduct to a new perspective to address important questions involving chromosome evolution in fishes.

A novel primer set for multilocus phylogenetic inference in East African cichlid fishes

The cichlid fishes in the East African Great Lakes are a prime model system for the study of adaptive radiation. Therefore, the availability of an elaborate phylogenetic framework is an important prerequisite. Previous phylogenetic hypotheses on East African cichlids are mainly based on mitochondrial and/or fragment-based markers, and, to date, no taxon-rich phylogeny exists that is based on multilocus DNA sequence data. Here, we present the design of an extensive new primer set (24 nuclear makers) for East African cichlids that will be used for multilocus phylogenetic analyses in the future. The primers are designed to work for both Sanger sequencing and next-generation sequencing with the 454 technology. As a proof of principle, we validate these primers in a phylogenetically representative set of 16 cichlid species from Lake Tanganyika and main river systems in the area and provide a basic evaluation of the markers with respect to marker length and diversity indices.

Genomics of adaptation and speciation in cichlid fishes: recent advances and analyses in African and Neotropical lineages

Cichlid fishes are remarkably phenotypically diverse and species-rich. Therefore, they provide an exciting opportunity for the study of the genetics of adaptation and speciation by natural and sexual selection. Here, we review advances in the genomics and transcriptomics of cichlids, particularly regarding ecologically relevant differences in body shape, trophic apparatus, coloration and patterning, and sex determination. Research conducted so far has focused almost exclusively on African cichlids. To analyse genomic diversity and selection in a Neotropical radiation, we conducted a comparative transcriptomic analysis between sympatric, ecologically divergent crater-lake Midas cichlids (Lake Xiloá Amphilophus amarillo and Amphilophus sagittae). We pyrosequenced (Roche 454) expressed sequence tag (EST) libraries and generated more than 178 000 000 ESTs and identified nine ESTs under positive selection between these sister species (Ka/Ks > 1). None of these ESTs were found to be under selection in African cichlids. Of 11 candidate genes for ecomorphological differentiation in African cichlids, none showed signs of selection between A. amarillo and A. sagittae. Although more population-level studies are now needed to thoroughly document patterns of divergence during speciation of cichlids, available information so far suggests that adaptive phenotypic diversification in Neotropical and African cichlids may be evolving through non-parallel genetic bases.

Social status predicts how sex steroid receptors regulate complex behavior across levels of biological organization

Social status strongly affects behavior and physiology, in part mediated by gonadal hormones, although how each sex steroid acts across levels of biological organization is not well understood. We examine the role of sex steroids in modulating social behavior in dominant (DOM) and subordinate (SUB) males of a highly social fish, Astatotilapia burtoni. We first used agonists and antagonists to each sex steroid receptor and found that androgens and progestins modulate courtship behavior only in DOM, whereas estrogens modulate aggressive behavior independent of social status. We then examined the hormonal and physiological responses to sex steroid receptor antagonist treatment and uncovered substantial changes in circulating steroid hormone levels and gonad size only in SUB, not in DOM. Consistent with status-based physiological sensitivities to drug manipulation, we found that neuropeptide and steroid receptor gene expression in the preoptic area was sensitive only in SUB. However, when we compared the transcriptomes of males that received either vehicle or an estrogen receptor antagonist, 8.25% of all genes examined changed expression in DOM in comparison with only 0.56% in SUB. Finally, we integrate behavior, physiology, and brain gene expression to infer functional modules that underlie steroid receptor regulation of behavior. Our work suggests that environmentally induced changes at one level of biological organization do not simply affect changes of similar magnitude at other levels, but that instead very few key pathways likely serve as conduits for executing plastic responses across multiple levels.

The function of anal fin egg-spots in the cichlid fish Astatotilapia burtoni

Color and pigmentation patterns of animals are often targets of sexual selection because of their role in communication. Although conspicuous male traits are typically implicated with intersexual selection, there are examples where sex-specific displays play a role in an intrasexual context, e.g. when they serve as signals for aggression level and/or status. Here, we focus on the function of a conspicuous male ornament in the most species-rich tribe of cichlid fishes, the haplochromines. A characteristic feature of these ca. 1500 species are so-called egg-spots in form of ovoid markings on the anal fins of males, which are made up of carotenoid based pigment cells. It has long been assumed that these yellow, orange or reddish egg-spots play an important role in the courtship and spawning behavior of these maternal mouth-brooding fishes by mimicking the eggs of a conspecific female. The exact function of egg-spots remains unknown, however, and there are several hypotheses about their mode of action. To uncover the function of this cichlid-specific male ornament, we used female mate choice experiments and a male aggression test in the haplochromine species Astatotilapia burtoni. We manipulated the number and arrangement of egg-spots on the anal fins of males, or removed them entirely, and tested (1) female preference with visual contact only using egg-traps, (2) female preference with free contact using paternity testing with microsatellites and (3) male aggression. We found that females did not prefer males with many egg-spots over males with fewer egg-spots and that females tended to prefer males without egg-spots over males with egg-spots. Importantly, males without egg-spots sired clutches with the same fertilization rate as males with egg-spots. In male aggression trials, however, males with fewer egg-spots received significantly more attacks, suggesting that egg-spots are an important signal in intrasexual communication.

Comparative gene expression profiles for highly similar aggressive phenotypes in male and female cichlid fishes (Julidochromis)

Julidochromis marlieri and Julidochromis transcriptus are two closely related Tanganyikan cichlids that have evolved different behavior and mating strategies since they diverged from their common ancestor. While J. transcriptus follows the ancestral pattern of male dominance, male-biased sexual size dimorphism and territoriality, the pattern is reversed in J. marlieri. In J. marlieri, females show all of these behavioral and morphological characteristics. This raises the question of whether female J. marlieri achieve the dominant phenotype by expressing the same genes as J. transcriptus males or whether novel brain gene expression patterns have evolved to produce a similar behavioral phenotype in the females of J. marlieri. This study used cDNA microarrays to investigate whether female J. marlieri and male J. transcriptus show conserved or divergent patterns of brain gene expression. Analysis of microarray data in both species showed certain gene expression patterns associated with sex role independent of gonadal sex and, to a lesser extent, gene expression patterns associated with sex independent of sex role. In general, these data suggest that while there has been substantial divergence in gene expression patterns between J. transcriptus and J. marlieri, we can detect a highly significant overlap for a core set of genes related to aggression in both species. These results suggest that the proximate mechanisms regulating aggressive behavior in J. transcriptus and J. marlieri may be shared.

Positive Darwinian selection drives the evolution of the morphology-related gene, EPCAM, in particularly species-rich lineages of African cichlid fishes

The study of genetic evolution within the context of adaptive radiations offers insights to genes and selection pressures that result in rapid morphological change. Cichlid fishes are very species-rich and variable in coloration, behavior, and morphology, and so provide a classical model system for studying the genetics of adaptive radiation. In this study, we researched the evolution of the epithelial cell adhesion molecule (EPCAM), a candidate gene for the adaptive evolution of morphology broadly, and skin development specifically, in fishes. We compared EPCAM gene sequences from a rapidly speciating African cichlid lineage (the haplochromines), a species-poor African lineage (Nile tilapia Oreochromis niloticus), and a very young adaptive radiation in the Neotropics (sympatric crater lake Midas cichlids, Amphilophus sp.). Our results, based on a hierarchy of evolutionary analyses of nucleotide substitution, demonstrate that there are different selection pressures on the EPCAM gene among the cichlid lineages. Several waves of positive natural selection were identified not only on the terminal branches, but also on ancestral branches. Interestingly, significant positive or directional selection was found in the haplochromine cichlids only but not the comparatively species-poor tilapia lineage. We hypothesize that the strong signal of selection in the ancestral African cichlid lineage coincided with the transition from riverine to lacustrine habitat. The two neotropical species for which we collected new sequence data were invariant in the EPCAM locus. Our results suggest that functional changes promoted by positive Darwinian selection are widespread in the EPCAM gene during African cichlid evolution.

Serotonergic modulation of startle-escape plasticity in an African cichlid fish: a single-cell molecular and physiological analysis of a vital neural circuit

Social life affects brain function at all levels, including gene expression, neurochemical balance, and neural circuits. We have previously shown that in the cichlid fish Astatotilapia burtoni brightly colored, socially dominant (DOM) males face a trade-off between reproductive opportunities and increased predation risk. Compared with camouflaged subordinate (SUB) males, DOMs exposed to a loud sound pip display higher startle responsiveness and increased excitability of the Mauthner cell (M-cell) circuit that governs this behavior. Using behavioral tests, intracellular recordings, and single-cell molecular analysis, we show here that serotonin (5-HT) modulates this socially regulated plasticity via the 5-HT receptor subtype 2 (5-HTR2). Specifically, SUBs display increased sensitivity to pharmacological manipulation of 5-HTR2 compared with DOMs in both startle-escape behavior and electrophysiological properties of the M-cell. Immunohistochemistry showed serotonergic varicosities around the M-cells, further suggesting that 5-HT impinges directly onto the startle-escape circuitry. To determine whether the effects of 5-HTR2 are pre- or postsynaptic, and whether other 5-HTR subtypes are involved, we harvested the mRNA from single M-cells via cytoplasmic aspiration and found that 5-HTR subtypes 5A and 6 are expressed in the M-cell. 5-HTR2, however, was absent, suggesting that it affects M-cell excitability through a presynaptic mechanism. These results are consistent with a role for 5-HT in modulating startle plasticity and increase our understanding of the neural and molecular basis of a trade-off between reproduction and predation.

Short read Illumina data for the de novo assembly of a non-model snail species transcriptome (Radix balthica, Basommatophora, Pulmonata), and a comparison of assembler performance

Background

Until recently, read lengths on the Solexa/Illumina system were too short to reliably assemble transcriptomes without a reference sequence, especially for non-model organisms. However, with read lengths up to 100 nucleotides available in the current version, an assembly without reference genome should be possible. For this study we created an EST data set for the common pond snail Radix balthica by Illumina sequencing of a normalized transcriptome. Performance of three different short read assemblers was compared with respect to: the number of contigs, their length, depth of coverage, their quality in various BLAST searches and the alignment to mitochondrial genes.

Results

A single sequencing run of a normalized RNA pool resulted in 16,923,850 paired end reads with median read length of 61 bases. The assemblies generated by VELVET, OASES, and SeqMan NGEN differed in the total number of contigs, contig length, the number and quality of gene hits obtained by BLAST searches against various databases, and contig performance in the mt genome comparison. While VELVET produced the highest overall number of contigs, a large fraction of these were of small size (< 200bp), and gave redundant hits in BLAST searches and the mt genome alignment. The best overall contig performance resulted from the NGEN assembly. It produced the second largest number of contigs, which on average were comparable to the OASES contigs but gave the highest number of gene hits in two out of four BLAST searches against different reference databases. A subsequent meta-assembly of the four contig sets resulted in larger contigs, less redundancy and a higher number of BLAST hits.

Conclusion

Our results document the first de novo transcriptome assembly of a non-model species using Illumina sequencing data. We show that de novo transcriptome assembly using this approach yields results useful for downstream applications, in particular if a meta-assembly of contig sets is used to increase contig quality. These results highlight the ongoing need for improvements in assembly methodology.

Comparative transcriptomics of Eastern African cichlid fishes shows signs of positive selection and a large contribution of untranslated regions to genetic diversity

The hundreds of endemic species of cichlid fishes in the East African Great Lakes Tanganyika, Malawi, and Victoria are a prime model system in evolutionary biology. With five genomes currently being sequenced, eastern African cichlids also represent a forthcoming genomic model for evolutionary studies of genotype-to-phenotype processes in adaptive radiations. Here we report the functional annotation and comparative analyses of transcriptome data sets for two eastern African cichlid species, Astatotilapia burtoni and Ophthalmotilapia ventralis, representatives of the modern haplochromines and ectodines, respectively. Nearly 647,000 expressed sequence tags were assembled in more than 46,000 contigs for each species using the 454 sequencing technology, largely expanding the current sequence data set publicly available for these cichlids. Total predicted coverage of their proteome diversity is approximately 50% for both species. Comparative qualitative and quantitative analyses show very similar transcriptome data for the two species in terms of both functional annotation and relative abundance of gene ontology terms expressed. Average genetic distance between species is 1.75% when all transcript types are considered including nonannotated sequences, 1.33% for annotated sequences only including untranslated regions, and decreases to nearly half, 0.95%, for coding sequences only, suggesting a large contribution of noncoding regions to their genetic diversity. Comparative analyses across the two species, tilapia and the outgroup medaka based on an overlapping data set of 1,216 genes (∼526 kb) demonstrate cichlid-specific signature of disruptive selection and provide a set of candidate genes that are putatively under positive selection. Overall, these data sets offer the genetic platform for future comparative analyses in light of the upcoming genomes for this taxonomic group.

The ontogenetic development of egg-spots in the haplochromine cichlid fish Astatotilapia burtoni

A series of Astatotilapia burtoni photographs were used to investigate the ontogenetic development of male egg-spots, a putative evolutionary key innovation of haplochromine cichlids. Four stages of egg-spot development were defined and all males had developed true egg-spots (stage 4) before reaching a standard length of 25 mm. Raising condition only slightly influenced the timing of the first appearance of true egg-spots.

Identification and characterization of gene expression involved in the coloration of cichlid fish using microarray and qRT-PCR approaches

It has been suggested that speciation on the basis of sexual selection is an important mechanism for the generation of new species for East African cichlids, where male body coloration is one of the major discriminatory factors used by females in mate choice. To gain insight into the molecular basis of cichlid coloration, we studied the Lake Malawi cichlid Pseudotropheus saulosi, comparing transcription in the bright blue skin of males to the yellow skin of females. Our cDNA microarray experiments identified 46 clones that exhibited expression differences between the two sexes, of which five were confirmed to be differentially expressed by relative quantitative real-time PCR (qRT-PCR). This gene list includes a representative from the endosomal-to-Golgi vesicle trafficking pathway, Coatomer protein complex, subunit zeta-1 (Copz-1), which is known to be a critical determinant of pigmentation in humans and zebrafish. With the support of microscopic images of the skin of these specimens, we interpret the transcriptional differences between the blue males and yellow females. Here, we provide insight into the putative functional diversification of genes involved in the coloration of cichlids and by extension, on the evolution of coloration in teleost fish.

Heritability and adaptive significance of the number of egg-dummies in the cichlid fish Astatotilapia burtoni

Cichlid fishes are a textbook example of rapid speciation and exuberant diversity--this applies especially to haplochromines, a lineage with approximately 1800 species. Haplochromine males uniquely possess oval, bright spots on their anal fin, called 'egg-spots' or 'egg-dummies'. These are presumed to be an evolutionary key innovation that contributed to the tribe's evolutionary success. Egg-spots have been proposed to mimic the ova of the mouthbrooding females of the corresponding species, contribute to fertilization success and even facilitate species recognition. Interestingly, egg-spot number varies extensively not only between species, but also within some populations. This high degree of intraspecific variation may appear to be counterintuitive since selection might be expected to act to stabilize traits that are correlated with fitness measures. We addressed this 'paradox' experimentally, and found that in the haplochromine cichlid Astatotilapia burtoni, the number of egg-spots was related to male age, body condition and dominance status. Intriguingly, the egg-spot number also had a high heritable component (narrow sense heritability of 0.5). These results suggest that the function of egg-spots might have less to do with fertilization success or species recognition, but rather relate to mate choice and/or male-male competition, helping to explain the high variability in this important trait.

Comparative physical maps derived from BAC end sequences of tilapia (Oreochromis niloticus)

Background

The Nile tilapia is the second most important fish in aquaculture. It is an excellent laboratory model, and is closely related to the African lake cichlids famous for their rapid rates of speciation. A suite of genomic resources has been developed for this species, including genetic maps and ESTs. Here we analyze BAC end-sequences to develop comparative physical maps, and estimate the number of genome rearrangements, between tilapia and other model fish species.

Results

We obtained sequence from one or both ends of 106,259 tilapia BACs. BLAST analysis against the genome assemblies of stickleback, medaka and pufferfish allowed identification of homologies for approximately 25,000 BACs for each species. We calculate that rearrangement breakpoints between tilapia and these species occur about every 3 Mb across the genome. Analysis of 35,000 clones previously assembled into contigs by restriction fingerprints allowed identification of longer-range syntenies.

Conclusions

Our data suggest that chromosomal evolution in recent teleosts is dominated by alternate loss of gene duplicates, and by intra-chromosomal rearrangements (~one per million years). These physical maps are a useful resource for comparative positional cloning of traits in cichlid fishes. The paired BAC end sequences from these clones will be an important resource for scaffolding forthcoming shotgun sequence assemblies of the tilapia genome.

Comparative genomics in teleost species: Knowledge transfer by linking the genomes of model and non-model fish species

Comparative genomics is a powerful tool to transfer knowledge coming from model fish species to non-model fish species of economic or/and evolutionary interest. Such transfer is of importance as functional studies either are difficult to perform with most non-model species. The first comparative map constructed using the human and the chimpanzee genome allowed the identification of putative orthologues. Although comparative mapping in teleosts is still in its infancy, five model teleost genomes from different orders have been fully sequenced to date and the sequencing of several commercially important species are also underway or near completion. The accessibility of these whole genome sequences and rapid developments in genomics of fish species are paving the way towards new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of non-model, but economically, physiologically or evolutionary important species is now feasible. Furthermore, comparison of low coverage gene maps of non-model fish species against fully sequenced fish species will enhance the efficiency of candidate gene identification projected for quantitative trait loci (QTL) scans for traits of special interest.

Rapid evolution and selection inferred from the transcriptomes of sympatric crater lake cichlid fishes

Crater lakes provide a natural laboratory to study speciation of cichlid fishes by ecological divergence. Up to now, there has been a dearth of transcriptomic and genomic information that would aid in understanding the molecular basis of the phenotypic differentiation between young species. We used next-generation sequencing (Roche 454 massively parallel pyrosequencing) to characterize the diversity of expressed sequence tags between ecologically divergent, endemic and sympatric species of cichlid fishes from crater lake Apoyo, Nicaragua: benthic Amphilophus astorquii and limnetic Amphilophus zaliosus. We obtained 24 174 A. astorquii and 21 382 A. zaliosus high-quality expressed sequence tag contigs, of which 13 106 pairs are orthologous between species. Based on the ratio of nonsynonymous to synonymous substitutions, we identified six sequences exhibiting signals of strong diversifying selection (K(a)/K(s) > 1). These included genes involved in biosynthesis, metabolic processes and development. This transcriptome sequence variation may be reflective of natural selection acting on the genomes of these young, sympatric sister species. Based on Ks ratios and p-distances between 3'-untranslated regions (UTRs) calibrated to previously published species divergence times, we estimated a neutral transcriptome-wide substitutional mutation rate of approximately 1.25 x 10(-6) per site per year. We conclude that next-generation sequencing technologies allow us to infer natural selection acting to diversify the genomes of young species, such as crater lake cichlids, with much greater scope than previously possible.

An EST resource for tilapia based on 17 normalized libraries and assembly of 116,899 sequence tags

BACKGROUND

Large collections of expressed sequence tags (ESTs) are a fundamental resource for analysis of gene expression and annotation of genome sequences. We generated 116,899 ESTs from 17 normalized and two non-normalized cDNA libraries representing 16 tissues from tilapia, a cichlid fish widely used in aquaculture and biological research.

RESULTS

The ESTs were assembled into 20,190 contigs and 36,028 singletons for a total of 56,218 unique sequences and a total assembled length of 35,168,415 bp. Over the whole project, a unique sequence was discovered for every 2.079 sequence reads. 17,722 (31.5%) of these unique sequences had significant BLAST hits (e-value < 10(-10)) to the UniProt database.

CONCLUSION

Normalization of the cDNA pools with double-stranded nuclease allowed us to efficiently sequence a large collection of ESTs. These sequences are an important resource for studies of gene expression, comparative mapping and annotation of the forthcoming tilapia genome sequence.

EST analysis and identification of gonad-related genes from the normalized cDNA library of large yellow croaker, Larimichthys crocea

On grounds of the especially limited numbers of identified gonad-specific or gonad-related genes of large yellow croaker Larimichthys crocea which may represent a major obstacle for the study of gonad development and sex differentiation, we initiated a sequencing program of Expressed Sequence Tags (ESTs) in large yellow croaker. In this study, we firstly constructed a normalized gonad cDNA library using the combination of SMART technique and DSN treatment. The titer of amplified cDNA library was 4.8x10(11) and the percentage of unique cDNA sequences of the library was 82.49%. 2916 unique cDNAs were clustered from the 3535 high quality ESTs. Among the 1785 ESTs which had significant homology with known genes in the NCBI database, about 64 significant gonad-related genes were found, accounting for 3.59% of the total unique cDNAs. Specifically, the testis-specific LRR gene and testis-specific chromodomain Y-like protein gene were identified from fish for the first time. Six gonad-related microsatellite-containing ESTs were identified from the 129 ESTs containing 149 microsatellites. Expression patterns of 10 of these gonad-related gene homologues in ovaries and testes were examined by qRT-PCR. The results will be powerful resources for our further investigation to establish the molecular mechanisms of gonad development and sex differentiation in large yellow croaker.

Ecological speciation in marine v. freshwater fishes

Absolute barriers to dispersal are not common in marine systems, and the prevalence of planktonic larvae in marine taxa provides potential for gene flow across large geographic distances. These observations raise the fundamental question in marine evolutionary biology as to whether geographic and oceanographic barriers alone can account for the high levels of species diversity observed in marine environments such as coral reefs, or whether marine speciation also operates in the presence of gene flow between diverging populations. In this respect, the ecological hypothesis of speciation, in which reproductive isolation results from divergent or disruptive natural selection, is of particular interest because it may operate in the presence of gene flow. Although important insights into the process of ecological speciation in aquatic environments have been provided by the study of freshwater fishes, comparatively little is known about the possibility of ecological speciation in marine teleosts. In this study, the evidence consistent with different aspects of the ecological hypothesis of speciation is evaluated in marine fishes. Molecular approaches have played a critical role in the development of speciation hypotheses in marine fishes, with a role of ecology suggested by the occurrence of sister clades separated by ecological factors, rapid cladogenesis or the persistence of genetically and ecologically differentiated species in the presence of gene flow. Yet, ecological speciation research in marine fishes is still largely at an exploratory stage. Cases where the major ingredients of ecological speciation, namely a source of natural divergent or disruptive selection, a mechanism of reproductive isolation and a link between the two have been explicitly documented are few. Even in these cases, specific predictions of the ecological hypothesis of speciation remain largely untested. Recent developments in the study of freshwater fishes illustrate the potential for molecular approaches to address specific questions related to the ecological hypothesis of speciation such as the nature of the genes underlying key ecological traits, the magnitude of their effect on phenotype and the mechanisms underlying their differential expression in different ecological contexts. The potential provided by molecular studies is fully realized when they are complemented with alternative (e.g. ecological, theoretical) approaches.

Genomic reaction norms: using integrative biology to understand molecular mechanisms of phenotypic plasticity

Phenotypic plasticity is the development of different phenotypes from a single genotype, depending on the environment. Such plasticity is a pervasive feature of life, is observed for various traits and is often argued to be the result of natural selection. A thorough study of phenotypic plasticity should thus include an ecological and an evolutionary perspective. Recent advances in large-scale gene expression technology make it possible to also study plasticity from a molecular perspective, and the addition of these data will help answer long-standing questions about this widespread phenomenon. In this review, we present examples of integrative studies that illustrate the molecular and cellular mechanisms underlying plastic traits, and show how new techniques will grow in importance in the study of these plastic molecular processes. These techniques include: (i) heterologous hybridization to DNA microarrays; (ii) next generation sequencing technologies applied to transcriptomics; (iii) techniques for studying the function of noncoding small RNAs; and (iv) proteomic tools. We also present recent studies on genetic model systems that uncover how environmental cues triggering different plastic responses are sensed and integrated by the organism. Finally, we describe recent work on changes in gene expression in response to an environmental cue that persist after the cue is removed. Such long-term responses are made possible by epigenetic molecular mechanisms, including DNA methylation. The results of these current studies help us outline future avenues for the study of plasticity.

Interspecific profiling of gene expression informed by comparative genomic hybridization: A review and a novel approach in African cichlid fishes

Modern genomic approaches have facilitated great progress in our understanding of the molecular and genetic underpinnings of ecological and evolutionary processes. Analysis of gene expression through heterologous hybridization in particular has enabled genome-scale studies in many ecologically and evolutionarily interesting species. However, these studies have been hampered by the difficulty of comparing-on a common array platform-gene-expression profiles across species due to sequence divergence altering the dynamics of hybridization. All too often, comparisons of expression profiles across species were limited to contrasting lists of gene or even of just functional categories. Here we review these issues and propose a novel solution. Exploiting the diverse cichlid lineages of East Africa as our model-system, we then present results from an experimental case study that compares the neural gene-expression profiles of males and females of two species that differ in mating system. Using a single microarray platform that contains genes from one species, Astatotilapia burtoni, we conducted a total of 16 direct comparisons for neural gene-expression level between individual males and females from a pair of sister species, the polygynous Enantiopus melanogenys and the monogamous Xenotilapia flavipinnis. Next, we conducted a meta-analysis with previously published data from two different intra-specific expression studies to determine whether sex-specific neural gene expression is more closely associated with behavioral phenotype than it is with gonadal sex. Our results indicate that the gene expression profiles are species-specific to a large extent, as relatively few genes show conserved expression patterns associated with either sex. Finally, we describe how competitive genomic DNA hybridizations between the two focal species allow us to assess the degree to which divergence of sequences biases the results. We propose a masking technique that correlates interspecific expression ratios obtained with cDNA with hybridization ratios obtained with genomic DNA for the same set of species and determines threshold sequence divergence to reduce false positives. Our approach should be applicable to a wide range of interesting questions related to the evolution and ecology of gene expression.

Genomic analysis of cichlid fish 'natural mutants'

In the lakes of East Africa, cichlid fishes have formed adaptive radiations that are each composed of hundreds of endemic, morphologically stunningly diverse, but genetically extremely similar species. In the past 20 years, it became clear that their extreme phenotypic diversity arose within very short time spans, and that phenotypically radically different species are exceptionally similar genetically; hence, they could be considered to be 'natural mutants'. Many species can be hybridized and, therefore, provide a unique opportunity to study the genetic underpinnings of phenotypic diversification. Comparative large-scale genomic analyses are beginning to unravel the patterns and processes that led to the formation of the cichlid species flocks. Cichlids are an emerging evolutionary genomic model system for fundamental questions on the origin of phenotypic diversity.

Fish and chips: functional genomics of social plasticity in an African cichlid fish

Behavior and physiology are regulated by both environment and social context. A central goal in the study of the social control of behavior is to determine the underlying physiological, cellular and molecular mechanisms in the brain. The African cichlid fish Astatotilapia burtoni has long been used as a model system to study how social interactions regulate neural and behavioral plasticity. In this species, males are either socially dominant and reproductively active or subordinate and reproductively suppressed. This phenotypic difference is reversible. Using an integrative approach that combines quantitative behavioral measurements, functional genomics and bioinformatic analyses, we examine neural gene expression in dominant and subordinate males as well as in brooding females. We confirm the role of numerous candidate genes that are part of neuroendocrine pathways and show that specific co-regulated gene sets (modules), as well as specific functional gene ontology categories, are significantly associated with either dominance or reproductive state. Finally, even though the dominant and subordinate phenotypes are robustly defined, we find a surprisingly high degree of individual variation in the transcript levels of the very genes that are differentially regulated between these phenotypes. The results of the present study demonstrate the molecular complexity in the brain underlying social behavior, identify novel targets for future studies, validate many candidate genes and exploit individual variation in order to gain biological insights.

A microsatellite-based genetic linkage map of the cichlid fish, Astatotilapia burtoni (Teleostei): a comparison of genomic architectures among rapidly speciating cichlids

Cichlid fishes compose an astonishingly large number of species and formed species flocks in record-breaking time. To facilitate efficient genome scans and comparisons of cichlid genomes, we constructed a medium-density genetic linkage map of microsatellite markers of Astatotilapia burtoni. The mapping cross was derived from two inbred laboratory lines to obtain F(2) progeny by intercrossing. The map revealed 25 linkage groups spanning 1249.3 cM of the genome (size approximately 950 Mb) with an average marker spacing of 6.12 cM. The seven Hox clusters, ParaHox C1, and two paralogs of Pdgfrbeta were mapped to different linkage groups, thus supporting the hypothesis of a teleost-specific genome duplication. The A. burtoni linkage map was compared to the other two available maps for cichlids using shared markers that showed conservation and synteny among East African cichlid genomes. Interesting candidate genes for cichlid speciation were mapped using SNP markers.

Adaptive sequence evolution in a color gene involved in the formation of the characteristic egg-dummies of male haplochromine cichlid fishes

Background

The exceptionally diverse species flocks of cichlid fishes in East Africa are prime examples of parallel adaptive radiations. About 80% of East Africa's more than 1 800 endemic cichlid species, and all species of the flocks of Lakes Victoria and Malawi, belong to a particularly rapidly evolving lineage, the haplochromines. One characteristic feature of the haplochromines is their possession of egg-dummies on the males' anal fins. These egg-spots mimic real eggs and play an important role in the mating system of these maternal mouthbrooding fish.

Results

Here, we show that the egg-spots of haplochromines are made up of yellow pigment cells, xanthophores, and that a gene coding for a type III receptor tyrosine kinase, colony-stimulating factor 1 receptor a (csf1ra), is expressed in egg-spot tissue. Molecular evolutionary analyses reveal that the extracellular ligand-binding and receptor-interacting domain of csf1ra underwent adaptive sequence evolution in the ancestral lineage of the haplochromines, coinciding with the emergence of egg-dummies. We also find that csf1ra is expressed in the egg-dummies of a distantly related cichlid species, the ectodine cichlid Ophthalmotilapia ventralis, in which markings with similar functions evolved on the pelvic fin in convergence to those of the haplochromines.

Conclusion

We conclude that modifications of existing signal transduction mechanisms might have evolved in the haplochromine lineage in association with the origination of anal fin egg-dummies. That positive selection has acted during the evolution of a color gene that seems to be involved in the morphogenesis of a sexually selected trait, the egg-dummies, highlights the importance of further investigations of the comparative genomic basis of the phenotypic diversification of cichlid fishes.