Rate of evolutionary change in cranial morphology of the marsupial genus Monodelphis is constrained by the availability of additive genetic variation

We tested the hypothesis that the rate of marsupial cranial evolution is dependent on the distribution of genetic variation in multivariate space. To do so, we carried out a genetic analysis of cranial morphological variation in laboratory strains of Monodelphis domestica and used estimates of genetic covariation to analyse the morphological diversification of the Monodelphis brevicaudata species group. We found that within-species genetic variation is concentrated in only a few axes of the morphospace and that this strong genetic covariation influenced the rate of morphological diversification of the brevicaudata group, with between-species divergence occurring fastest when occurring along the genetic line of least resistance. Accounting for the geometric distribution of genetic variation also increased our ability to detect the selective regimen underlying species diversification, with several instances of selection only being detected when genetic covariances were taken into account. Therefore, this work directly links patterns of genetic covariation among traits to macroevolutionary patterns of morphological divergence. Our findings also suggest that the limited distribution of Monodelphis species in morphospace is the result of a complex interplay between the limited dimensionality of available genetic variation and strong stabilizing selection along two major axes of genetic variation.

Relationship of age and body mass index to the expression of obesity and osteoarthritis-related genes in human meniscus

OBJECTIVE

Aging and obesity contribute to the initiation and progression of osteoarthritis with little information on their relation to gene expression in joint tissues, particularly the meniscus. Here, we test the hypothesis that patient age and body mass index (BMI) correlate with the expression of osteoarthritis- and obesity-related gene signatures in the meniscus.

DESIGN

Meniscus was obtained from patients (N=68) undergoing arthroscopic partial meniscectomy. The mRNA expression of 24 osteoarthritis-related and 4 obesity-related genes in meniscus was assessed by quantitative real-time PCR. The relationship between gene expression and patient age and BMI was analyzed using Spearman's rank-order correlation. Hierarchical cluster dendrogram and heat map were generated to study inter-gene associations.

RESULTS

Age was negatively correlated (P<0.05) with the expression of MMP-1 (r=-0.447), NFκB2 (r=-0.361), NFκBIA (r=-0.312), IκBA (r=-0.308), IL-8 (r=-0.305), ADAMTS-4 (r=-0.294), APLN (apelin) (r=-0.250) and IL-6 (r=-0.244). Similarly, BMI was negatively correlated with the expression of APLN (r=-0.328), ACAN (r=-0.268) and MMP-1 (r=-0.261). After adjusting for the correlation between age and BMI (r=0.310; P=0.008), the only independent effect of BMI on gene expression was for APLN (r=-0.272). However, age had an independent effect on the expression on ADAMTS-4 (r=-0.253), MMP-1 (r=-0.399), IL-8 (r=-0.327), COL1A1 (r=-0.287), NFκBIA (r=-0.278), NFκB2 (r=-0.312) and IκBA (r=-0.299). The gene correlation analysis identified four clusters of potentially relevant genes: transcription factors, matrix-degrading enzymes, cytokines and chemokines, and obesity genes.

CONCLUSION

Age and BMI were negatively correlated with several osteoarthritis- and obesity-related genes. Although the bulk of these changes appeared to be driven by age, expression of APLN was related to BMI. Inter-gene correlation analysis implicated a common role for strongly correlated genes. Although age-related variations in gene expression appear to be more relevant than obesity-related differences for the role of the meniscus in osteoarthritis development, further investigation into the role of APLN in meniscus and joint health is warranted.

Healing quantitative trait loci in a combined cross analysis using related mouse strain crosses

Inbred mouse strains MRL and LG share the ability to fully heal ear hole punches with the full range of appropriate tissues without scarring. They also share a common ancestry, MRL being formed from a multi-strain cross with two final backcrosses to LG before being inbred by brother-sister mating. Many gene-mapping studies for healing ability have been performed using these two strains, resulting in the location of about 20 quantitative trait loci (QTLs). Here, we combine two of these crosses (N = 638), MRL/lpr × C57BL/6NTac and LG/J × SM/J, in a single combined cross analysis to increase the mapping power, decrease QTL support intervals, separate multiple QTLs and establish allelic states at individual QTL. The combined cross analysis located 11 QTLs, 6 affecting only one cross (5 LG × SM and 1 MRL × B6) and 5 affecting both crosses, approximately the number of common QTLs expected given strain SNP similarity. Amongst the five QTLs mapped in both crosses, three had significantly different genetic effects, additive in one cross and over or underdominant in the other. It is possible that allelic states at these three loci are different in SM and B6 because they lead to differences in dominance interactions with the LG and MRL alleles. QTL support intervals are 40% smaller in the combined cross analysis than in either of the single crosses. Combined cross analysis was successful in enhancing the interpretation of earlier QTL results for these strains.

Genotype-dependent responses to levels of sibling competition over maternal resources in mice

Research on phenotypic plasticity has often focused on how a given genotype responds to the changing physical environments such as temperature or diet. However, for many species the social environment has an equally important role because of competition for resources. During early development, the level of competition for limited (maternally provided) resources will often depend critically on the number of siblings. Therefore, competition among siblings should drive the evolution of genes that allow flexible responses to realized levels of competition and maternal resource availability. However, it is unknown whether genetically based differences between individuals exist in their response to the social environment that affect their future development. Using a quantitative trait locus approach in an experimental population of mice we demonstrate that effects of sibling number on body weight depend on individual genotype at seven loci, over and above the general negative litter size effect. Overall, these litter size-by-genotype interactions considerably modified the degree to which increasing litter size caused reduced weight. For example at one locus this effect leads to a 7% difference in body weight at week 7 between individuals experiencing the extremes of the normal range of litter sizes in our population (five to nine litter mates). The observed interaction between genotype and the competitive environment can produce differences in body weight that are similar in magnitude to the main effect of litter size on weight. Our results show that different genotypes respond to the social environment differentially and that interaction effects of genotype with litter size can be as important as genotype-independent effects of litter size.

The effect of a population bottleneck on the evolution of genetic variance/covariance structure

It is well known that standard population genetic theory predicts decreased additive genetic variance (V(a) ) following a population bottleneck and that theoretical models including interallelic and intergenic interactions indicate such loss may be avoided. However, few empirical data from multicellular model systems are available, especially regarding variance/covariance (V/CV) relationships. Here, we compare the V/CV structure of seventeen traits related to body size and composition between control (60 mating pairs/generation) and bottlenecked (2 mating pairs/generation; average F = 0.39) strains of mice. Although results for individual traits vary considerably, multivariate analysis indicates that V(a) in the bottlenecked populations is greater than expected. Traits with patterns and amounts of epistasis predictive of enhanced V(a) also show the largest deviations from additive expectations. Finally, the correlation structure of weekly weights is not significantly different between control and experimental lines but correlations between necropsy traits do differ, especially those involving the heart, kidney and tail length.

