Between a Rock and a Hard Polytomy: Phylogenomics of the Rock-Dwelling Mbuna Cichlids of Lake Malawi.

Whole genome sequences are beginning to revolutionize our understanding of phylogenetic relationships. Yet, even whole genome sequences can fail to resolve the evolutionary history of the most rapidly radiating lineages, where incomplete lineage sorting, standing genetic variation, introgression, and other factors obscure the phylogenetic history of the group. To overcome such challenges, one emerging strategy is to integrate results across different methods. Most such approaches have been implemented on reduced representation genomic data sets, but whole genomes should provide the maximum possible evidence approach. Here, we test the ability of single nucleotide polymorphisms extracted from whole genome resequencing data, implemented in an integrative genomic approach, to resolve key nodes in the phylogeny of the mbuna, rock-dwelling cichlid fishes of Lake Malawi, which epitomize the phylogenetic intractability that often accompanies explosive lineage diversification. This monophyletic radiation has diversified at an unparalleled rate into several hundred species in less than 2 million years. Using an array of phylogenomic methods, we consistently recovered four major clades of mbuna, but a large basal polytomy among them. Although introgression between clades apparently contributed to the challenge of phylogenetic reconstruction, reduction of the data set to nonintrogressed sites still did not help to resolve the basal polytomy. On the other hand, relationships among six congeneric species pairs were resolved without ambiguity, even in one case where existing data led us to predict that resolution would be difficult. We conclude that the bursts of diversification at the earliest stages of the mbuna radiation may be phylogenetically unresolvable, but other regions of the tree are phylogenetically clearly supported. Integration of multiple phylogenomic approaches will continue to increase confidence in relationships inferred from these and other whole-genome data sets. [Incomplete lineage sorting; introgression; linkage disequilibrium; multispecies coalescence; rapid radiation; soft polytomy.].

The repeated evolution of stripe patterns is correlated with body morphology in the adaptive radiations of East African cichlid fishes.

Color patterns are often linked to the behavioral and morphological characteristics of an animal, contributing to the effectiveness of such patterns as antipredatory strategies. Species-rich adaptive radiations, such as the freshwater fish family Cichlidae, provide an exciting opportunity to study trait correlations at a macroevolutionary scale. Cichlids are also well known for their diversity and repeated evolution of color patterns and body morphology. To study the evolutionary dynamics between color patterns and body morphology, we used an extensive dataset of 461 species. A phylogenetic supertree of these species shows that stripe patterns evolved ~70 times independently and were lost again ~30 times. Moreover, stripe patterns show strong signs of correlated evolution with body elongation, suggesting that the stripes' effectiveness as antipredatory strategy might differ depending on the body shape. Using pedigree-based analyses, we show that stripes and body elongation segregate independently, indicating that the two traits are not genetically linked. Their correlation in nature is therefore likely maintained by correlational selection. Lastly, by performing a mate preference assay using a striped CRISPR-Cas9 mutant of a nonstriped species, we show that females do not differentiate between striped CRISPR mutant males and nonstriped wild-type males, suggesting that these patterns might be less important for species recognition and mate choice. In summary, our study suggests that the massive rates of repeated evolution of stripe patterns are shaped by correlational selection with body elongation, but not by sexual selection.

The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes.

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Contrasting signatures of genomic divergence during sympatric speciation.

The transition from 'well-marked varieties' of a single species into 'well-defined species'-especially in the absence of geographic barriers to gene flow (sympatric speciation)-has puzzled evolutionary biologists ever since Darwin1,2. Gene flow counteracts the buildup of genome-wide differentiation, which is a hallmark of speciation and increases the likelihood of the evolution of irreversible reproductive barriers (incompatibilities) that complete the speciation process3. Theory predicts that the genetic architecture of divergently selected traits can influence whether sympatric speciation occurs4, but empirical tests of this theory are scant because comprehensive data are difficult to collect and synthesize across species, owing to their unique biologies and evolutionary histories5. Here, within a young species complex of neotropical cichlid fishes (Amphilophus spp.), we analysed genomic divergence among populations and species. By generating a new genome assembly and re-sequencing 453 genomes, we uncovered the genetic architecture of traits that have been suggested to be important for divergence. Species that differ in monogenic or oligogenic traits that affect ecological performance and/or mate choice show remarkably localized genomic differentiation. By contrast, differentiation among species that have diverged in polygenic traits is genomically widespread and much higher overall, consistent with the evolution of effective and stable genome-wide barriers to gene flow. Thus, we conclude that simple trait architectures are not always as conducive to speciation with gene flow as previously suggested, whereas polygenic architectures can promote rapid and stable speciation in sympatry.

