Convergent Evolution of Cichlid Fish Pharyngeal Jaw Dentitions in Mollusk-Crushing Predators: Comparative X-Ray Computed Tomography of Tooth Sizes, Numbers, and Replacement

Dental convergence is a hallmark of cichlid fish adaptive radiations. This type of repeated evolution characterizes both the oral jaws of these fishes as well as their pharyngeal jaws that are modified gill arches used to functionally process prey like hard-shelled mollusks. To test several hypotheses regarding the evolution of cichlid crushing pharyngeal dentitions, we used X-ray computed tomography scans to comparatively examine dental evolution in the pharyngeal jaw of a diversity of New World Heroine cichlid lineages. The substantial variation in erupted tooth sizes and numbers as well as replacement teeth found in these fishes showed several general patterns. Larger toothed species tended to have fewer teeth suggesting a potential role of spatial constraints in cichlid dental divergence. Species with larger numbers of erupted pharyngeal teeth also had larger numbers of replacement teeth. Replacement tooth size is almost exactly predicted (r = 0.99) from the size of erupted teeth across all of the species. Mollusk crushing was, therefore, highly associated with not only larger pharyngeal teeth, but also larger replacement teeth. Whether dental divergence arises as a result of environmental induced plasticity or originates via trophic polymorphism as found in the species Herichthys minckleyi, there appear to be general rules that structure interspecific divergence in cichlid pharyngeal erupted and replacement dentitions.

Ecotype differences in aggression, neural activity and behaviorally relevant gene expression in cichlid fish

In Lake Malawi, two ecologically distinct lineages of cichlid fishes (rock- vs sand-dwelling ecotypes, each comprised of over 200 species) evolved within the last million years. The rock-dwelling species (Mbuna) are aggressively territorial year-round and males court and spawn with females over rocky substrate. In contrast, males of sand-dwelling species are not territorial and instead aggregate on seasonal breeding leks in which males construct courtship "bowers" in the sand. However, little is known about how phenotypic variation in aggression is produced by the genome. In this study, we first quantify and compare behavior in seven cichlid species, demonstrating substantial ecotype and species differences in unconditioned mirror-elicited aggression. Second, we compare neural activity in mirror-elicited aggression in two representative species, Mchenga conophoros (sand-dwelling) and Petrotilapia chitimba (rock-dwelling). Finally, we compare gene expression patterns between these two species, specifically within neurons activated during mirror aggression. We identified a large number of genes showing differential expression in mirror-elicited aggression, as well as many genes that differ between ecotypes. These genes, which may underly species differences in behavior, include several neuropeptides, genes involved in the synthesis of steroid hormones and neurotransmitter activity. This work lays the foundation for future experiments using this emerging genetic model system to investigate the genomic basis of evolved species differences in both brain and behavior.

Genome-wide protein phylogenies for four African cichlid species

Background

The thousands of species of closely related cichlid fishes in the great lakes of East Africa are a powerful model for understanding speciation and the genetic basis of trait variation. Recently, the genomes of five species of African cichlids representing five distinct lineages were sequenced and used to predict protein products at a genome-wide level. Here we characterize the evolutionary relationship of each cichlid protein to previously sequenced animal species.

Results

We used the Treefam database, a set of preexisting protein phylogenies built using 109 previously sequenced genomes, to identify Treefam families for each protein annotated from four cichlid species: Metriaclima zebra, Astatotilapia burtoni, Pundamilia nyererei and Neolamporologus brichardi. For each of these Treefam families, we built new protein phylogenies containing each of the cichlid protein hits. Using these new phylogenies we identified the evolutionary relationship of each cichlid protein to its nearest human and zebrafish protein. This data is available either through download or through a webserver we have implemented.

Conclusion

These phylogenies will be useful for any cichlid researchers trying to predict biological and protein function for a given cichlid gene, understanding the evolutionary history of a given cichlid gene, identifying recently duplicated cichlid genes, or performing genome-wide analysis in cichlids that relies on using databases generated from other species.

Electronic supplementary material

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

A compendium of developmental gene expression in Lake Malawi cichlid fishes

BACKGROUND

Lake Malawi cichlids represent one of a growing number of vertebrate models used to uncover the genetic and developmental basis of trait diversity. Rapid evolutionary radiation has resulted in species that share similar genomes but differ markedly in phenotypes including brains and behavior, nuptial coloration and the craniofacial skeleton. Research has begun to identify the genes, as well as the molecular and developmental pathways that underlie trait divergence.

RESULTS

