ANALYSIS OF GENOTYPIC DIFFERENCES IN DEVELOPMENTAL STABILITY IN ANNONA CHERIMOLA

The genetic basis of developmental stability, measured as asymmetry (fluctuating asymmetry in leaves), was analyzed in leaves and flowers of cherimoya (Annona cherimola Mill) and atemoya (A. cherimola × A. squamosa). The individuals analyzed belonged to a controlled collection of cultivars (clones) that had previously been characterized by means of isozymes. We used a nested design to analyze the differences in asymmetry at several sampling levels: individual leaves and flowers, individual trees, and genotypes. The clonal repeatability of developmental stability was not significantly different from zero, thus suggesting the absence of heritability of the asymmetry for leaves and flowers under these environmental conditions. No relationship between asymmetry and individual heterozygosity was found, but leaf fluctuating asymmetry was significantly related to particular isozymic genes. Petal and leaf size showed a phenotypically plastic response to the exposure zone of the tree (mainly due to light). Leaf fluctuating asymmetry also showed such a plastic response. No significant correlation was found between asymmetry and any pomological characters (some of these being fitness related). Finally, the hybrid species (atemoya) did not show larger developmental instability than did the parental species (cherimoya). All these data show that cherimoya asymmetry reveals the random nature of developmental noise, with developmental stability for leaves being possibly related to specific chromosome regions, but with weak evidence for genotypic differences in developmental stability.     

DEVELOPMENTAL STABILITY OF RAINBOW TROUT HYBRIDS: GENOMIC COADAPTATION OR HETEROZYGOSITY?

I compare the developmental stability of first generation hybrids between hatchery strains of rainbow trout (Salmo gairdneri) to that of the three pure parental strains raised in a common environment. Two of three reciprocal hybrid pairs show significantly less fluctuating asymmetry of four meristic characters than is found in parental strains. In contrast, the third hybrid pair shows reduced but not significantly lower developmental stability compared to pure strains. These hybrids had previously been found to develop slower than their maternal parental strains, indicating divergence of parental regulatory mechanisms controlling early ontogeny. A significant positive association between the degree of relative delay in hybrid developmental rate and their degree of developmental instability supports this view. For example, the only hybrid pair with decreased developmental stability also had the largest relative delay in development time of all hybrids. Neither absolute developmental rate nor enzyme heterozygosity at 42 loci alone can explain the degree of fluctuating asymmetry in these hybrids. The developmental stability of hybrids is apparently a result of the interaction between the developmental divergence between parental strains and their genomic heterozygosity due to hybridization.     

A meta-analysis of the heritability of developmental stability

The existence of additive genetic variance in developmental stability has important implications for our understanding of morphological variation. The heritability of individual fluctuating asymmetry and other measures of developmental stability have frequently been estimated from parent-offspring regressions, sib analyses, or from selection experiments. Here we review by meta-analysis published estimates of the heritability of developmental stability, mainly the degree of individual fluctuating asymmetry in morphological characters. The overall mean effect size of heritabilities of individual fluctuating asymmetry was 0.19 from 34 studies of 17 species differing highly significantly from zero (P < 0.0001). The mean heritability for 14 species was 0.27. This indicates that there is a significant additive genetic component to developmental stability. Effect size was larger for selection experiments than for studies based on parent-offspring regression or sib analyses, implying that genetic estimates were unbiased by maternal or common environment effects. Additive genetic coefficients of variation for individual fluctuating asymmetry were considerably higher than those for character size per se. Developmental stability may be significantly heritable either because of strong directional selection, or fluctuating selection regimes which prevent populations from achieving a high degree of developmental stability to current environmental and genetic conditions.     

THE GENETIC BASIS OF DEVELOPMENTAL STABILITY IN APIS MELLIFERA: HETEROZYGOSITY VERSUS GENIC BALANCE

The genetic basis for developmental stability in the haplo-diploid honeybee Apis mellifera was determined by comparing the level of asymmetry between diploid females and haploid males both among and within inbreeding levels. There was no significant relationship between the level of inbreeding and the level of fluctuating asymmetry for both females and males. It is therefore argued that the general level of genomic heterozygosity is not an important factor for the determination and maintenance of developmental stability in this system, but rather that the balance of genes within chromosomes plays the major role. The observation that males were generally more asymmetric than females suggests that developmental stability in females may also be influenced by additional factors such as gene dosage, sex-limited genes or cytoplasmic elements.     

FLUCTUATING ASYMMETRY IN CENTRAL AND MARGINAL POPULATIONS OF LYCHNIS VISCARIA IN RELATION TO GENETIC AND ENVIRONMENTAL FACTORS

