Postmarital residence and within-sex genetic diversity among the Urubu-Ka'apor Indians, Brazilian Amazon

The analysis of biologic variation in prehistoric human populations separately by sex has been used as a tool to recover post-marital residential rules. These studies, which focus on the sexual distribution of skeletal traits, assume that the degree of intragroup or intergroup biologic diversity is higher in one sex with regard to unilocality (uxori- or virilocality). Despite a recent attempt to interpret this phenomenon in terms of population genetics (Konigsberg 1988), the main assumption has never been tested in situations in which the real residential practice of an indigenous population is known and in which genetic rather than phenotypic data are available. We investigated the within-group and between-group genetic variability among males and females from 4 villages of an uxorilocal Amazonian tribe, the Urubu-Ka'apor, on the basis of 20 polymorphic loci. The results were only partly concordant with the expected. Individual mean per locus heterozygosities were not different between the sexes, and the analysis of genetic heterogeneity showed similar gene frequencies for males and females in all villages. On the other hand, the intergroup approach detected a level of variation significantly greater among females than among males. The ethnographic evidence shows that three of the four subgroups studied belong to the same gamic unity, with the fourth subgroup belonging to another gamic network. Within-sex differences in intergroup analysis turned out to be more evident; yet, when those 3 villages were investigated separately, the female FST (0.0609) proved to be significantly higher than the male FST (0.0218). Such results suggest that the intergroup analysis is more sensitive to the genetic effects of differential migration rates between the sexes. In prehistoric contexts, therefore, an intergroup genetic approach can provide more reliable grounds for sociocultural inferences.     

Intraspecific Trait Variation Driven by Plasticity and Ontogeny in Hypochaeris radicata

The importance of intraspecific variation in plant functional traits for structuring communities and driving ecosystem processes is increasingly recognized, but mechanisms governing this variation are less studied. Variation could be due to adaptation to local conditions, plasticity in observed traits, or ontogeny. We investigated 1) whether abiotic stress caused individuals, maternal lines, and populations to exhibit trait convergence, 2) whether trait variation was primarily due to ecotypic differences or trait plasticity, and 3) whether traits varied with ontogeny. We sampled three populations of Hypochaeris radicata that differed significantly in rosette diameter and specific leaf area (SLA). We grew nine maternal lines from each population (27 lines total) under three greenhouse conditions: ambient conditions (control), 50% drought, or 80% shade. Plant diameter and relative chlorophyll content were measured throughout the experiment, and leaf shape, root∶shoot ratio, and SLA were measured after five weeks. We used hierarchical mixed-models and variance component analysis to quantify differences in treatment effects and the contributions of population of origin and maternal line to observed variation. Observed variation in plant traits was driven primarily by plasticity. Shade significantly influenced all measured traits. Plant diameter was the only trait that had a sizable proportion of trait variation (30%) explained by population of origin. There were significant ontogenetic differences for both plant diameter and relative chlorophyll content. When subjected to abiotic stress in the form of light or water limitation, Hypochaeris radicata exhibited significant trait variability. This variation was due primarily to trait plasticity, rather than to adaptation to local conditions, and also differed with ontogeny.     

Craniometric variation,genetic theory,and modern human origins

