Developmental integration in a complex morphological structure: how distinct are the modules in the mouse mandible?

The mouse mandible has long served as a model system for studying the development and evolution of complex morphological structures. We used the methods of geometric morphometrics to reassess the hypothesis that the mandible consists of two separate modules: an anterior part bearing the teeth and a posterior part with muscle attachment surfaces and articulating with the skull. The analyses particularly focused on covariation of fluctuating asymmetry, because such covariation is due exclusively to direct interactions between the developmental processes that produce the traits of interest, whereas variation of traits among individuals also reflects other factors. The patterns of fluctuating asymmetry and individual variation were only partly consistent, indicating that developmental processes contribute differentially to variation at different levels. The results were in agreement with the hypothesis that the anterior and posterior parts of the mandible are separate develop-mental modules. Comparison of all alternative partitions of the landmarks into two contiguous subsets confirmed the hypothesis for the location of the boundary between modules but also underscored that the separation between them is not complete. Modularity is therefore manifest as the relative independence of parts within the framework of overall integration of the mandible as a whole-it is a matter of degrees, not all or nothing.     

Drosophila wing modularity revisited through a quantitative genetic approach

To predict the response of complex morphological structures to selection it is necessary to know how the covariation among its different parts is organized. Two key features of covariation are modularity and integration. The Drosophila wing is currently considered a fully integrated structure. Here, we study the patterns of integration of the Drosophila wing and test the hypothesis of the wing being divided into two modules along the proximo‐distal axis, as suggested by developmental, biomechanical, and evolutionary evidence. To achieve these goals we perform a multilevel analysis of covariation combining the techniques of geometric morphometrics and quantitative genetics. Our results indicate that the Drosophila wing is indeed organized into two main modules, the wing base and the wing blade. The patterns of integration and modularity were highly concordant at the phenotypic, genetic, environmental, and developmental levels. Besides, we found that modularity at the developmental level was considerably higher than modularity at other levels, suggesting that in the Drosophila wing direct developmental interactions are major contributors to total phenotypic shape variation. We propose that the precise time at which covariance‐generating developmental processes occur and/or the magnitude of variation that they produce favor proximo‐distal, rather than anterior‐posterior, modularity in the Drosophila wing.     

PLASTICITY, CANALIZATION, AND DEVELOPMENTAL STABILITY OF THE DROSOPHILA WING: JOINT EFFECTS OF MUTATIONS AND DEVELOPMENTAL TEMPERATURE

The phenotypic effects of genetic and environmental manipulations have been rarely investigated simultaneously. In addition to phenotypic plasticity, their effect on the amount and directions of genetic and phenotypic variation is of particular evolutionary importance because these constitute the material for natural selection. Here, we used heterozygous insertional mutations of 16 genes involved in the formation of the Drosophila wing. The flies were raised at two developmental temperatures (18°C and 28°C). Landmark-based geometric morphometrics was used to analyze the variation of the wing size and shape at different hierarchical levels: among genotypes and temperatures; among individuals within group; and fluctuating asymmetry (FA). Our results show that (1) the phenotypic effects of the mutations depend on temperature; (2) reciprocally, most mutations affect wing plasticity; (3) both temperature and mutations modify the levels of FA and of among individuals variation within lines. Remarkably, the patterns of shape FA seem unaffected by temperature whereas those associated with individual variation are systematically altered. By modifying the direction of available phenotypic variation, temperature might thus directly affect the potential for further evolution. It suggests as well that the developmental processes responsible for developmental stability and environmental canalization might be partially distinct.     

Natural genetic variation in fluctuating asymmetry of wing shape in Drosophila melanogaster

Fluctuating asymmetry (FA), defined as random deviation from perfect symmetry, has been used to assay the inability of individuals to buffer their developmental processes from environmental perturbations (i.e., developmental instability). In this study, we aimed to characterize the natural genetic variation in FA of wing shape in Drosophila melanogaster, collected from across the Japanese archipelago. We quantified wing shapes at whole wing and partial wing component levels and evaluated their mean and FA. We also estimated the heritability of the mean and FA of these traits. We found significant natural genetic variation in all the mean wing traits and in FA of one of the partial wing components. Heritability estimates for mean wing shapes were significant in two and four out of five wing traits in males and females, respectively. On the contrary, heritability estimates for FA were low and not significant. This is a novel study of natural genetic variation in FA of wing shape. Our findings suggest that partial wing components behave as distinct units of selection for FA, and local adaptation of the mechanisms to stabilize developmental processes occur in nature.     

Inferring developmental modularity from morphological integration: analysis of individual variation and asymmetry in bumblebee wings

