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.     

Functional diversification of kir7.1 in cichlids accelerated by gene duplication

Mutation in the inward rectifier potassium channel gene, kir7.1, was previously identified as being responsible for the broader stripe zebrafish skin pattern mutant, jaguar/obelix. An amino acid substitution in this channel causes a broader stripe pattern than that of wild type zebrafish. In this study we analyzed cichlid homologs of the zebrafish kir7.1 gene. We identified two kinds of homologous genes in cichlids and named them cikir7.1 and cikir7.2. Southern hybridization using cichlid genome revealed that cichlids from the African Great Lakes, South America and Madagascar have two copies of the gene. Cichlids from Sri Lanka, however, showed only one band in this experiment. Database analysis revealed that only one copy of the kir7.1 gene exists in the genomes of the teleosts zebrafish, tetraodon, takifugu, medaka and stickleback. The deduced amino acid sequence of cikir7.1 is highly conserved among African cichlids, whereas that of cikir7.2 has several amino acid substitutions even in conserved transmembrane domains. Gene expression analysis revealed that cikir7.1 is expressed specifically in brain and eye, and cikir7.2 in testis and ovary; zebrafish kir7.1, however, is expressed in brain, eye, skin, caudal fin, testis and ovary. These results suggest that gene duplication of the cichlid kir7.1 occurred in a common ancestor of the family Cichlidae, that the function of parental kir7.1 was then divided into two genes, cikir7.1 and cikir7.2, and that the evolutionary rate of cikir7.2 might have been accelerated, thereby effecting functional diversification in the cichlid lineage. Thus, the evolution of kir7.1 genes in cichlids provides a typical example of gene duplication--one gene is conserved while the other becomes specialized for a novel function.     

Molecular phylogeny and speciation patterns in host-specific monogeneans (Cichlidogyrus,Dactylogyridae) parasitizing cichlid fishes (Cichliformes,Cichlidae) in Lake Tanganyika

Cichlidogyrus (including Scutogyrus) is the most speciose dactylogyridean monogenean genus known from African and Levantine cichlid fishes (Cichlidae). While its taxonomy is well established, little is known about the phylogenetic relationships and evolutionary history of this ectoparasite, especially from hosts belonging to one of the most impressive vertebrate radiations, the cichlid fishes from the East African Great Lakes and surrounding hydrological systems. Phylogenetic inference based on DNA sequences of the nuclear 18S, internal transcribed spacer 1 and 28S rDNA genes revealed that Cichlidogyrus parasitizing mainly West African cichlid tribes is paraphyletic with respect to species parasitizing hosts belonging to the East African cichlid radiation, which constitute a well-supported monophylum. Members of Cichlidogyrus from tylochromine and oreochromine hosts that colonised Lake Tanganyika only recently, cluster with their non-Lake Tanganyika relatives, indicating that they colonised Lake Tanganyika with their current host species, and did not jump over from any of the many cichlid species already present in the lake. The diversification of Cichlidogyrus in Lake Tanganyika seems to be driven by failure to diverge in old lineages of cichlids, cospeciation in more recently evolved ones, and host switching followed by parasite duplication at the level of the various host tribes. Evaluation of host specificity and structural evolution of haptoral and reproductive organs in Lake Tanganyika Cichlidogyrus revealed that strict specialist species together with larval hook size represent the ancestral state of haptor configuration, suggesting that members of Cichlidogyrus in this system evolved from a very simple form to a more complex one similarly to their West African congeners. Generalist species among Cichlidogyrus with a sclerotized vagina parasitizing ancient Lake Tanganyika lineages seem to have developed a different hook configuration, most probably to ensure successful colonisation of new, phylogenetically unrelated hosts.     

Cone opsin genes of african cichlid fishes: tuning spectral sensitivity by differential gene expression