Recent controversies surrounding models of modern human origins have focused on among-group variation, particularly the reconstruction of phylogenetic trees from mitochondrial DNA (mtDNA) and, the dating of population divergence. Problems in tree estimation have been seen as weakening the case for a replacement model and favoring a multiregional evolution model. There has been less discussion of patterns of within-group variation, although the mtDNA evidence has consistently shown the greatest diversity within African populations. Problems of interpretation abound given the numerous factors that can influence within-group variation, including the possibility of earlier divergence, differences in population size, patterns of population expansion, and variation in migration rates. We present a model of within-group phenotypic variation and apply it to a large set of craniometric data representing major Old World geographic regions (57 measurements for 1,159 cases in four regions: Europe, Sub-Saharan Africa, Australasia, and the Far East). The model predicts a linear relationship between variation within populations (the average within-group variance) and variation between populations (the genetic distance of populations to pooled phenotypic means). On a global level this relationship should hold if the long-term effective population sizes of each region are correctly specified. Other potential effects on withingroup variation are accounted for by the model. Comparison of observed and expected variances under the assumption of equal effective sizes for four regions indicates significantly greater within-group variation in Africa and significantly less within-group variation in Europe. These results suggest that the long-term effective population size was greatest in Africa. Closer examination of the model suggests that the long-term African effective size was roughly three times that of any other geographic region. Using these estimates of relative population size, we present a method for analyzing ancient population structure, which provides estimates of ancient migration. This method allows us to reconstruct migration history between geographic regions after adjustment for the effect of genetic drift on interpopulational distances. Our results show a clear isolation of Africa from other regions. We then present a method that allows direct estimation of the ancient migration matrix, thus providing us with information on the actual extent of interregional migration. These methods also provide estimates of time frames necessary to reach genetic equilibrium. The ultimate goal is extracting as much information from present-day patterns of human variation relevannt to issues of human origins. Our results are in agreement with mismatch distribution analysis of mtDNA, and they support a “weak Garden o Eden” model. In this model, modern-day variation can be explained by divergence from an initial source (perhaps Africa) into a number o small isolated populations, followed by later population expansion throughout our species. The major populationn expansions of Homo sapiens during and after the late Pleistocene have had the effect of “freezing” ancient patterns of population structure. While this is not the only possible scenario, we do note the close agreement with ecent analyses of mtDNA mismatch distibutions. © 1994 Wiley-Liss, Inc.     

Intraspecific variation in fruit size and shape in Corema album (Ericaceae) along a latitudinal gradient: from fruits to populations

Plant trait variation can be the result of environmental variability, developmental instability, and plasticity, although it can also arise from previous selective pressures on fruit traits themselves or directly on their variation. We aimed to quantify fruit size and shape variability at within‐plant, among‐plant, among‐population, and among‐year levels. For 2 years, we measured fruit size and shape along the geographical range of Corema album, aiming to determine whether (1) the population level among‐plant variation is lower than within‐plant and among‐population variation and (2) the ratio of within‐plant to among‐plant variation follows a latitudinal pattern. Levels of fruit variation were in accordance with the reported mean levels for reproductive organs. Most variance concentrated on within‐ and among‐individual levels for size, showing higher values for among‐individual variation in fruit shape. Although fruit size retained important variation among populations, this source of variance was negligible for fruit shape. This difference could arise from contrasting mechanical or developmental constraints. Despite the marked climatic differences along the latitudinal range of the species, latitude did not affect the ratio of within‐ to among‐plant variation. We show that an analysis of the often‐disregarded sources of intraspecific variation can produce unexpected results that deserve further investigation.     

Morphological variation and population structure on the island of Korčula,Yugoslavia

A comparison of head and body morphology on the ecologically uniform island of Korčula provides a good example of the balance that occurs in microevolution between selective pressures toward homogenization and selective inertia toward heterogeneity. thirty-eight measurements were made from a sample of 471 males and 526 females. Head variables, being more eco-stable than body variables, remain relatively more different between two distinct populations (eastern and western villages) than do the more eco-labile body variables, although both do vary significantly between east and west. The differences apply to both men and women. These east-west contrasts reflect the differential migration of Slavs to the island over the past three centuries, with a new wave of immigrants settling mainly in the east and introducing a new gene pool to the pre-existing Slavic hybrid population which had settled the island in the sixth to the eighth centuries. In addition to the predominant east-west differences in morphology, we also find significant variation between all villages in both head and body variables. We conclude that this is an indication of the considerable reproductive isolation that has persisted between all villages until the post-World War II period. This paper is a part of a study conducted with financial suport from the joint Yugoslav-United States Board for Bilateral Research, as a common Institute for Medical Research and Occupational Health (Zagreb) and Smithsonian Institution (Washington, DC) research, projects JFP-429 and JFP-674.     

