Origins of shared genetic variation in African cichlids

The genomic substrate for adaptive radiation in African cichlid fish

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Genetic mapping of horizontal stripes in Lake Victoria cichlid fishes: benefits and pitfalls of using RAD markers for dense linkage mapping

The genetic dissection of naturally occurring phenotypes sheds light on many fundamental and longstanding questions in speciation and adaptation and is a central research topic in evolutionary biology. Until recently, forward-genetic approaches were virtually impossible to apply to nonmodel organisms, but the development of next-generation sequencing techniques eases this difficulty. Here, we use the ddRAD-seq method to map a colour trait with a known adaptive function in cichlid fishes, well-known textbook examples for rapid rates of speciation and astonishing phenotypic diversification. A suite of phenotypic key innovations is related to speciation and adaptation in cichlids, among which body coloration features prominently. The focal trait of this study, horizontal stripes, evolved in parallel in several cichlid radiations and is associated with piscivorous foraging behaviour. We conducted interspecific crosses between Haplochromis sauvagei and H. nyererei and constructed a linkage map with 867 SNP markers distributed on 22 linkage groups and total size of 1130.63 cM. Lateral stripes are inherited as a Mendelian trait and map to a single genomic interval that harbours a paralog of a gene with known function in stripe patterning. Dorsolateral and mid-lateral stripes were always coinherited and are thus under the same genetic control. Additionally, we directly quantify the genotyping error rates in RAD markers and offer guidelines for identifying and dealing with errors. Uncritical marker selection was found to severely impact linkage map construction. Fortunately, by applying appropriate quality control steps, a genotyping accuracy of >99.9% can be reached, thus allowing for efficient linkage mapping of evolutionarily relevant traits.

Phylogeography of Heliconius cydno and its closest relatives: disentangling their origin and diversification

The origins of the extraordinary diversity within the Neotropics have long fascinated biologists and naturalists. Yet, the underlying factors that have given rise to this diversity remain controversial. To test the relative importance of Quaternary climatic change and Neogene tectonic and paleogeographic reorganizations in the generation of biodiversity, we examine intraspecific variation across the Heliconius cydno radiation and compare this variation to that within the closely related Heliconius melpomene and Heliconius timareta radiations. Our data, which consist of both mtDNA and genome-scan data from nearly 2250 amplified fragment length polymorphism (AFLP) loci, reveal a complex history of differentiation and admixture at different geographic scales. Both mtDNA and AFLP phylogenies suggest that H. timareta and H. cydno are probably geographic extremes of the same radiation that probably diverged from H. melpomene prior to the Pliocene-Pleistocene boundary, consistent with hypotheses of diversification that rely on geological events in the Pliocene. The mtDNA suggests that this radiation originated in Central America or the northwestern region of South America, with a subsequent colonization of the eastern and western slopes of the Andes. Our genome-scan data indicate significant admixture among sympatric H. cydno/H. timareta and H. melpomene populations across the extensive geographic ranges of the two radiations. Within H. cydno, both mtDNA and AFLP data indicate significant population structure at local scales, with strong genetic differences even among adjacent H. cydno colour pattern races. These genetic patterns highlight the importance of past geoclimatic events, intraspecific gene flow, and local population differentiation in the origin and establishment of new adaptive forms.

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.

Interspecific variation in Rx1 expression controls opsin expression and causes visual system diversity in African cichlid fishes

The mechanisms underlying natural phenotypic diversity are key to understanding evolution and speciation. Cichlid fishes are among the most speciose vertebrates and an ideal model for identifying genes controlling species differences. Cichlids have diverse visual sensitivities that result from species expressing subsets of seven cichlid cone opsin genes. We previously identified a quantitative trait locus (QTL) that tunes visual sensitivity by varying SWS2A (short wavelength sensitive 2A) opsin expression in a genetic cross between two Lake Malawi cichlid species. Here, we identify Rx1 (retinal and anterior neural fold homeobox) as the causative gene for the QTL using fine mapping and RNAseq in retinal transcriptomes. Rx1 is differentially expressed between the parental species and correlated with SWS2A expression in the F2 progeny. Expression of Rx1 and SWS2A is also correlated in a panel of 16 Lake Malawi cichlid species. Association mapping in this panel identified a 413-bp deletion located 2.5-kb upstream of the Rx1 translation start site that is correlated with decreased Rx1 expression. This deletion explains 62% of the variance in SWS2A expression across 53 cichlid species in 29 genera. The deletion occurs in both the sand and rock-dwelling cichlid clades, suggesting that it is an ancestral polymorphism. Our finding supports the hypothesis that mixing and matching of ancestral polymorphisms can explain the diversity of present day cichlid phenotypes.

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

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

The complex hybrid origins of the root knot nematodes revealed through comparative genomics

