Search

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.

Abstract Background Divergence within cis-regulatory sequences may contribute to the adaptive evolution of gene expression, but functional alleles in these regions are difficult to identify without abundant genomic resources. Among African cichlid fishes, the differential expression of seven opsin genes has produced adaptive differences in visual sensitivity. Quantitative genetic analysis suggests that cis-regulatory alleles near the SWS2-LWS opsins may contribute to this variation. Here, we sequence BACs containing the opsin genes of two cichlids, Oreochromis niloticus and Metriaclima zebra. We use phylogenetic footprinting and shadowing to examine divergence in conserved non-coding elements, promoter sequences, and 3'-UTRs surrounding each opsin in search of candidate cis-regulatory sequences that influence cichlid opsin expression. Results We identified 20 conserved non-coding elements surrounding the opsins of cichlids and other teleosts, including one known enhancer and a retinal microRNA. Most conserved elements contained computationally-predicted binding sites that correspond to transcription factors that function in vertebrate opsin expression; O. niloticus and M. zebra were significantly divergent in two of these. Similarly, we found a large number of relevant transcription factor binding sites within each opsin's proximal promoter, and identified five opsins that were considerably divergent in both expression and the number of transcription factor binding sites shared between O. niloticus and M. zebra. We also found several microRNA target sites within the 3'-UTR of each opsin, including two 3'-UTRs that differ significantly between O. niloticus and M. zebra. Finally, we examined interspecific divergence among 18 phenotypically diverse cichlids from Lake Malawi for one conserved non-coding element, two 3'-UTRs, and five opsin proximal promoters. We found that all regions were highly conserved with some evidence of CRX transcription factor binding site turnover. We also found three SNPs within two opsin promoters and one non-coding element that had weak association with cichlid opsin expression. Conclusions This study is the first to systematically search the opsins of cichlids for putative cis-regulatory sequences. Although many putative regulatory regions are highly conserved across a large number of phenotypically diverse cichlids, we found at least nine divergent sequences that could contribute to opsin expression differences in cis and stand out as candidates for future functional analyses.

We constructed a second-generation linkage map of tilapia from the [F.sub.2] progeny of an interspecific cross between Oreochromis niloticus and Oreochromis aureus. The map reported here contains 525 microsatellite and 21 gene-based markers. It spans 1311 cM in 24 linkage groups, for an average marker spacing of 2.4 cM. We detected associations of sex and red color with markers on linkage group 3. This map will enable mapping and selective breeding of quantitative traits important to the economic culture of tilapia as a food fish and will contribute to the study of closely related cichlids that have undergone explosive adaptive radiation in the lakes of East Africa.

Abstract Background Divergence within cis-regulatory sequences may contribute to the adaptive evolution of gene expression, but functional alleles in these regions are difficult to identify without abundant genomic resources. Among African cichlid fishes, the differential expression of seven opsin genes has produced adaptive differences in visual sensitivity. Quantitative genetic analysis suggests that cis-regulatory alleles near the SWS2-LWS opsins may contribute to this variation. Here, we sequence BACs containing the opsin genes of two cichlids, Oreochromis niloticus and Metriaclima zebra. We use phylogenetic footprinting and shadowing to examine divergence in conserved non-coding elements, promoter sequences, and 3'-UTRs surrounding each opsin in search of candidate cis-regulatory sequences that influence cichlid opsin expression. Results We identified 20 conserved non-coding elements surrounding the opsins of cichlids and other teleosts, including one known enhancer and a retinal microRNA. Most conserved elements contained computationally-predicted binding sites that correspond to transcription factors that function in vertebrate opsin expression; O. niloticus and M. zebra were significantly divergent in two of these. Similarly, we found a large number of relevant transcription factor binding sites within each opsin's proximal promoter, and identified five opsins that were considerably divergent in both expression and the number of transcription factor binding sites shared between O. niloticus and M. zebra. We also found several microRNA target sites within the 3'-UTR of each opsin, including two 3'-UTRs that differ significantly between O. niloticus and M. zebra. Finally, we examined interspecific divergence among 18 phenotypically diverse cichlids from Lake Malawi for one conserved non-coding element, two 3'-UTRs, and five opsin proximal promoters. We found that all regions were highly conserved with some evidence of CRX transcription factor binding site turnover. We also found three SNPs within two opsin promoters and one non-coding element that had weak association with cichlid opsin expression. Conclusions This study is the first to systematically search the opsins of cichlids for putative cis-regulatory sequences. Although many putative regulatory regions are highly conserved across a large number of phenotypically diverse cichlids, we found at least nine divergent sequences that could contribute to opsin expression differences in cis and stand out as candidates for future functional analyses.