Organisms are built from distinct modules, which are internally coherent but flexible in their relationships among one another. We examined morphological variation within and between two candidate modules: the fore- and hindwings of bumblebees (Hymenoptera: Apidae: Bombus empatiens). We used the techniques of geometric morphometrics (Procrustes superimposition) to analyze the variation of landmark configurations in fore- and hindwings. Regression was used to correct for size-related shape variation (allometry). Principal component analysis revealed patterns of variation that were remarkably similar for individual variation and fluctuating asymmetry (FA). Because covariation of FA among parts must be due to direct transmission of the developmental perturbations causing FA, this agreement of patterns suggests that much of individual variation is also due to direct developmental interactions within each developing wing. Moreover, partial least squares analysis indicated that the patterns of shape covariation between fore- and hindwings were nearly the same as the patterns of within-wing variation. Shape covariation of FA was only found in bees that had been reared under elevated CO(2) concentration but not in bees from the control treatment, suggesting that the mechanisms of developmental interactions between fore- and hindwings are related to gas exchange. We conclude that the fore- and hindwings are developmental modules that maintain internal coherence through direct developmental interactions and are connected to each other only by relatively few links that use the system of interactions within modules.     

Morphological integration and modularity in Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae) hind wings

The morphological integration of the hind wings of the western corn rootworm Diabrotica virgifera virgifera LeConte was investigated to get a better insight of the undergone by this invasive species. Geometric morphometric methods were used to test two modularity hypotheses associated with the wing development and function (hypothesis H1: anterior/posterior or H2: distal/proximal wing parts). Both hypotheses were rejected and the results showed the integrated behavior of the hind wings of D. v. virgifera. The hypothesized modules do not represent separate units of variation, so in a similar fashion as exhibited by the model species Drosophila melanogaster, the hind wings of D. v. virgifera act as a single functional unit. The moderate covariation strength found between anterior and posterior and distal and proximal parts of the hind wing of D. v. virgifera confirms its integrated behavior. We conclude that the wing shape shows internal integration, which could enable flexibility and thus enhance flight maneuverability. This study contributes to the understanding of morphological integration and modularity on a non-model organism. Additionally, these findings lay the groundwork for future flight performance and biogeographical studies on how wing shape and size vary across the endemic and expanded/invaded range in the USA and Europe infested with D. v. virgifera.     

Modularity of a Cambrian ptychoparioid trilobite cranidium

SUMMARY Modularity of the cranidium of Crassifimbra? metalaspis, a Cambrian ptychoparioid trilobite, is investigated using landmark‐based geometric morphometric methods to gain insight into the integration among morphogenetic processes responsible for shaping the head of an ancient arthropod. Of particular interest is the extent to which the structure of phenotypic integration was governed by direct interactions among developmental pathways, because these interactions may generate long‐term constraints on evolutionary innovation. A modified two‐way ANOVA decomposes cranidial shape variation into components representing symmetric variation among individuals and fluctuating asymmetry (FA). The structure of integration of each of these components is inferred from correlated deviations in shape among nine partitions of the cranidium. Significant correlation among partitions in FA indicates direct interactions among their respective developmental pathways. An a priori hypothesis that modularity was determined by functional association among partitions is not well supported by the among‐partition correlation structure for either component of variation. Instead, exploratory analyses reveal that phenotypic integration was strongly influenced by spatially localized morphogenetic controls. Comparison of the structures of the Individuals and FA components of variation reveals that the two share relatively few commonalities: the structure of phenotypic integration was only weakly influenced by direct interactions. The large contribution of parallel variation to phenotypic integration suggests that modularity was unlikely to have imposed a long‐term constraint on evolutionary innovation in these early trilobites.     

Morpho morphometrics: Shared ancestry and selection drive the evolution of wing size and shape in Morpho butterflies

Butterfly wings harbor highly diverse phenotypes and are involved in many functions. Wing size and shape result from interactions between adaptive processes, phylogenetic history, and developmental constraints, which are complex to disentangle. Here, we focus on the genus Morpho (Nymphalidae: Satyrinae, 30 species), which presents a high diversity of sizes, shapes, and color patterns. First, we generate a comprehensive molecular phylogeny of these 30 species. Next, using 911 collection specimens, we quantify the variation of wing size and shape across species, to assess the importance of shared ancestry, microhabitat use, and sexual selection in the evolution of the wings. While accounting for phylogenetic and allometric effects, we detect a significant difference in wing shape but not size among microhabitats. Fore and hindwings covary at the individual and species levels, and the covariation differs among microhabitats. However, the microhabitat structure in covariation disappears when phylogenetic relationships are taken into account. Our results demonstrate that microhabitat has driven wing shape evolution, although it has not strongly affected forewing and hindwing integration. We also found that sexual dimorphism of forewing shape and color pattern are coupled, suggesting a common selective force.     