Root knot nematodes (RKN) can infect most of the world's agricultural crop species and are among the most important of all plant pathogens. As yet however we have little understanding of their origins or the genomic basis of their extreme polyphagy. The most damaging pathogens reproduce by obligatory mitotic parthenogenesis and it has been suggested that these species originated from interspecific hybridizations between unknown parental taxa. We have sequenced the genome of the diploid meiotic parthenogen Meloidogyne floridensis, and use a comparative genomic approach to test the hypothesis that this species was involved in the hybrid origin of the tropical mitotic parthenogen Meloidogyne incognita. Phylogenomic analysis of gene families from M. floridensis, M. incognita and an outgroup species Meloidogyne hapla was carried out to trace the evolutionary history of these species' genomes, and we demonstrate that M. floridensis was one of the parental species in the hybrid origins of M. incognita. Analysis of the M. floridensis genome itself revealed many gene loci present in divergent copies, as they are in M. incognita, indicating that it too had a hybrid origin. The triploid M. incognita is shown to be a complex double-hybrid between M. floridensis and a third, unidentified, parent. The agriculturally important RKN have very complex origins involving the mixing of several parental genomes by hybridization and their extreme polyphagy and success in agricultural environments may be related to this hybridization, producing transgressive variation on which natural selection can act. It is now clear that studying RKN variation via individual marker loci may fail due to the species' convoluted origins, and multi-species population genomics is essential to understand the hybrid diversity and adaptive variation of this important species complex. This comparative genomic analysis provides a compelling example of the importance and complexity of hybridization in generating animal species diversity more generally.

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.

Identification of cichlid fishes from Lake Malawi using computer vision

BACKGROUND

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

METHODOLOGY/PRINCIPAL FINDING

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

CONCLUSIONS

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

Speciation in fishes

The field of speciation has seen much renewed interest in the past few years, with theoretical and empirical advances that have moved it from a descriptive field to a predictive and testable one. The goal of this review is to provide a general background on research on speciation as it pertains to fishes. Three major components to the question are first discussed: the spatial, ecological and sexual factors that influence speciation mechanisms. We then move to the latest developments in the field of speciation genomics. Affordable and rapidly available, massively parallel sequencing data allow speciation studies to converge into a single comprehensive line of investigation, where the focus has shifted to the search for speciation genes and genomic islands of speciation. We argue that fish present a very diverse array of scenarios, making them an ideal model to study speciation processes.

Temporal patterns of diversification across global cichlid biodiversity (Acanthomorpha: Cichlidae)

The contrasting distribution of species diversity across the major lineages of cichlids makes them an ideal group for investigating macroevolutionary processes. In this study, we investigate whether different rates of diversification may explain the disparity in species richness across cichlid lineages globally. We present the most taxonomically robust time-calibrated hypothesis of cichlid evolutionary relationships to date. We then utilize this temporal framework to investigate whether both species-rich and depauperate lineages are associated with rapid shifts in diversification rates and if exceptional species richness can be explained by clade age alone. A single significant rapid rate shift increase is detected within the evolutionary history of the African subfamily Pseudocrenilabrinae, which includes the haplochromins of the East African Great Lakes. Several lineages from the subfamilies Pseudocrenilabrinae (Australotilapiini, Oreochromini) and Cichlinae (Heroini) exhibit exceptional species richness given their clade age, a net rate of diversification, and relative rates of extinction, indicating that clade age alone is not a sufficient explanation for their increased diversity. Our results indicate that the Neotropical Cichlinae includes lineages that have not experienced a significant rapid burst in diversification when compared to certain African lineages (rift lake). Neotropical cichlids have remained comparatively understudied with regard to macroevolutionary patterns relative to African lineages, and our results indicate that of Neotropical lineages, the tribe Heroini may have an elevated rate of diversification in contrast to other Neotropical cichlids. These findings provide insight into our understanding of the diversification patterns across taxonomically disparate lineages in this diverse clade of freshwater fishes and one of the most species-rich families of vertebrates.

Colour variation in cichlid fish: developmental mechanisms, selective pressures and evolutionary consequences

Cichlid fishes constitute one of the most species-rich families of vertebrates. In addition to complex social behaviour and morphological versatility, they are characterised by extensive diversity in colouration, both within and between species. Here, we review the cellular and molecular mechanisms underlying colour variation in this group and the selective pressures responsible for the observed variation. We specifically address the evidence for the hypothesis that divergence in colouration is associated with the evolution of reproductive isolation between lineages. While we conclude that cichlid colours are excellent models for understanding the role of animal communication in species divergence, we also identify taxonomic and methodological biases in the current research effort. We suggest that the integration of genomic approaches with ecological and behavioural studies, across the entire cichlid family and beyond it, will contribute to the utility of the cichlid model system for understanding the evolution of biological diversity.