Developmental instability in the form of increased fluctuating asymmetry can be caused by either genetic or environmental stress. Because extinctions can be attributed broadly to these factors, fluctuating asymmetry may provide a sensitive tool for detecting such stresses. We studied the level of fluctuating asymmetry of flowers of a perennial outcrossing plant species, Lychnis viscaria, both in natural and common-garden populations. The degree of flower asymmetry was higher in small, isolated, and marginal populations of the species range. These marginal populations also were the most homozygous. In the core area of the species' range, flowers were more symmetrical The level of asymmetry was correlated with both population size and heterozygosity. However, a partial correlation analysis revealed that when the impact of population size was controlled for, there was a negative relationship between fluctuating asymmetry and heterozygosity, whereas when controlling for heterozygosity, no relationship between population size and fluctuating asymmetry was found. This indicates that genetic consequences of small population size probably underlie the relationship between the level of asymmetry and population size. Results from a transplantation experiment showed that individuals subjected to a higher environmental stress had an increased level of asymmetry compared to control plants. In the common-garden conditions the level of fluctuating asymmetry did not differ between the central and marginal populations. This suggests that presumably both genetic and environmental factors affected to the higher level of asymmetry among marginal populations compared to central ones. In all we conclude that even though fluctuating asymmetry seems to be a sensitive tool for detecting stresses, results from studies focusing on only one factor should be interpreted with caution.     

Developmental stability in brown trout: are there any effects of heterozygosity or environmental stress?

We studied the developmental stability of brown trout, Salmo trutta L., in 10 populations (five acidified, five control) in Norway, measured as fluctuating asymmetry (FA) and departure from the morphological norm. We measured four meristic and four morphometric characters, and scored the level of biochemical heterozygosity at 49 loci (20 polymorphic). We reared eggs of a single population in a hatchery using four different water qualities (three replicates of each treatment) to test the effect of acidification stress on developmental instability. There were no significant differences in the level of FA, in departure from the morphological norm between brown trout sampled from lakes with acidified or control water qualities, or in brown trout hatched at different water qualities. There was no correlation between level of heterozygosity and FA or departure from the morphological norm, either when tested within populations or among populations. There were no single-locus effects on developmental stability tested for 11 loci. We conclude that measures of developmental stability or morphological variability are not useful for detecting acidification stress in brown trout. Furthermore, we conclude that developmental stability in our material varies independently of heterozygosity.     

A study of wing morphology and fluctuating asymmetry in interspecific hybrids between Drosophila buzzatii and D. koepferae

In this work we investigate the effect of interspecific hybridization on wing morphology using geometric morphometrics in the cactophilic sibling species D. buzzatii and D. koepferae. Wing morphology in F₁ hybrids exhibited an important degree of phenotypic plasticity and differs significantly from both parental species. However, the pattern of morphological variation between hybrids and the parental strains varied between wing size and wing shape, across rearing media, sexes, and crosses, suggesting a complex genetic architecture underlying divergence in wing morphology. Even though there was significant fluctuating asymmetry for both, wing size and shape in F₁ hybrids and both parental species, there was no evidence of an increased degree of fluctuating asymmetry in hybrids as compared to parental species. These results are interpreted in terms of developmental stability as a function of a balance between levels of heterozygosity and the disruption of coadaptation as an indirect consequence of genomic divergence.     

Developmental stability and canalization of limb bones of brown haresLepus europaeus with varying levels of heterozygosity

We have studied fluctuating asymmetry (FA), as indicator of developmental stability, and between-individual variation, as surrogate of developmental canalization (DC), in long bones (humerus, ulna, radius, femur, tibia) of 72 wild-living adult-sized brown haresLepus europaeus Pallas, 1778 with variable individual heterozygosity (H).H was calculated from 13 polymorphic allozyme loci. According to the “over-dominance hypothesis”, we expected increased developmental stability and canalization at higherH-levels. But at the individual level we did not find any significant correlation between overall FA (FA_I) andH. Also, standard deviations (SD) of mean length (over both body sides) of bones did not differ between individuals from two intentionally created groups of hares, namely one with high and one with lowH. FA-indices and variances of FA-indices of bone lengths did not differ significantly when compared between two intentionally created groups of hares with high and low SD of bone lengths, respectively. These latter findings suggest that developmental stability and DC are two separate or partly separate mechanisms of developmental homeostasis in the studied appendicular skeleton, and thatH has no traceable effect on develop-mental homeostasis. If there is still such an effect, it should be clearly smaller than a possibly combined effect of (presently uncontrolled) environmental stressors.     

