Phenotypic plasticity and modularity allow for the production of novel mosaic phenotypes in ants
© Londe et al. 2015
Received: 23 August 2015
Accepted: 12 November 2015
Published: 1 December 2015
The origin of discrete novelties remains unclear. Some authors suggest that qualitative phenotypic changes may result from the reorganization of preexisting phenotypic traits during development (i.e., developmental recombination) following genetic or environmental changes. Because ants combine high modularity with extreme phenotypic plasticity (queen and worker castes), their diversified castes could have evolved by developmental recombination. We performed a quantitative morphometric study to investigate the developmental origins of novel phenotypes in the ant Mystrium rogeri, which occasionally produces anomalous ‘intercastes.’ Our analysis compared the variation of six morphological modules with body size using a large sample of intercastes.
We confirmed that intercastes are conspicuous mosaics that recombine queen and worker modules. In addition, we found that many other individuals traditionally classified as workers or queens also exhibit some level of mosaicism. The six modules had distinct profiles of variation suggesting that each module responds differentially to factors that control body size and polyphenism. Mosaicism appears to result from each module responding differently yet in an ordered and predictable manner to intermediate levels of inducing factors that control polyphenism. The order of module response determines which mosaic combinations are produced.
Because the frequency of mosaics and their canalization around a particular phenotype may evolve by selection on standing genetic variation that affects the plastic response (i.e., genetic accommodation), developmental recombination is likely to play an important role in the evolution of novel castes in ants. Indeed, we found that most mosaics have queen-like head and gaster but a worker-like thorax congruent with the morphology of ergatoid queens and soldiers, respectively. Ergatoid queens of M. oberthueri, a sister species of M. rogeri, could have evolved from intercastes produced ancestrally through such a process.
KeywordsQueen Intercaste Developmental recombination Caste evolution
Evolution by reorganization of preexisting elements. Different subunits are rearranged, and can be deleted, duplicated, and moved in various ways. This occurs not only at molecular level  but also at higher phenotypic levels. Typical examples of such process are heterochrony (changes in the timing of a developmental process) and heterotopy (changes in spatial location of a developmental process) . Such evolutionary changes are likely to result from selection of mutations in cis-regulatory elements within particular gene regulatory networks 
Reorganization of ancestral phenotypic traits in a particular individual, before genetic accommodation has fixed the phenotype in the population. The process has also been termed ‘chimeric’ or ‘somatic’ recombination
Selection on standing genetic variation that molds the plastic response of a phenotypic trait. This occurs when the developmental-genetic system is sensitized, because genetic variation becomes exposed as phenotypic variation when the organism encounters a different environment
Phenotype recombining within single individual traits that are normally found in distinct individuals
Some authors suggest that evolutionary novelty may result from the recombination of ancestral phenotypes following changes in the timing or the location of preexisting developmental processes [1, 5, 9–13]. This mechanism is referred to as developmental recombination ([9, 10, 14–16]; Table 1). For example, the recurrent evolution of limnetic vs. benthic forms in stickleback (Gasterosteus aculeatus) could have resulted from the altered expression of adaptive traits evolved in ancestral populations, emerging when these traits were expressed alternatively at different times during the life cycle [9, 17]. It has been proposed that evolution by developmental recombination may occur following three steps : (1) Initially, a population consists of individuals that respond differently to some environmental and genomic inputs, because of genetic variation that either encodes distinct fixed traits, epistatic relationships, or distinct sensitivity levels to environmental factors (plastic response). (2) Some individuals are affected by a new genetic or environmental input that causes a reorganization of their phenotype because some phenotypic subunits differ in their responsiveness to the new input (developmental recombination). (3) If the resulting change in phenotype has a positive effect on fitness, it will be favored through selection of genetic factors involved in the plastic response. Selection therefore can lead to further changes in the regulation of development (frequency, timing, and circumstances of the new response) and the characteristics of new traits expressed as differences in morphology, behavior, physiology, etc.  (See notion of “genetic accommodation,” [9, 18–22]; Table 1). This model of evolution by developmental recombination is based on the association of two widespread characteristics of living organisms: plasticity and modularity.
Phenotypic plasticity is the ability of a single genotype to produce alternative phenotypes in response to different environmental factors [9, 23, 24]. It results from the fact that the environment not only filters phenotypic variation but may also induce it [9, 13]. Some environmental cues can regulate gene expression and determine which structures will develop. When there is selection for the genetic ability to respond to environmental cues, phenotypic plasticity becomes adaptive . Phenotypic plasticity can range from a subtle adjustment in growth rate to complete polyphenism involving the production of discrete alternative phenotypes . The plastic response to environmental factors can be continuous (reaction norms) or discrete (polyphenisms) [27, 28]. Classical examples of phenotypic plasticity include the adjustment of biomass allocation to leaf tissue in plants in response to light intensity , the development of a more robust head capsule, and a stronger bite in grasshoppers that have experienced hard food during ontogeny [30, 31], and the development of morphological castes in social hymenoptera . Phenotypic plasticity may play an important role in evolution because responsiveness to the environment affects the strength of selection under alternative environmental circumstances and therefore the intensity of drift and the accumulation of standing genetic variation . Such standing variation is available to be acted upon by genetic accommodation and the canalization of an initially environmentally induced trait can evolve rapidly in response to selection . Experimental approaches confirm that some adaptations derive from characters that were environmentally induced in ancestral populations [10, 11, 35, 36]. Because phenotypic plasticity involves alternative phenotypes, it may also facilitate combinatorial evolution by increasing the number of distinct versions of each module, and thus the number of building blocks available for developmental recombination (e.g., [9, 37, 38]).
Modularity is a universal property of organisms resulting from the branching nature of development . Modules are units, the subparts of which are strongly integrated by numerous genetic, developmental, or functional interactions , but are more weakly connected to other modules . Modularity at the anatomical scale is underpinned by the modularity of developmental processes including interactions among proteins and gene regulation [41–44]. Indeed, regulatory gene interactions form a network with multilevel modular subcircuits. Although most proteins regulating development have pleiotropic effects, the timing of their expression is specifically regulated in different spatial domains by cis-regulatory elements within particular gene regulatory networks (GRNs) . Genetic input affecting particular cis-regulatory modules can thus affect some phenotypic subparts independently of others. As a result of this modular architecture, some mutations can cause gains, losses, or redeployments in other contexts of a function or a morphological structure controlled downstream of the mutations, without affecting other components of the phenotype . A deep level of anatomical modularity is governed by Hox genes that control body plan morphogenesis along the anterior–posterior axis. For instance, the control of forelimbs and hind-limbs by different phases of Hox gene expression regulated by separate cis-regulatory enhancer elements has allowed birds and bats to evolve wings as well as legs . At the phenotypic level, because the subparts that form one module experience the same selective pressure, they co-evolve semi-independently from subparts of other modules [9, 41]. A consequence of this independence is that modules such as arthropod imaginal disks can develop numerous times along the body plan or be deleted independently of other traits [2, 48, 49]. For instance, the Ultrabithorax gene prevents wing development from the imaginal disk of the third thoracic segment by repressing cis-regulatory modules at multiple nodes of the wing patterning GRN . This property allows for the reorganization of different phenotypic building blocks, i.e., evolution by developmental recombination .
Ants combine an extreme polyphenism (queen and worker castes) with a high modularity typical of arthropods. They are a highly diversified clade with various morphological castes [51–53 for an overview of this diversity]. The ancestral and most common caste system features winged queens that perform aerial dispersal and reproduction, and wingless workers that are involved in brood care, foraging, defense, cleaning, and construction of the nest. Accordingly, queens have a fully functional reproductive system with a spermatheca, a large articulated thorax with flight sclerites, large wing muscles and wings, and three ocelli. In contrast, workers have an atrophied reproductive system, a reduced fused thorax lacking wings and flight muscles, and generally no ocelli [51, 52]. In many species, workers may vary in size but do not differ in shape (minor vs. major workers). However, many species have evolved new castes, such as permanently wingless queens called ergatoid queens [54–59] and soldiers that differ from workers in the relative growth of some body parts .
Typically, queen and worker castes are determined mainly by the environment [61, 62]. Environmental factors (e.g., food quality and quantity, temperature and queen pheromones) affect hormonal secretions (mainly juvenile hormone) that induce physiological and cellular responses, ultimately resulting in a developmental switch toward either caste [61, 63]. Switches are determined by several factors such as hormone titers and timing of tissue sensitivity to hormones, and result in caste-specific patterns of gene expression possibly mediated by DNA methylation [64–66]. Regarding wing polyphenism, Abouheif and Wray  have shown that wingless worker and soldier phenotypes are produced through interruptions in the gene regulatory networks yielding wings in queens in the imaginal wing disks of larvae. Interruptions occur at different points in the network depending on species. Caste determination in some species also involves various degrees of genetic control .
Molet et al.  proposed that some new castes in ants may have evolved via developmental recombination of existing castes. Indeed, wingless queens and soldiers appear to be mosaics (Table 1) of winged queen and worker morphologies [60, 69]. This hypothesis implies that environmental or genetic inputs can generate new phenotypes, a known phenomenon in ants. In addition to the discrete queen-worker dichotomy, some authors have described individuals that are neither queens nor workers, called intercastes [54, 70–75]. Intercastes are rare, anomalous adults with various morphologies, visually ranging from almost similar to winged queens to almost similar to workers. Although generally not winged, intercastes often have simplified flight sclerites and on occasion wing stubs. They also may have a spermatheca and one to several ocelli. Thus, intercastes seem to be mosaic phenotypes recombining queen and worker traits. Importantly, they survive as adults, although their behaviors are unstudied. They can be produced when unusual genetic and environmental inputs exceed the buffering capacities of development. This could disturb signaling pathways upstream of some caste-specific GRNs, leading to anomalous gene expression during the ontogeny of the corresponding morphological structures. As a result, modules of the developing larvae do not consistently follow queen and worker pathways. More specifically, departure from normal developmental processes could be the consequence of changes (gains, losses, or modifications) in linkages within GRNs caused by the evolution of cis-regulatory elements (CREs). This includes the co-option of new transcription factor inputs by mutations in existing CREs [76, 77], the co-option of transposable elements as new CREs , the loss of transcriptional inputs in existing CREs , and the remodeling of CREs .
Intercastes have been described in about 20 species  and are likely to be taxonomically widespread. These mosaic individuals probably either go unnoticed due to their rarity or are discarded by researchers because of their abnormal features. Because some intercastes look morphologically similar to ergatoid queens, we suggest that they represent an early step in the evolution of ergatoid queens, before the selection of genetic factors involved in the induction of their phenotypes has fixed a particular phenotype (i.e., genetic accommodation). The emergence of a new caste from environmentally induced anomalies followed by genetic accommodation has also been proposed for the evolution of super soldiers in Pheidole . Accordingly, studying intercastes, and more generally, developmental mechanisms allowing for the production of mosaic phenotypes, will contribute to our understanding of caste evolution.
The intuitive concept of mosaicism has allowed for the description of intercastes based on striking, discrete traits such as the presence or absence of wings and ocelli and of a broad or narrow thorax [73, 74, 81, 82]. However, no quantitative measure of mosaicism has been performed, and consequently phenotypes with less obvious mosaicism have likely been ignored. This means intercastes as described in the literature probably only represent a fraction of the existing range of mosaicism. Indeed, a continuous range of mosaic phenotypes, ranging from worker-like to queen-like, probably exists. Intercastes following the classical definition may only be highly striking cases of mosaic phenotypes (i.e., clearly intermediate between workers and queens), and less distinctive individuals at the extreme of this continuum (i.e., more worker-like or more queen-like) may remain undetected by researchers. Therefore, we propose a new procedure based on morphometric data to quantify the degree of mosaicism and precisely describe the range of combinations among queen and worker modules. We test whether individuals initially identified as intercastes based on discrete characters are effectively mosaics for quantitative morphometric traits, and whether additional mosaic individuals have previously been overlooked.
We used the ant Mystrium rogeri Forel, 1899 (Amblyoponinae) as it is an ideal model to study the production of intercastes and the evolution of new castes. This species erratically produces intercastes [83, 84], and the genus Mystrium includes both species with winged queens (M. rogeri and M. camillae) and species with ergatoid queens (M. oberthueri and M. voeltzkowi) . Accordingly, evolutionary transitions from winged queens to ergatoid queens have occurred in this genus in the past, and may be ongoing in extant M. rogeri populations.
Colony collection and choice of samples
Number of queens, intercastes, and workers from each colony in our sample
2D shape analysis
We used geometric morphometrics to extract shape information from the four modules analyzed in 2D. Head, thorax, and legs were separated and laid flat. Photographs were taken under a stereomicroscope and each module was photographed twice to take into account small deviations caused by optical biases such as lighting variations or position relative to camera lens . We then digitalized the photographs using TPSDig2 software . Measurement error was computed using Procrustes ANOVA  based on the two photograph sessions. Polyphenism induces large shape differences between castes, so identifying numerous homologous anatomical landmarks for shape analysis is not possible. Therefore, we studied outlines using the sliding semi-landmarks method . For head shape, we defined one set of semi-landmarks along the edges of the head capsule and another set along the edges of the clypeus in dorsal view. For thorax shape, we digitalized one set of semi-landmarks along the edges in dorsal view. The thoracic configuration was then split into three partial configurations: pronotum, mesonotum, and propodeum (Fig. 2c). Semi-landmarks were aligned using the minimum Procrustes distance criterion [89, 90]. For each individual and module, we calculated an average configuration based on the two photographs. We focused our analysis on symmetrical variance. In this way, we reflected lateral semi-landmark configurations across the symmetry axes using the “object symmetry” procedure . The shape component of each configuration was extracted using the Procrustes superimposition method [85, 87, 92, 93]. All modules were first rescaled with their own unitary centroid size in order to extract shape information. Rescaled configurations were then superimposed and rotated around their centroid so as to minimize the sum of squared distances between corresponding semi-landmarks. All analyses were performed in R 3.0.1 using routines from the geomorph library  and function from Claude .
The second step was to generate a queen-worker axis for each module. This axis was computed using queen and worker affiliations as defined by the PCA. The data from the two modules measured in 1D could be used directly to generate two queen-worker axes. The data from the four modules measured in 2D required a reduction of dimensionality: we used Between Group Analyses (BGA) [95, 96] that maximize between-groups variance, and projected the results on the first component to generate four additional queen-worker axes.
The final step was to project our full dataset (the three complete colonies plus the additional samples) onto these axes. Then, in order to standardize the range of variation of each module from 0 for the most worker-like morphology to 1 for the most queen-like morphology, we shifted and rescaled the projected data. Each queen-likeness value V was, therefore, transformed as V′ = (V − Vmin)/(Vmax − Vmin).
Two alternative methods to compute the queen-worker axes were rejected. The first axis of a principal component analysis is not suitable because it does not maximize variance between groups. The first component of a linear or quadratic discriminant analysis not only maximizes the discrimination between groups, but also dramatically distorts intra-group variances by transforming relations among initial variables . This effect could have biased the among-groups measure of mosaicism because our mosaicism estimator is the variance among queen-likeness indexes of each module within individuals (see “Mosaicism” section below).
Profiles of variation
The hypothesis that mosaicism originates from differential responses to the caste-determining factors among modules implies that each module shows a specific pattern of variation along the worker-queen axis that we call a ‘profile of variation.’ Profiles of variation were obtained by plotting the queen-likeness of each module against a reference variable. We chose to use individual body length (sum of head and thorax length extracted from 2D configurations) as a reference variable because it is both methodologically independent from queen-likeness and illustrates queen-worker polyphenism well because queens are larger than workers . Gaster length was not included in this reference variable because the gaster is an articulated structure that can be contracted or extended considerably depending on nutrition, reproduction, and preservation. Profiles of variation may be considered allometric relationships between an index of shape and body size.
In order to test whether profiles of variation differed among modules, we fitted a parametric model with queen-likeness data for each module using the R package ‘grofit.’ Given this model, we then extracted the individual body length for which queen-likeness of this module equals 0.5. We called this individual body length the ‘transition point.’ Differences in transition points among modules were tested using pairwise comparisons of mean transition points computed from 1000 bootstrap samplings .
Mosaicism is the degree to which a phenotype recombines modules that normally occur in alternative phenotypes. For instance, an individual with a queen-like head and a strongly worker-like pronotum has a higher level of mosaicism than another individual with a queen-like head and a slightly worker-like pronotum. We computed mosaicism as the standard error of queen-likeness indexes for the six modules, i.e., the dispersion of queen-likenesses among modules. We compared the average mosaicism in queens, workers, and intercastes using Kruskal–Wallis tests. In order to visualize the link between mosaicism and the different profiles of variation among modules, mosaicism was also plotted against individual body length. Slope significance for queens, intercastes, and workers was tested against zero using a bootstrap procedure (1000 samples). In order to verify that mosaicism patterns were not artifacts caused by the queen-likeness computation method, we re-computed mosaicism following inter-individual permutation of queen-likeness for each module. We then tested whether mosaicism patterns disappeared; if so, we concluded there was no methodological artifact.
A preliminary analysis showed that shape of head, pronotum, mesonotum, and propodeum presented significant allometric relationships with centroid sizes (bootstrap procedure as advised by : P < 0.01), except in queens where a significant allometric effect was only found for pronotum (P = 0.049) (we used standardized major axis regression in order to take into account errors in both queen-likeness indexes and size measurements, ). These allometric relationships mean that a proportion of the variance observed in our shape data (and consequently in queen-likeness) is explained by the size of the modules. We did not correct queen-likeness for these allometric effects because this would have affected the profiles of variation, which are themselves allometric relationships between the queen-likeness of each module and individual body size.
Procrustes ANOVA showed that measurement errors were almost 10 times smaller than inter-individual variation (11.4 % for head shape, 0.4 % for pronotum shape, 2.1 % for mesonotum shape, 4.5 % for propodeum shape, 10.1 % for leg length, and 2.2 % for gaster width), confirming that our measurement methods were appropriate and that estimates of shape variances were not compromised by measurement error. Although significant, colony effect ranged between 2.2 and 7.5 times smaller than the effect of caste on the queen-likeness of different modules.
Profiles of variation
Individual mosaicism was quantified as the standard error of queen-likeness for the six modules. Mosaicism was higher for intercastes than for queens and workers (respectively, Chi2 = 40.74, P < 10−9; Chi2 = 36.47, P < 10−8), and higher for workers than for queens (Chi2 = 16.98, P < 10−4) (Fig. 6b).
Plotting mosaicism against individual body length showed that mosaicism is continuous between workers and intercastes (Fig. 6b). Large workers were more mosaic than small workers (slope: 0.08; bootstrap P < 0.001), and as mosaic as intercastes. Mosaicism also increased with body length in intercastes (slope: 0.11; bootstrap P = 0.049) but not in queens (slope: −0.007; bootstrap P = 0.32). Inter-individual permutations of queen-likeness indexes erased this general pattern of mosaicism, proving that it reflected a biological reality and was not a methodological artifact.
The pattern of mosaicism against body length resulted from the variance among the profiles of variation of the different modules. Mosaicism was low for low body length, where all modules had a low queen-likeness. As body length increased, the parametric fits diverged among modules and mosaicism increased. Finally, at high body length, all modules had high queen-likeness and mosaicism became low again (Fig. 6).
Although they are rare (1.6 % of adult females), we obtained a reasonable number of intercastes by collecting 60 colonies of M. rogeri. Our geometric morphometric study of 29 queens, 37 intercastes, and 124 workers showed that the six modules (head, pronotum, mesonotum propodeum, legs, and gaster) grow differently as body length increases; i.e., they have distinct profiles of variation. This is in agreement with the hypothesis that the production of mosaics results from the differential responses of modules to the same inducing factors. Our results suggest that, for intermediate levels of inducing factors, some modules develop as in queens, while others develop as in workers, resulting in the production of mosaic phenotypes.
In our analyses, queen-likeness of each module was plotted against individual body length (Fig. 5). However, the significant correlation between queen-likeness and body length probably does not reflect a causal relationship, but rather a common response of both modules and body length to others factors that control polyphenic development. Body length is also a component of caste polyphenism, queens being larger than workers in many social insect species  including M. rogeri . We used body length as a reference to compare the profiles of variation among several developmental modules, but the same divergence in scaling relationships among modules is expected with any other reference variable as long as this variable is also correlated with factors determining caste polyphenism. Indeed, the scaling relationships linking determining factors and body length necessarily affect the profiles of variation of the modules in the same way regardless of the reference used and thus do not qualitatively affect their divergence. The robustness of our results was confirmed using average queen-likeness as reference instead of body length (data not shown). Therefore, the divergence in scaling relationships among modules suggests that they respond differently to the factors driving caste polyphenism.
In ants, the factors driving development toward a worker or queen phenotype are primarily environmental [61, 63]. Yet, increasing evidence indicates that caste polyphenism can be controlled by genetic, genomic, and maternal factors to various extents [68, 101]; depending on the species, polyphenisms may range from fully environmentally induced to fully genetically determined. In the first case, profiles of variation as established in our procedure would be interpreted as reaction norms [27, 28] with an unknown rescaling function in the abscissa depending on the link between environmental factors and body length. In the second case, divergences in profiles of variation would represent pleiotropic effects among modules in response to various genetic inputs. Our experimental design cannot distinguish genetic components from environmental components. Both might be responsible for the observed effects. The results, therefore, suggest that mosaicism is generated by intermediate levels of determining factors, regardless of their genetic or environmental nature, inducing differential responses among modules.
Mosaicism increases gradually between workers and intercastes and abruptly falls in queens (Fig. 6b). This is in accordance with the progressive divergence among profiles of variation from low to intermediate body lengths, and their rapid convergence from intermediate to high body lengths (Fig. 6a). This suggests that the absence of mosaic phenotypes intermediate between intercastes and queens is not due to an incomplete sample, from which queen-like intercastes would be missing by chance. Rather, it results from the late and sharp increase of queen-likeness for the mesonotum (Fig. 6a), which dramatically reduces the window of sensitivity and leads to very queen-like intercastes. An obvious explanation for the late increase of queen-likeness for the mesonotum is that a queen-like mesonotum is an adaptation to flight that is energetically costly to produce . Selection could allow small variations in thorax morphologies for workers and intermediate phenotypes, but canalize queen morphology by allowing the expression of an expensive, queen-like mesonotum only in individuals expressing queen-like modules, i.e., in perfect queen phenotypes. Without a markedly dimorphic mesonotum, mosaicism would probably decrease more continuously between intercastes and queens, following a reverse U-shape (i.e., increase then decrease) gradually connecting workers to queens. The gradual increase in mosaicism from low to intermediate body length shows that intercastes as described in the literature are only particularly noticeable mosaic individuals with striking qualitative characters among a continuous range of mosaics. Our quantitative approach leads to a better understanding of the developmental mechanisms underlying queen/worker polyphenism.
Using the profiles of variation of the different modules, we can predict which mosaic phenotypes can or cannot be produced in M. rogeri. Among modules, queen-likeness increased at different rates with body length; as a result, the transition points (body length at which queen-likeness exceeds 0.5) differed among modules (Fig. 7). As body size increased, legs were the first module to pass the transition point, followed successively by gaster and head, then propodeum, and finally mesonotum and pronotum. This sequential switch from worker-like to queen-like modules restrains the range of possible combinations. For instance, an individual with a queen-like pronotum is unlikely to have worker-like legs, gaster, or head. Okada et al.  found a similar sequential limitation in possible phenotypes in intercastes of the ant Temnothorax nylanderi. This sequence is likely shared among ants, but further investigation in other species is required.
Previous works suggested that plasticity and modularity allow for combinatorial evolution by reorganizing phenotypes through developmental recombination [9–11, 13, 69]. In social hymenoptera, it is known that the reorganization of phenotypes (i.e., production of intercastes) can be induced by artificially manipulating the environment or physiology of larvae [102–105]. However, for evolution by developmental recombination to occur, two steps are necessary. First, new combinations of characters must occur spontaneously in nature and be exposed to selection. The present study shows that this is the case in M. rogeri. Second, the production of new combinations must be heritable. Importantly, although polyphenisms and reaction norms are plastic responses to environmental factors (i.e., nonheritable factors), their features are genetically controlled, so they can evolve under selection [106, 107]. In addition, phenotypes that are environmentally determined can become genetically determined following adjustments in the sensitivity of development , an evolutionary process known as genetic accommodation [9, 19–22]. Since any mosaic phenotype is visible to natural selection through its direct fitness when fertile or its contribution to colony fitness when sterile, its production could increase as a result of selection on additive genetic factors that change the shape of the profiles of variation for some modules and expand the window of sensitivity allowing induction of this phenotype. At the molecular level, evolution by developmental recombination can rely on the selection of cis-regulatory mutations (i.e., co-option of external transcription factors by mutations in existing CREs, co-option of transposable elements as new CREs, loss of transcriptional factors inputs in existing CREs, and remodeling of CREs) that modify timing and threshold responses in caste-specific signalization within a particular GRN but not in others. This results in a change in pattern of differential gene expression among caste-specific modules. In cases where polyphenism is controlled by the social environment, the evolution of the frequencies of intercaste production and the range of mosaics produced could also occur via the selection of genetic determinants of social behaviors, i.e., brood care and food supply .
Yang and Abouheif  described asymmetric male–female mosaics (gynandromorphs) as monsters with no evolutionary significance, but one possible outcome of symmetric queen-worker mosaics is the evolution of novel castes in ants. For instance, we found that most mosaics have a queen-like head and gaster but worker-like thorax. These phenotypes are congruent with the morphology of ergatoid queens and soldiers, two castes that have repetitively evolved across ants and are suspected to have evolved by phenotypic recombination . Indeed, soldiers recombine a worker thorax (winglessness) with a queen head (defense or seed milling) or a queen gaster (trophic eggs). Similarly, ergatoid queens recombine a worker thorax (winglessness) with a queen gaster (reproductive organs). Therefore, ergatoid queens of M. oberthueri [110, 111], a sister species of M. rogeri, could have evolved from intercastes produced ancestrally through a process of genetic accommodation [9, 19–22].
Overall, our results show that the production of reorganized phenotypes can occur as a consequence of modularity and developmental plasticity, due to differential plastic responses among modules. This provides a parsimonious explanation for the propensity of ants to evolve new ergatoid queen and soldier castes because most mosaics phenotypes are congruent with the morphology of these castes. This scenario still needs to be refined by comparing the regulatory gene networks underlying development of both intercastes and novel castes (e.g., ), and by studying the behavior of developmental anomalies and quantifying their contribution to colony fitness. More generally, our work underlines the need to take into account developmental plasticity in modern evolutionary thought because it determines which phenotypes can or cannot be produced and thus significantly affects the evolutionary potential of populations.
between group analyses
gene regulatory networks
principal component analysis
SL participated in the collection of ant colonies, carried out the data acquisition, performed morphometric analyses, and wrote the manuscript. TM participated in the design of the study and in the writing of the manuscript. RC brought technical expertise to data acquisition and morphometric analyses. VD assisted in morphometric analyses and interpretation of the results. BF participated in the collection of ant colonies and brought expertise in the field of ant evolution. MM participated in the design of the study, in the writing of the manuscript, and supervised the project. All authors read and approved the final manuscript.
This work was funded by Agence Nationale de la Recherche grant ANTEVO ANR-12-JSV7-0003-01 to M. Molet and T. Monnin. We thank Claudie Doums for fruitful discussions and Romain Péronnet for help with ant rearing. Field work for this study was partially supported by the National Science Foundation under Grant Nos. DEB-0072713 and DEB-0842395 and could not have been completed without the gracious support of the Malagasy people and the Arthropod Inventory Team (Balsama Rajemison, Jean Claude Rakotonirina, Jean-Jacques Rafanomezantsoa, Chrislain Ranaivo, Hanitriniana Rasoazanamavo, Nicole Rasoamanana, Manoa Razafimamonjy, and Clavier Randrianandrasana).
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Moczek AP. On the origins of novelty in development and evolution. BioEssays. 2008;30:432–47.View ArticlePubMedGoogle Scholar
- Lewis EB. Genes and developmental pathways. Am Zool. 1963;3:33–56.View ArticleGoogle Scholar
- Gehring WJ, Hiromi Y. Homeotic genes and the homeobox. Annu Rev Genet. 1986;20:147–73.View ArticlePubMedGoogle Scholar
- Shubin N, Tabin C, Carroll S. Deep homology and the origins of evolutionary novelty. Nature. 2009;457:818–23.View ArticlePubMedGoogle Scholar
- Rajakumar R, San Mauro D, Dijkstra MB, Huang MH, Wheeler DE, Hiou-Tim F, Khila A, Cournoyea M, Abouheif E. Ancestral developmental potential facilitates parallel evolution in ants. Science. 2012;335:79–82.View ArticlePubMedGoogle Scholar
- Jacob F. Evolution and tinkering. Science. 1977;196:1161–6.View ArticlePubMedGoogle Scholar
- Wake DB, Roth G: The linkage between ontogeny and phylogeny in the evolution of complex systems. NewYork: Wiley, 1989.
- Duboule D, Wilkins AS. The evolution of ‘bricolage’. Trends Genet. 1998;14:54–9.View ArticlePubMedGoogle Scholar
- West-Eberhard MJ: Developmental plasticity and evolution. Oxford University Press; 2003.
- West-Eberhard MJ. Developmental plasticity and the origin of species differences. Proc Natl Acad Sci. 2005;102(suppl 1):6543–9.PubMed CentralView ArticlePubMedGoogle Scholar
- West-Eberhard MJ. Phenotypic accommodation: adaptive innovation due to developmental plasticity. J Exp Zoolog B Mol Dev Evol. 2005;304B:610–8.View ArticleGoogle Scholar
- Prud’homme B, Minervino C, Hocine M, Cande JD, Aouane A, Dufour HD, Kassner VA, Gompel N. Body plan innovation in treehoppers through the evolution of an extra wing-like appendage. Nature. 2011;473:83–6.View ArticlePubMedGoogle Scholar
- Abouheif E, Favé M-J, Ibarrarán-Viniegra AS, Lesoway MP, Rafiqi AM, Rajakumar R. Eco-Evo-Devo: the time has come. In: Landry CR, Aubin-Horth N (eds) Ecological Genomics. Springer Netherlands; 2014:107–125. [Advances in Experimental Medicine and Biology, vol. 781].
- Davidson EH. Lineage-specific gene expression and the regulative capacities of the sea urchin embryo: a proposed mechanism. Development. 1989;105:421–45.PubMedGoogle Scholar
- Ray TS. Metamorphosis in the Araceae. Am J Bot. 1990;77:1599–609.View ArticleGoogle Scholar
- Raff RA, Kaufman TC: Embryos, genes, and evolution: the developmental—genetic basis of evolutionary change. 1991(XXVI):395.
- Foster SA, Baker JA. Evolution in parallel: new insights from a classic system. Trends Ecol Evol. 2004;19:456–9.View ArticlePubMedGoogle Scholar
- Waddington CH. Genetic assimilation of an acquired character. Evolution. 1953;7:118.View ArticleGoogle Scholar
- Suzuki Y, Nijhout HF. Evolution of a polyphenism by genetic accommodation. Science. 2006;311:650–2.View ArticlePubMedGoogle Scholar
- Moczek AP. Developmental capacitance, genetic accommodation, and adaptive evolution. Evol Dev. 2007;9:299–305.View ArticlePubMedGoogle Scholar
- Nijhout HF, Suzuki Y. Environment and genetic accommodation. Biol Theory. 2008;3:204–12.View ArticleGoogle Scholar
- Suzuki Y, Nijhout HF. Genetic basis of adaptive evolution of a polyphenism by genetic accommodation. J Evol Biol. 2008;21:57–66.PubMedGoogle Scholar
- Debat V, David P. Mapping phenotypes: canalization, plasticity and developmental stability. Trends Ecol Evol. 2001;16:555–61.View ArticleGoogle Scholar
- Whitman D, Agrawal A. What is phenotypic plasticity and why is it important. Phenotypic Plastic Insects. 2009; 1–63.
- Gilbert SF, Epel D: Ecological developmental biology. Edition 1. Sunderland: Sinauer Associates; 2009.
- Nijhout HF. Development and evolution of adaptive polyphenisms. Evol Dev. 2003;5:9–18.View ArticlePubMedGoogle Scholar
- Woltereck R. Weitere experimentelle untersuchungen über Artänderung, speziell über das Wesen quantitativer Artunterschiede bei Daphniden. Z Für Indukt Abstamm-Vererbungslehre. 1909;9:146–146.Google Scholar
- Stearns SC. The evolutionary significance of phenotypic plasticity. Bioscience. 1989;39:436–45.View ArticleGoogle Scholar
- Sultan SE, Bazzaz FA. Phenotypic plasticity in Polygonum persicaria. I. Diversity and uniformity in genotypic norms of reaction to light. Evolution. 1009;1993:47.Google Scholar
- Thompson DB. Consumption rates and the evolution of diet-induced plasticity in the head morphology of Melanoplus femurrubrum (Orthoptera: Acrididae). Oecologia. 1992;89:204–13.View ArticleGoogle Scholar
- Thompson. Genotype–environment interaction and the ontogeny of diet-induced phenotypic plasticity in size and shape of Melanoplus femurrubrum (Orthoptera: Acrididae). J Evol Biol. 1999;12:38–48.View ArticleGoogle Scholar
- Oster GF, Wilson EO. Caste and ecology in the social insects. Princeton University Press; 1978.
- Snell-Rood EC, Van Dyken JD, Cruickshank T, Wade MJ, Moczek AP. Toward a population genetic framework of developmental evolution: the costs, limits, and consequences of phenotypic plasticity. BioEssays. 2010;32:71–81.PubMed CentralView ArticlePubMedGoogle Scholar
- Scoville AG, Pfrender ME. Phenotypic plasticity facilitates recurrent rapid adaptation to introduced predators. Proc Natl Acad Sci. 2010;107:4260–3.PubMed CentralView ArticlePubMedGoogle Scholar
- Shaw KA, Scotti ML, Foster SA. Ancestral plasticity and the evolutionary diversification of courtship behaviour in threespine sticklebacks. Anim Behav. 2007;73:415–22.View ArticleGoogle Scholar
- Wund MA, Baker JA, Clancy B, Golub JL, Foster SA. A test of the “Flexible Stem” model of evolution: ancestral plasticity, genetic accommodation, and morphological divergence in the threespine stickleback radiation. Am Nat. 2008;172:449–62.View ArticlePubMedGoogle Scholar
- Ghalambor CK, McKAY JK, Carroll SP, Reznick DN. Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments. Funct Ecol. 2007;21:394–407.View ArticleGoogle Scholar
- Moczek AP, Sultan S, Foster S, Ledón-Rettig C, Dworkin I, Nijhout HF, Abouheif E, Pfennig DW: The role of developmental plasticity in evolutionary innovation. Proc R Soc Lond B Biol Sci 2011; rspb20110971.
- Wagner GP. Adaptation and the modular design of organisms. In: Morán F, Moreno A, Merelo JJ, Chacón P, editors. Advances in Artificial Life.. Springer Berlin Heidelberg; 1995. p. 315–328. [Lecture Notes in Computer Science, vol. 929].
- Breuker CJ, Debat V, Klingenberg CP. Functional evo-devo. Trends Ecol Evol. 2006;21:488–92.View ArticlePubMedGoogle Scholar
- Klingenberg CP. Morphological integration and developmental modularity. Annu Rev Ecol Evol Syst. 2008;39:115–32.View ArticleGoogle Scholar
- Klingenberg CP, Leamy LJ, Cheverud JM. Integration and modularity of quantitative trait locus effects on geometric shape in the mouse mandible. Genetics. 2004;166:1909–21.PubMed CentralView ArticlePubMedGoogle Scholar
- Wagner GP, Pavlicev M, Cheverud JM. The road to modularity. Nat Rev Genet. 2007;8:921–31.View ArticlePubMedGoogle Scholar
- Burgio G, Baylac M, Heyer E, Montagutelli X. Exploration of the genetic organization of morphological modularity on the mouse mandible using a set of interspecific recombinant congenic strains between C57BL/6 and mice of the Mus spretus species. G3 Genes Genomes Genetics. 2012;2:1257–68.PubMed CentralPubMedGoogle Scholar
- Carroll SB. Evo-devo and an expanding evolutionary synthesis: a genetic theory of morphological evolution. Cell. 2008;134:25–36.View ArticlePubMedGoogle Scholar
- Erwin DH, Davidson EH. The evolution of hierarchical gene regulatory networks. Nat Rev Genet. 2009;10:141–8.View ArticlePubMedGoogle Scholar
- Shubin N, Tabin C, Carroll S. Fossils, genes and the evolution of animal limbs. Nature. 1997;388:639–48.View ArticlePubMedGoogle Scholar
- Morata G, Garcia-Bellido A. Developmental analysis of some mutants of the bithorax system of Drosophila. Wilhelm Rouxs Arch Dev Biol. 1976;179:125–43.View ArticleGoogle Scholar
- Morata G, Lawrence PA. Development of the eye-antenna imaginal disc of Drosophila. Dev Biol. 1979;70:355–71.View ArticlePubMedGoogle Scholar
- Galant R, Carroll SB. Evolution of a transcriptional repression domain in an insect Hox protein. Nature. 2002;415:910–3.View ArticlePubMedGoogle Scholar
- Wilson EO: The insect societies. 1971. p. 548.
- Hölldobler B, Wilson EO. The Ants. Cambridge: Harvard University Press; 1990.View ArticleGoogle Scholar
- Fisher BL. AntWeb-Ants of the world. 2012. http://www.antweb.org/ (2015). Accessed 05 Nov 2015.
- Heinze J, Tsuji K. Ant reproductive strategies. Res Popul Ecol. 1995;37:135–49.View ArticleGoogle Scholar
- Heinze J, Keller L. Alternative reproductive strategies: a queen perspective in ants. Trends Ecol Evol. 2000;15:508–12.View ArticlePubMedGoogle Scholar
- Peeters C, Ito F. Colony dispersal and the evolution of queen morphology in social hymenoptera. Annu Rev Entomol. 2001;46:601–30.View ArticlePubMedGoogle Scholar
- Peeters C, Molet M: Chapter 9: colonial reproduction and life histories. In: Lori Lach, Catherine L. Parr, Kirsti L, editors. Abbott, Ant ecology. Oxford University Press; 2010.
- Peeters C. Convergent evolution of wingless reproductives across all subfamilies of ants, and sporadic loss of winged queens (Hymenoptera: Formicidae). Myrmecol News. 2012;16:75–91.Google Scholar
- Cronin AL, Molet M, Doums C, Monnin T, Peeters C. Recurrent evolution of dependent colony foundation across eusocial insects. Annu Rev Entomol. 2013;58:37–55.View ArticlePubMedGoogle Scholar
- Molet M, Maicher V, Peeters C. Bigger helpers in the Ant Cataglyphis bombycina: increased worker polymorphism or novel soldier caste? PLoS One. 2014;9:e84929.PubMed CentralView ArticlePubMedGoogle Scholar
- Wheeler DE. Developmental and physiological determinants of caste in social hymenoptera: evolutionary implications. Am Nat. 1986;128:13–34.View ArticleGoogle Scholar
- Wheeler DE. The developmental basis of worker caste polymorphism in ants. Am Nat. 1991;138:1218–38.View ArticleGoogle Scholar
- Hartfelder K, Engels W. Current topics in developmental biology. Elsevier; 1998.
- Evans JD, Wheeler DE. Gene expression and the evolution of insect polyphenisms. BioEssays. 2001;23:62–8.View ArticlePubMedGoogle Scholar
- Moczek AP, Snell-Rood EC. The basis of bee-ing different: the role of gene silencing in plasticity. Evol Dev. 2008;10:511–3.View ArticlePubMedGoogle Scholar
- Elango N, Hunt BG, Goodisman MAD, Yi SV. DNA methylation is widespread and associated with differential gene expression in castes of the honeybee, Apis mellifera. Proc Natl Acad Sci. 2009;106:11206–11.PubMed CentralView ArticlePubMedGoogle Scholar
- Abouheif E, Wray GA. Evolution of the gene network underlying wing polyphenism in ants. Science. 2002;297:249–52.View ArticlePubMedGoogle Scholar
- Schwander T, Lo N, Beekman M, Oldroyd BP, Keller L. Nature versus nurture in social insect caste differentiation. Trends Ecol Evol. 2010;25:275–82.View ArticlePubMedGoogle Scholar
- Molet M, Wheeler DE, Peeters C. Evolution of novel mosaic castes in ants: modularity, phenotypic plasticity, and colonial buffering. Am Nat. 2012;180:328–41.View ArticlePubMedGoogle Scholar
- Wheeler WM. Worker ants with vestiges of wings. Order Trutees Am Mus Nat Hist 1905.
- Wheeler WM, Weber NA. Mosaics and Other Anomalies among Ants. 1937.
- Brian MV. Studies of caste differentiation in «Myrmica Rubra» L. Insectes Soc. 1955;2:1–34.View ArticleGoogle Scholar
- Plateaux L: Sur le polymorphisme social de la fourrni Leptothorax nylanderi (Förster). I. Morphologie et biologie comparee des castes. 1970.
- Peeters CP. Ergatoid queens and intercastes in ants: two distinct adult forms which look morphologically intermediate between workers and winged queens. Insectes Soc. 1991;38:1–15.View ArticleGoogle Scholar
- Heinze J. Intercastes, intermorphs, and ergatoid queens: who is who in ant reproduction? Insectes Soc. 1998;45:113–24.View ArticleGoogle Scholar
- Wang X, Chamberlin HM. Multiple regulatory changes contribute to the evolution of the Caenorhabditis lin-48 ovo gene. Genes Dev. 2002;16:2345–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Gompel N, Prud’homme B, Wittkopp PJ, Kassner VA, Carroll SB. Chance caught on the wing: cis-regulatory evolution and the origin of pigment patterns in Drosophila. Nature. 2005;433:481–7.View ArticlePubMedGoogle Scholar
- Bejerano G, Lowe CB, Ahituv N, King B, Siepel A, Salama SR, Rubin EM, James Kent W, Haussler D. A distal enhancer and an ultraconserved exon are derived from a novel retroposon. Nature. 2006;441:87–90.View ArticlePubMedGoogle Scholar
- Hinman VF, Davidson EH. Evolutionary plasticity of developmental gene regulatory network architecture. Proc Natl Acad Sci. 2007;104:19404–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Zinzen RP, Cande J, Ronshaugen M, Papatsenko D, Levine M. Evolution of the Ventral Midline in Insect Embryos. Dev Cell. 2006;11:895–902.View ArticlePubMedGoogle Scholar
- Düssmann O, Peeters C, Hölldobler B. Morphology and reproductive behaviour of intercastes in the ponerine ant Pachycondyla obscuricornis. Insectes Soc. 1996;43:421–5.View ArticleGoogle Scholar
- Okada Y, Plateaux L, Peeters C. Morphological variability of intercastes in the ant Temnothorax nylanderi: pattern of trait expression and modularity. Insectes Soc. 2013;60:319–28.View ArticleGoogle Scholar
- Molet M, Fisher BL, Ito F, Peeters C. Shift from independent to dependent colony foundation and evolution of “multi-purpose” ergatoid queens in Mystrium ants (subfamily Amblyoponinae). Biol J Linn Soc. 2009;98:198–207.View ArticleGoogle Scholar
- Yoshimura M, Fisher BL. A revision of the ant genus Mystrium in the Malagasy region with description of six new species and remarks on Amblyopone and Stigmatomma (Hymenoptera, Formicidae, Amblyoponinae). ZooKeys. 2014;1–99.
- Claude J. Morphometrics with R. Springer; 2008.
- Rohlf J: tpsDig2. 2009.
- Klingenberg CP, McIntyre GS. Geometric morphometrics of developmental instability: analyzing patterns of fluctuating asymmetry with procrustes methods. Evolution. 1998;52:1363–75.View ArticleGoogle Scholar
- Bookstein FL. Landmark methods for forms without landmarks: morphometrics of group differences in outline shape. Med Image Anal. 1997;1:225–43.View ArticlePubMedGoogle Scholar
- Andresen PR, Bookstein FL, Couradsen K, Ersboll BK, Marsh JL, Kreiborg S. Surface-bounded growth modeling applied to human mandibles. IEEE Trans Med Imaging. 2000;19:1053–63.View ArticlePubMedGoogle Scholar
- Sheets HD, Kim K, Mitchell CE. A combined landmark and outline-based approach to ontogenetic shape change in the Ordovician trilobite Triarthrus becki. In: Elewa DAMT, editor. Morphometrics. Berlin Heidelberg: Springer; 2004. p. 67–82.Google Scholar
- Klingenberg CP, Barluenga M, Meyer A. Shape analysis of symmetric structures: quantifying variation among individuals and asymmetry. Evolution. 2002;56:1909–20.View ArticlePubMedGoogle Scholar
- Goodall C. Procrustes methods in the statistical analysis of shape. J R Stat Soc Ser B Methodol. 1991;53:285–339.Google Scholar
- Smith DR, Crespi BJ, Bookstein FL. Fluctuating asymmetry in the honey bee, Apis mellifera: effects of ploidy and hybridization. J Evol Biol. 1997;10:551–74.View ArticleGoogle Scholar
- Adams DC, Otárola-Castillo E. Geomorph: an r package for the collection and analysis of geometric morphometric shape data. Methods Ecol Evol. 2013;4:393–9.View ArticleGoogle Scholar
- Culhane AC, Perrière G, Considine EC, Cotter TG, Higgins DG. Between-group analysis of microarray data. Bioinformatics. 2002;18:1600–8.View ArticlePubMedGoogle Scholar
- Doledec S, Chessel D. Rythmes saisonniers et composantes stationnelles en milieu aquatique. I: Description d’un plan d’observation complet par projection de variables. Acta Oecologica Oecologia Gen. 1987;8:403–26.Google Scholar
- Fukunaga K. Introduction to statistical pattern recognition. Academic Press; 1990.
- Efron B, Tibshirani R. Bootstrap methods for standard errors, confidence intervals, and other measures of statistical accuracy. Stat Sci. 1986;1:54–75.View ArticleGoogle Scholar
- Klingenberg CP: Multivariate allometry. In Marcus LF, Corti M, Loy A, Naylor GJP, Slice DE, editors. Advances in Morphometrics. USA: Springer; 1996. p. 23–49. [NATO ASI Series, vol. 284].
- Fairbairn DJ. Allometry for sexual size dimorphism: pattern and process in the coevolution of body size in males and females. Annu Rev Ecol Syst. 1997;28:659–87.View ArticleGoogle Scholar
- Schwander T, Humbert J-Y, Brent CS, Cahan SH, Chapuis L, Renai E, Keller L. Maternal effect on female caste determination in a social insect. Curr Biol. 2008;18:265–9.View ArticlePubMedGoogle Scholar
- Brian MV. Caste determination in a Myrmicine ant. Experientia. 1951;7:182–3.View ArticlePubMedGoogle Scholar
- Windig JJ. Reaction norms and the genetic basis of phenotypic plasticity in the wing pattern of the butterfly Bicyclus anynana. J Evol Biol. 1994;7:665–95.View ArticleGoogle Scholar
- Murakami T, Ohkawara K, Higashi S. Morphology and developmental plasticity of reproductive females in Myrmecina nipponica (Hymenoptera: Formicidae). Ann Entomol Soc Am. 2002;95:577–82.View ArticleGoogle Scholar
- Linksvayer TA, Kaftanoglu O, Akyol E, Blatch S, Amdam GV, Page RE. Larval and nurse worker control of developmental plasticity and the evolution of honey bee queen–worker dimorphism. J Evol Biol. 2011;24:1939–48.PubMed CentralView ArticlePubMedGoogle Scholar
- Scheiner SM. Genetics and evolution of phenotypic plasticity. Annu Rev Ecol Syst. 1993;24:35–68.View ArticleGoogle Scholar
- Nicoglou A. The evolution of phenotypic plasticity: genealogy of a debate in genetics. Stud Hist Philos Sci Part C Stud Hist Philos Biol Biomed Sci. 2015;50:67–76.View ArticleGoogle Scholar
- Moczek AP. The behavioral ecology of threshold evolution in a polyphenic beetle. Behav Ecol. 2003;14:841–54.View ArticleGoogle Scholar
- Yang AS, Abouheif E. Gynandromorphs as indicators of modularity and evolvability in ants. J Exp Zoolog B Mol Dev Evol. 2011;316B:313–8.View ArticleGoogle Scholar
- Molet M, Peeters C, Fisher BL. Winged queens replaced by reproductives smaller than workers in Mystrium ants. Naturwissenschaften. 2007;94:280–7.View ArticlePubMedGoogle Scholar
- Molet M, Peeters C, Follin I, Fisher BL. Reproductive caste performs intranidal tasks instead of workers in the Ant Mystrium oberthueri. Ethology. 2007;113:721–9.View ArticleGoogle Scholar
- McGinnis W, Krumlauf R. Homeobox genes and axial patterning. Cell. 1992;68:283–302.View ArticlePubMedGoogle Scholar
- Hall BK. Evo-Devo: evolutionary developmental mechanisms. Int J Dev Biol. 2003;47:491–6.PubMedGoogle Scholar
- Maeshiro T, Kimura M. The role of robustness and changeability on the origin and evolution of genetic codes. Proc Natl Acad Sci. 1998;95:5088–93.PubMed CentralView ArticlePubMedGoogle Scholar