Lack of response to artificial selection on developmental stability of partial wing shape components in Drosophila melanogaster

Developmental stability, the ability of organisms to buffer their developmental processes against developmental noise is often evaluated with fluctuating asymmetry (FA). Natural genetic variation in FA has been investigated using Drosophila wings as a model system and the recent estimation of the heritability of wing shape FA was as large as 20 %. Because natural genetic variation in wing shape FA was found to localize in a partial component of the wings, heritable variation in specific parts of the wings might be responsible for FA estimation based on the whole wing shape. In this study, we quantified the shape of three partial components of the wings, and estimated the heritability of the wing shape FA based on artificial selections. As a result, FA values for the partial wing shape components did not respond to artificial selections and the heritability scores estimated were very small. These results indicate that natural additive genetic variation in FA of partial wing components was very small compared with that in a complex wing trait.     

Compartments and the control of growth in the Drosophila wing imaginal disc

The mechanisms that control organ growth are among the least known in development. This is particularly the case for the process in which growth is arrested once final size is reached. We have studied this problem in the wing disc of Drosophila, the developmental and growth parameters of which are well known. We have devised a method to generate entire fast-growing Minute(+) (M(+)) discs or compartments in slow developing Minute/+ (M/+) larvae. Under these conditions, a M(+) wing disc gains at least 20 hours of additional development time. Yet it grows to the same size of Minute/+ discs developing in M/+ larvae. We have also generated wing discs in which all the cells in either the anterior (A) or the posterior (P) compartment are transformed from M/+ to M(+). We find that the difference in the cell division rate of their cells is reflected in autonomous differences in the developmental progression of these compartments: each grows at its own rate and manifests autonomous regulation in the expression of the developmental genes wingless and vestigial. In spite of these differences, ;mosaic' discs comprising fast and slow compartments differentiate into adult wings of the correct size and shape. Our results demonstrate that imaginal discs possess an autonomous mechanism with which to arrest growth in anterior and posterior compartments, which behave as independent developmental units. We propose that this mechanism does not act by preventing cell divisions, but by lengthening the division cycle.     

Multidimensional analysis of Drosophila wing variation in Evolution Canyon

Environmental stress has been suggested to be a major evolutionary force, both through inducing strong selection and because of its direct impact on developmental buffering processes that alter the evolvability of organisms. In particular, temperature has attracted much attention because of its importance as an ecological feature and the relative ease with which it can be experimentally manipulated in the lab. Evolution Canyon, Lower Nahal Oren, Israel, is a well studied natural site where ecological parameters are suspected to drive evolutionary differentiation. In this study, using Drosophila melanogaster isofemale lines derived from wild flies collected on both slopes of the canyon, we investigated the effect of developmental temperature upon the different components of phenotypic variation of a complex trait: the wing. Combining geometric and traditional morphometrics, we find only limited evidence for a differentiation among slopes. Investigating simultaneously phenotypic plasticity, genetic variation among isofemale lines, variation among individuals and fluctuating asymmetry, we could not identify a consistent effect of the stressful conditions encountered on the south facing slope. The prevailing structuring effect is that of the experimentally manipulated temperature which clearly influences wing mean size and shape. Variability, in contrast, is not consistently affected by temperature. Finally, we investigated the specific relationship between individual variation and fluctuating asymmetry. Using metric multi-dimensional scaling we show that the related patterns of wing shape variation are not identical, supporting the view that the underlying developmental processes are to a certain extent different.     

Temperature and genotypic effects on life history and fluctuating asymmetry in a field strain of Culex pipiens

Fluctuating asymmetry (FA) has been proposed as a tool to measure levels of stress experienced by populations of organisms during development. To be of value as a bio-marker to highlight conditions at particular sites, it is important that variation in FA is due to environmental (eg pollution) variation and not genetic variation among populations and families, in other words heritability for FA should be very close to zero. A full-sib design was set up in which families of Culex pipiens mosquitoes collected from the field were reared at three different developmental temperatures. The effects of temperature and family on developmental rate, egg to adult survival and four wing morphological measures were assessed. There was both a temperature and a family effect on development rate and survival. Temperature affected all four wing traits, but an influence of family was only evident in two of the wing traits. Two separate measures of FA for each of the wing traits were obtained. The mean estimates of FA were mainly around 1% of the value of the character measured. There was evidence of an increase in FA with increase in temperature stress. Heritability was estimated for the wing traits and wing trait FA's using restricted estimation maximum likelihood. The estimates of heritability for the wing traits were small and, individually, did not differ significantly from zero. There was also no evidence of heritable genetic variation for any of the wing trait FA's. The results are discussed in relation to other studies where FA heritabilities have been estimated and in relation to the use of FA as an indicator of environmental stress.     