META-ANALYSES OF THE ASSOCIATION BETWEEN MULTILOCUS HETEROZYGOSITY AND FITNESS

Meta-analyses of published correlation coefficients between multilocus heterozygosity (MLH) and two fitness surrogates, growth rate and fluctuating asymmetry, suggested that the strength of these correlations are generally weak. A variety of plants and animals was included in the meta-analyses. A statistically homogeneous group of MLH–growth rate correlation coefficients that included both plants and animals yielded a common correlation of r_z = 0.133. A common correlation of r_z = –0.170 was estimated for correlations between MLH and fluctuating asymmetry in three species of salmonid fishes. These results suggest that selection, including overdominance, has at most a weak effect at allozyme loci and cast some doubt on the widely held notion that heterozygosity and individual fitness are strongly correlated.     

The genetic basis of developmental stability in Apis mellifera II. Relationships between character size, asymmetry and single-locus heterozygosity

Haploid male and diploid female honey bees, Apis mellifera, from colonies headed by queens polymorphic at the malate dehydrogenase (MDH) locus were examined for the influence of MDH genotype and heterozygosity on the size and asymmetry of six morphological characters. Although there were significant differences among MDH genotypes for mean character size within colonies, these effects were inconsistent between colonies. There were no significant relationships between MDH genotype or heterozygosity and asymmetry, indicating that genetic variation at this locus has no impact on developmental stability in this species.     

Fluctuating asymmetry for specific bristle characters in Notch mutants of Drosophila melanogaster

Asymmetry has been used as a measure of developmental stability for bilaterally symmetrical organisms. Most studies have failed to partition the genetic and environmental contributions to the asymmetry phenotype due to the limitations of the systems used or the shortcomings in experimental design. The Notch mutants of Drosophila melanogaster were used to study the genetic contribution to asymmetry for six different bristle characters. Asymmetry response was character specific for the mutants examined. For N ^spl, N ^Co, N ^264–47, Ax ^71d, Ax ^9B2, Ax ^E2, l(1)N ^B and nd ² significant asymmetry responses, relative to wildtype Canton‐S, were observed for some characters. N ^60g11 and nd ¹ did not exhibit significant asymmetry for any of the characters examined. All of the mutants except N ^60g11 and nd ¹ showed thoracic bristle asymmetry. However, when asymmetry scores were pooled over the five bristle characters which individually exhibited fluctuating asymmetry, no significant differences were found between any genotypes. Therefore pooling asymmetry values across characters obscures the significant character specific asymmetry values observed. Thus caution is necessary when using the asymmetry phenotype of specific characters to draw organism wide conclusions about developmental stability. This revised version was published online in July 2006 with corrections to the Cover Date.     

DEVELOPMENTAL STABILITY IN LEAVES OF CLARKIA TEMBLORIENSIS (ONAGRACEAE) AS RELATED TO POPULATION OUTCROSSING RATES AND HETEROZYGOSITY

Four natural populations of Clarkia tembloriensis, whose levels of heterozygosity and rates of outcrossing were previously found to be correlated, are examined for developmental instability in their leaves. From the northern end of the species range, we compare a predominantly selfing population (t? = 0.26) with a more outcrossed population (t? = 0.84), which is genetically similar. From the southern end of the range, we compare a highly selfing population (t? = 0.03) with a more outcrossed population (t? = 0.58). We measured developmental stability in the populations using two measures of within-plant variation in leaf length as well as calculations of fluctuating asymmetry (FA) for several leaf traits. Growth-chamber experiments show that selfing populations are significantly more variable in leaf length than more outcrossed populations. Developmental instability can contribute to this difference in population-level variance. Plants from more homozygous populations tend to have greater within-plant variance over developmentally comparable nodes than plants from more heterozygous populations, but the difference is not significant. At the upper nodes of the plant, mature leaf length declines steadily with plant age, allowing for a regression of leaf length on node. On average, the plants from more homozygous populations showed higher variance about the regression (MSE) and lower R² values, suggesting that the decline in leaf length with plant age is less stable in plants from selfing populations than in plants from outcrossing populations. Fluctuating asymmetry (FA) was calculated for four traits within single leaves at up to five nodes per plant. At the early nodes of the plant where leaf arrangement is opposite, FA was also calculated for the same traits between opposite leaves at a node. Fluctuating asymmetry is significantly greater in the southern selfing population than in the neighboring outcrossed population. Northern populations do not differ in FA. Fluctuating asymmetry can vary significantly between nodes. The FA values of different leaf traits were not correlated. We show that developmental stability can be measured in plants using FA and within-plant variance. Our data suggest that large differences in breeding system are associated with differences in stability, with more inbred populations being the least stable.     

