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.

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.].

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.

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.

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.

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.

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.

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.

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.

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.

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

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

Pelagic larval duration predicts extinction risk in a freshwater fish clade.

Pelagic larval duration (PLD) can influence evolutionary processes ranging from dispersal to extinction in aquatic organisms. Using estimates of PLD obtained from species of North American darters (Percidae: Etheostomatinae), we demonstrate that this freshwater fish clade exhibits surprising variation in PLD. Comparative analyses provide some evidence that higher stream gradients favour the evolution of shorter PLD. Additionally, similar to patterns in the marine fossil record in which lower PLD is associated with greater extinction probability, we found a reduced PLD in darter lineages was evolutionarily associated with extinction risk. Understanding the causes and consequences of PLD length could lead to better management and conservation of organisms in our increasingly imperiled aquatic environments.

Biomechanical trade-offs bias rates of evolution in the feeding apparatus of fishes.

Morphological diversification does not proceed evenly across the organism. Some body parts tend to evolve at higher rates than others, and these rate biases are often attributed to sexual and natural selection or to genetic constraints. We hypothesized that variation in the rates of morphological evolution among body parts could also be related to the performance consequences of the functional systems that make up the body. Specifically, we tested the widely held expectation that the rate of evolution for a trait is negatively correlated with the strength of biomechanical trade-offs to which it is exposed. We quantified the magnitude of trade-offs acting on the morphological components of three feeding-related functional systems in four radiations of teleost fishes. After accounting for differences in the rates of morphological evolution between radiations, we found that traits that contribute more to performance trade-offs tend to evolve more rapidly, contrary to the prediction. While ecological and genetic factors are known to have strong effects on rates of phenotypic evolution, this study highlights the role of the biomechanical architecture of functional systems in biasing the rates and direction of trait evolution.

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.

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.

Species tree estimation and the historical biogeography of heroine cichlids.

Heroine cichlids are major components of the fish faunas in both Central America and the Caribbean. To examine the evolutionary patterns of how cichlids colonized both of these regions, we reconstructed the phylogenetic relationships among 23 cichlid lineages. We used three phylogenetically novel nuclear markers (Dystropin b, Myomesin1, and Wnt7b) in combination with sequence data from seven other gene regions (Nd2, Rag1, Enc1, Sreb2, Ptr, Plagl2, and Zic1) to elucidate the species tree of these cichlids. The species examined represent major heroine lineages in South America, Central America, and the Greater Antilles. The individual gene trees of these groups were topologically quite discordant. Therefore, we combined the genetic partitions and inferred the species tree using both concatenation and a coalescent-based Bayesian method. The two resulting phylogenetic topologies were largely concordant but differed in two fundamental ways. First, more nodes in the concatenated tree were supported with substantial or 100% Bayesian posterior support than in the coalescent-based tree. Second, there was a minor, but biogeographically critical, topological difference between the concatenated and coalescent-based trees. Nevertheless, both analyses recovered topologies consistent with the Greater Antillean heroines being phylogenetically nested within the largely Central American heroine radiation. This study suggests that reconstructions of cichlid phylogeny and historical biogeography should account for the vagaries of individual gene histories.

Temporal diversification of Central American cichlids.

BACKGROUND: Cichlid fishes are classic examples of adaptive radiation because of their putative tendency to explosively diversify after invading novel environments. To examine whether ecological opportunity increased diversification (speciation minus extinction) early in a species-rich cichlid radiation, we determined if Heroine cichlids experienced a burst of diversification following their invasion of Central America. RESULTS: We first reconstructed the Heroine phylogeny and determined the basal node to use as the root of Central American Heroine diversification. We then examined the influence of incomplete taxon sampling on this group's diversification patterns. First, we added missing species randomly to the phylogeny and assessed deviations from a constant rate of lineage accumulation. Using a range of species numbers, we failed to recover significant deviations from a pure-birth process and found little support for an early burst of diversification. Then, we examined patterns of lineage accumulation as nodes were increasingly truncated. We assumed that as we removed more recently diverged lineages that sampling would become more complete thereby increasing the power to detect deviations from a pure-birth model. However, truncation of nodes provided even less support for an early burst of diversification. CONCLUSIONS: Contrary to expectations, our analyses suggest Heroine cichlids did not undergo a burst of diversification when they invaded from South America. Throughout their history in Central America, Heroine cichlids appear to have diversified at a constant rate.

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.

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.