We assemble a compendium of gene expression for Lake Malawi cichlids, across pharyngula (the phylotypic stage) and larval stages of development, encompassing hundreds of gene transcripts. We chart patterns of expression in Bone morphogenetic protein (BMP), Fibroblast growth factor (FGF), Hedgehog (Hh), Notch and Wingless (Wnt) signaling pathways, as well as genes involved in neurogenesis, calcium and endocrine signaling, stem cell biology, and numerous homeobox (Hox) factors-in three planes using whole-mount in situ hybridization. Because of low sequence divergence across the Malawi cichlid assemblage, the probes we employ are broadly applicable in hundreds of species. We tabulate gene expression across general tissue domains, and highlight examples of unexpected expression patterns.

CONCLUSIONS

On the heels of recently published genomes, this compendium of developmental gene expression in Lake Malawi cichlids provides a valuable resource for those interested in the relationship between evolution and development.

Quantitative multivariate analysis of dynamic multicellular morphogenic trajectories

Interrogating fundamental cell biology principles that govern tissue morphogenesis is critical to better understanding of developmental biology and engineering novel multicellular systems. Recently, functional micro-tissues derived from pluripotent embryonic stem cell (ESC) aggregates have provided novel platforms for experimental investigation; however elucidating the factors directing emergent spatial phenotypic patterns remains a significant challenge. Computational modelling techniques offer a unique complementary approach to probe mechanisms regulating morphogenic processes and provide a wealth of spatio-temporal data, but quantitative analysis of simulations and comparison to experimental data is extremely difficult. Quantitative descriptions of spatial phenomena across multiple systems and scales would enable unprecedented comparisons of computational simulations with experimental systems, thereby leveraging the inherent power of computational methods to interrogate the mechanisms governing emergent properties of multicellular biology. To address these challenges, we developed a portable pattern recognition pipeline consisting of: the conversion of cellular images into networks, extraction of novel features via network analysis, and generation of morphogenic trajectories. This novel methodology enabled the quantitative description of morphogenic pattern trajectories that could be compared across diverse systems: computational modelling of multicellular structures, differentiation of stem cell aggregates, and gastrulation of cichlid fish. Moreover, this method identified novel spatio-temporal features associated with different stages of embryo gastrulation, and elucidated a complex paracrine mechanism capable of explaining spatiotemporal pattern kinetic differences in ESC aggregates of different sizes.

Developmental Plasticity of Patterned and Regenerating Oral Organs

In many aquatic vertebrates, including bony and cartilaginous fishes, teeth and taste buds colocalize on jaw elements. In these animals, taste buds are renewed continuously throughout life, whereas teeth undergo cycled whole-organ replacement by various means. Recently, studies of cichlid fishes have yielded new insights into the development and regeneration of these dental and sensory oral organs. Tooth and taste bud densities covary positively across species with different feeding strategies, controlled by common regions of the genome and integrated molecular signals. Developing teeth and taste buds share a bipotent epithelium during early patterning stages, from which dental and taste fields are specified. Moreover, these organs share a common epithelial ribbon that supports label-retaining cells during later stages of regeneration. During both patterning and regeneration stages, dental organs can be converted to taste bud fate by manipulation of BMP signaling. These observations highlight a surprising long-term plasticity between dental and sensory organ types. Here, we review these findings and discuss the implications of developmental plasticity that spans the continuum of craniofacial organ patterning and regeneration.

Lake Malawi cichlid evolution along a benthic/limnetic axis

Divergence along a benthic to limnetic habitat axis is ubiquitous in aquatic systems. However, this type of habitat divergence has largely been examined in low diversity, high latitude lake systems. In this study, we examined the importance of benthic and limnetic divergence within the incredibly species-rich radiation of Lake Malawi cichlid fishes. Using novel phylogenetic reconstructions, we provided a series of hypotheses regarding the evolutionary relationships among 24 benthic and limnetic species that suggests divergence along this axis has occurred multiple times within Lake Malawi cichlids. Because pectoral fin morphology is often associated with divergence along this habitat axis in other fish groups, we investigated divergence in pectoral fin muscles in these benthic and limnetic cichlid species. We showed that the eight pectoral fin muscles and fin area generally tended to evolve in a tightly correlated manner in the Lake Malawi cichlids. Additionally, we found that larger pectoral fin muscles are strongly associated with the independent evolution of the benthic feeding habit across this group of fish. Evolutionary specialization along a benthic/limnetic axis has occurred multiple times within this tropical lake radiation and has produced repeated convergent matching between exploitation of water column habitats and locomotory morphology.

Common developmental pathways link tooth shape to regeneration