Null alleles at two lactate dehydrogenase loci in rainbow trout are associated with decreased developmental stability

Previous studies with rainbow trout (Oncorhynchus mykiss) have shown that allozymic heterozygotes have increased developmental stability, as measured by reduced fluctuating bilateral asymmetry. In this paper, we examine the phenotypic effects of null alleles at two lactate dehydrogenase (LDH) loci. If the association between allozymic heterozygosity and developmental stability is due largely to linked chromosomal segments, then we would expect null allele heterozygotes to have increased developmental stability. In contrast, heterozygotes for LDH null alleles in three populations have reduced developmental stability. This suggests that the reduction in enzyme activity at these loci is having a deleterious effect on development that is strong enough to mask any beneficial effects that may be associated with heterozygosity for these chromosomal segments. The LDH loci examined in this study are members of two different paralogous pairs of duplicate genes produced by the polyploidization of the ancestral salmonid genome. The apparent deleterious effects of these null alleles in heterozygotes could retard the possible loss of duplicate gene expression.     

MORPHOLOGICAL VARIABILITY AND ASYMMETRY IN THE CHEETAH (ACINONYX JUBATUS), A GENETICALLY UNIFORM SPECIES

The African cheetah (Acinonyx jubatus) is an unusual species because of its extremely low amount of biochemical genetic variation. A comparative analysis of morphological variation of 16 cranial characters from four species of Felidae (ocelot, Leopardus pardalus; margay, L. wiedii; leopard, Panthera pardus; and cheetah) was undertaken to evaluate the consequence of biochemical monomorphism on morphological variation. The species were selected because the cheetah has been shown previously to possess extremely low amounts of biochemical genetic variation as opposed to the other three species which retain comparatively high levels of allozyme heterozygosity. The cheetah sample showed dramatically greater fluctuating asymmetry but was not outstanding in morphological variability. Elevated levels of fluctuating asymmetry have been interpreted as a reflection of developmental instability, which is a common consequence of inbreeding. The inverse correlation of genetic variation and developmental stability (homeostasis) observed here fulfills prior expectations and further emphasizes the genetic invariability of the cheetah species.     

Relationship between fluctuating asymmetry,morphological modality and heterozygosity in an elderly Israeli population

We have been checking the following working hypotheses: 1) There is a negative correlation between genetic heterozygosity and fluctuating asymmetry (FA); 2) FA is a measure of developmental stability/instability of the whole organism, i.e. we expect negative correlation between FA and morphological proximity of a set of mass-size variables of an individual to a population centroid; and 3) FA is a measure of character-specific stability in a population, i.e. we expect correlation between magnitude of FA and deviation of an individual from the population centroid of the bilateral characters themselves. For this purpose each individual in a sample of about 200 elderly individuals was assessed for 11 polymorphic blood systems (14 genetic loci) as well as for a set of 26 anthropometric traits: 1) a set of ten mass size variables; and 2) a set of eight pairs of bilateral measurements. Four multivariate measures of morphological centrality were computed, two measures for size and two measures of shape distances from the i^th individual to the population centroid for mass-size variables and also for the bilateral variables. A multivariate measure of FA for 8 bilateral pairs was also computed. No relationship was detected between FA and heterozygosity, or between FA and any of the four multivariate deviations. Thus, we concluded that our data do not support the listed hypotheses.     

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.     

GENOMIC COADAPTATION AND DEVELOPMENTAL STABILITY WITHIN INTROGRESSED POPULATIONS OF ENNEACANTHUS GLORIOSUS AND E. OBESUS (PISCES,CENTRARCHIDAE)

Using fluctuating bilateral asymmetry as a measure of developmental stability, we tested the hypothesis that genomic coadaptation mediates developmental stability in natural populations. Hybrid populations were more asymmetrical than populations of the parental species, and ranks of overall developmental instability were positively correlated with ranks of mean heterozygosity in these populations. The failure to find increased asymmetry in previous studies of natural hybrid populations (Jackson, 1973a, 1973b; Felley, 1980) suggests that such populations may have re-evolved coadapted genomes. Increased asymmetry in hybrid Enneacanthus populations may reflect the youthfulness of these populations.     

Fluctuating asymmetry and genetic variability in the roe deer (Capreolus capreolus): a test of the developmental stability hypothesis in mammals using neutral molecular markers