Multiple group principal component analysis and population differentiation

This paper explores the requirements and advantages of multiple group principal component analysis (MGPCA) when it is used to investigate population differentiation. A distinction is drawn between equality of orientation of the within-group axes and equality of variance along these axes. Several examples of the use of MGPCA are discussed and it is shown that MGPCA per se does not require equality of variance along the axes although it may be a requirement of some of the techniques subsequently used to analyse the component scores. MGPCA is simple and direct, being based on the mathematically well defined eigenvector analysis of a symmetric positive definite (pooled within-group covariance) matrix and it can be thought of as a step in the computation of canonical variate analysis (CVA). It can be used with CVA (which is the most popular method of biometrically assessing population affinities) to assess the contribution of within-group components to among-group discrimination. It is also one of a range of appropriate techniques that can be used to define (and delete if required) within-group growth effects and is particularly suitable when CVA is being used to assess the population affinities. When used in this way it has the advantage of being more influenced by the groups with the greatest growth range.     

Prey‐associated head‐size variation in an invasive lizard in the Hawaiian Islands

Biological invasions are recognized as a primary driver of large‐scale changes in global ecosystems. This study addresses ecomorphological variation in head size within and among populations of an ecologically destructive invasive predator, and evaluates the potential roles of environmental components in phenotypic differentiation. We used four size‐corrected measurements of head morphology in Jackson's chameleons, Trioceros jacksonii xantholophus (N = 319), collected from multiple Hawaiian Islands to assess phenotypic variation among and within islands. Results of analysis of variance (ANOVA) comparing chameleon head size (PC1) among islands revealed significant differences (mean difference > 5%) associated with variation in both rainfall and diet composition using Mann–Whitney U‐tests and chi‐squared analyses. These results suggest that morphological differentiation among populations from different islands has occurred over a relatively short ecological timescale, and is likely the result of ecomorphological adaptation to differences in exploited prey hardness. Intra‐island allopatric population variation, however, was also detected in this study. Although we might expect that genetic change is the more likely explanation for differences between islands than within, and that plasticity may be more likely an explanation for the within‐ than the between‐island differences, it is also possible that both within‐ and between‐island patterns are the results of genetic change, or of plasticity.     

Statistical assessment of population variability: a methodological approach

We have implemented a multivariate statistical methodology to assess the degree and pattern of cranial variability in skeletal samples. Specifically, the method is designed to test whether variability in a skeletal sample exceeds "normal within-group variability" defined in the present instance as variability present among crania from a cemetery associated with a village. It involves comparing a covariance matrix derived from a sample of unknown composition to one representing "normal within-group variability." The method has been applied to two Plains Indian craniometric samples. The Leavenworth site (39CO9) represents the remnants of previously autonomous Arikara Indian villages devastated by epidemic diseases in the late 1700s. The Bad River 2 Phase is an archaeological designation grouping together closely related sites in the Bad-Cheyenne region of South Dakota dating from 1740-1795 AD. We were able to show substantial heterogeneity among crania from Leavenworth. District burial areas at Leavenworth account for some of the heterogeneity, supporting the notion that they represent an attempt to maintain former social distinctions. We were unable to differentiate among sites within the Bad River 2 Phase, suggesting that it is a valid biological unit.     

Intra-individual variation in sperm tail length in murine rodents