Spectral tuning of visual pigments is typically accomplished through changes in opsin amino acid sequence. Within a given opsin class, changes at a few key sites control wavelength specificity. To investigate known differences in the visual pigment spectral sensitivity of the Lake Malawi cichlids, Metriaclima zebra (368, 488, and 533 nm) and Dimidiochromis compressiceps (447, 536, and 569 nm), we sequenced cone opsin genes from these species as well as Labeotropheus fuelleborni and Oreochromis niloticus. These cichlids have five distinct classes of cone opsin genes, including two unique SWS-2 genes. Comparisons of the inferred amino acid sequences from the five cone opsin genes of M. zebra, D. compressiceps, and L. fuelleborni show the sequences to be nearly identical. Therefore, evolution of key opsin sites cannot explain the differences in visual pigment sensitivities. Real-time PCR demonstrates that different cichlid species express different subsets of the available cone opsin genes. Metriaclima zebra and L. fuelleborni express a complement of genes which give them UV-shifted visual pigments, while D. compressiceps expresses a different set to produce a red-shifted visual system. Thus, variations in cichlid spectral sensitivity have arisen through evolution of gene regulation, rather than through changes in opsin amino acid sequence.     

Multilocus phylogeny of cichlid fishes (Pisces: Perciformes): evolutionary comparison of microsatellite and single-copy nuclear loci

Among vertebrates, cichlid fishes are the paradigmatic example of adaptive radiation and ecological specialization. In turn, molecular genetic studies of cichlids have focused primarily on more recently diverged groups. Here, we present an evolutionary hypothesis of the major lineages of cichlid fishes based on DNA sequence data from two nuclear loci. One marker, Tmo-4C4, is a single-copy locus containing a region of amino acid similarity to the muscle protein TITIN. Flanking sequence from a second, microsatellite, locus Tmo-M27, shows similarity to mammalian RAS guanine nucleotide-releasing factor. We compare and combine data from these loci to evaluate phylogenetic performance. In separate and combined analyses, the sequence data support and clarify previous morphological hypotheses of cichlid major-group relationships. Indian and Malagasy cichlids form a basal, paraphyletic group. Neotropical cichlids are the sister clade to an African assemblage composed of the paraphyletic west and Pan-African lineages and a group of east African rift lake taxa. We use a consensus phylogeny of the Cichlidae to trace evolutionary changes in the microsatellite repeat motif at Tmo-M27. Analysis reveals that the repeat region was nearly lost in the ancestor to cichlids and then amplified extensively in African taxa. Results demonstrate that the two new DNA markers could be widely applied in perciform systematics. Furthermore, the comparative approach can unveil mutational dynamics of simple-sequence repeat loci over long periods of fish evolution. Simple-sequence repeat regions are increasingly being found in introns of important regulatory genes. We address issues involving their function and suggest caution in making assumptions of strict neutrality.     

Complete Mitochondrial DNA Sequences of the Threadfin Cichlid (Petrochromis trewavasae) and the Blunthead Cichlid (Tropheus moorii) and Patterns of Mitochondrial Genome Evolution in Cichlid Fishes

The cichlid fishes of the East African Great Lakes represent a model especially suited to study adaptive radiation and speciation. With several African cichlid genome projects being in progress, a promising set of closely related genomes is emerging, which is expected to serve as a valuable data base to solve questions on genotype-phenotype relations. The mitochondrial (mt) genomes presented here are the first results of the assembly and annotation process for two closely related but eco-morphologically highly distinct Lake Tanganyika cichlids, Petrochromis trewavasae and Tropheus moorii. The genomic sequences comprise 16,588 bp (P. trewavasae) and 16,590 bp (T. moorii), and exhibit the typical mitochondrial structure, with 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a non-coding control region. Analyses confirmed that the two species are very closely related with an overall sequence similarity of 96%. We analyzed the newly generated sequences in the phylogenetic context of 21 published labroid fish mitochondrial genomes. Consistent with other vertebrates, the D-loop region was found to evolve faster than protein-coding genes, which in turn are followed by the rRNAs; the tRNAs vary greatly in the rate of sequence evolution, but on average evolve the slowest. Within the group of coding genes, ND6 evolves most rapidly. Codon usage is similar among examined cichlid tribes and labroid families; although a slight shift in usage patterns down the gene tree could be observed. Despite having a clearly different nucleotide composition, ND6 showed a similar codon usage. C-terminal ends of Cox1 exhibit variations, where the varying number of amino acids is related to the structure of the obtained phylogenetic tree. This variation may be of functional relevance for Cox1 synthesis.     

Age and spread of the haplochromine cichlid fishes in Africa

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

