Building Developmental Integration into Functional Systems: Function-Induced Integration of Mandibular Shape

The mammalian mandible is a developmentally modular but functionally integrated system. Whether morphological integration can evolve to match the optimal pattern of functional integration may depend on the developmental origin of integration, specifically, on the role that direct epigenetic interactions play in shaping integration. These interactions are hypothesized to integrate modules and also to be highly conservative, potentially constraining the evolution of integration. Using the fox squirrel (Sciurus niger) mandible as a model system, we test five a priori developmental hypotheses that predict mandibular integration and we also explore for correlations between shapes of mandibular regions not anticipated by any of the developmental models. To determine whether direct epigenetic interactions are highly conserved in rodents, we examine the correlation structure of fluctuating asymmetry, and compare integration patterns between fox squirrels and prairie deer mice (Peromyscus maniculatus bairdii). In fox squirrels, we find a correlation structure unanticipated by all a priori developmental models: adjacent parts along the proximodistal jaw axis are correlated whereas more distant ones are not. The most notable exception is that the shape of the anterior incisor alveolus is correlated with the shape of the ramus (FA component) or coronoid (symmetric component). Those exceptions differ between species; in prairie deer mice, the molar alveolus is connected to more parts, and the incisor alveolus to fewer, than in fox squirrels. The structure of integration suggests that the mandible cannot be decomposed into parts but rather is a single connected unit, a result consistent with its functional integration. That match between functional and developmental integration may arise, at least in part, from function-induced growth, building developmental integration into the functional system and enabling direct epigenetic interactions to evolve when function does.     

Fluctuating asymmetry and morphometric variation of hand bones

The major aim of this study was to test three hypotheses: 1) more complex traits of the hand are less prone to developmental insults and therefore show lower fluctuating asymmetry (FA) as compared with simple traits; 2) the manifestation of FA correlates with the variability of the trait (i.e., CV); and 3) FA is an organ-wide property, and therefore a concordance exists between the FA measures of different traits in hand bones. Seventy-two bilateral measurements of hand bones, were made from plain-film radiographs of 365 cadavers. A complex trait was considered as the total length of the three phalanges of a finger and their contiguous metacarpals. Simple traits were considered to be the lengths of individual bone that made up the complex trait. The following results were obtained: 1) on the average simple traits, composing the complex trait, show much higher FA than the corresponding complex trait, but this result is expected if there is no correlation (or low correlation) between FA of simple traits within the complex trait, due to random direction of right-left differences; 2) strong and highly significant correlation was observed between FA and CV of studied traits, regardless of sex and age of individuals; and 3) the majority of FA measurements of hand bones showed no correlation. However, correlations between some sets of FA traits were highly significant. They were interpreted, although not specifically tested, as the result of a tight relationship between traits related not only developmentally but also by active performance of the same function. Am J Phys Anthropol 107:125–136, 1998. © 1998 Wiley-Liss, Inc.     

GENETICS OF FLUCTUATING ASYMMETRY: A DEVELOPMENTAL MODEL OF DEVELOPMENTAL INSTABILITY

Although numerous studies have found that fluctuating asymmetry (FA) can have a heritable component, the genetic and developmental basis of FA is poorly understood. We used a developmental model of a trait, according to a diffusion-threshold process, whose parameters are under genetic control. We added a small amount of random variation to the parameter values of this model to simulate developmental noise. As a result of the nonlinearity of the model, different genotypes differed in their sensitivity to developmental noise, even though the noise is completely random and independent of the genotype. The heritable component of FA can thus be understood as genetically modulated expression of variation that is itself entirely nongenetic. The loci responsible for this genetic variation of FA are the same that affect the left/right mean of the trait, showing that genetic variation for FA does not require genes that specifically control FA. Furthermore, the model offers alternative explanations for phenomena widely discussed in the literature on FA, for instance, the correlations between FA and heterozygosity and between FA and trait size. The model underscores the importance of dominance and epistasis, and therefore unites the study of FA with the classical theory of quantitative genetics.     

Morphological variability and developmental instability in subpopulations of the Eurasian badger (Meles meles) in Denmark