Non-trophic Functional Ecology of Vertebrate Teeth: A Review.

Teeth are critical to the functional ecology of vertebrate trophic abilities, but are also used for a diversity of other non-trophic tasks. Teeth can play a substantial role in how animals move, manipulate their environment, positively interact with conspecifics, antagonistically interact with other organisms, and sense the environment. We review these non-trophic functions in an attempt to place the utility of human and all other vertebrate dentitions in a more diverse framework that emphasizes an expanded view of the functional importance and ecological diversity of teeth. In light of the extensive understanding of the developmental genetics, trophic functions, and evolutionary history of teeth, comparative studies of vertebrate dentitions will continue to provide unique insights into multi-functionality, many-to-one mapping, and the evolution of novel abilities.

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.

A Genomic Cluster Containing Novel and Conserved Genes is Associated with Cichlid Fish Dental Developmental Convergence.

The two toothed jaws of cichlid fishes provide textbook examples of convergent evolution. Tooth phenotypes such as enlarged molar-like teeth used to process hard-shelled mollusks have evolved numerous times independently during cichlid diversification. Although the ecological benefit of molar-like teeth to crush prey is known, it is unclear whether the same molecular mechanisms underlie these convergent traits. To identify genes involved in the evolution and development of enlarged cichlid teeth, we performed RNA-seq on the serially homologous-toothed oral and pharyngeal jaws as well as the fourth toothless gill arch of Astatoreochromis alluaudi. We identified 27 genes that are highly upregulated on both tooth-bearing jaws compared with the toothless gill arch. Most of these genes have never been reported to play a role in tooth formation. Two of these genes (unk, rpfA) are not found in other vertebrate genomes but are present in all cichlid genomes. They also cluster genomically with two other highly expressed tooth genes (odam, scpp5) that exhibit conserved expression during vertebrate odontogenesis. Unk and rpfA were confirmed via in situ hybridization to be expressed in developing teeth of Astatotilapia burtoni. We then examined expression of the cluster's four genes in six evolutionarily independent and phylogenetically disparate cichlid species pairs each with a large- and a small-toothed species. Odam and unk commonly and scpp5 and rpfA always showed higher expression in larger toothed cichlid jaws. Convergent trophic adaptations across cichlid diversity are associated with the repeated developmental deployment of this genomic cluster containing conserved and novel cichlid-specific genes.

Phenotypic Plasticity in Vertebrate Dentitions.

Vertebrates interact directly with food items through their dentition, and these interactions with trophic resources could often feedback to influence tooth structure. Although dentitions are often considered to be a fixed phenotype, there is the potential for environmentally induced phenotypic plasticity in teeth to extensively influence their diversity. Here, we review the literature concerning phenotypic plasticity of vertebrate teeth. Even though only a few taxonomically disparate studies have focused on phenotypic plasticity in teeth, there are a number of ways teeth can change their size, shape, or patterns of replacement as a response to the environment. Elucidating the underlying physiological, developmental, and genetic mechanisms that generate phenotypic plasticity can clarify its potential role in the evolution of dental phenotypes.

Grand Challenges in Comparative Tooth Biology.

Teeth are a model system for integrating developmental genomics, functional morphology, and evolution. We are at the cusp of being able to address many open issues in comparative tooth biology and we outline several of these newly tractable and exciting research directions. Like never before, technological advances and methodological approaches are allowing us to investigate the developmental machinery of vertebrates and discover both conserved and excitingly novel mechanisms of diversification. Additionally, studies of the great diversity of soft tissues, replacement teeth, and non-trophic functions of teeth are providing new insights into dental diversity. Finally, we highlight several emerging model groups of organisms that are at the forefront of increasing our appreciation of the mechanisms underlying tooth diversification.

Conservation and novelty in the microRNA genomic landscape of hyperdiverse cichlid fishes.