Abstract Background Periodic patterning of iterative structures is a fundamental process during embryonic organization and development. Studies have shown how gene networks are employed to pattern butterfly eyespots, fly bristles and vertebrate epithelial appendages such as teeth, feathers, hair and mammary glands. Despite knowledge of how these features are organized, little is known about how diversity in periodic patterning is generated in nature. We address this problem through the molecular analysis of oral jaw dental diversity in Lake Malawi cichlids, where closely related species exhibit from 1 to 20 rows of teeth, with total teeth counts ranging from around 10 to 700. Results We investigate the expression of conserved gene networks (involving bmp2, bmp4, eda, edar, fgf8, pax9, pitx2, runx2, shh and wnt7b) known to pattern iterative structures and teeth in other vertebrates. We show that spatiotemporal variation in expression pattern reflects adult morphological diversity among three closely related Malawi cichlid species. Combinatorial epithelial expression of pitx2 and shh appears to govern the competence both of initial tooth sites and future tooth rows. Epithelial wnt7b and mesenchymal eda are expressed in the inter-germ and inter-row regions, and likely regulate the spacing of these shh-positive units. Finally, we used chemical knockdown to demonstrate the fundamental role of hedgehog signalling and initial placode formation in the organization of the periodically patterned cichlid dental programme. Conclusion Coordinated patterns of gene expression differ among Malawi species and prefigure the future-ordered distribution of functional teeth of specific size and spacing. This variation in gene expression among species occurs early in the developmental programme for dental patterning. These data show how a complex multi-rowed vertebrate dentition is organized and how developmental tinkering of conserved gene networks during iterative pattern formation can impact upon the evolution of trophic novelty.

Abstract Background Cichlid fishes have radiated into hundreds of species in the Great Lakes of Africa. Brightly colored males display on leks and vie to be chosen by females as mates. Strong discrimination by females causes differential male mating success, rapid evolution of male color patterns and, possibly, speciation. In addition to differences in color pattern, Lake Malawi cichlids also show some of the largest known shifts in visual sensitivity among closely related species. These shifts result from modulated expression of seven cone opsin genes. However, the mechanisms for this modulated expression are unknown. Results In this work, we ask whether these differences might result from changes in developmental patterning of cone opsin genes. To test this, we compared the developmental pattern of cone opsin gene expression of the Nile tilapia, Oreochromis niloticus, with that of several cichlid species from Lake Malawi. In tilapia, quantitative polymerase chain reaction showed that opsin gene expression changes dynamically from a larval gene set through a juvenile set to a final adult set. In contrast, Lake Malawi species showed one of two developmental patterns. In some species, the expressed gene set changes slowly, either retaining the larval pattern or progressing only from larval to juvenile gene sets (neoteny). In the other species, the same genes are expressed in both larvae and adults but correspond to the tilapia adult genes (direct development). Conclusion Differences in visual sensitivities among species of Lake Malawi cichlids arise through heterochronic shifts relative to the ontogenetic pattern of the tilapia outgroup. Heterochrony has previously been shown to be a powerful mechanism for change in morphological evolution. We found that altering developmental expression patterns is also an important mechanism for altering sensory systems. These resulting sensory shifts will have major impacts on visual communication and could help drive cichlid speciation.

Abstract Background Evolutionary biologists want to explain the origin of novel features and functions. Two recent but separate lines of research address this question. The first describes one possible outcome of hybridization, called transgressive segregation, where hybrid offspring exhibit trait distributions outside of the parental range. The second considers the explicit mapping of form to function and illustrates manifold paths to similar function (called many to one mapping, MTOM) when the relationship between the two is complex. Under this scenario, functional novelty may be a product of the number of ways to elicit a functional outcome (i.e., the degree of MTOM). We fuse these research themes by considering the influence of MTOM on the production of transgressive jaw biomechanics in simulated hybrids between Lake Malawi cichlid species. Results We characterized the component links and functional output (kinematic transmission, KT) of the 4-bar mechanism in the oral jaws of Lake Malawi cichlids. We demonstrated that the input and output links, the length of the lower jaw and the length of the maxilla respectively, have consistent but opposing relationships with KT. Based on these data, we predicted scenarios in which species with different morphologies but similar KT (MTOM species) would produce transgressive function in hybrids. We used a simple but realistic genetic model to show that transgressive function is a likely outcome of hybridization among Malawi species exhibiting MTOM. Notably, F2 hybrids are transgressive for function (KT), but not the component links that contribute to function. In our model, transgression is a consequence of recombination and assortment among alleles specifying the lengths of the lower jaw and maxilla. Conclusion We have described a general and likely pervasive mechanism that generates functional novelty. Simulations of hybrid offspring among Lake Malawi cichlids exhibiting MTOM produce transgressive function in the majority of cases, and at appreciable frequency. Functional transgression (i) is a product of recombination and assortment between alleles controlling the lengths of the lower jaw and the maxilla, (ii) occurs in the absence of transgressive morphology, and (iii) can be predicted from the morphology of parents. Our genetic model can be tested by breeding Malawi cichlid hybrids in the laboratory and examining the resulting range of forms and functions.