Aim Local populations from different geographical regions may differ in the selection regimes to which they are exposed. Differences in environmental factors and population density may affect the relative importance of different selective forces (e.g. natural vs. sexual selection). We suggest a direction of investigation concerned with the developmental instability of morphological traits. The goal is to disclose putative small‐scale geographical differences in the evolutionary forces, which may be hard to detect. Location Craniometrical investigations were carried out on ninety‐eight skulls and teeth of the Eurasian badger (Meles meles) collected during the period 1995–97 from three different populations in Denmark. One of these thrives at low population density, whereas the two others are characterized by high local density. Methods The skulls were investigated for developmental instability (DI) using fluctuating asymmetry (FA) as its estimator. FA was measured on canines, molars, premolar teeth and other skull and mandible traits. For the statistical analyses, we applied nonparametric permutation tests. Results Evidence was found suggesting differentiation among populations in mean degree of FA, and the FA values measured on canines were higher in the high‐density populations. FA of the canines was significantly higher in males than females, in contrast to FA of the other traits. Evidence of a negative relationship between canine size and their FA was found, whereas no significant correlations were found between the molar and premolar teeth measures and their FA. Main conclusions Our results suggest that canines could be under directional selection stemming from intrasexual competition, which may be stronger in high‐density zones. The other teeth investigated seem to be under a stabilizing regime hence their FA is mainly affected by environmental stresses. The negative relationship between canine size and FA found in males suggests the capacity of badgers to respond in an evolutionary way to environmental changes, despite the low genetic variability previously found at the molecular level.     

Correlated development, organism‐wide asymmetry and patterns of asymmetry in two moth species

Understanding the mechanisms that determine the development of a bilaterally symmetrical trait is crucial to the interpretation of patterns of fluctuating asymmetry (FA). Experimental and theoretical studies have indicated that feedback mechanisms both within and between developing traits, may participate in the developmental control of asymmetry. This study provides evidence that naturally occurring patterns of FA are affected by interactions between different traits. We found positive between‐trait correlations in signed FA values for tibia lengths on different legs, but not between wing and tibia FA in two moth species. Further research should investigate if trait functionality is related to this presumed correlated development. An extension of the Rashevsky–Turing model of morphogenesis further showed that correlations between the signed FA values can be generated by feedback mechanisms that regulate growth patterns between traits. We argue that such feedback mechanisms can be expected to be widespread and show that between‐trait correlations in the unsigned FA then become confounded with correlations in the signed FA. In addition, correlated development appeared to invalidate the use of the hypothetical repeatability to translate correlations between the unsigned FA values into correlations in the presumed underlying developmental instability. In conclusion, the presence of an organism‐wide asymmetry, which are most frequently found in morphologically integrated traits, may be even less common than previously thought. This revised version was published online in July 2006 with corrections to the Cover Date.     

Fluctuating asymmetry in fetuses of diabetic rhesus macaques

This study examines dental fluctuating asymmetry (FA) in two samples of fetal rhesus monkeys, one composed of 19 fetuses from diabetic mothers (FDM) and the other of 20 fetuses from nondiabetic mothers. Seventeen measurements were taken on the deciduous dentition of right and left mandibles. The degree of FA was assessed by comparing FDM to fetuses of normal mothers by correlation between right and left sides, and analysis of variation differences between right and left sides. Significant FA was found for three traits based on the correlation between right and left sides and for seven traits by the between-treatment ratio of variance between sides. Distal teeth, both within and outside of a morphologic field, exhibit significantly greater FA than mesial teeth. Our results support the hypothesis that developmental instability is detectable by dental FA.     

A latent variable model of developmental instability in relation to men's sexual behaviour

A single trait's fluctuating asymmetry (FA) is expected to be a poor measure of developmental instability. Hence, studies that examine associations between FA and outcomes expected to covary with developmental instability often have little power in detecting meaningful relationships. One way of increasing the power of detecting relationships between developmental instability and outcomes is through the use of multiple traits' FA. The way multiple traits have typically been used is in trait aggregates. Here, we illustrate another way of examining relationships with developmental instability using multiple traits' FA: through structural equation modelling. Covariances between measures of FA and an outcome variable are interpreted within the context of an explicit model of associations between variables, which is tested for fit and the parameters specified within the model are estimated. We used nine traits' FA as markers of a latent variable of men's developmental instability, which was associated with the number of sexual partners. The results indicate a sizeable correlation between developmental instability and men's sexual history, despite small correlations between individual traits' FA and sexual history.     