Fluctuating asymmetry (FA), used as an indicator of developmental stability, has long been hypothesized to be negatively correlated with genetic variability as a consequence of more variable organisms being better suited to buffer developmental pathways against environmental stress. However, it is still a matter of debate if this is due to metabolic properties of enzymes encoded by certain key loci or rather to overall genomic heterozygosity. Previous analyses suggest that there might be a general difference between homeo- and poikilotherms in that only the latter tend to exhibit the negative correlation predicted by theory. In the present study, we addressed these questions by analysing roe deer (Capreolus capreolus) from five German populations with regard to FA in metric and non-metric skull and mandible traits as well as variability at eight microsatellite loci. Genetic variability was quantified by heterozygosity and mean d2 parameters, and although the latter did not show any relationship with FA, we found for the first time a statistically significant negative correlation of microsatellite heterozygosity and non-metric FA among populations. Because microsatellites are non-coding markers, this may be interpreted as evidence for the role of overall genomic heterozygosity in determining developmental stability. To test if the threshold character of non-metric traits is responsible for the metric vs non-metric difference we also carried out calculations where we treated our metric traits as threshold values. This, however, did not yield significant correlations between FA and genetic variability either.     

Enzyme heterozygosity,metabolism, and developmental stability

Developmental homeostasis, measured as either fluctuating asymmetry or variance of morphological characters, increases with enzyme heterozygosity in many, but not all, natural populations. These results have been reported forDrosophila, monarch butterflies, honeybees, blue mussels, side-blotched lizards, killifish, salmonid fishes, guppies, Sonoran topminnows, herring, rufous-collared sparrows, house sparrows, brown hares, white-tailed deer, and humans. Because heterozygosity at a few loci can not predict heterozygosity of the entiry genome, these loci must be detecting localized zones that influence the developmental environment. Studies of malate dehydrogenase in honeybees,Apis mellifera, and lactate dehydrogenase in killifish,Fundulus heteroclitus, revealed that developmental homeostasis varied with heterozygosity of individual loci. Heterozygotes differed from homozygotes in fluctuating asymmetry, morphological variance, and in correlations between morphological characters. The protein loci in these studies code for enzymes, and therefore do not directly influence morphological characters. However, some enzymatic loci substantially influence metabolism, and contribute to variation in the amount of energy available for development and growth. This argument can be made most convincingly for the LDH polymorphism in killifish. LDH genotypes differ in enzyme kinetic properties that measure differences in physiological efficiency, and these differences produce measurable and predictable differences in physiology and development. Under environmental conditions which impose a stress upon development, genotypes at these loci may have different amounts of energy available for development, and consequently exhibit different levels of developmental homeostasis.     

Bilateral variation and the evolutionary origin of macroscopic asymmetries

Given that characters exhibiting macroscopic asymmetry have evolved in a wide variety of taxa, heritable variation for bilateral asymmetry must have arisen at some point in their history. The recognition that heritable variation may underlie some statistical asymmetries not only raises concerns about the incautious use of statistical estimates of FA in studies of developmental stability, but it suggests some intriguing questions about the possible evolutionary origins of macroscopic asymmetries. First, we developed an additive model of bilateral variation based on some simple assumptions about the developmental control of bilateral variation. Second, using a new approach for studying statistical asymmetries, we conducted an analysis of bilateral variation in eight metrical traits of a harpacticoid copepod (Tigriopus californicus) to search for novel forms of statistical asymmetries. The model we developed revealed three independent statistical asymmetries of potential evolutionary significance:a) a previously unrecognized form of asymmetry (referred to here asnormal covariant asymmetry),b) antisymmetry, andc) directional asymmetry. Because each pattern of variation would seem to require different amounts and kinds of developmental-genetic information [a- only negative feedback between sides (bilateral inhibition),b- bothbilateral inhibition and average departure from symmetry (bilateral offset),c- bilateral inhibition, bilateral offset, and a consistent overdevelopment of one side or the other (side-bias control)], those requiring less information would seem more likely to represent earlier stages in the evolution of macroscopic asymmetries. Our analysis of bilateral variation inTigriopus revealed no evidence for any form of statistical asymmetry other than fluctuating asymmetry. However, a significant positive covariation between sides, even after correction for body size variation, suggested that factors influencing relative limb length (whether genetic or environmental) affected both sides equally rather than one side at the expense of the other. Finally, we note that certain statistical asymmetries (directional asymmetry, any form of covariant asymmetry) may render characters unreliable for estimating developmental stability because, unlike pure fluctuating asymmetry, they may signal a genetic component to asymmetry variation.