In many non-mammalian vertebrates, adult dentitions result from cyclical rounds of tooth regeneration wherein simple unicuspid teeth are replaced by more complex forms. Therefore and by contrast to mammalian models, the numerical majority of vertebrate teeth develop shape during the process of replacement. Here, we exploit the dental diversity of Lake Malawi cichlid fishes to ask how vertebrates generally replace their dentition and in turn how this process acts to influence resulting tooth morphologies. First, we used immunohistochemistry to chart organogenesis of continually replacing cichlid teeth and discovered an epithelial down-growth that initiates the replacement cycle via a labial proliferation bias. Next, we identified sets of co-expressed genes from common pathways active during de novo, lifelong tooth replacement and tooth morphogenesis. Of note, we found two distinct epithelial cell populations, expressing markers of dental competence and cell potency, which may be responsible for tooth regeneration. Related gene sets were simultaneously active in putative signaling centers associated with the differentiation of replacement teeth with complex shapes. Finally, we manipulated targeted pathways (BMP, FGF, Hh, Notch, Wnt/β-catenin) in vivo with small molecules and demonstrated dose-dependent effects on both tooth replacement and tooth shape. Our data suggest that the processes of tooth regeneration and tooth shape morphogenesis are integrated via a common set of molecular signals. This linkage has subsequently been lost or decoupled in mammalian dentitions where complex tooth shapes develop in first generation dentitions that lack the capacity for lifelong replacement. Our dissection of the molecular mechanics of vertebrate tooth replacement coupled to complex shape pinpoints aspects of odontogenesis that might be re-evolved in the lab to solve problems in regenerative dentistry.

Origins of shared genetic variation in African cichlids

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

The genetic basis of a complex functional system

The relationship between form and function can have profound effects on evolutionary dynamics and such effects may differ for simple versus complex systems. In particular, functions produced by multiple structural configurations (many-to-one mapping, MTOM) may dampen constituent trade-offs and promote diversification. Unfortunately, we lack information about the genetic architecture of MTOM functional systems. The skulls of teleost fishes contain both simple (lower jaw levers) as well as more complex (jaws modeled as 4-bar linkages) functional systems within the same craniofacial unit. We examined the mapping of form to function and the genetic basis of these systems by identifying quantitative trait loci (QTL) in hybrids of two Lake Malawi cichlid species. Hybrid individuals exhibited novelty (transgressive segregation) in morphological components and function of the simple and complex jaw systems. Functional novelty was proportional to the prevalence of extreme morphologies in the simple levers; by contrast, recombination of parental morphologies produced transgression in the MTOM 4-bar linkage. We found multiple loci of moderate effect and epistasis controlling jaw phenotypes in both the simple and complex systems, with less phenotypic variance explained by QTL for the 4-bar. Genetic linkage between components of the simple and complex systems partly explains phenotypic correlations and may constrain functional evolution.

Divergence in cis-regulatory sequences surrounding the opsin gene arrays of African cichlid fishes

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.

The macroecology of rapid evolutionary radiation

A long-standing debate in ecology addresses whether community composition is the result of stochastic factors or assembly rules. Non-random, over-dispersed patterns of species co-occurrence have commonly been attributed to competition--a particularly important force in adaptive radiation. We thus examined the macroecology of the recently radiated cichlid rock-fish assemblage in Lake Malawi, Africa at a spectrum of increasingly fine spatial scales (entire lake to depth within rock-reef sites). Along this range of spatial scales, we observed a signal of community structure (decreased co-occurrence of species) at the largest and smallest scales, but not in between. Evidence suggests that the lakewide signature of structure is driven by extreme endemism and micro-allopatric speciation, while patterns of reduced co-occurrence with depth are indicative of species interactions. We identified a 'core' set of rock-reef species, found in combination throughout the lake, whose depth profiles exhibited replicated positive and negative correlation. Our results provide insight into how ecological communities may be structured differently at distinct spatial scales, re-emphasize the importance of local species interactions in community assembly, and further elucidate the processes shaping speciation in this model adaptive radiation.

Evolution of microRNAs and the diversification of species

MicroRNAs (miRNAs) are ancient, short noncoding RNA molecules that regulate the transcriptome through post-transcriptional mechanisms. miRNA riboregulation is involved in a diverse range of biological processes, and misregulation is implicated in disease. It is generally thought that miRNAs function to canalize cellular outputs, for instance as “fail-safe” repressors of gene misexpression. Genomic surveys in humans have revealed reduced genetic polymorphism and the signature of negative selection for both miRNAs themselves and the target sequences to which they are predicted to bind. We investigated the evolution of miRNAs and their binding sites across cichlid fishes from Lake Malawi (East Africa), where hundreds of diverse species have evolved in the last million years. Using low-coverage genome sequence data, we identified 100 cichlid miRNA genes with mature regions that are highly conserved in other animal species. We computationally predicted target sites on the 3′-untranslated regions (3′-UTRs) of cichlid genes to which miRNAs may bind and found that these sites possessed elevated single nucleotide polymorphism (SNP) densities. Furthermore, polymorphic sites in predicted miRNA targets showed higher minor allele frequencies on average and greater genetic differentiation between Malawi lineages when compared with a neutral expectation and nontarget 3′-UTR SNPs. Our data suggest that divergent selection on miRNA riboregulation may have contributed to the diversification of cichlid species and may similarly play a role in rapid phenotypic evolution of other natural systems.