QUANTITATIVE GENETIC MODELS FOR DEVELOPMENT,EPIGENETIC SELECTION,AND PHENOTYPIC EVOLUTION

Herein we describe a general multivariate quantitative genetic model that incorporates two potentially important developmental phenomena, maternal effects and epigenetic effects. Maternal and epigenetic effects are defined as partial regression coefficients and phenotypic variances are derived in terms of age-specific genetic and environmental variances. As a starting point, the traditional quantitative genetic model of additive gene effects and random environmental effects is cast in a developmental time framework. From this framework, we first extend a maternal effects model to include multiple developmental ages for the occurrence of maternal effects. An example of maternal effects occurring at multiple developmental ages is prenatal and postnatal maternal effects in mammals. Subsequently, a model of intrinsic and epigenetic effects in the absence of maternal effects is described. It is shown that genetic correlations can arise through epigenetic effects, and in the absence of other developmental effects, epigenetic effects are in general confounded with age-specific intrinsic genetic effects. Finally, the two effects are incorporated into the basic quantitative genetic model. For this more biologically realistic model combining maternal and epigenetic effects, it is shown that the phenotypic regressions of offspring on mother and offspring on father can be used in some cases to estimate simultaneously maternal effects and epigenetic effects.     

Parasitism, developmental plasticity and bilateral asymmetry of wing feathers in alpine swift, Apus melba, nestlings

The hypothesis that developmental instability is a cost of developmental plasticity is explored using the alpine swift ( Apus melba ) as a model organism. In a previous study, experimentally parasitized nestlings showed a reduced wing growth rate in the first half of the rearing period when parasites were abundant (i.e. peak infestation) and an accelerated growth rate (i.e. compensatory growth) in the second half when parasites decreased in number. This suggests that alpine swifts are able to adjust growth rate in relation to variation in parasite loads. Because developmental plasticity may entail fitness costs, the energy required to sustain compensatory growth may be invested at the expense of developmental stability, potentially resulting in larger deviations from symmetry in paired, bilateral traits (i.e. fluctuating asymmetry, FA). This hypothesis predicts higher FA in parasitized than deparasitized nestlings because of compensatory growth, and hence individuals sustaining the highest level of compensatory growth rate should exhibit the highest FA levels. Another non-mutually exclusive hypothesis argues that parasites directly cause FA by diverting energy required by host for maintenance and growth, and predicts that individuals suffering the most from parasitism during peak infestation should exhibit the highest FA levels. The present study shows that wing feathers of experimentally parasitized nestlings were more asymmetrical than those of experimentally deparasitized ones 50 days after hatching. Furthermore, in parasitized individuals FA was negatively correlated with wing growth rate during the period of peak infestation but not during the period of compensatory growth. These findings suggest that developmental homeostasis is more sensitive to parasites than to compensatory growth.     

Developmental quantitative genetic models of evolutionary change

Discussions about evolutionary change in developmental processes or morphological structures are predicated on specific quantitative genetic models whose parameters predict whether evolutionary change can occur, its relative rate and direction, and if correlated change will occur in other related and unrelated structures. The appropriate genetic model should reflect the relevant genetical and developmental biology of the organisms, yet be simple enough in its parameters so that deductions can be made and hypotheses tested. As a consequence, the choice of the most appropriate genetic model for polygenically controlled traits is a complex tissue and the eventual choice of model is often a compromise between completeness of the model and computational expediency. Herein, we discuss several developmental quantitative genetic models for the evolution of development and morphology. The models range from the classical direct effects model to complex epigenetic models. Further, we demonstrate the algebraic equivalency of the Cowley and Atchley epigenetic model and Wagner's developmental mapping model. Finally, we propose a new multivariate model for continuous growth trajectories. The relative efficacy of these various models for understanding evolutionary change in developmental and morphological traits is discussed. © 1994 Wiley-Liss, Inc.     

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.     