Neuroanatomic and behavioral traits for autistic disorders in age-specific restricted index selection mice

The pathogenesis of neurodevelopmental disorders such as autism is believed to be influenced by interactions between genetic and environmental factors, and appropriate animal models are needed to assess the influence of such factors on relevant neurodevelopmental phenotypes. A set of inbred mouse strains (Atchley strains) including A12 (E+L0) and A22 (E-L0) were generated by age-specific restricted index selection from a baseline random-bred ICR mouse population obtained from Harlan Sprague-Dawley [Atchley et al. (1997) Genetics 146(2):629-640; Indianapolis, IN, USA). As compared with the A22 strain, A12 mice had significantly increased early (P0-P10) body weight gain with minimal changes in late (P28-P56) body weight gain. We found that these strains also differed in brain weight, brain volume, cell proliferation, and FGF-2 levels in certain brain regions. Specifically, brain weight and volume were significantly greater in A12 mice than that in A22 mice at P10 and P28. Quantitative analysis of bromodeoxyuridine (BrdU) labeling of proliferating cells showed that the number of BrdU-positive cells in the A12 strain were significantly greater in the frontal cortex and lesser in the dentate gyrus than that in the A22 strain at P28. Western blot revealed that fibroblast growth factors-2 (FGF-2), but not brain-derived neurotrophic factor (BDNF), expression was significantly increased in the frontal cortex of A12 strain at P28. Also, A12 mice exhibited decreased intra-strain social interaction and increased repetitive stereotyped behaviors at P28. Our study suggests that A12 mice may partially mimic the anatomic and behavioral traits of patients with neurodevelopmental disorders such as autism spectrum disorders, and therefore may yield insights into the developmental mechanisms involved in their pathogenesis.

Reply to Lawler: feeding competition, cooperation, and the causes of primate sociality