In eutherian mammals, there are marked interspecific differences in sperm head shape and tail length. In a few species, sperm head variability occurs but intra-individual variation in sperm tail length has rarely been investigated or commented upon. Here, we ask the question: Do murine rodent species that have variable sperm head shapes exhibit greater intra-individual variation in sperm midpiece and total tail lengths than closely related species where little, or no, sperm head variability occurs? From three separate lineages, we selected three pairs of murine rodents, one of which has monomorphic, and the other variable, sperm head shape. These were from southern Asia the bandicoot rats Bandicota bengalensis and Bandicota indica , from southern Africa the veld rats, Aethomys chrysophilus and Aethomys ineptus and from Australia the fawn hopping mouse Notomys cervinus and the spinifex hopping mouse Notomys alexis . Cauda epididymal sperm smears were prepared and sperm midpiece and total tail lengths were determined. A linear mixed-effects model was used to estimate intra-individual variance. The results showed that in all three species where there are variable sperm head shapes ( B. indica , A. ineptus and N. alexis ), statistically significantly greater intra-individual variability of sperm midpiece and total tail lengths occurs ( P <0.0001 in all cases). These species all have relatively smaller testes mass compared with the closely related species with monomorphic sperm populations. This suggests that depressed levels of intermale sperm competition may result in the occurrence of variability in not only the divergent sperm head shape but also in the length of the midpiece as well as that of the total length of the sperm tail.     

Melanesian mtDNA complexity

Melanesian populations are known for their diversity, but it has been hard to grasp the pattern of the variation or its underlying dynamic. Using 1,223 mitochondrial DNA (mtDNA) sequences from hypervariable regions 1 and 2 (HVR1 and HVR2) from 32 populations, we found the among-group variation is structured by island, island size, and also by language affiliation. The more isolated inland Papuan-speaking groups on the largest islands have the greatest distinctions, while shore dwelling populations are considerably less diverse (at the same time, within-group haplotype diversity is less in the most isolated groups). Persistent differences between shore and inland groups in effective population sizes and marital migration rates probably cause these differences. We also add 16 whole sequences to the Melanesian mtDNA phylogenies. We identify the likely origins of a number of the haplogroups and ancient branches in specific islands, point to some ancient mtDNA connections between Near Oceania and Australia, and show additional Holocene connections between Island Southeast Asia/Taiwan and Island Melanesia with branches of haplogroup E. Coalescence estimates based on synonymous transitions in the coding region suggest an initial settlement and expansion in the region at approximately 30-50,000 years before present (YBP), and a second important expansion from Island Southeast Asia/Taiwan during the interval approximately 3,500-8,000 YBP. However, there are some important variance components in molecular dating that have been overlooked, and the specific nature of ancestral (maternal) Austronesian influence in this region remains unresolved.     

Variation in short tandem repeats is deeply structured by genetic background on the human Y chromosome

Eleven biallelic polymorphisms and seven short-tandem-repeat (STR) loci mapping on the nonrecombining portion of the human Y chromosome have been typed in men from northwestern Africa. Analysis of the biallelic markers, which represent probable unique events in human evolution, allowed us to characterize the stable backgrounds or haplogroups of Y chromosomes that prevail in this geographic region. Variation in the more rapidly mutating genetic markers (STRs) has been used both to estimate the time to the most recent common ancestor for STR variability within these stable backgrounds and to explore whether STR differentiation among haplogroups still retains information about their phylogeny. When analysis of molecular variance was used to study the apportionment of STR variation among both genetic backgrounds (i.e., those defined by haplogroups) and population backgrounds, we found STR variability to be clearly structured by haplogroups. More than 80% of the genetic variance was found among haplogroups, whereas only 3.72% of the genetic variation could be attributed to differences among populations-that is, genetic variability appears to be much more structured by lineage than by population. This was confirmed when two population samples from the Iberian Peninsula were added to the analysis. The deep structure of the genetic variation in old genealogical units (haplogroups) challenges a population-based perspective in the comprehension of human genome diversity. A population may be better understood as an association of lineages from a deep and population-independent gene genealogy, rather than as a complete evolutionary unit.     

Variation in estimated recombination rates across human populations