Learning the Formation Mechanism of Domain-Level Chromatin States with Epigenomics Data

Epigenetic modifications can extend over long genomic regions to form domain-level chromatin states that play critical roles in gene regulation. The molecular mechanism for the establishment and maintenance of these states is not fully understood and remains challenging to study with existing experimental techniques. Here, we took a data-driven approach and parameterized an information-theoretic model to infer the formation mechanism of domain-level chromatin states from genome-wide epigenetic modification profiles. This model reproduces statistical correlations among histone modifications and identifies well-known states. Importantly, it predicts drastically different mechanisms and kinetic pathways for the formation of euchromatin and heterochromatin. In particular, long, strong enhancer and promoter states grow gradually from short but stable regulatory elements via a multistep process. On the other hand, the formation of heterochromatin states is highly cooperative, and no intermediate states are found along the transition path. This cooperativity can arise from a chromatin looping-mediated spreading of histone methylation mark and supports collapsed, globular three-dimensional conformations rather than regular fibril structures for heterochromatin. We further validated these predictions using changes of epigenetic profiles along cell differentiation. Our study demonstrates that information-theoretic models can go beyond statistical analysis to derive insightful kinetic information that is otherwise difficult to access.     

Evidence for organism-wide asymmetry in five bird species of a fragmented afrotropical forest

Empirical evidence for between-trait correlation in fluctuating asymmetry (FA) at the individual level is generally lacking or contradictory. Yet the assumption of organism-wide asymmetry, estimated by the asymmetry of any given trait, is inherent to most studies that use FA as a measure of developmental stability (DS). A commonly assumed reason for this weak between-trait correlation is the low repeatability of individual, single-trait asymmetry. In this paper we describe high repeatability and significant between-trait correlation in population- and individual-level FA in five afrotropical bird species inhabiting a fragmented cloud forest. Absence of anti-symmetry and of between-trait correlation in signed FA levels permits us to translate the observed patterns into the presumed underlying DS, using the concept of hypothetical repeatability. This correction, which has not been applied before in this context, proved adequate as it yielded correlations comparable to those found at the population level.     

Dental fluctuating asymmetry in the Gullah: tests of hypotheses regarding developmental stability in deciduous vs. permanent and male vs. female teeth

In this investigation, deciduous teeth (canines, c; first molars, m1; second molars, m2) and their permanent successors (canines, C; first premolars, P1; second premolars, P2) were used to test two related hypotheses about fluctuating asymmetry (FA). First, based on the biology of the developing dentition, it was predicted that deciduous teeth would be more developmentally stable and thus exhibit less dimensional FA than their permanent successors. Second, based on sex differences in tooth development, it was predicted that female canines would have greater developmental stability (less FA) than male canines. Bucco-lingual measurements were made on dental casts from a single Gullah population. Using a repeated-measures study design (n = 3 repeated measures), we tested these hypotheses on sample sizes ranging from 63-82 antimeric pairs. Neither hypothesis was supported by our data. In most cases, Gullah deciduous teeth did not exhibit statistically significantly less FA than their permanent successors; indeed, statistically significant differences were found for only 3 of 12 deciduous vs. permanent contrasts, and in two cases, the deciduous tooth had greater FA. Female mandibular canines exhibited statistically significantly greater FA than those of males, while there was no statistically significant sex difference in FA for the maxillary canine. FA in these Gullah samples is high when compared to Archaic and late prehistoric Ohio Valley Native Americans, consistent with historical and archaeological evidence that environmental stress was relatively higher in the Gullah population. We suggest that when environmental stress in a population is high, the impact of differences in tooth formation time spans and developmental buffering upon FA may be minor relative to the effect of developmental noise.     

SOME MODELS FOR DEVELOPMENT,GROWTH, AND MORPHOMETRIC CORRELATION