This is a reply to Richard Lawler's commentary on our previous work [Lawler, 2011; this issue] in which he develops a set of operational models to test socioecological theories of the evolutionary importance of feeding competition. We strongly agree that we need to critically re-evaluate the basic assumptions of all models of primate sociality, and to verify the explanatory power of alternative models. We also feel Lawler's commentary provides an important opportunity to broaden the debate concerning the fundamental roles of cooperation, competition, and aggression in understanding primate social systems. Lawler provides a number of suggestions as to how models developed in primate socioecology might be tested. We agree with these suggestions, make further suggestions, and call for specific operational definitions so that researchers might begin to develop and test various methodologies. However, we also call for testing alternative theories. Current socioecological theory is based on the assumption that competition and positive selection is always in operation and has driven the evolution of living organisms. We believe that this "explanation of choice" often is treated as an assumed truth to which data are forced to fit, rather than being seen as a theory to be tested. Furthermore, we agree with Weiss and Buchanan [2009. The Mermaid's Tale: Four Billion Years of Cooperation in the Making of Living Things] that on ecological and developmental scales, where organisms actually live out their lives, cooperation may play a more fundamental role than competition.

Sex-specific quantitative trait loci linked to autoresuscitation failure in SWR/J mice

Autoresuscitation (AR) is a highly conserved response among mammals, which allows survival from transient extreme hypoxia. During hypoxia, bradycardia, and hypoxic gasping develop after a brief period of hyperactivity. Normally, AR occurs if oxygen is restored during the gasping period where an initial heart rate increase is rapidly followed resumption or eupneic breathing. Humans and other mammals can survive multiple immediately repeated AR. A defective AR capacity has been implicated in Sudden Infant Death Syndrome. We had reported earlier that inbred strains of mice such as BALB/cJ could survive a characteristic number of immediately repeated AR trials, but that SWR/J mice failed to AR from a single hypoxic episode. We now report that strains closely related to SWR/J, FVB/N and SJL/J exhibit partial resuscitation defects relative to BALB/cJ or other mouse strains, establishing a genetic basis for variation in AR failure. The AR trial phenotype of BALB/cJ x SWR/J intercross F(1) and F(2) mice was consistent with BALB/cJ dominance and a discrete number of loci. Genome-wide mapping conducted with 60 intercross F(2) animals linked two loci to the number of AR trials survived, including one sex-specific locus with male expression, consistent with the observed 50% male bias for Sudden Infant Death Syndrome in humans. A locus carried on SWR/J chromosome 10 seems to be particularly important in AR failure and was confirmed in a partial consomic line. These results establish a genetic basis for AR failure phenotype in mice, with relevance to Sudden Infant Death Syndrome.

Epistasis and the evolutionary dynamics of measured genotypic values during simulated serial bottlenecks

The evolutionary effects of epistasis have been primarily explored analytically and most empirical studies have utilized yeast, viral and bacterial populations. Empirical analyses in multi-cellular organisms are rare because of experimental constraints. Here, we report the results of a genome-wide scan for two-way epistasis in 16 traits related to body size and composition in F(2) mice from the LG/J by SM/J intercross. We analyze two-locus genotypic values at quantitative trait loci (QTL), which provides an especially detailed view of epistatic architectures, to evaluate their predicted evolutionary consequences via Monte Carlo simulations. Epistatic profiles vary, but all traits show complicated genetic architectures which are largely hidden in single locus QTL scans. On average, detected epistatic effects are comparable in size to marginal effects. Simulations demonstrate an expected preservation, and often inflation, of heritable variance across several generations of small effective population size for many identified epistatic pairs over a range of starting allele frequencies.

A search for quantitative trait loci exhibiting imprinting effects on mouse mandible size and shape

Genomic imprinting refers to the pattern of monoallelic parent-of-origin-dependent gene expression where one of the two alleles at a locus is expressed and the other silenced. Although some genes in mice are known to be imprinted, the true scope of imprinting and its impact on the genetic architecture of a wide range of morphometric traits is mostly unknown. We therefore searched for quantitative trait loci (QTL) exhibiting imprinting effects on mandible size and shape traits in a large F(3) population of mice originating from an intercross of the LG/J (Large) and SM/J (Small) inbred strains. We discovered a total of 51 QTL affecting mandible size and shape, 6 of which exhibited differences between reciprocal heterozygotes, the usual signature of imprinting effects. However, our analysis showed that only one of these QTL (affecting mandible size) exhibited a pattern consistent with true imprinting effects, whereas reciprocal heterozygote differences in the other five all were due to maternal genetic effects. We concluded that genomic imprinting has a negligible effect on these specific morphometric traits, and that maternal genetic effects may account for many of the previously reported instances of apparent genomic imprinting.

Quantitative trait loci linked to thalamus and cortex gray matter volumes in BXD recombinant inbred mice

To investigate whether there are separate or shared genetic influences on the development of the thalamus and cerebral cortex, we identified quantitative trait loci (QTLs) for relevant structural volumes in BXD recombinant inbred (RI) strains of mice. In 34 BXD RI strains and two parental strains (C57BL/6J and DBA/2J), we measured the volumes of the entire thalamus and cortex gray matter using point counting and Cavalieri's rule. Heritability was calculated using analysis of variance (ANOVA), and QTL analysis was carried out using WebQTL (http://www.genenetwork.org). The heritability of thalamus volume was 36%, and three suggestive QTLs for thalamus volume were identified on chromosomes 10, 11 and 16. The heritability of cortical gray matter was 43%, and four suggestive QTLs for cortex gray matter volume were identified on chromosomes 2, 8, 16 and 19. The genetic correlation between thalamus and cortex gray matter volumes was 0.64. Also, a single QTL on chromosome 16 (D16Mit100) was identified for thalamus volume, cortex gray matter volume and Morris water maze search-time preference (r=0.71). These results suggest that there are separate and shared genetic influences on the development of the thalamus and cerebral cortex.

Genetic architecture of Arabidopsis thaliana response to infection by Pseudomonas syringae

Plant pathogens can severely reduce host yield and fitness. Thus, investigating the genetic basis of plant response to pathogens is important to further understand plant-pathogen coevolution and to improve crop production. The interaction between Arabidopsis thaliana and Pseudomonas syringae is an important model for studying the genetic basis of plant-pathogen interactions. Studies in this model have led to the discovery of many genes that differentiate a resistant from a susceptible plant. However, little is known about the genetic basis of quantitative variation in response to P. syringae. In this study, we investigate the genetic basis of three aspects of A. thaliana's response to P. syringae: symptom severity, bacterial population size and fruit production using a quantitative trait loci (QTL) analysis. We found two QTL for symptom severity and two for fruit production (possible candidate genes for observed QTL are discussed). We also found significant two-locus epistatic effect on symptom severity and fruit production. Although bacterial population size and symptom severity were strongly phenotypically correlated, we did not detect any QTL for bacterial population size. Despite the detected genetic variation observed for susceptibility, we found only a weak overall relationship between susceptibility traits and fitness, suggesting that these traits may not respond to selection.

An epistatic genetic basis for fluctuating asymmetry of tooth size and shape in mice

Although there typically is little additive genetic variation for fluctuating asymmetry (FA), or variation in nondirectional differences between left and right sides of bilateral characters, several investigators have hypothesized that FA may have an epistatic genetic basis. We tested this hypothesis by conducting a whole genome scan of FA of size and shape of the mandibular molars in house mice from an F2 intercross population generated from crossing the Large (LG/J) and Small (SM/J) inbred strains. Although no individual genes (QTLs=quantitative trait loci) on any of the 19 autosomes significantly affected FA for centroid size, and only two affected shape FA, a number of pairwise combinations of QTLs exhibited significant epistasis for FA in both molar size and shape. The QTLs involved in these interactions differed for FA in molar size versus FA in molar shape, but their epistatic contributions to the total variance was nearly the same (about 20%) for FA in both molar characters. It was noted that the genetic architecture of FA in the molar characters, consisting of little or no additive genetic variance but an abundance of epistatic genetic variance, is consistent with that of other typical fitness components such as litter size.

Epistasis affecting litter size in mice

Litter size is an important reproductive trait as it makes a major contribution to fitness. Generally, traits closely related to fitness show low heritability perhaps because of the corrosive effects of directional natural selection on the additive genetic variance. Nonetheless, low heritability does not imply, necessarily, a complete absence of genetic variation because genetic interactions (epistasis and dominance) contribute to variation in traits displaying strong heterosis in crosses, such as litter size. In our study, we investigated the genetic architecture of litter size in 166 females from an F2 intercross of the SM/J and LG/J inbred mouse strains. Litter size had a low heritability (h2 = 12%) and a low repeatability (r = 33%). Using interval-mapping methods, we located two quantitative trait loci (QTL) affecting litter size at locations D7Mit21 + 0 cM and D12Mit6 + 8 cM, on chromosomes 7 and 12 respectively. These QTL accounted for 12.6% of the variance in litter size. In a two-way genome-wide epistasis scan we found eight QTL interacting epistatically involving chromosomes 2, 4, 5, 11, 14, 15 and 18. Taken together, the QTL and their interactions explain nearly 49% (39.5% adjusted multiple r2) of the phenotypic variation for litter size in this cross, an increase of 36% over the direct effects of the QTL. This indicates the importance of epistasis as a component of the genetic architecture of litter size and fitness in our intercross population.

Cranial evolution in sakis (Pithecia, Platyrrhini) II: evolutionary processes and morphological integration

Patterns of interspecific differentiation in saki monkeys (Pithecia) are quantitatively described and possible evolutionary processes producing them are examined. The comparison of species correlation matrices to expected patterns of morphological integration reveal significant and similar patterns of development-based cranial integration among species. Aspects of the facial region are more heavily influenced by general size variation than features of the neural region. The comparison of pooled within- and between-groups V/CV matrices suggests that genetic drift might be a sufficient explanation for saki cranial evolution. Differential natural selection gradients are also reconstructed because selection may also have caused population differentiation through evolutionary time. These gradients illustrate the inherent multivariate nature of selection, being a consequence of the interaction between existing morphological integration (correlation) among traits and the action of natural selection. Yet, our attempt to interpret selection gradients in terms of their functional significance did not result in any clear association between selection and function. Perhaps this is also an indication that morphological evolution in sakis was mostly neutral.

Contribution of maternal effect QTL to genetic architecture of early growth in mice

Existing approaches to characterizing quantitative trait loci (QTL) utilize a paradigm explicitly focused on the direct effects of genes, where phenotypic variation among individuals is mapped onto genetic variation of those individuals. For many characters, however, the genotype of the mother via its maternal effect accounts for a considerable portion of the genetically based variation in progeny phenotypes. Thus the focus on direct effect QTL may result in an insufficient or misleading characterization of genetic architecture due to the omission of the potentially important source of genetic variance contributed by maternal effects. We analyze the relative contribution of direct and maternal effect (ME) QTL to early growth in mice using a three-generation intercross of the Small (SM/J) and Large (LG/J) inbred mouse lineages. Using interval mapping and composite interval mapping, direct effect (DE) QTL for early growth (change in body mass during the interval from week 1 to 2) were detected in the F(2) generation of the intercross (n = 510), where no maternal genetic effect variance is present (all individuals are progeny of genetically identical F(1) mothers). ME QTL were detected by treating the phenotypes of cross-fostered F(3) pups as a characteristic of their nurse-dam (n = 168 dams with cross-fostered progeny). Five DE QTL, significant at a chromosome wide level (alpha = 0.05), were detected, with two significant at a genome wide level. FourME QTL significant at the chromosome wide level were detected, with three significant at the genome wide level. A model containing only DE QTL accounted for 11.8% of phenotypic variance, while a model containing only ME QTL accounted for 31.5% of the among litter variance in growth. There was no evidence for pleiotropy of DE and ME loci since there was no overlap between loci detected in these two analyses. Epistasis between all pairs of loci was analyzed for both DEs and MEs. Ten pairs of loci showed significant epistasis for MEs (alpha = 0.05 corrected for multiple comparisons) while four pairs showed significant epistasis for DEs on early growth.

A comparison of phenotypic variation and covariation patterns and the role of phylogeny, ecology, and ontogeny during cranial evolution of new world monkeys

Similarity of genetic and phenotypic variation patterns among populations is important for making quantitative inferences about past evolutionary forces acting to differentiate populations and for evaluating the evolution of relationships among traits in response to new functional and developmental relationships. Here, phenotypic co variance and correlation structure is compared among Platyrrhine Neotropical primates. Comparisons range from among species within a genus to the superfamily level. Matrix correlation followed by Mantel's test and vector correlation among responses to random natural selection vectors (random skewers) were used to compare correlation and variance/covariance matrices of 39 skull traits. Sampling errors involved in matrix estimates were taken into account in comparisons using matrix repeatability to set upper limits for each pairwise comparison. Results indicate that covariance structure is not strictly constant but that the amount of variance pattern divergence observed among taxa is generally low and not associated with taxonomic distance. Specific instances of divergence are identified. There is no correlation between the amount of divergence in covariance patterns among the 16 genera and their phylogenetic distance derived from a conjoint analysis of four already published nuclear gene datasets. In contrast, there is a significant correlation between phylogenetic distance and morphological distance (Mahalanobis distance among genus centroids). This result indicates that while the phenotypic means were evolving during the last 30 millions years of New World monkey evolution, phenotypic covariance structures of Neotropical primate skulls have remained relatively consistent. Neotropical primates can be divided into four major groups based on their feeding habits (fruit-leaves, seed-fruits, insect-fruits, and gum-insect-fruits). Differences in phenotypic covariance structure are correlated with differences in feeding habits, indicating that to some extent changes in interrelationships among skull traits are associated with changes in feeding habits. Finally, common patterns and levels of morphological integration are found among Platyrrhine primates, suggesting that functional/developmental integration could be one major factor keeping covariance structure relatively stable during evolutionary diversification of South American monkeys.

A simple correction for multiple comparisons in interval mapping genome scans

Several approaches have been proposed to correct point-wise significance thresholds used in interval-mapping genome scans. A method for significance threshold correction based on the Bonferroni test is presented. This test involves calculating the effective number of independent comparisons performed in a genome scan from the variance of the eigenvalues of the observed marker correlation matrix. The more highly correlated the markers, the higher the variance of the eigenvalues and the lower the number of independent tests performed on a chromosome. This approach was evaluated by mapping 1000 normally distributed phenotypes along chromosomes of varying length and marker density in a population size of 500. Experiment-wise significance thresholds obtained from the simulation are compared to those calculated using the Bonferroni criterion and the newly developed measure of the effective number of independent tests in a genome scan. The Bonferroni calculation produced significance thresholds very similar to those obtained by simulation. The threshold levels for both Bonferroni and simulation analysis depended strongly on the marker density and size of chromosomes. There was a slight bias of about 1% in the thresholds obtained at the 5% and 10% point-wise significance levels. The method introduced here provides a relatively simple correction for multiple comparisons that can be easily calculated using standard statistics packages.

Subspecific genetic differences in the saddle-back tamarin (Saguinus fuscicollis) postcranial skeleton

Saddle-back tamarins (Saguinus fuscicollis) have been differentiated into subspecies based on coat color and facial morphology. Morphology of first-generation genetic hybrids can provide an understanding of population differences and evolution that may otherwise be difficult to obtain. This study compares postcranial morphology in two hybrids of saddle-back tamarins (S. f. illigeri x S. f. lagonotus and S. f. illigeri x S. f. leucogenys) to the purebred subspecies. We examine whether postcranial morphology differs between the purebred subspecies and their hybrids, and whether hybrids exhibit heterosis indicative of genetic differences between the parental subspecies. Fifty postcranial dimensions, or traits, were measured on 287 S. fuscicollis skeletons. Measurements were corrected for sex differences and environmental differences (wild-derived or laboratory-born). Females were larger than males for most postcranial dimensions, especially in the pelvis. Significant environmental differences were found throughout the postcranial skeleton, with wild-derived animals generally larger than animals born in the laboratory. Significant heterosis was found for six postcranial traits in the S. f. illigeri x S. f. leucogenys cross, and heterosis was found for 15 postcranial traits in the S. f. illigeri x S. f. lagonotus cross. These results suggest that subspecies differences in skeletal morphology are generally due to genetic rather than environmental differences. The subspecies appear to represent geographic isolates that differ in gene frequencies and directional dominance for skeletal traits.

Genetic architecture of mandible shape in mice: effects of quantitative trait loci analyzed by geometric morphometrics

This study introduces a new multivariate approach for analyzing the effects of quantitative trait loci (QTL) on shape and demonstrates this method for the mouse mandible. We quantified size and shape with the methods of geometric morphometrics, based on Procrustes superimposition of five morphological landmarks recorded on each mandible. Interval mapping for F(2) mice originating from an intercross of the LG/J and SM/J inbred strains revealed 12 QTL for size, 25 QTL for shape, and 5 QTL for left-right asymmetry. Multivariate ordination of QTL effects by principal component analysis identified two recurrent features of shape variation, which involved the positions of the coronoid and angular processes relative to each other and to the rest of the mandible. These patterns are reminiscent of the knockout phenotypes of a number of genes involved in mandible development, although only a few of these are possible candidates for QTL in our study. The variation of shape effects among the QTL showed no evidence of clustering into distinct groups, as would be expected from theories of morphological integration. Further, for most QTL, additive and dominance effects on shape were markedly different, implying overdominance for specific features of shape. We conclude that geometric morphometrics offers a promising new approach to address problems at the interface of evolutionary and developmental genetics.

Genetic architecture of adiposity in the cross of LG/J and SM/J inbred mice

The genetic basis of variation in obesity in human populations is thought to be owing to many genes of relatively small effect and their interactions. The LG/J by SM/J intercross of mouse inbred strains provides an excellent model system in which to investigate multigenic obesity. We previously mapped a large number of quantitative trait loci (QTLs) affecting adult body weight in this cross. We map body composition traits, adiposity, and skeletal size, in a replicate F2 intercross of the same two strains containing 510 individuals. Using interval-mapping methods, we located eight QTLs affecting adiposity (Adip1-8). Two of these adiposity loci also affected tail length (Adip4 and Adip6) along with seven additional tail length QTLs (Skl1-7). A further four QTLs (Wt1-4) affect adult weight but not body composition. These QTLs have relatively small effects, typically about 0.2-0.4 standard deviation units, and account for between 3% and 10% of the variance in individual characters. All QTLs participated in epistatic interactions with other QTLs. Most of these interactions were due to additive-by-additive epistasis, which can nullify the apparent effects of single loci in our population. Adip8 interacts with all the other adiposity QTLs and seems to play a central role in the genetic system affecting obesity in this cross. Only two adiposity QTLs, Adip4 and Adip6, also affect tail length, indicating largely separate genetic control of variation in adiposity and skeletal size. Body size and obesity QTLs in the same locations as those discovered here are commonly found in mapping experiments with other mouse strains.

Is the genotype-phenotype map modular? A statistical approach using mouse quantitative trait loci data

Various theories about the evolution of complex characters make predictions about the statistical distribution of genetic effects on phenotypic characters, also called the genotype-phenotype map. With the advent of QTL technology, data about these distributions are becoming available. In this article, we propose simple tests for the prediction that functionally integrated characters have a modular genotype-phenotype map. The test is applied to QTL data on the mouse mandible. The results provide statistical support for the notion that the ascending ramus region of the mandible is modularized. A data set comprising the effects of QTL on a more extensive portion of the phenotype is required to determine if the alveolar region of the mandible is also modularized.

Interactions between infant growth and survival: evidence for selection on age-specific body weight in captive common marmosets (Callithrix jacchus)

The objective of this study is to investigate factors influencing infant survival in captive common marmosets. We investigated the influence of age-specific weight, litter size, caging, and the presence of helpers on survival to 6 months of age in 189 Callithrix jacchus infants. Infant survival was analyzed using Cox Proportional Hazards regression, and fitness functions were plotted to explore the relationship between survival and growth. Results indicate that weights at birth and 120 days significantly affect future survival probability. Litter size significantly influences survival prior to 60 days of age with larger litters having poorer survival. Males and females did not have significantly different survival and the presence of helpers in the group did not influence survival probability. Patterns of survival with respect to age-specific weights suggest stabilizing selection on birth weight and directional selection on weight at 120 days of age.

Age of epiphyseal closure in tamarins and marmosets

Estimates of the chronological age for animals of unknown age provide useful information for medical, demographic, and evolutionary studies. Skeletal development, as indicated by epiphyseal closure, can be used to estimate an animal's chronological age or specify its stage of development. Many studies of Primate skeletal development have used animals of unknown age, with the order of epiphyseal closure providing a relative age for each animal. This study examines the age of epiphyseal closure at 22 epiphyseal sites using animals of known age at death in three calitrichid species (Saguinus fuscicollis, Saguinus oedipus, and Callithrix jacchus). The observed average age of epiphyseal closure is similar in these tamarins and marmosets. There is a significant difference in rate of development between the species. Regression equations can predict the age of unknown animals to within 4.8 months for S. fuscicollis, 8.6 months for S. oedipus, and 7.6 months for C. jacchus (twice the standard error of the estimate). These age estimates allow us to determine if an animal is relatively mature or immature, but are largely unacceptable for studies in which precise age estimates are necessary. The order of epiphyseal closure is similar across 11 monkey species (using additional data from published literature) and supports the suggestion of a general pattern in Primate skeletal development.

Quantitative trait loci for directional but not fluctuating asymmetry of mandible characters in mice

Non-directional variation in right minus left differences in bilateral characters, referred to as fluctuating asymmetry (FA), often has been assumed to be largely or entirely environmental in origin. FA increasingly has been used as a measure of developmental stability, and its presumed environmental origin has facilitated the comparisons of populations believed to differ in their levels of stability. Directional asymmetry (DA), in which one side is consistently larger than the other, has been assumed to be at least partially heritable. Both these assumptions were tested with interval mapping techniques designed to detect any quantitative trait loci (QTLs) affecting FA or DA in 15 bilateral mandible characters in house mice resulting from a cross of the F1 between CAST/Ei (wild strain) and M16i (selected for rapid growth rate) back to M16i. For purposes of the analysis, all mandibles were triply measured and 92 microsatellite markers were scored in a total of 350 mice. No significant QTLs were found for FA, but three QTLs significantly affected DA in several characters, confirming both assumptions. The QTLs for DA were similar in location to those affecting the size of several of the mandible characters, although they accounted for an average of only 1% of the total phenotypic variation in DA.

Phenotypic covariance structure in tamarins (genus Saguinus): a comparison of variation patterns using matrix correlation and common principal component analysis

Constancy of variation/covariation structure among populations is frequently assumed in order to measure the differential selective forces which have caused population differentiation through evolutionary time. Following Steppan ([1997] Evolution 51:571-594), this assumption is examined among closely related tamarin species (genus Saguinus), using two distinct approaches applied to the task of evaluating similarity in patterns of morphological variation: common principal component analysis and matrix correlations. While the results of these analyses may appear contradictory, closer examination reveals them as complementary, highlighting the wisdom of combined methodologies. Overall, the results reveal a close relationship among the morphologically based variance structures of the tamarin species a relationship whose pattern is consistent with the pattern of phylogenetic relatedness as found via a molecular genetic study. More specifically, both methodological approaches provide some support for divergence of S. geoffroyi and S. oedipus (with regards to their patterns of morphological variation) from other tamarin species. This suggests that variance/covariance structure may have diverged through evolutionary time in the tamarin lineage, placing assumptions of constancy in doubt.

Mapping quantitative trait loci for murine growth: a closer look at genetic architecture

Over 20 years ago, D. S. Falconer and others launched an important avenue of research into the quantitative of body size growth in mice. This study continues in that tradition by locating quantitative trait loci (QTLs) responsible for murine growth, such as age-specific weights and growth periods, and examining the genetic architecture for body weight. We identified a large number of potential QTLs in an earlier F2 intercross (Intercross I) of the SM/J and LG/J inbred mouse strains. Many of these QTLs are replicated in a second F2 intercross (Intercross II) between the same two strains. These replicated regions provide candidate regions for future fine-mapping studies. We also examined body size and growth QTLs using the combined data set from these two intercrosses, resulting in 96 microsatellite markers being scored for 1045 individuals. An examination of the genetic architecture for age-specific weight and growth periods resulted in locating 20 separate QTLs, which were mainly additive in nature, although dominance was found to affect early growth and body size. QTLs affecting early and late growth were generally distinct, mapping to separate chromosome locations. This QTL pattern indicates largely separate genetic and physiological systems for early and later murine growth, as Falconer suggested. We also found sex-specific QTLs for body size with implications for the evolution of sexual dimorphism.

Differential response to dietary fat in large (LG/J) and small (SM/J) inbred mouse strains

The "large" (LG/J) and "small" (SM/J) inbred mouse strains differ for a wide variety of traits related to body size and obesity. Ninety-three LG/J and SM/J mice were divided into two treatment categories and fed a moderately high-fat diet (21% kcal fat) or a low-fat diet (12% kcal fat) from weaning to necropsy. Strain differences in obesity-related traits and differential response to dietary fat increases were analyzed using ANOVA. LG/J animals grow faster from 3 to 10 wk, have longer tails, and have heavier body weight, liver weight, and fat pad weight than SM/J animals. SM/J animals grow faster after 10 wk of age and have higher fasting glucose levels than LG/J animals. SM/J mice were more responsive to increased dietary fat than LG/J mice for growth after 10 wk, necropsy weight, liver weight, fat pad weights, and fasting glucose levels (in males). The growth from 3 to 10 wk had a much greater response in the LG/J strain, whereas tail length had no response. This pattern of dietary response is similar to that expected under the "thrifty" phenotype hypothesis. Genes affecting strain differences and the differential response of the strains to dietary fat can be successfully mapped in the intercross of the LG/J and SM/J strains. This intercross provides an excellent multigenic model for the genetic basis of complex traits and diseases related to body size and obesity.

Historical biogeography of tamarins, genus Saguinus: the molecular phylogenetic evidence

Hypotheses of the historical biogeography of tamarins (genus Saguinus) based on variation in coat colors and body size are tested using phylogenetic relationships inferred from mitochondrial DNA (mtDNA) sequence data. Samples from all 12 species of Saguinus and several subspecies are included in the analysis. Approximately 1,200 bases of mtDNA sequence from the cytochrome b and D-loop regions are reported for the tamarins and several outgroup taxa. Parsimony analysis of the mtDNA sequence data reveals Saguinus to be a monophyletic taxon composed of two major clades: one, the Small-bodied clade, contains S. nigricollis, S. tripartitus, and S. fuscicollis, and the other, the Large-bodied clade, contains the other nine species. The phylogenetic relationships among tamarins inferred from the mtDNA sequence data reject previous hypotheses for the historical biogeography of tamarins and suggest different dispersal routes for this group of New World monkeys. The molecular data suggest that tamarins dispersed across South America in two major waves from an origin somewhere south of the Amazon. One wave moved in a westerly direction, whereas the other moved in a northeastern direction toward the Amazon delta and then west along the northern portion of the continent into northern Colombia and Panama.

Quantitative trait loci for fluctuating asymmetry of discrete skeletal characters in mice

Levels of fluctuating asymmetry (FA) are often taken as indicators of the degree to which genotypes differ in their ability to buffer genetic and environmental sources of variation. Interval mapping techniques were used to search for quantitative trait loci (QTLs) affecting directional asymmetry (DA) and fluctuating asymmetry (FA) in six bilateral discrete skeletal traits in house mice. These six characters as well as 76 microsatellite markers were scored in over 500 mice that resulted from crosses of F1 mice originally produced from matings of the Large (LG/J) and Small (SM/J) inbred strains. The number of QTLs affecting DA in each of the characters was no more than expected by chance alone so it was concluded that there was little evidence for individual genes affecting DA. There appeared to be a genetical basis for FA in these characters, however, because the number of QTLs significantly affecting FA (10 at the 5% level, three at the 1% level) was greater than expected by chance alone. The 10 QTLs significantly affecting FA in any given character were located on eight different chromosomes, mostly at locations for QTLs affecting other characters or DA in other characters. Their cumulative contribution to the total phenotypic variance was small, averaging only 3.9% per locus. Dominance genotypic values for these QTLs were more extreme than additive genotypic values, suggesting that heterozygotes at many loci are better buffered than homozygotes and that allelic interactions (dominance) may play an important role in the production of FA.

Precision of surface measurements for below-knee residua

OBJECTIVE

To determine the absolute and relative precision of geometric measurements made of below knee (BK) residua and their BK plaster positive casts using calipers, electromagnetic digitizer, optical surface scanner (OSS), and spiral x-ray computed tomography (SXCT).

DESIGN

The experimental measurement protocol for a single measurement session was as follows: Dot markers were placed on the residuum, and volume and distances were measured using water displacement and calipers; residuum was measured using electromagnetic digitizer; residuum was scanned using three-dimensional (3D) OSS; a negative plaster cast of subject's residuum was made; the residuum was scanned using SXCT scanner. These steps were repeated at a second measurement session. Plaster positive casts were constructed and subsequently measured using the same protocol.

PARTICIPANTS

Thirteen adult below-knee amputee volunteers (subjects) participated in the study, and nine subjects returned for a second measurement session. The study group consisted of 9 men and 4 women; 10 Caucasians and 3 African Americans.

RESULTS

Distance measurements for all measurement devices were repeatable within 1% in vivo and within 0.5% on plaster casts; and volumes were within 1% in vivo and within 0.1% on plaster casts. Distance measurements for each device were precise within 3% in vivo and within 1% on plaster casts; and volumes were within 5% in vivo and within 6% on plaster casts when compared with caliper and water displacement measures.

CONCLUSION

These measurement systems were found to be substantially equivalent in terms of repeatability and precision for measurement of lower extremity residua.

Quantitative trait loci for murine growth

Body size is an archetypal quantitative trait with variation due to the segregation of many gene loci, each of relatively minor effect, and the environment. We examine the effects of quantitative trait loci (QTLs) on age-specific body weights and growth in the F2 intercross of the LG/J and SM/J strains of inbred mice. Weekly weights (1-10 wk) and 75 microsatellite genotypes were obtained for 535 mice. Interval mapping was used to locate and measure the genotypic effects of QTLs on body weight and growth. QTL effects were detected on 16 of the 19 autosomes with several chromosomes carrying more than one QTL. The number of QTLs for age-specific weights varied from seven at 1 week to 17 at 10 wk. The QTLs were each of relatively minor, subequal effect. QTLs affecting early and late growth were generally distinct, mapping to different chromosomal locations indicating separate genetic and physiological systems for early and later murine growth.

Genetic and environmental impacts on litter size and early infant survival in three species of callitrichids

Callitrichids are unusual among anthropoid primates in that they can deliver one to four offspring per litter in captivity. Factors underlying intraspecific variation in litter size are unclear. Data from the Oak Ridge Associated Universities (ORAU) Marmoset Research Center provided an opportunity to investigate determinants of variation in (1) average litter size at birth and (2) average litter size surviving to 2 weeks of age in Saguinus fuscicollis, S. oedipus, and Callithrix jacchus. The objectives were to (1) investigate the impact of parity, husbandry, origin (wild or captive), and subspecific hybridization on phenotypic variance in litter size and (2) to estimate the heritability of litter size at birth and at 14 days. Husbandry changes were the only significant covariate and increased litter size at birth in C. jacchus only. Heritability estimates for average litter size at birth were significant only in S. fuscicollis (h2 = 0.306). Heritability of average litter size at 2 weeks of age was not significantly greater than zero. Estimates of genetic variability were relatively high, however both traits had very high levels of residual variance. The results suggest that litter size is a malleable trait that may respond to environmental factors.

Asymmetric vault modification in Hopi crania

Cradleboarding was practiced by numerous prehistoric and historic populations, including the Hopi. In this group, one result of cradleboarding was bilateral or asymmetric flattening of the posterior occipital. We test whether cradleboarding had significant effects on the morphology of the cranial vault, cranial base, and face. Additionally, we examine associations between direction of flattening and asymmetric craniofacial growth. A skeletal sample of Hopi from the Old Walpi site includes both nonmodified (N = 43) and modified individuals (N = 39). Three-dimensional coordinates of 53 landmarks were obtained using a diagraph. Thirty-six landmarks were used to define nine finite elements in the cranial vault, cranial base, and face. Finite element scaling was used to compare average nonmodified individuals, with averages of bilaterally, right, and left modified individuals. The significance of variation among "treatment" groups was evaluated using a bootstrap test. Pearson product-moment correlations test the association of asymmetry with direction of modification. Hopi cradleboarding has a significant effect on growth of the cranial vault, but does not affect morphology of the cranial base or face. Bilateral flattening of the cranial vault leads to decreased length and increased width of the cranial vault. Flattening of the right or left cranial vault results in ipsilaterally decreased length and width coupled with a corresponding increased length and width on the contralateral side of the cranial vault. There is a significant correlation of size asymmetry with direction of modification in the cranial vault, but not with size or shape change in the cranial base or face.

Fluctuating asymmetry in tamarin (Saguinus) cranial morphology: intra- and interspecific comparisons between taxa with varying levels of genetic heterozygosity

Fluctuating asymmetry has been proposed as a measure of developmental homeostasis and an indicator of populations under stress. However, controversy surrounds not only the validity of an association between fluctuating asymmetry and levels of protein heterozygosity, but also whether fluctuating asymmetry can be used to identify populations under genetic and environmental stress. The relationship between levels of heterozygosity and developmental homeostasis is considered by comparing levels of cranial fluctuating asymmetry in three tamarin samples with contrasting levels of heterozygosity: (1) low heterozygosity cotton-top tamarins (N = 324), (2) presumably normally heterozygous Illiger's saddle-back tamarins (N = 208), and (3) relatively highly heterozygous hybrids between saddle-back tamarin subspecies (N = 31). All specimens originated at the Oak Ridge Associated Universities' Marmoset Research Center. A nested ANOVA design was used to separate out variation due to individual differences, side-to-side differences (fluctuating asymmetry), and measurement error. We found statistically significant levels of fluctuating asymmetry in nearly all of the traits surveyed and a negative correlation between levels of fluctuating asymmetry and genetic heterozygosity. Efforts to use fluctuating asymmetry to identify populations endangered by reduced genetic variability and/or under stress may be inhibited by small sample sizes, neglect of repeated measures, and lack of appropriate reference populations.

Polymorphism for PCR-analyzed microsatellites between the inbred mouse strains LG and SM

Using agarose gel electrophoresis, we surveyed four strains of inbred mice (AKR/J, C57BL/J, LG/J, and SM/J) for 472 microsatellite loci. Agarose electrophoresis proved to be extremely efficient in separating alleles differing by six or more base pairs and detected a majority of allelic differences of between two and six base pairs. Overall, 64.4% of loci showed polymorphism among the four strains, and pairwise comparisons ranged from 42.1% to 48.1%. Microsatellite polymorphism for strains LG/J and SM/J has not been previously described and was sufficiently high (47.1%) to make these size-divergent strains excellent candidates for quantitative trait loci (QTL) analysis of normal growth.

Anthropometric optical surface imaging system repeatability, precision, and validation

Disciplines using human body surface dimensions require accurate, repeatable measurements. This study presents a design for the analysis of repeatability, precision, and validation of a new anthropometric device. This model enables estimation of the proportion of the total variation attributable to each level of data collection. This model is applied to an analysis of repeatability, precision, and validation of the Cencit Imaging System, a new optical surface scanner. Twenty-seven facial landmarks were marked on 10 men and 10 women at two measurement sessions. Two images were scanned during each session, and each image was digitized twice. The Cencit Imaging System results were compared with a previously validated digitizer. The Cencit Imaging System was found to produce accurate, highly repeatable images. Much of the error in this study is attributable to human error in marking landmarks on the subjects. The new imaging system will prove useful in a variety of anthropometric applications.

Epistasis and its contribution to genetic variance components

We present a new parameterization of physiological epistasis that allows the measurement of epistasis separate from its effects on the interaction (epistatic) genetic variance component. Epistasis is the deviation of two-locus genotypic values from the sum of the contributing single-locus genotypic values. This parameterization leads to statistical tests for epistasis given estimates of two-locus genotypic values such as can be obtained from quantitative trait locus studies. The contributions of epistasis to the additive, dominance and interaction genetic variances are specified. Epistasis can make substantial contributions to each of these variance components. This parameterization of epistasis allows general consideration of the role of epistasis in evolution by defining its contribution to the additive genetic variance.

Effect of premature sagittal suture closure on craniofacial morphology in a prehistoric male Hopi

Scaphocephaly is caused by premature sagittal suture closure. It restricts medial-lateral growth of the cranial vault which appears longer and narrower than normal. We examined how this natural malformation affected morphology of the cranial base and face. A prehistoric adult male Hopi with scaphocephaly was compared to 19 normal male Hopi (Old Walpi Series, Field Museum of Natural History, Chicago, IL). Three-dimensional coordinates of 53 landmarks were recorded on the cranial vault, cranial base, and face. Finite element scaling analysis was used to test for significant differences between the scaphocephalic and normal individuals in craniofacial morphology. Finite element scaling results indicated that scaphocephaly has a profound effect on cranial vault morphology. However, morphology of the cranial base and face were only slightly affected and, typically, are within the normal range of variation for the sample. An abnormality in the cranial base is not associated with scaphocephaly in this individual.

Individual genes underlying quantitative traits: molecular and analytical methods

Because of innovations in quantitative and molecular genetics, it is becoming possible to study the individual loci affecting quantitative traits, quantitative trait loci (QTL). We review a number of recently developed molecular and analytical techniques that allow the identification and localization of QTL and the quantification of their effects, and present an example. We suggest that this area of investigation, already in widespread use in the agricultural and medical sciences, will become increasingly important in evolutionary biology as techniques continue to improve.

Genetic epidemiology of colon cancer in the cotton-top tamarin (Saguinus oedipus)

Colon cancer is known to be heritable in humans, but the opportunity to investigate the genetic epidemiology of cancer in nonhuman primates has been limited by the size of available populations. The cotton-top tamarin (Saguinus oedipus) colony at the Oak Ridge Associated Universities (ORAU) Marmoset Research Center is a large population with a high rate of spontaneous colon cancer that has been monitored over several years, thus allowing investigation of the genetic basis for colon cancer in this colony. The presence of colon cancer at death was scored in 392 necropsies at the colony. Genealogical and demographic data for these animals were obtained from colony records. The heritability of the liability to colon cancer was estimated using maximum-likelihood-based pedigree analyses after evaluating the effects of gender, origin (wild-born or laboratory-born), and age at death on cancer experience. Cancer rates were not significantly different between males and females or between wild- and laboratory-born animals. Differences in age at death were also statistically insignificant for both laboratory- and wild-born animals. The heritability estimate for the liability to contract colon cancer is 17% for the ORAU cotton-top tamarins. This heritability estimate is not significantly different from zero, indicating no evidence for heritable variation in cancer experience in this population. If genetic factors affect cotton-top tamarin colon cancer, they are fixed or nearly fixed in this population.

Effects of annular cranial vault modification on the cranial base and face

Artificial modification of the cranial vault was practiced by a number of prehistoric and protohistoric populations, frequently during an infant's first year of life. We test the hypothesis that, in addition to its direct effects on the cranial vault, annular cranial vault modification has a significant indirect effect on cranial base and facial morphology. Two skeletal series from the Pacific Northwest Coast, which include both nonmodified and modified crania, were used: the Kwakiutl (62 nonmodified, 45 modified) and Nootka (28 nonmodified, 20 modified). Three-dimensional coordinates of 53 landmarks were obtained using a diagraph, and 36 landmarks were used to define nine finite elements in the cranial vault, cranial base, and face. Finite element scaling was used to compare average nonmodified and average modified crania, and the significance of the results were evaluated using a bootstrap test. Annular modification of the cranial vault produces significant effects on the morphology of the cranial base and face. Annular modification in the Kwakiutl resulted in restrictions of the cranial vault in the medial-lateral and superior-inferior dimensions and an increase in anterior-posterior growth. Similar dimensional changes are observed in the cranial base. The Kwakiutl face is increased anterior-posteriorly and reduced anterior-laterally to posterior-medially. Similar effects of modification are observed in the Nootka cranial vault and cranial base, though not in the face. These results demonstrate the developmental interdependence of the cranial vault, cranial base, and face.

Cranial deformation and nonmetric trait variation

Cranial deformation is known to influence many traditional craniometric variables, but its effects on nonmetric trait variation are not well characterized. In this study, we examine the effects of three types of deformation (annular, lambdoid flattening, and fronto-occipital) on nonmetric traits, using a large sample of protohistoric and prehistoric crania. Our results indicate that a few traits are increased or decreased in relative frequency by particular types of deformation, but that these effects have little impact on the calculation of biological distances between groups.

Systematics of the Saguinus oedipus group of the bare-face tamarins: evidence from facial morphology

The systematics of the Saguinus oedipus group within the bare-face tamarins remains open to question. Hershkovitz (Living New World Monkeys (Platyrrhini), Vol.1. Chicago: University of Chicago Press, 1977) places the cotton-top and rufus-naped tamarins as subspecies of Saguinus oedipus (S. o. oedipus and S. o. geoffroyi, respectively). In contrast, several other authors have argued that these two taxa should be considered separate species (S. oedipus and S. geoffroyi). Phylogenetic relationships within the group are also disputed. Resolving these different interpretations has been difficult in part because no study of this group has included an objective measure of expected levels of specific vs. subspecific variation. We used facial measurements from 179 adult crania to address the systematics of this group and included a related species that is known to include multiple subspecies. Our sample included three taxa from the S. oedipus group of the bare-face tamarins (S. oedipus, S. geoffroyi, and S. leucopus) and six subspecies from the related hairy-face tamarin species S. fuscicollis. Comparisons to S. leucopus provided a relative measure of species-level differences. Analyses that included S. fuscicollis provided a measure of subspecific variation. There was no evidence of facial sexual dimorphism in any of these taxa. A variety of multivariate statistical analyses including discriminant function and cluster analysis suggest that S. oedipus and S. geoffroyi differ morphologically at a level consistent with species-level distinctions. The extent of differences between these taxa is large. The differences in their facial morphology was on the order of differences between S. oedipus or S. geoffroyi and S. leucopus rather than the extent of variation among S. fuscicollis subspecies. Furthermore, a comparison of collecting localities revealed that the variation we observed among S. oedipus and S. geoffroyi was not clinical but presented a large morphological discontinuity at the boundary between taxa. Our analyses also suggested that S. leucopus is more similar to S. oedipus than is either to S. geoffroyi. Finally, it may be that there are some distinct species within the S. fuscicollis group. However, this hypothesis, along with other phylogenetic relationships suggested by this study, will require more data and further study.

Effects of fronto-occipital artificial cranial vault modification on the cranial base and face

Artificial reshaping of the cranial vault has been practiced by many human groups and provides a natural experiment in which the relationships of neurocranial, cranial base, and facial growth can be investigated. We test the hypothesis that fronto-occipital artificial reshaping of the neurocranial vault results in specific changes in the cranial base and face. Fronto-occipital reshaping results from the application of pads or a cradle board which constrains cranial vault growth, limiting growth between the frontal and occipital and allowing compensatory growth of the parietals in a mediolateral direction. Two skeletal series including both normal and artificially modified crania are analyzed, a prehistoric Peruvian Ancon sample (47 normal, 64 modified crania) and a Songish Indian sample from British Columbia (6 normal, 4 modified). Three-dimensional coordinates of 53 landmarks were measured with a diagraph and used to form 9 finite elements as a prelude to finite element scaling analysis. Finite element scaling was used to compare average normal and modified crania and the results were evaluated for statistical significance using a bootstrap test. Fronto-occipitally reshaped Ancon crania are significantly different from normal in the vault, cranial base, and face. The vault is compressed along an anterior-superior to posterior-inferior axis and expanded along a mediolateral axis in modified individuals. The cranial base is wider and shallower in the modified crania and the face is foreshortened and wider with the anterior orbital rim moving inferior and posterior towards the cranial base. The Songish crania display a different modification of the vault and face, indicating that important differences may exist in the morphological effects of fronto-occipital reshaping from one group to another.

Effects of fronto-occipital cranial reshaping on mandibular form

Cultural reshaping (artificial deformation or modification) of the neurocranial vault provides an artificially increased range of morphological variability within which the relationship between the growing neurocranium and face can be investigated. We analyze crania which have been fronto-occipitally compressed to ascertain possible morphological effects on the mandible. We collected measures of mandibular breadth, length, and height from 82 modified (N = 48) and unmodified (N = 34) crania from a Peruvian Ancon series. Angle classification was also scored in order to investigate whether or not occlusal relationships were affected by neurocranial reshaping. Only intercondylar distance (posterior mandibular breadth) exhibited significant differences between unmodified and modified groups, though this difference was relatively small compared with vault deformation. The modified crania had a higher frequency of normal occlusion (Class I) than the unmodified crania. Increased intercondylar breadth in modified skulls is due to a cascade of effects which begin with a direct effect of the fronto-occipital deforming device on neurocranial shape (increased neurocranial width). The increase in mandibular breadth may be a compensatory response to increased cranial base breadth and maintains articulation between the cranial base and mandible. The increased posterior breadth, coupled with a slight decrease in mandibular depth, may contribute to the change in occlusal relationships suggested for this sample.