Molecular and fossil evidence place the origin of cichlid fishes long after Gondwanan rifting

Cichlid fishes are a key model system in the study of adaptive radiation, speciation and evolutionary developmental biology. More than 1600 cichlid species inhabit freshwater and marginal marine environments across several southern landmasses. This distributional pattern, combined with parallels between cichlid phylogeny and sequences of Mesozoic continental rifting, has led to the widely accepted hypothesis that cichlids are an ancient group whose major biogeographic patterns arose from Gondwanan vicariance. Although the Early Cretaceous (ca 135 Ma) divergence of living cichlids demanded by the vicariance model now represents a key calibration for teleost molecular clocks, this putative split pre-dates the oldest cichlid fossils by nearly 90 Myr. Here, we provide independent palaeontological and relaxed-molecular-clock estimates for the time of cichlid origin that collectively reject the antiquity of the group required by the Gondwanan vicariance scenario. The distribution of cichlid fossil horizons, the age of stratigraphically consistent outgroup lineages to cichlids and relaxed-clock analysis of a DNA sequence dataset consisting of 10 nuclear genes all deliver overlapping estimates for crown cichlid origin centred on the Palaeocene (ca 65–57 Ma), substantially post-dating the tectonic fragmentation of Gondwana. Our results provide a revised macroevolutionary time scale for cichlids, imply a role for dispersal in generating the observed geographical distribution of this important model clade and add to a growing debate that questions the dominance of the vicariance paradigm of historical biogeography.     

Temporal diversification of Mesoamerican cichlid fishes across a major biogeographic boundary

The Mexican Neovolcanic Plateau sharply divides the vertebrate fauna of Mesoamerica where the climate of both the neotropics and temperate North America gradually blend. Only a few vertebrate groups such as the Heroine cichlids, distributed from South America to the Rio Grande in North America, are found both north and south of the Neovolcanic Plateau. To better understand the geography and temporal diversification of cichlids at this geologic boundary, we used mitochondrial DNA sequences of the cytochrome b (cyt b) gene to reconstruct the relationships of 52 of the approximately 80 species of Heroine cichlids in Mesoamerica. Our analysis suggests several cichlids in South America should be considered as part of the Mesoamerican Heroine clade because they and the cichlids north of the Isthmus of Panama are clearly supported as monophyletic with respect to all other Neotropical cichlids. We also recovered a group containing species in Paratheraps+Paraneetroplus+Vieja as the sister clade to Herichthys. Herichthys is the only cichlid clade north of the Mexican Plateau and it is monophyletic. Non-parametric rate smoothing of cichlid cyt b sequence resulted in an estimated divergence time of approximately 6 million years for Herichthys. This temporal diversification is concordant with divergence times estimated for anurans in the genus Bufo, a group that exhibits a similar geographic distribution. Our results indicate the 5-million-year-old extension of the Mexican Neovolcanic Plateau to the Gulf Coast of Mexico has strongly influenced the current transition between the vertebrate faunas of the Neotropics and Nearctic.     

Molecular Characterization of Two Endothelin Pathways in East African Cichlid Fishes

The adaptive radiations of cichlid fishes in East Africa have been associated with the acquisition of evolutionary novelties as well as the ecological opportunities existing in the East African Great lakes. Two remarkable evolutionary innovations are the pharyngeal jaw apparatus, found in all cichlid species, and the anal fin egg-spots of mouthbrooding cichlids. Based on their conserved functions during the development of both the jaw apparatus and pigmentation, the endothelin ligands and receptors form a putative link between these naturally and sexually selected traits. Here we study the evolutionary history of four members of two endothelin pathways (Edn1/EdnrAa and Edn3b/EdnrB1a) to elucidate their possible roles during the evolution and development of key innovations in East African cichlids species. The analyses performed on partial sequences (ca. 6,000 bp per taxon) show that all four endothelin family members evolved under purifying selection, although both ligands are characterized by an accelerated rate of protein evolution in comparison to the receptors. In accordance with earlier findings, we show that the mature protein sequence of Edn1 and Edn3 are highly conserved, also in cichlids, whereas the preproendothelin parts are variable indicating relaxed selective constraints. In the receptors, nonsynonymous substitutions were mainly found in the ligand-binding domains suggesting functional divergence. Gene expression assays with Real-Time PCR indeed reveal that the two studied endothelin pathways are expressed in the cichlid pharyngeal jaw and in the haplochromine egg-spot (among other pigment-cell containing tissues), suggesting their involvement during morphogenesis of naturally and sexually selected traits in cichlids.     