Genetic and phenotypic correlations between morphometric traits can be a direct consequence of shared developmental history and common systems of growth regulation. Correlation between traits, therefore, need not imply direct functional or adaptive constraints on those traits. Useful models of the developmental origins of correlations will consider mechanisms that can reduce initially high correlation of traits that arise from a single developmental precursor. Several models presented here predict such correlations for different modes of fission of a precursor. Timing of developmental events may also affect correlations and respond to selection on adult traits. The models may apply to development of the tetrapod limb bud, including variance and covariance induced by known developmental mutants.     

Mononuclear phagocyte system in SLE. II. A kinetic model of immune complex handling in systemic lupus erythematosus

Using the principles of reaction kinetics, we constructed a model for the handling of immune complexes and the pathogenesis of SLE immune complex disease. The model incorporates rate constants for complement- and Fc-mediated clearance, parameters for autoantibody, complement and immune complex levels, and scores for clinical disease activity. The model assumes that complement fixation by immune complexes is a prerequisite for complement-mediated clearance and that disease activity results from immune complex deposition. To test the relationships derived, data from 32 lupus patients were analyzed and the predictions were compared with actual findings. The model predicts a low correlation coefficient between disease activity and immune complex levels (found, r = 0.25, p greater than 0.1). The model also predicts a poor correlation between disease activity and impaired Fc-mediated clearance in patients with normal complement levels (found, r = 0.10, p greater than 0.1), but a high correlation coefficient between disease activity and impaired Fc-mediated clearance in patients with hypocomplementemia (found, r = 0.61, p less than 0.001). In patients with normal complement levels, the model predicts a good correlation between anti-DNA antibody and immune complex levels (found, r = 0.71, p less than 0.001), whereas hypocomplementemic patients should have a good correlation between anti-DNA to CH50 ratios and immune complex levels (found, r = 0.73, p less than 0.001). The model predicts that disease activity should correlate better with the product of the anti-DNA to CH50 ratio and the rate constant for Fc-mediated clearance than with any single parameter (found, r = 0.85, p less than 0.0001). These significant correlations, which were predicted by the model, suggest that complement-mediated mechanisms are the first line of host defense against immune complex-induced injury, that the efficiency of complement opsonization plays a central role, and that both abnormal complement- and Fc-receptor function leads to active renal disease in SLE.     

Mixture analysis of asymmetry: modelling directional asymmetry, antisymmetry and heterogeneity in fluctuating asymmetry

The occurrence of different forms of asymmetry complicates the analysis and interpretation of patterns in asymmetry. Furthermore, between-individual heterogeneity in developmental stability (DS) and thus fluctuating asymmetry (FA), is required to find relationships between DS and other factors. Separating directional asymmetry (DA) and antisymmetry (AS) from real FA and understanding between-individual heterogeneity in FA is therefore crucial in the analysis and interpretation of patterns in asymmetry. In this paper we introduce and explore mixture analysis to (i) identify FA, DA and AS from the distribution of the signed asymmetry, and (ii) to model and quantify between-individual heterogeneity in developmental stability and FA. In addition, we expand mixtures to the estimation of the proportion of variation in the unsigned FA that can be attributed to between-individual heterogeneity in the presumed underlying developmental stability (the so-called hypothetical repeatability). Finally, we construct weighted normal probability plots to investigate the assumption of underlying normality of the different components. We specifically show that (i) model selection based on the likelihood ratio test has the potential to yield models that incorporate nearly all heterogeneity in FA; (ii) mixtures appear to be a powerful and sensitive statistical technique to identify the different forms of asymmetry; (iii) restricted measurement accuracy and the occurrence of many zero observations results in an overestimation of the hypothetical repeatability on the basis of the model parameters; and (iv) as judged from the high correlation coefficients of the normal probability plots, the underlying normality assumption appears to hold for the empirical data we analysed. In conclusion, mixtures provide a useful statistical tool to study patterns in asymmetry.     

Associations between leaf developmental stability,variability, canalization,and phenotypic plasticity in Abutilon theophrasti