Recently it has been reported that recombination hotspots appear to be highly variable between humans and chimpanzees, and there is evidence for between-person variability in hotspots, and evolutionary transience. To understand the nature of variation in human recombination rates, it is important to describe patterns of variability across populations. Direct measurement of recombination rates remains infeasible on a large scale, and population-genetic approaches can be imprecise, and are affected by demographic history. Reports to date have suggested broad similarity in recombination rates at large genomic scales and across human populations. Here, we examine recombination rate estimates at a finer population and genomic scale: 28 worldwide populations and 107 SNPs in a 1 Mb stretch of chromosome 22q. We employ analysis of variance of recombination rate estimates, corrected for differences in effective population size using genome-wide microsatellite mutation rate estimates. We find substantial variation in fine-scale rates between populations, but reduced variation within continental groups. All effects examined (SNP-pair, region, population and interactions) were highly significant. Adjustment for effective population size made little difference to the conclusions. Observed hotspots tended to be conserved across populations, albeit at varying intensities. This holds particularly for populations from the same region, and also to a considerable degree across geographical regions. However, some hotspots appear to be population-specific. Several results from studies on the population history of humans are in accordance with our analysis. Our results suggest that between-population variation in DNA sequences may underly recombination rate variation.     

A Method for Detecting Positive Growth Autocorrelation without Marking Individuals

In most ecological studies, within-group variation is a nuisance that obscures patterns of interest and reduces statistical power. However, patterns of within-group variability often contain information about ecological processes. In particular, such patterns can be used to detect positive growth autocorrelation (consistent variation in growth rates among individuals in a cohort across time), even in samples of unmarked individuals. Previous methods for detecting autocorrelated growth required data from marked individuals. We propose a method that requires only estimates of within-cohort variance through time, using maximum likelihood methods to obtain point estimates and confidence intervals of the correlation parameter. We test our method on simulated data sets and determine the loss in statistical power due to the inability to identify individuals. We show how to accommodate nonlinear growth trajectories and test the effects of size-dependent mortality on our method''s accuracy. The method can detect significant growth autocorrelation at moderate levels of autocorrelation with moderate-sized cohorts (for example, statistical power of 80% to detect growth autocorrelation ρ ² = 0.5 in a cohort of 100 individuals measured on 16 occasions). We present a case study of growth in the red-eyed tree frog. Better quantification of the processes driving size variation will help ecologists improve predictions of population dynamics. This work will help researchers to detect growth autocorrelation in cases where marking is logistically infeasible or causes unacceptable decreases in the fitness of marked individuals.     

Microsatellite variation among divergent populations of stalk-eyed flies, genus Cyrtodiopsis

Microsatellite primers are often developed in one species and used to assess neutral variability in related species. Such analyses may be confounded by ascertainment bias (i.e. a decline in amplification success and allelic variability with increasing genetic distance from the source of the microsatellites). In addition, other factors, such as the size of the microsatellite, whether it consists of perfect or interrupted tandem repeats, and whether it is autosomal or X-linked, can affect variation. To test the relative importance of these factors on microsatellite variation, we examine patterns of amplification and allelic diversity in 52 microsatellite loci amplified from five individuals in each of six populations of Cyrtodiopsis stalk-eyed flies that range from 2.2 % to 11.2% mitochondrial DNA sequence divergence from the population used for microsatellite development. We find that amplification success and most measures of allelic diversity declined with genetic distance from the source population, in some cases an order of magnitude faster than in birds or mammals. The median and range of the repeat array length did not decline with genetic distance. In addition, for loci on the X chromosome, we find evidence of lower observed heterozygosity compared with loci on autosomes. The differences in variability between X-linked and autosomal loci are not adequately explained by differences in effective population sizes of the chromosomes. We suggest, instead, that periodic selection events associated with X-chromosome meiotic drive, which is present in many of these populations, reduces X-linked variation.     

Geographical gradients in the population dynamics of North American prairie ducks