Phylogenetic utility of the tyrosine kinase gene X-src for assessing relationships among representative cichlid fishes

The nuclear gene X-src is a member of the tyrosine-kinase class of proto-oncogenes whose normal product is localized within the cytoplasm of the cell. The X-src gene has been used in only a few phylogenetic studies, each focusing on systematics of killifishes (Cyprinodontiformes). The present study is an attempt to examine the phylogenetic utility of X-src for uncovering relationships of representative cichlid fishes, especially the cichlids of Middle America. The family Cichlidae is a species-rich group of tropical freshwater fish made up of more than 1000 species which show a Gondwanan pattern of distribution. Cichlid fish have been the focus of numerous studies ranging from behavioral to biogeographical to systematic in nature. Particular emphasis has focused on the cichlids of the African Great Lakes and the explosive adaptive radiation of this group. However, Neotropical cichlids have received considerably less attention than their African counterparts. Our findings regarding the utility of X-src concur with those of previous phylogenetic analyses showing the exons of X-src to be highly conserved and useful mostly for revealing deep relationships among taxa. Like previous X-src studies, we also found the intron sequences of the gene to be variable in length and difficult to align across distantly related taxa but they provided useful information for resolving relationships among more closely related taxa. The X-src phylogeny supports the monophyly of Neotropical cichlids and cichlasomines groups A (=heroines) + B (=cichlasomines). A highly resolved tree is obtained within the heroines but little support is evident for most nodes based on the low number of unambiguous substitutions. The X-src gene is likely to be quite useful for resolving deep phylogenetic relationships such as those among major groups of actinopterygian fishes.     

Phylogenetic analysis of cichlid fishes using nuclear DNA markers

The recent explosive adaptive radiation of cichlids in the great lakes of Africa has attracted the attention of both morphologists and molecular biologists. To decipher the phylogenetic relationships among the various taxa within the family Cichlidae is a prerequisite for answering some fundamental questions about the nature of the speciation process. In the present study, we used the random amplification of polymorphic DNA (RAPD) technique to obtain sequence differences between selected cichlid species. We then designed specific primers based on these sequences and used them to amplify template DNA from a large number of species by the polymerase chain reaction (PCR). We sequenced the amplified products and searched the sequences for indels and shared substitutions. We identified a number of such characters at three loci-- DXTU1, DXTU2, and DXTU3--and used them for phylogenetic and cladistic analysis of the relationships among the various cichlid groups. Our studies assign an outgroup position to Neotropical cichlids in relation to African cichlids, provide evidence for a sister-group relationship of tilapiines to the haplochromines, group Cyphotilapia frontosa with the lamprologines of Lake Tanganyika, place Astatoreochromis alluaudi to an outgroup position with respect to other haplochromines of Lakes Victoria and Malawi, and provide additional support for the monophyly of the remaining Lake Victoria haplochromines and the Lake Malawi haplochromines. The described approach holds great promise for further resolution of cichlid phylogeny.      

Cichlid genomics and phenotypic diversity in a comparative context

Cichlid fishes provide an excellent natural system for integrating studies of genomics and adaptive radiation. Cichlids are unique in comprising a substantial fraction of all vertebrate species, possessing unique jaw structures, displaying an exceptional range of breeding systems, and exhibiting rampant phenotypic convergence. The rate of divergence in cichlid jaws, teeth, color patterns, visual systems, reproductive biology, and mating behaviors is unparalleled among vertebrates. I discuss ways rapid divergence in cichlids and other adaptive radiations make understanding the genomic basis of adaptive divergence more tractable. Then, I briefly overview some major findings and insights into vertebrate adaptation that have been gained through cichlid genetic studies. Finally, I discuss the extensive evolutionary replication provided by cichlid adaptive radiations and their potential for studies of genotype-to-phenotype mapping.     

Mix and match color vision: tuning spectral sensitivity by differential opsin gene expression in Lake Malawi cichlids