The odontode explosion: the origin of tooth-like structures in vertebrates

Essentially we show recent data to shed new light on the thorny controversy of how teeth arose in evolution. Essentially we show (a) how teeth can form equally from any epithelium, be it endoderm, ectoderm or a combination of the two and (b) that the gene expression programs of oral versus pharyngeal teeth are remarkably similar. Classic theories suggest that (i) skin denticles evolved first and odontode-inductive surface ectoderm merged inside the oral cavity to form teeth (the 'outside-in' hypothesis) or that (ii) patterned odontodes evolved first from endoderm deep inside the pharyngeal cavity (the 'inside-out' hypothesis). We propose a new perspective that views odontodes as structures sharing a deep molecular homology, united by sets of co-expressed genes defining a competent thickened epithelium and a collaborative neural crest-derived ectomesenchyme. Simply put, odontodes develop 'inside and out', wherever and whenever these co-expressed gene sets signal to one another. Our perspective complements the classic theories and highlights an agenda for specific experimental manipulations in model and non-model organisms.

Comparative rates of lower jaw diversification in cichlid adaptive radiations

The lower jaw (LJ) provides an ideal trophic phenotype to compare rates and patterns of macroevolution among cichlid radiations. Using a novel phylogeny of four genes (ND2, dlx2, mitfb, and s7), we examined the evolutionary relationships among two of the most phylogenetically disparate cichlid radiations: (i) the Central America Heroines; and (ii) the East African Lake Malawi flock. To quantify jaw morphology, we measured two LJ lever systems in approximately 40 species from each lineage. Using geologic calibrations, we generated a chronogram for both groups and examined the rates of jaw evolution in the two radiations. The most rapidly evolving components of the LJ differed between the two radiations. However, the Lake Malawi flock exhibited a much faster rate of evolution in several components of the LJ. This rapid rate of divergence is consistent with natural selection, promoting unparalleled trophic diversification in Lake Malawi cichlids.

Evolvable network architectures

There is significant research interest recently to understand the evolution of the current Internet, as well as to design clean-slate Future Internet architectures. Clearly, even when network architectures are designed from scratch, they have to evolve as their environment (i.e., technological constraints, service requirements, applications, economic conditions, etc) always changes. A key question then is: what makes a network architecture evolvable? What determines the ability of a network architecture to evolve as its environment changes? In this paper, we review some relevant ideas about evolvability from the biological literature. We examine the role of robustness and modularity in evolution, and their relation with evolvability. We also discuss evolutionary kernels and punctuated equilibria, two important concepts that may be relevant to the so-called ossification of the core Internet protocols. Finally, we examine optimality, a design objective that is often of primary interest in engineering but that does not seem to be abundant in biology.

An ancient gene network is co-opted for teeth on old and new jaws

Vertebrate dentitions originated in the posterior pharynx of jawless fishes more than half a billion years ago. As gnathostomes (jawed vertebrates) evolved, teeth developed on oral jaws and helped to establish the dominance of this lineage on land and in the sea. The advent of oral jaws was facilitated, in part, by absence of hox gene expression in the first, most anterior, pharyngeal arch. Much later in evolutionary time, teleost fishes evolved a novel toothed jaw in the pharynx, the location of the first vertebrate teeth. To examine the evolutionary modularity of dentitions, we asked whether oral and pharyngeal teeth develop using common or independent gene regulatory pathways. First, we showed that tooth number is correlated on oral and pharyngeal jaws across species of cichlid fishes from Lake Malawi (East Africa), suggestive of common regulatory mechanisms for tooth initiation. Surprisingly, we found that cichlid pharyngeal dentitions develop in a region of dense hox gene expression. Thus, regulation of tooth number is conserved, despite distinct developmental environments of oral and pharyngeal jaws; pharyngeal jaws occupy hox-positive, endodermal sites, and oral jaws develop in hox-negative regions with ectodermal cell contributions. Next, we studied the expression of a dental gene network for tooth initiation, most genes of which are similarly deployed across the two disparate jaw sites. This collection of genes includes members of the ectodysplasin pathway, eda and edar, expressed identically during the patterning of oral and pharyngeal teeth. Taken together, these data suggest that pharyngeal teeth of jawless vertebrates utilized an ancient gene network before the origin of oral jaws, oral teeth, and ectodermal appendages. The first vertebrate dentition likely appeared in a hox-positive, endodermal environment and expressed a genetic program including ectodysplasin pathway genes. This ancient regulatory circuit was co-opted and modified for teeth in oral jaws of the first jawed vertebrate, and subsequently deployed as jaws enveloped teeth on novel pharyngeal jaws. Our data highlight an amazing modularity of jaws and teeth as they coevolved during the history of vertebrates. We exploit this diversity to infer a core dental gene network, common to the first tooth and all of its descendants.