Developmental stability, canalization, and phenotypic plasticity are the most common sources of phenotypic variation, yet comparative studies investigating the relationships between these sources, specifically in plants, are lacking. To investigate the relationships among developmental stability or instability, developmental variability, canalization, and plasticity in plants, we conducted a field experiment with Abutilon theophrasti, by subjecting plants to three densities under infertile vs. fertile soil conditions. We measured the leaf width (leaf size) and calculated fluctuating asymmetry (FA), coefficient of variation within and among individuals (CV_intra and CV_inter), and plasticity (PI_rel) in leaf size at days 30, 50, and 70 of plant growth, to analyze the correlations among these variables in response to density and soil conditions, at each of or across all growth stages. Results showed increased density led to lower leaf FA, CV_intra, and PI_rel and higher CV_inter in fertile soil. A positive correlation between FA and PI_rel occurred in infertile soil, while correlations between CV_inter and PI_rel and between CV_inter and CV_intra were negative at high density and/or in fertile soil, with nonsignificant correlations among them in other cases. Results suggested the complexity of responses of developmental instability, variability, and canalization in leaf size, as well as their relationships, which depend on the strength of stresses. Intense aboveground competition that accelerates the decrease in leaf size (leading to lower plasticity) will be more likely to reduce developmental instability, variability, and canalization in leaf size. Increased developmental instability and intra‐ and interindividual variability should be advantageous and facilitate adaptive plasticity in less stressful conditions; thus, they are more likely to positively correlate with plasticity, whereas developmental stability and canalization with lower developmental variability should be beneficial for stabilizing plant performance in more stressful conditions, where they tend to have more negative correlations with plasticity.     

数量性状发育遗传模型及其分析方法的研究进展

发育遗传模型是同时反映性状遗传和发育本质、提供影响遗传变异及调整发育进程的有关因素的信息的模型。建立在群体遗传学基础上的直接效应模型适用于单一基因控制的简单性状。渐成模型将遗传变异分解成直接分量和渐成分量(母体效应和互作效应),能更好地反映有机体遗传和发育的生物学机制。生长轨迹模型有效地综合了复杂性状各分量的发育动态,可获得连续的、综合的、详细的、动态的发育信息。条件遗传分析方法不仅可以估算特定时间段的净效应,且可将净效应分解为不同遗传分量,了解各效应分量的相对贡献。 Abstract：Developmental genetic models and analysis methods for quantitative traits are presented.Developmental genetic models should reflect the genetic and developmental essence,and provide the information of the factors influencing the genetic variation and the developmental process.Direct effect models,which based on the population genetics,may be suitable to analyze simple traits with single gene.Epigenetic models can decompose the whole genetic variation into direct and epigenetic components (maternal effects and epigenetic interaction effects),so that biological mechanism can be better understood.Growth trace models effectively synthesize the developmental dynamics of components of complex traits.With them,continuous,compositive,detailed,and dynamic information of development is available.Conditional analysis method can not only estimate the net effects in a specific time interval,but also depose them into genetic components and help to appreciate the contributions of different effects.     

The effect of birthweight on tooth-size variability in twins.

Studies indicate that low birthweight (LBW) children display reduced deciduous tooth size but there is little information about permanent tooth size. It has also been shown that dental fluctuating asymmetry (FA) increases in response to various environmental influences, but the relationship between birthweight and FA remains unclear. The aim of this study was to compare tooth size and asymmetry, according to birthweight, in the deciduous and permanent dentitions of a sample of Australian twins. The study sample comprised 436 twins, classified into 2 groups: normal birthweight (NBW > 2500 g) and low birthweight (LBW < or = 2500 g). For each individual it was generally possible to measure maximum mesiodistal crown diameters of both deciduous and permanent central incisors from serial dental models. Correlations were calculated between tooth-size variables and birthweight; subsequently comparisons of tooth size and FA were made between the LBW and NBW samples using Student's t tests. Small positive correlations (around .1) were noted between birthweight and tooth-size variables. There was no evidence of tooth-size reduction in the LBW male sample, but the LBW females displayed tooth-size reduction of approximately 2-3% for both deciduous and permanent incisors, compared to the NBW females. There was no evidence of increased FA in the LBW individuals of either sex. These findings indicate that developing teeth are generally well-protected from developmental disturbances during prenatal and perinatal periods. Further research is needed to clarify the biological basis of an apparently true but weak association between tooth size and birthweight.