Epigenetic interactions and the structure of phenotypic variation in the cranium

Understanding the developmental and genetic basis for evolutionarily significant morphological variation in complex phenotypes such as the mammalian skull is a challenge because of the sheer complexity of the factors involved. We hypothesize that even in this complex system, the expression of phenotypic variation is structured by the interaction of a few key developmental processes. To test this hypothesis, we created a highly variable sample of crania using four mouse mutants and their wild-type controls from similar genetic backgrounds with developmental perturbations to particular cranial regions. Using geometric morphometric methods we compared patterns of size, shape, and integration in the sample within and between the basicranium, neurocranium, and face. The results highlight regular and predictable patterns of covariation among regions of the skull that presumably reflect the epigenetic influences of the genetic perturbations in the sample. Covariation between relative widths of adjoining regions is the most dominant factor, but there are other significant axes of covariation such as the relationship between neurocranial size and basicranial flexion. Although there are other sources of variation related to developmental perturbations not analyzed in this study, the patterns of covariation created by the epigenetic interactions evident in this sample may underlie larger scale evolutionary patterns in mammalian craniofacial form.     

GEOMETRIC MORPHOMETRICS OF DEVELOPMENTAL INSTABILITY: ANALYZING PATTERNS OF FLUCTUATING ASYMMETRY WITH PROCRUSTES METHODS

Although fluctuating asymmetry has become popular as a measure of developmental instability, few studies have examined its developmental basis. We propose an approach to investigate the role of development for morphological asymmetry by means of morphometric methods. Our approach combines geometric morphometrics with the two-way ANOVA customary for conventional analyses of fluctuating asymmetry and can discover localized features of shape variation by examining the patterns of covariance among landmarks. This approach extends the notion of form used in studies of fluctuating asymmetry from collections of distances between morphological landmarks to an explicitly geometric concept of shape characterized by the configuration of landmarks. We demonstrate this approach with a study of asymmetry in the wings of tsetse flies (Glossina palpalis gambiensis). The analysis revealed significant fluctuating and directional asymmetry for shape as well as ample shape variation among individuals and between the offspring of young and old females. The morphological landmarks differed markedly in their degree of variability but multivariate patterns of landmark covariation identified by principal component analysis were generally similar between fluctuating asymmetry (within-individual variability) and variation among individuals. Therefore there is no evidence that special developmental processes control fluctuating asymmetry. We relate some of the morphometric patterns to processes known to be involved in the development of fly wings.     

Multilevel assessment of the Lacertid lizard cranial modularity

Different factors and processes that produce phenotypic variation at the individual, population, or interspecific level can influence or alter the covariance structure among morphological traits. Therefore, studies of the patterns of integration and modularity at multiple levels—static, ontogenetic, and evolutionary, can provide invaluable data on underlying factors and processes that structured morphological variation, directed, or constrained evolutionary changes. Our dataset, consisting of cranium shape data for 14 lizard species from the family Lacertidae, with substantial samples of hatchlings and adults along with their inferred evolutionary relationships, enabled us to assess modularity and morphological integration at all three levels. Five, not mutually exclusive modularity hypotheses of lizard cranium, were tested, and the effects of allometry on intensity and the pattern of integration and modularity were estimated. We used geometric morphometrics to extract symmetric and asymmetric, as well as allometric and nonallometric, components of shape variation. At the static level, firm confirmation of cranial modularity was found for hypotheses which separate anterior and posterior functional compartments of the skull. At the ontogenetic level, two alternative hypotheses (the “anteroposterior” and “neurodermatocranial” hypotheses) of ventral cranial modularity were confirmed. At the evolutionary level, the “neurodermatocranial” hypothesis was confirmed for the ventral cranium, which is in accordance with the pattern observed at the ontogenetic level. The observed pattern of static modularity could be driven by functional demands and can be regarded as adaptive. Ontogenetic modularity and evolutionary modularity show the same developmental origin, indicating conservatism of modularity patterns driven by developmental constraints.     

Developmental plasticity associated with early structural integration and evolutionary patterns: Examples of developmental bias and developmental facilitation in the skeletal system