Natural selection governs local, but not global, evolutionary gene coexpression networks in Caenorhabditis elegans

Background

Large-scale evaluation of gene expression variation among Caenorhabditis elegans lines that have diverged from a common ancestor allows for the analysis of a novel class of biological networks – evolutionary gene coexpression networks. Comparative analysis of these evolutionary networks has the potential to uncover the effects of natural selection in shaping coexpression network topologies since C. elegans mutation accumulation (MA) lines evolve essentially free from the effects of natural selection, whereas natural isolate (NI) populations are subject to selective constraints.

Results

We compared evolutionary gene coexpression networks for C. elegans MA lines versus NI populations to evaluate the role that natural selection plays in shaping the evolution of network topologies. MA and NI evolutionary gene coexpression networks were found to have very similar global topological properties as measured by a number of network topological parameters. Observed MA and NI networks show node degree distributions and average values for node degree, clustering coefficient, path length, eccentricity and betweeness that are statistically indistinguishable from one another yet highly distinct from randomly simulated networks. On the other hand, at the local level the MA and NI coexpression networks are highly divergent; pairs of genes coexpressed in the MA versus NI lines are almost entirely different as are the connectivity and clustering properties of individual genes.

Conclusion

It appears that selective forces shape how local patterns of coexpression change over time but do not control the global topology of C. elegans evolutionary gene coexpression networks. These results have implications for the evolutionary significance of global network topologies, which are known to be conserved across disparate complex systems.

A periodic pattern generator for dental diversity

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.

Comparative analysis reveals signatures of differentiation amid genomic polymorphism in Lake Malawi cichlids

Low coverage survey sequencing shows that although Lake Malawi cichlids are phenotypically and behaviorally diverse, they appear genetically like a subdivided population.

Background

Cichlid fish from East Africa are remarkable for phenotypic and behavioral diversity on a backdrop of genomic similarity. In 2006, the Joint Genome Institute completed low coverage survey sequencing of the genomes of five phenotypically and ecologically diverse Lake Malawi species. We report a computational and comparative analysis of these data that provides insight into the mechanisms that make closely related species different from one another.

Results

We produced assemblies for the five species ranging in aggregate length from 68 to 79 megabase pairs, identified putative orthologs for more than 12,000 human genes, and predicted more than 32,000 cross-species single nucleotide polymorphisms (SNPs). Nucleotide diversity was lower than that found among laboratory strains of the zebrafish. We collected around 36,000 genotypes to validate a subset of SNPs within and among populations and across multiple individuals of about 75 Lake Malawi species. Notably, there were no fixed differences observed between focal species nor between major lineages. Roughly 3% to 5% of loci surveyed are statistical outliers for genetic differentiation (F_ST) within species, between species, and between major lineages. Outliers for F_ST are candidate genes that may have experienced a history of natural selection in the Malawi lineage.

Conclusion

We present a novel genome sequencing strategy, which is useful when evolutionary diversity is the question of interest. Lake Malawi cichlids are phenotypically and behaviorally diverse, but they appear genetically like a subdivided population. The unique structure of Lake Malawl cichlid genomes should facilitate conceptually new experiments, employing SNPs to identity genotype-phenotype association, using the entire species flock as a mapping panel.

Visual sensitivities tuned by heterochronic shifts in opsin gene expression

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.

Hybridization produces novelty when the mapping of form to function is many to one

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, F₂ 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.

Evolution and development of complex biomechanical systems: 300 million years of fish jaws

The jaws of teleost fishes are diverse and complex musculoskeletal systems. The focus in this review is on the major biomechanical systems in the teleost head, and the range and interplay of functional, developmental, and genetic influences that shape the modular and integrated evolution of elements. Insights possible from comparative studies are discussed in the context of traditional and new models for studies of craniofacial evolution and development.

Do constructional constraints influence cichlid craniofacial diversification?