1. Geographic gradients in population dynamics may occur because of spatial variation in resources that affect the deterministic components of the dynamics (i.e. carrying capacity, the specific growth rate at small densities or the strength of density regulation) or because of spatial variation in the effects of environmental stochasticity. To evaluate these, we used a hierarchical Bayesian approach to estimate parameters characterizing deterministic components and stochastic influences on population dynamics of eight species of ducks (mallard, northern pintail, blue-winged teal, gadwall, northern shoveler, American wigeon, canvasback and redhead (Anas platyrhynchos, A. acuta, A. discors, A. strepera, A. clypeata, A. americana, Aythya valisineria and Ay. americana, respectively) breeding in the North American prairies, and then tested whether these parameters varied latitudinally. 2. We also examined the influence of temporal variation in the availability of wetlands, spring temperature and winter precipitation on population dynamics to determine whether geographical gradients in population dynamics were related to large-scale variation in environmental effects. Population variability, as measured by the variance of the population fluctuations around the carrying capacity K, decreased with latitude for all species except canvasback. This decrease in population variability was caused by a combination of latitudinal gradients in the strength of density dependence, carrying capacity and process variance, for which details varied by species. 3. The effects of environmental covariates on population dynamics also varied latitudinally, particularly for mallard, northern pintail and northern shoveler. However, the proportion of the process variance explained by environmental covariates, with the exception of mallard, tended to be small. 4. Thus, geographical gradients in population dynamics of prairie ducks resulted from latitudinal gradients in both deterministic and stochastic components, and were likely influenced by spatial differences in the distribution of wetland types and shapes, agricultural practices and dispersal processes. 5. These results suggest that future management of these species could be improved by implementing harvest models that account explicitly for spatial variation in density effects and environmental stochasticity on population abundance.     

Within-group human variation in the Asian Pleistocene: the three Upper Cave crania

Numerous studies on Pleistocene samples have shown that within-group cranial variation was greater than that seen today. The three anatomically modern Upper Cave crania (UC 101, UC 102, and UC 103) from Zhoukoudian, China provide one of the best samples available for addressing the issue of the antiquity of the modern pattern of within-group cranial variability because archaeological evidence indicates that they are spatially and temporally restricted. Research on the Upper Cave fossils usually only includes UC 101 and UC 103 because of postmortem damage to UC 102's cranial vault. However, the face of UC 102 is undamaged, allowing for most facial measurements to be performed accurately. In this study we use facial dimensions to compare all three Upper Cave specimens, and we evaluate whether the variation seen among them is larger than that observed in extant populations.Using a worldwide sample of modern populations to establish a baseline, the three Upper Cave crania were compared to each other. Since there is disagreement over the sex of UC 102, this specimen is treated alternately as a female and as a male. Results show that the Upper Cave specimens exhibit significantly more variation than do individuals within more recent human populations, especially if UC 102 is considered male. Furthermore, results indicate that the fossils never fall into the same modern human group, and that each specimen is significantly atypical of its nearest modern neighbor in multivariate space.We conclude that the three Upper Cave crania do not represent a family group but are representative of the larger contemporaneous heterogeneous Asian Pleistocene population. Our results support the contention that today's within-group homogeneity is a relatively recent phenomenon, and is likely the result of a Neolithic population expansion and its many effects.     

Frog population viability under present and future climate conditions: a Bayesian state-space approach

1.?World-wide extinctions of amphibians are at the forefront of the biodiversity crisis, with climate change figuring prominently as a potential driver of continued amphibian decline. As in other taxa, changes in both the mean and variability of climate conditions may affect amphibian populations in complex, unpredictable ways. In western North America, climate models predict a reduced duration and extent of mountain snowpack and increased variability in precipitation, which may have consequences for amphibians inhabiting montane ecosystems. 2.?We used Bayesian capture-recapture methods to estimate survival and transition probabilities in a high-elevation population of the Columbia spotted frog (Rana luteiventris) over 10?years and related these rates to interannual variation in peak snowpack. Then, we forecasted frog population growth and viability under a range of scenarios with varying levels of change in mean and variance in snowpack. 3.?Over a range of future scenarios, changes in mean snowpack had a greater effect on viability than changes in the variance of snowpack, with forecasts largely predicting an increase in population viability. Population models based on snowpack during our study period predicted a declining population. 4.?Although mean conditions were more important for viability than variance, for a given mean snowpack depth, increases in variability could change a population from increasing to decreasing. Therefore, the influence of changing climate variability on populations should be accounted for in predictive models. The Bayesian modelling framework allows for the explicit characterization of uncertainty in parameter estimates and ecological forecasts, and thus provides a natural approach for examining relative contributions of mean and variability in climatic variables to population dynamics. 5.?Longevity and heterogeneous habitat may contribute to the potential for this amphibian species to be resilient to increased climatic variation, and shorter-lived species inhabiting homogeneous ecosystems may be more susceptible to increased variability in climate conditions.     

Genetic basis of between‐individual and within‐individual variance of docility

Between‐individual variation in phenotypes within a population is the basis of evolution. However, evolutionary and behavioural ecologists have mainly focused on estimating between‐individual variance in mean trait and neglected variation in within‐individual variance, or predictability of a trait. In fact, an important assumption of mixed‐effects models used to estimate between‐individual variance in mean traits is that within‐individual residual variance (predictability) is identical across individuals. Individual heterogeneity in the predictability of behaviours is a potentially important effect but rarely estimated and accounted for. We used 11 389 measures of docility behaviour from 1576 yellow‐bellied marmots (Marmota flaviventris) to estimate between‐individual variation in both mean docility and its predictability. We then implemented a double hierarchical animal model to decompose the variances of both mean trait and predictability into their environmental and genetic components. We found that individuals differed both in their docility and in their predictability of docility with a negative phenotypic covariance. We also found significant genetic variance for both mean docility and its predictability but no genetic covariance between the two. This analysis is one of the first to estimate the genetic basis of both mean trait and within‐individual variance in a wild population. Our results indicate that equal within‐individual variance should not be assumed. We demonstrate the evolutionary importance of the variation in the predictability of docility and illustrate potential bias in models ignoring variation in predictability. We conclude that the variability in the predictability of a trait should not be ignored, and present a coherent approach for its quantification.     

Developmental variability during early embryonic development of zebra fish, Danio rerio

Early vertebrate embryos pass through a period of remarkable morphological similarity. Possible causes for such similarity of early embryos include modularity, developmental constraints, stabilizing selection, canalization, and exhausted genetic variability. Supposedly, each process creates different patterns of variation and covariation of embryonic traits. We study the patterns of variation of the embryonic phenotype to test ideas about possible evolutionary mechanisms shaping the early embryonic development. We use the zebra fish, Danio rerio, as a model organism and apply repeated measures of individual embryos to study temporal changes of phenotypic variability during development. In particular, we are looking at the embryonic development from 12 hours post fertilization until 27 hours post fertilization. During this time period, the development of individual embryos is documented at hourly intervals. We measured maximum diameter of the eye, length of embryo, number of somites, inclination of somites, and the yolk size (as a maternal effect). The coefficient of variation (CV) was used as a measure of variability that was independent of size. We used a principal component analysis for analysis of morphological integration. The experimental setup kept environment x genotype interactions constant. Nongenetic parental contributions had no significant effects on interindividual variability. Thus all observed phenotypic variation was based on additive genetic variance and error variance. The average CV declined from 14% to 7.7%. The decline of the CV was in particular expressed during 15-19 h post fertilization and occurred in association with multiple correlations among embryonic traits and a relatively high degree of morphological integration. We suggest that internal constraints determine the patterns of variability during early embryonic development of zebra fish.