MicroRNAs (miRNAs) play crucial roles in the post-transcriptional control of messenger RNA (mRNA). These miRNA-mRNA regulatory networks are present in nearly all organisms and contribute to development, phenotypic divergence, and speciation. To examine the miRNA landscape of cichlid fishes, one of the most species-rich families of vertebrates, we profiled the expression of both miRNA and mRNA in a diverse set of cichlid lineages. Among these, we found that conserved miRNAs differ from recently arisen miRNAs (i.e. lineage specific) in average expression levels, number of target sites, sequence variability, and physical clustering patterns in the genome. Furthermore, conserved miRNA target sites tend to be enriched at the 5' end of protein-coding gene 3' UTRs. Consistent with the presumed regulatory role of miRNAs, we detected more negative correlations between the expression of miRNA-mRNA functional pairs than in random pairings. Finally, we provide evidence that novel miRNA targets sites are enriched in genes involved in protein synthesis pathways. Our results show how conserved and evolutionarily novel miRNAs differ in their contribution to the genomic landscape and highlight their particular evolutionary roles in the adaptive diversification of cichlids.

Pleiotropic jaw morphology links the evolution of mechanical modularity and functional feeding convergence in Lake Malawi cichlids.

Complexity in how mechanistic variation translates into ecological novelty could be critical to organismal diversification. For instance, when multiple distinct morphologies can generate the same mechanical or functional phenotype, this could mitigate trade-offs and/or provide alternative ways to meet the same ecological challenge. To investigate how this type of complexity shapes diversity in a classic adaptive radiation, we tested several evolutionary consequences of the anterior jaw four-bar linkage for Lake Malawi cichlid trophic diversification. Using a novel phylogenetic framework, we demonstrated that different mechanical outputs of the same four jaw elements are evolutionarily associated with both jaw protrusion distance and jaw protrusion angle. However, these two functional aspects of jaw protrusion have evolved independently. Additionally, although four-bar morphology showed little evidence for attraction to optima, there was substantial evidence of adaptive peaks for emergent four-bar linkage mechanics and jaw protrusion abilities among Malawi feeding guilds. Finally, we highlighted a clear case of two cichlid species that have -independently evolved to graze algae in less than 2 Myr and have converged on similar jaw protrusion abilities as well as four-bar linkage mechanics, but have evolved these similarities via non-convergent four-bar morphologies.

Agouti-related peptide 2 facilitates convergent evolution of stripe patterns across cichlid fish radiations.

The color patterns of African cichlid fishes provide notable examples of phenotypic convergence. Across the more than 1200 East African rift lake species, melanic horizontal stripes have evolved numerous times. We discovered that regulatory changes of the gene agouti-related peptide 2 (agrp2) act as molecular switches controlling this evolutionarily labile phenotype. Reduced agrp2 expression is convergently associated with the presence of stripe patterns across species flocks. However, cis-regulatory mutations are not predictive of stripes across radiations, suggesting independent regulatory mechanisms. Genetic mapping confirms the link between the agrp2 locus and stripe patterns. The crucial role of agrp2 is further supported by a CRISPR-Cas9 knockout that reconstitutes stripes in a nonstriped cichlid. Thus, we unveil how a single gene affects the convergent evolution of a complex color pattern.

Evolutionary divergence of 3' UTRs in cichlid fishes.

BACKGROUND: Post-transcriptional regulation is crucial for the control of eukaryotic gene expression and might contribute to adaptive divergence. The three prime untranslated regions (3' UTRs), that are located downstream of protein-coding sequences, play important roles in post-transcriptional regulation. These regions contain functional elements that influence the fate of mRNAs and could be exceptionally important in groups such as rapidly evolving cichlid fishes. RESULTS: To examine cichlid 3' UTR evolution, we 1) identified gene features in nine teleost genomes and 2) performed comparative analyses to assess evolutionary variation in length, functional motifs, and evolutionary rates of 3' UTRs. In all nine teleost genomes, we found a smaller proportion of repetitive elements in 3' UTRs than in the whole genome. We found that the 3' UTRs in cichlids tend to be longer than those in non-cichlids, and this was associated, on average, with one more miRNA target per gene in cichlids. Moreover, we provided evidence that 3' UTRs on average have evolved faster in cichlids than in non-cichlids. Finally, analyses of gene function suggested that both the top 5% longest and 5% most rapidly evolving 3' UTRs in cichlids tended to be involved in ribosome-associated pathways and translation. CONCLUSIONS: Our results reveal novel patterns of evolution in the 3' UTRs of teleosts in general and cichlids in particular. The data suggest that 3' UTRs might serve as important meta-regulators, regulators of other mechanisms governing post-transcriptional regulation, especially in groups like cichlids that have undergone extremely fast rates of phenotypic diversification and speciation.