Constraints on form should determine how organisms diversify. Owing to competition for the limited space within the body, investment in adjacent structures may frequently represent an evolutionary compromise. For example, evolutionary trade-offs between eye size and jaw muscles in cichlid fish of the African great lakes are thought to represent a constructional constraint that influenced the diversification of these assemblages. To test the evolutionary independence of these structures in Lake Malawi cichlid fish, we measured the mass of the three major adductor mandibulae (AM) muscles and determined the eye volume in 41 species. Using both traditional and novel methodologies to control for resolved and unresolved phylogenetic relationships, we tested the evolutionary independence of these four structures. We found that evolutionary change in the AM muscles was positively correlated, suggesting that competition for space in the head has not influenced diversification among these jaw muscles. Furthermore, there was no negative relationship between change in total AM muscle mass and eye volume, indicating that there has been little effect of the evolution of eye size on AM evolution in Lake Malawi cichlids. The comparative approach used here should provide a robust method to test whether constructional constraints frequently limit phenotypic change in adaptive radiations.

Genome size is negatively correlated with effective population size in ray-finned fish

A recent theory suggesting that genome size and complexity can increase as a passive consequence of small effective population size has generated much controversy. In this article, we demonstrate that freshwater fish species, which have smaller effective population sizes than marine fish species, have larger genomes. We show that genome size is negatively correlated with genetic variability, independent of phylogeny, body size and generation time. Genome duplication is also observed predominantly in freshwater fish. These results suggest that the raw materials of complexity originate under conditions of reduced selection efficiency.

Integration and evolution of the cichlid mandible: the molecular basis of alternate feeding strategies

African cichlid fishes have repeatedly evolved highly specialized modes of feeding through adaptations of their oral jaws. Here, we explore the molecular genetic basis of the opening and closing lever mechanisms of the cichlid lower jaw, which have traditionally been used to describe the mechanics of feeding behavior in bony fishes. Quantitative genetic analyses demonstrate that the opening and closing mechanisms are genetically modular and therefore free to evolve independently. Bmp4 (bone morphogenetic protein 4) is one of two loci that segregate with the mechanical advantage of closing and that together account for >30% of the phenotypic variance in this trait. Species-specific differences in jaw shape are obvious early in cichlid larval development and are correlated with patterns of bmp4 expression in the mandibular primordium. When bmp4 is overexpressed in the obligate suction feeder Danio rerio, mandibular morphology exhibits specific transformations of opening and closing lever ratios. We conclude that patterns of morphological integration of the cichlid jaw reflect a balance among conflicting functional demands. Further, we demonstrate that bmp4 has the potential to alter mandibular morphology in a way that mimics adaptive variation among fish species.

The transcriptional consequences of mutation and natural selection in Caenorhabditis elegans

The evolutionary importance of gene-expression divergence is unclear: some studies suggest that it is an important mechanism for evolution by natural selection, whereas others claim that most between-species regulatory changes are neutral or nearly neutral. We examined global transcriptional divergence patterns in a set of Caenorhabditis elegans mutation-accumulation lines and natural isolate lines to provide insights into the evolutionary importance of transcriptional variation and to discriminate between the forces of mutation and natural selection in shaping the evolution of gene expression. We detected the effects of selection on transcriptional divergence patterns and characterized them with respect to coexpressed gene sets, chromosomal clustering of expression changes and functional gene categories. We directly compared observed transcriptional variation patterns in the mutation-accumulation and natural isolate lines to a neutral model of transcriptome evolution to show that strong stabilizing selection dominates the evolution of transcriptional change for thousands of C. elegans expressed sequences.

A second-generation genetic linkage map of tilapia (Oreochromis spp.)