The relation of developmental plasticity to evolutionary diversification is a key component of evolutionary theory involving developmental bias, but the basis of the relationship varies among traits and among taxa. Here I review some scenarios of how structural integration during early organogenesis could influence this relationship. When condensations are highly integrated and dependent on each other during early organogenesis, both plasticity and evolution are restricted, for example size proportions in molar tooth rows and phalanges within a digit. When similar condensations develop and remain separate (in tracheal cartilages and feather buds), they show high levels of variation and diversity in number but not in shape and size, at least at early stages. When non‐similar structures form separately and then integrate while still undergoing patterning, high levels of plasticity (in number, size, shape; in rib uncinate processes) or new dimensions of ecologically‐significant variation (cusp offset, in mammal teeth) are seen. Although each of these structural integration scenarios is unique, the modulation of evolvability is detectable and informative. Parsing the influence of structural integration at these developmental levels, rather than later‐stage structural correlations or only through genetic covariation, may be necessary to advance understanding of evolvability of the phenotype.     

The "eyespot module" and eyespots as modules: development, evolution, and integration of a complex phenotype

Organisms are inherently modular, yet modules also evolve in response to selection for functional integration or functional specialization of traits. For serially repeated homologous traits, there is a clear expectation that selection on the function of individual traits will reduce the integration between traits and subdivide a single ancestral module. The eyespots on butterfly wings are one example of serially repeated morphological traits that share a common developmental mechanism but are subject to natural and sexual selection for divergent functions. Here, I test two hypotheses about the organization of the eyespot pattern into independent dorsal-ventral and anterior-posterior modules, using a graphical modeling technique to examine patterns of eyespot covariation among and within wing surfaces in the butterfly Bicyclus anynana. Although there is a hierarchical and complex pattern of integration among eyespots, the results show a surprising mismatch between patterns of eyespot integration and the developmental and evolutionary eyespot units identified in previous empirical studies. These results are discussed in light of the relationships between developmental, functional, and evolutionary modules, and they suggest that developmental sources of independent trait variation are often masked by developmental sources of trait integration.     

Morphological integration and serial homology: A case study of the cranium and anterior vertebrae in salamanders

Serial homology or the repetition of equivalent developmental units and their derivatives is a phenomenon encountered in a variety of organisms, with the vertebrate axial skeleton as one of the most notable examples. Serially homologous structures can be viewed as an appropriate model system for studying morphological integration and modularity, due to the strong impact of development on their covariation. Here, we explored the pattern of morphological integration of the cranium and the first three serially homologous structures (atlas, first, and second trunk vertebrae) in salamandrid salamanders, using micro-CT scanning and three-dimensional geometric morphometrics. We explored the integration between structures at static and evolutionary levels. Effects of allometry on patterns of modularity were also taken into account. At the static level (within species), we analyzed inter-individual variation in shape to detect functional modules and intra-individual variation to detect developmental modules. Significant integration (based on inter-individual variation) among all structures was detected and allometry is shown to be an important integrating factor. The pattern of intra-individual, asymmetric variation indicates statistically significant developmental integration between the cranium and the atlas and between the first two trunk vertebrae. At the evolutionary level (among species), the cranium, atlas, and trunk vertebrae separate as different modules. Our results show that morphological integration at the evolutionary level coincides with morphological and functional differentiation of the axial skeleton, allowing the more or less independent evolutionary changes of the cranial skeleton and the vertebral column, regardless of the relatively strong integration at the static level. The observed patterns of morphological integration differ across levels, indicating different impacts of developmental and phylogenetic constraints and functional demands.     

Condition dependence of developmental stability in the sexually dimorphic fly Telostylinus angusticollis (Diptera: Neriidae)

Developmental stability is widely regarded as a condition‐dependent trait, but its relation to genotype and environment, and extent of developmental integration, remain contentious. In Telostylinus angusticollis, the dorsocentral bristles exhibit striking variation in developmental stability, manifested as fluctuating asymmetry (FA) in bristle position (‘positional FA’) and failure to develop some bristles (‘bristle loss’), in natural and laboratory populations. To determine whether this variation reflects condition, I tested for effects of genotype and environment (larval diet quality), and examined covariation with condition‐dependent traits. Positional FA was not affected by genotype or environment. However, positional FA covaried negatively with secondary sexual trait expression in males, and with sexual dimorphism in body shape, but covaried positively with body size in females. Bristle loss reflected both genotype and larval diet. Flies reared on poor‐quality diet exhibited a similar rate of bristle loss as wild flies. Both positional FA and bristle loss were greater in males. These results suggest that the relation between developmental stability and condition is complex and sex dependent.