The Integrated Genomic Architecture and Evolution of Dental Divergence in East African Cichlid Fishes (Haplochromis chilotes x H. nyererei).

The independent evolution of the two toothed jaws of cichlid fishes is thought to have promoted their unparalleled ecological divergence and species richness. However, dental divergence in cichlids could exhibit substantial genetic covariance and this could dictate how traits like tooth numbers evolve in different African Lakes and on their two jaws. To test this hypothesis, we used a hybrid mapping cross of two trophically divergent Lake Victoria species (Haplochromis chilotes × Haplochromis nyererei) to examine genomic regions associated with cichlid tooth diversity. Surprisingly, a similar genomic region was found to be associated with oral jaw tooth numbers in cichlids from both Lake Malawi and Lake Victoria. Likewise, this same genomic location was associated with variation in pharyngeal jaw tooth numbers. Similar relationships between tooth numbers on the two jaws in both our Victoria hybrid population and across the phylogenetic diversity of Malawi cichlids additionally suggests that tooth numbers on the two jaws of haplochromine cichlids might generally coevolve owing to shared genetic underpinnings. Integrated, rather than independent, genomic architectures could be key to the incomparable evolutionary divergence and convergence in cichlid tooth numbers.

Mechanical Transgressive Segregation and the Rapid Origin of Trophic Novelty.

Hybrid phenotypes are often intermediate between those of parental species. However, hybridization can generate novel phenotypes when traits are complex. For instance, even when the morphologies of individual musculo-skeletal components do not segregate outside the parental range in hybrid offspring, complex functional systems can exhibit emergent phenotypes whose mechanics exceed the parental values. To determine if transgression in mechanics could facilitate divergence during an adaptive radiation, we examined three functional systems in the trophic apparatus of Lake Malawi cichlid fishes. We conducted a simulation study of hybridization between species pairs whose morphology for three functional systems was empirically measured, to determine how the evolutionary divergence of parental species influences the frequency that hybridization could produce mechanics that transgress the parental range. Our simulations suggest that the complex mechanical systems of the cichlid trophic apparatus commonly exhibit greater transgression between more recently diverged cichlid species. Because (1) all three mechanical systems produce hybrids with transgressive mechanics in Lake Malawi cichlids, (2) hybridization is common, and (3) single hybrid crosses often recapitulate a substantial diversity of mechanics, we conclude that mechanical transgressive segregation could play an important role in the rapid accumulation of phenotypic variation in adaptive radiations.

Biting into the Genome to Phenome Map: Developmental Genetic Modularity of Cichlid Fish Dentitions.

Within vertebrates, teleost fishes provide a rich evolutionary context for studying the mechanisms of dental divergence because of the numerous axes along which their teeth have diverged phenotypically and presumably developmentally. Using both a review of teleost in situ hybridization and de novo transcriptome sequencing in a cichlid fish, we examined whether 341 gene homologs thought to play a role in developing mice teeth are expressed in the tooth-bearing jaws of teleosts. The similarities and putative differences in gene expression documented between the two most commonly used models, zebrafish and cichlids, highlight what can be learned from using a greater diversity of teleost model systems in studies of tooth development. Both types of gene expression analysis also provide substantial evidence for conservation of tooth gene expression from teleosts to mammals as well as between initial and replacement teeth. Additionally, we found that the cichlid oral and pharyngeal jaws share expression for a large percentage of genes that influence tooth development. Our transcriptome analyses also suggest sub-functionalization between gene paralogs expressed in teeth and paralogs expressed in other structures is likely a common pattern across teleost diversity. Teleost dentitions will continue to provide a potent system in which to examine the importance of both gene duplication as well as the conservation of gene expression for phenotypic diversification.

Adaptive simplification and the evolution of gecko locomotion: morphological and biomechanical consequences of losing adhesion.