We constructed a second-generation linkage map of tilapia from the F(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.

The cusp of evolution and development: a model of cichlid tooth shape diversity

Tooth shape is a hallmark of repeated evolutionary radiations among cichlid fishes from East Africa. Cusp shape and number vary both within populations and among closely related species with different feeding behaviors and ecologies. Here, we use histology and scanning electron microscopy to chart the developmental trajectory of tooth shape differences in fishes from Lake Malawi. We demonstrate that species with bi- or tricuspid adult (replacement) teeth initially possess a first-generation unicuspid dentition. Notably, the timing of turnover from first-generation to replacement teeth differs among species and is correlated with feeding ecology. Next, we use field data for cichlid species with adult unicuspid, bicuspid, and tricuspid teeth to demonstrate a strong and positive relationship between the number of teeth in a row and tooth shape. We discuss cichlid tooth ontogeny in the context of morphogenetic models designed to explain the developmental basis of tooth shape variation in mammals. We suggest that the dramatic differences in cichlid dentitions can be explained by variation in the expression of common activators and inhibitors acting at multiple stages of odontogenesis.

Hybridization and contemporary evolution in an introduced cichlid fish from Lake Malawi National Park

Rapidly evolving systems offer the chance to observe genetic and phenotypic change in real time. We exploit a well-characterized introduction of cichlid fish into Lake Malawi National Park to document a short history of habitat colonization and the evolution of genes and colour pattern. In the early 1960s, a fish exporter introduced individuals of Cynotilapia afra to a single site (Mitande Point) of Thumbi West Island and, as late as 1983, the species was confined to this location. In 2001, C. afra had colonized the entire perimeter of Thumbi West. In July of that year, we sampled C. afra individuals from six sites around the island and scored variation in dorsal fin colour as well as allelic diversity at six microsatellite loci. We found that, in two decades, C. afra had diverged into genetically distinct, phenotypically different northern and southern populations. We observed a high proportion of hybrids between the introduced C. afra and the native Metriaclima zebra on the southern coast of Thumbi West, and speculate that hybridization is facilitated by low water clarity at these windward sites. The short history of C. afra at Thumbi West is a microcosm of contemporary evolutionary divergence and may provide the opportunity to study the process from start to finish in genetic detail.

Genome mapping of the orange blotch colour pattern in cichlid fishes

The dramatic variation of cichlid fish colour pattern is thought to function in mate choice, evolve by sexual selection, and contribute to explosive speciation. Here, we combine linkage mapping and population genetic analyses to identify a single region of the cichlid genome responsible for the orange blotch (OB) colour phenotype. In each analysis, OB is tightly linked to the c-ski1 gene. Additionally, we use comparative mapping information from the Takifugu rubripes and human genomes to suggest positional candidate loci for OB. Our work should engender a more comprehensive understanding of the molecular ecology of OB and its role in cichlid speciation. Moreover, we have assembled the components of a method to focus upon the genetic basis of evolutionarily and ecologically significant phenotypes.

Genetic basis of adaptive shape differences in the cichlid head

East African cichlids exhibit an extraordinary level of morphological diversity. Key to their success has been a dramatic radiation in trophic biology, which has occurred rapidly and repeatedly in different lakes. In this report we take the first step in understanding the genetic basis of differences in cichlid oral jaw design. We estimate the effective number of genetic factors that control differences in the cichlid head through a comprehensive morphological assessment of two Lake Malawi cichlid species and their F(1) and F(2) hybrid progeny. We estimate that between one and 11 factors underlie shape difference of individual bony elements. We show that many of the skeletal differences in the head and oral jaw apparatus are inherited together, suggesting a degree of pleiotropy in the genetic architecture of this character complex. Moreover, we find that cosegregation of shape differences in different elements corresponds to developmental, rather than functional, units.

Directional selection has shaped the oral jaws of Lake Malawi cichlid fishes

East African cichlid fishes represent one of the most striking examples of rapid and convergent evolutionary radiation among vertebrates. Models of ecological speciation would suggest that functional divergence in feeding morphology has contributed to the origin and maintenance of cichlid species diversity. However, definitive evidence for the action of natural selection has been missing. Here we use quantitative genetics to identify regions of the cichlid genome responsible for functionally important shape differences in the oral jaw apparatus. The consistent direction of effects for individual quantitative trait loci suggest that cichlid jaws and teeth evolved in response to strong, divergent selection. Moreover, several chromosomal regions contain a disproportionate number of quantitative trait loci, indicating a prominent role for pleiotropy or genetic linkage in the divergence of this character complex. Of particular interest are genomic intervals with concerted effects on both the length and height of the lower jaw. Coordinated changes in this area of the oral jaw apparatus are predicted to have direct consequences for the speed and strength of jaw movement. Taken together, our results imply that the rapid and replicative nature of cichlid trophic evolution is the result of directional selection on chromosomal packages that encode functionally linked aspects of the craniofacial skeleton.

Evolutionary history of the parrotfishes: biogeography, ecomorphology, and comparative diversity

The family Scaridae comprises about 90 species of herbivorous coral reef, rock reef, and seagrass fishes. Parrotfishes are important agents of marine bioerosion who rework the substrate with their beaklike oral jaws. Many scarid populations are characterized by complex social systems including highly differentiated sexual stages, territoriality, and the defense of harems. Here, we test a hypothesis of relationships among parrotfish genera derived from nearly 2 kb of nuclear and mitochondrial DNA sequence. The DNA tree is different than a phylogeny based on comparative morphology and leads to important reinterpretations of scarid evolution. The molecular data suggest a split among seagrass and coral reef associated genera with nearly 80% of all species in the coral reef clade. Our phylogenetic results imply an East Tethyan origin of the family and the recurrent evolution of excavating and scraping feeding modes. It is likely that ecomorphological differences played a significant role in the initial divergence of major scarid lineages, but that variation in color and breeding behavior has triggered subsequent diversification. We present a two-phase model of parrotfish evolution to explain patterns of comparative diversity. Finally, we discuss the application of this model to other adaptively radiating clades.

Microsatellite variation associated with prolactin expression and growth of salt-challenged tilapia

Biologists have long argued that runs of alternating purines and pyrimidines could form alternative DNA structures, which might regulate transcription. Here, we report that simple sequence repeat polymorphisms in the tilapia prolactin 1 (prl 1) promoter are associated with differences in prl 1 gene expression and the growth response of salt-challenged fishes. Individuals homozygous for long microsatellite alleles express less prl 1 in fresh water but more prl 1 in half-seawater than fishes with other genotypes. Our work provides the first in vivo evidence that differences in microsatellite length among individuals may indeed affect gene expression and that variance in expression has concomitant physiological consequences. These results suggest that dinucleotide microsatellites represent an under-appreciated source of genetic variation for regulatory evolution.

Rapid isolation of CA microsatellites from the tilapia genome

We have developed (CA)n microsatellite markers for the cichlid fish, Oreochromis niloticus using a variation of the hybrid capture method. The resulting genomic library was highly enriched in repetitive DNA with 96% of clones containing CA repeats. The number of repeats ranged from four to 45 with an average of 19. Two-thirds of the sequenced clones had 12 or more repeats and sufficient flanking sequence to design primers. The resulting markers were tested in an F2 cross of O. niloticus x O. aureus. Nearly 90% of the markers amplified in this cross and 74% of these were informative. This work demonstrates the importance of minimizing the number of polymerase chain reaction (PCR) amplification cycles before and after the enrichment steps to reduce PCR recombination and the generation of chimaeric clones.

From phenotype to genotype

Advances in the field of genomics have made it possible to decipher the genetic and transcriptional changes that underlie differences among organisms. Here we discuss the merits and drawbacks of these strategies as they relate to evolutionary biology. We suggest that the molecular basis of natural variation can best be interpreted via the intersection of genomic and transcriptome approaches. We outline how this might be accomplished in practice and assemble the components of a method to comprehend the evolutionary interplay between genotype and phenotype.

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.

Reconstructing labroid evolution with single-copy nuclear DNA

Fifteen per cent of all living fishes are united in a single suborder (Labroidei) and display a dazzling array of behavioural and ecological traits. The labroids are considered monophyletic and members share a pharyngeal jaw apparatus (PJA) modified for crushing and processing prey. Outside of the explicitly functional PJA, there is no corroborative evidence for a monophyletic Labroidei. Here, we report the first molecular phylogenetic analysis of the suborder. Contrary to morphology-based phylogenies, our single-copy nuclear DNA data do not support labroid families as a natural group. Our data indicate that pharyngognathy has evolved independently among labroid families and that characters of the PJA are not reliable markers of perciform evolution. This work 'crushes' conventional views of fish phylogeny and should engender novel concepts of piscine life history evolution.

Evolutionary conservation of microsatellite flanking regions and their use in resolving the phylogeny of cichlid fishes (Pisces: Perciformes)

A phylogeny of the principal lineages of cichlid fishes and two other fish families of the suborder Labroidei was based on phylogenetic information from DNA sequences of the flanking region of a (CA)n microsatellite locus. Microsatellite (CA)n containing clones from a genomic library of an African cichlid fish from Lake Tanganyika, Tropheus moorii, were sequenced and primers for the polymerase chain reaction designed. All primers amplified the homologous microsatellite loci in many more than the source species and one microsatellite flanking locus (TmoM27) was particularly conserved and amplified in several lineages of perciform fishes that diverged more than 80-100 million years ago. Despite the extensive level of evolutionary conservation of this microsatellite flanking region (MFR), this nuclear region contained reliable phylogenetic information in the form of both point and length mutations. A phylogeny of cichlids based on this MFR agrees with other phylogenetic hypotheses based on morphological, mitochondrial, and anonymous nuclear DNA. Madagascan and Indian cichlids are found to be paraphyletic and the most basal group in the family Cichlidae. African and Neotropical cichlids are both monophyletic and sistergroups. Within African lineages, the East African cichlids are most likely to be monophyletic and the West African cichlids are probably paraphyletic and basal to all African species. The focal microsatellite locus contained much variation in (CA)n repeats in African cichlids and in surfperches (up to 64 repeats), but was short (with only 2-4 repeats) and almost invariant in Neotropical cichlids. The design of phylogenetically highly versatile MFR-primers will be of use not only for phylogeny reconstruction among families of perciform fishes, but also for population-level work in the thousands of species belonging to this highly species-rich suborder of fishes.