Cichlid fish of the East African Rift Lakes are renowned for their diversity and offer a unique opportunity to study adaptive changes in the visual system in rapidly evolving species flocks. Since color plays a significant role in mate choice, differences in visual sensitivities could greatly influence and even drive speciation of cichlids. Lake Malawi cichlids inhabiting rock and sand habitats have significantly different cone spectral sensitivities. By combining microspectrophotometry (MSP) of isolated cones, sequencing of opsin genes, and spectral analysis of recombinant pigments, we have established the cone complements of four species of Malawi cichlids. MSP demonstrated that each of these species predominately expresses three cone pigments, although these differ between species to give three spectrally different cone complements. In addition, rare populations of spectrally distinct cones were found. In total, seven spectral classes were identified. This was confirmed by opsin gene sequencing, expression, and in vitro reconstitution. The genes represent the four major classes of cone opsin genes that diverged early in vertebrate evolution. All four species possess a long-wave-sensitive (LWS), three spectrally distinct green-sensitive (RH2), a blue-sensitive (SWS2A), a violet-sensitive (SWS2B), and an ultraviolet-sensitive (SWS1) opsin. However, African cichlids determine their spectral sensitivity by differential expression of primarily only three of the seven available cone opsin genes. Phylogenetic analysis suggests that all percomorph fish have similar potential.     

Genomics of speciation and introgression in Princess cichlid fishes from Lake Tanganyika

How variation in the genome translates into biological diversity and new species originate has endured as the mystery of mysteries in evolutionary biology. African cichlid fishes are prime model systems to address speciation‐related questions for their remarkable taxonomic and phenotypic diversity, and the possible role of gene flow in this process. Here, we capitalize on genome sequencing and phylogenomic analyses to address the relative impacts of incomplete lineage sorting, introgression and hybrid speciation in the Neolamprologus savoryi‐complex (the ‘Princess cichlids’) from Lake Tanganyika. We present a time‐calibrated species tree based on whole‐genome sequences and provide strong evidence for incomplete lineage sorting in the early phases of diversification and multiple introgression events affecting different stages. Importantly, we find that the Neolamprologus chromosomes show centre‐to‐periphery biases in nucleotide diversity, sequence divergence, GC content, incomplete lineage sorting and rates of introgression, which are likely modulated by recombination density and linked selection. The detection of heterogeneous genomic landscapes has strong implications on the genomic mechanisms involved in speciation. Collinear chromosomal regions can be protected from gene flow and harbour incompatibility genes if they reside in lowly recombining regions, and coupling can evolve between nonphysically linked genomic regions (chromosome centres in particular). Simultaneously, higher recombination towards chromosome peripheries makes these more dynamic, evolvable regions where adaptation polymorphisms have a fertile ground. Hence, differences in genome architecture could explain the levels of taxonomic and phenotypic diversity seen in taxa with collinear genomes and might have contributed to the spectacular cichlid diversity observed today.     

Using SINEs to probe ancient explosive speciation: "hidden" radiation of African cichlids?

Cichlid fishes of the east African Great Lakes represent a paradigm of adaptive radiation. We conducted a phylogenetic analysis of cichlids including pan-African and west African species by using insertion patterns of short interspersed elements (SINEs) at orthologous loci. The monophyly of the east African cichlids was consistently supported by seven independent insertions of SINE sequences that are uniquely shared by these species. In addition, data from four other loci indicated that the genera Tilapia (pan-African) and Steatocranus (west African) are the closest relatives to east African cichlids. However, relationships among Tilapia, Steatocranus, and the east African clade were ambiguous because of incongruencies among topologies suggested by insertion patterns of SINEs at six other loci. One plausible explanation for this phenomenon is incomplete lineage sorting of alleles containing or missing a SINE insertion at these loci during ancestral speciation. Such incomplete sorting may have taken place earlier than 14 MYA, followed by random and stochastic fixation of the alleles in subsequent lineages. These observations prompted us to consider the possibility that cichlid speciation occurred at an accelerated rate during this period when the African Great Lakes did not exist. The SINE method could be useful for detecting ancient exclusive speciation events that tend to remain hidden during conventional sequence analyses because of accumulated point mutations.