Innovations permit the diversification of lineages, but they may also impose functional constraints on behaviors such as locomotion. Thus, it is not surprising that secondary simplification of novel locomotory traits has occurred several times among vertebrates and could potentially lead to exceptional divergence when constraints are relaxed. For example, the gecko adhesive system is a remarkable innovation that permits locomotion on surfaces unavailable to other animals, but has been lost or simplified in species that have reverted to a terrestrial lifestyle. We examined the functional and morphological consequences of this adaptive simplification in the Pachydactylus radiation of geckos, which exhibits multiple unambiguous losses or bouts of simplification of the adhesive system. We found that the rates of morphological and 3D locomotor kinematic evolution are elevated in those species that have simplified or lost adhesive capabilities. This finding suggests that the constraints associated with adhesion have been circumvented, permitting these species to either run faster or burrow. The association between a terrestrial lifestyle and the loss/reduction of adhesion suggests a direct link between morphology, biomechanics, and ecology.

Influence of substrate orientation on feeding kinematics and performance of algae-grazing Lake Malawi cichlid fishes.

Lake Malawi cichlids have been studied extensively in an effort to elucidate the mechanisms underlying their adaptive radiation. Both habitat partitioning and trophic specialization have been suggested to be critical ecological axes underlying the exceptional diversification of these fishes, but the mechanisms facilitating this divergence are often unclear. For instance, in the rock-dwelling mbuna of Lake Malawi, coexistence is likely tightly linked to how and where species feed on the algae coating all the surfaces of the rocky reefs they exclusively inhabit. Yet, although mbuna species often preferentially graze from particular substrate orientations, we understand very little about how substrate orientation influences feeding kinematics or feeding rates in any group of organisms. Therefore, for three species of mbuna, we quantified feeding kinematics and inferred the rates that algae could be ingested on substrates that mimicked the top, side and bottom of the algae-covered boulders these species utilize in Lake Malawi. A number of differences in feeding kinematics were found among species, and several of the kinematic variables were found to differ even within species when the fish grazed from different surface orientations. However, despite their preferences for particular microhabitats, we found no evidence for clear trade-offs in the rates that the three species were inferred to be able to obtain algae from different substrate orientations. Nevertheless, our results indicate microhabitat divergence linked to differences in feeding kinematics could have played a role in the origin and maintenance of the vast diversity of co-occurring Lake Malawi mbuna species.

Continental monophyly of cichlid fishes and the phylogenetic position of Heterochromis multidens.

The incredibly species-rich cichlid fish faunas of both the Neotropics and Africa are generally thought to be reciprocally monophyletic. However, the phylogenetic affinity of the African cichlid Heterochromis multidens is ambiguous, and this distinct lineage could make African cichlids paraphyletic. In past studies, Heterochromis has been variously suggested to be one of the earliest diverging lineages within either the Neotropical or the African cichlid radiations, and it has even been hypothesized to be the sister lineage to a clade containing all Neotropical and African cichlids. We examined the phylogenetic relationships among a representative sample of cichlids with a dataset of 29 nuclear loci to assess the support for the different hypotheses of the phylogenetic position of Heterochromis. Although individual gene trees in some instances supported alternative relationships, a majority of gene trees, integration of genes into species trees, and hypothesis testing of putative topologies all supported Heterochromis as belonging to the clade of African cichlids.

Explosive diversification following a benthic to pelagic shift in freshwater fishes.

BACKGROUND: Interspecific divergence along a benthic to pelagic habitat axis is ubiquitous in freshwater fishes inhabiting lentic environments. In this study, we examined the influence of this habitat axis on the macroevolution of a diverse, lotic radiation using mtDNA and nDNA phylogenies for eastern North America's most species-rich freshwater fish clade, the open posterior myodome (OPM) cyprinids. We used ancestral state reconstruction to identify the earliest benthic to pelagic transition in this group and generated fossil-calibrated estimates of when this shift occurred. This transition could have represented evolution into a novel adaptive zone, and therefore, we tested for a period of accelerated lineage accumulation after this historical habitat shift. RESULTS: Ancestral state reconstructions inferred a similar and concordant region of our mtDNA and nDNA based gene trees as representing the shift from benthic to pelagic habitats in the OPM clade. Two independent tests conducted on each gene tree suggested an increased diversification rate after this inferred habitat transition. Furthermore, lineage through time analyses indicated rapid early cladogenesis in the clade arising after the benthic to pelagic shift. CONCLUSIONS: A burst of diversification followed the earliest benthic to pelagic transition during the radiation of OPM cyprinids in eastern North America. As such, the benthic/pelagic habitat axis has likely influenced the generation of biodiversity across disparate freshwater ecosystems.