Ocular Deformations in Spaceflight-Associated Neuro-Ocular Syndrome and Idiopathic Intracranial Hypertension

Purpose

Spaceflight-associated neuro-ocular syndrome (SANS) shares several clinical features with idiopathic intracranial-hypertension (IIH), namely disc edema, globe-flattening, hyperopia, and choroidal folds. Globe-flattening is caused by increased intracranial pressure (ICP) in IIH, but the cause in SANS is uncertain. If increased ICP alone causes SANS, then the ocular deformations should be similar to IIH; if not, alternative mechanisms would be implicated.

Methods

Using optical coherence tomography (OCT) axial images of the optic nerve head, we compared "pre to post" ocular deformations in 22 patients with IIH to 25 crewmembers with SANS. We used two metrics to assess ocular deformations: displacements of Bruch's membrane opening (BMO-displacements) and Geometric Morphometrics to analyze peripapillary shape changes of Bruch's membrane layer (BML-shape).

Results

We found a large disparity in the mean retinal nerve-fiber layer thickness between SANS (108 um; 95% confidence interval [CI] = 105-111 um) and IIH (300 um; 95% CI = 251-350.1 um). The pattern of BML-shape and BMO-displacements in SANS were significantly different from IIH (P < 0.0001). Deformations in IIH were large and preponderantly anterior, whereas the deformations in SANS were small and bidirectional. The degree of disc edema did not explain the differences in ocular deformations.

Conclusions

This study showed substantial differences in the degree of disc edema and the pattern of ocular deformations between IIH and SANS. The precise cause for these differences is unknown but suggests that there may be fundamental differences in the underlying biomechanics of each consistent with the prevailing hypothesis that SANS is consequent to multiple factors beyond ICP alone. We propose a hypothetical model to explain the differences between IIH and SANS based on the pattern of indentation loads.

The evolution of mammalian brain size

Relative brain size has long been considered a reflection of cognitive capacities and has played a fundamental role in developing core theories in the life sciences. Yet, the notion that relative brain size validly represents selection on brain size relies on the untested assumptions that brain-body allometry is restrained to a stable scaling relationship across species and that any deviation from this slope is due to selection on brain size. Using the largest fossil and extant dataset yet assembled, we find that shifts in allometric slope underpin major transitions in mammalian evolution and are often primarily characterized by marked changes in body size. Our results reveal that the largest-brained mammals achieved large relative brain sizes by highly divergent paths. These findings prompt a reevaluation of the traditional paradigm of relative brain size and open new opportunities to improve our understanding of the genetic and developmental mechanisms that influence brain size.

Efficacy of diffeomorphic surface matching and 3D geometric morphometrics for taxonomic discrimination of Early Pleistocene hominin mandibular molars

Morphometric assessments of the dentition have played significant roles in hypotheses relating to taxonomic diversity among extinct hominins. In this regard, emphasis has been placed on the statistical appraisal of intraspecific variation to identify morphological criteria that convey maximum discriminatory power. Three-dimensional geometric morphometric (3D GM) approaches that utilize landmarks and semi-landmarks to quantify shape variation have enjoyed increasingly popular use over the past twenty-five years in assessments of the outer enamel surface (OES) and enamel-dentine junction (EDJ) of fossil molars. Recently developed diffeomorphic surface matching (DSM) methods that model the deformation between shapes have drastically reduced if not altogether eliminated potential methodological inconsistencies associated with the a priori identification of landmarks and delineation of semi-landmarks. As such, DSM has the potential to better capture the geometric details that describe tooth shape by accounting for both homologous and non-homologous (i.e., discrete) features, and permitting the statistical determination of geometric correspondence. We compare the discriminatory power of 3D GM and DSM in the evaluation of the OES and EDJ of mandibular permanent molars attributed to Australopithecus africanus, Paranthropus robustus and early Homo sp. from the sites of Sterkfontein and Swartkrans. For all three molars, classification and clustering scores demonstrate that DSM performs better at separating the A. africanus and P. robustus samples than does 3D GM. The EDJ provided the best results. P. robustus evinces greater morphological variability than A. africanus. The DSM assessment of the early Homo molar from Swartkrans reveals its distinctiveness from either australopith sample, and the "unknown" specimen from Sterkfontein (Stw 151) is notably more similar to Homo than to A. africanus.

THE ACCURACY OF PHYLOGENETIC ESTIMATION USING THE NEIGHBOR-JOINING METHOD

We studied the factors affecting the accuracy of the neighbor-joining (NJ) method for estimating phylogenies by simulating character change under different evolutionary models applied to twenty different 8-OTU tree topologies that varied widely with respect to tree imbalance and stemminess. The models incorporated three evolutionary rates-constant, varying among lineages, varying among characters-and three evolutionary contexts concerning patterns of character change relative to speciation events-phyletic, speciational, and punctuational. All combinations of the rate and context models were studied. In addition, three different absolute rates of change were investigated. To measure the accuracy, the strict consensus index was computed between the estimated tree and the tree topology along which the data had been generated. The results were analyzed by analysis of variance and compared to a previous study that evaluated UPGMA clustering and maximum parsimony (MP) as phylogenetic estimation techniques. We found evolutionary context and tree imbalance to be the most important factors affecting the accuracy of the NJ method. NJ was more accurate than UPGMA or MP in terms of the average strict consensus index over all treatments. However, no one method was more accurate than the other two for all combinations of treatments. Higher absolute rate of change generally resulted in higher accuracy for all three methods.

ACCURACY OF ESTIMATED PHYLOGENIES: EFFECTS OF TREE TOPOLOGY AND EVOLUTIONARY MODEL

A simulation study was carried out to investigate the relative importance of tree topology (both balance and stemminess), evolutionary rates (constant, varying among characters, and varying among lineages), and evolutionary models in determining the accuracy with which phylogenetic trees can be estimated. The three evolutionary context models were phyletic (characters can change at each simulated time step), speciational (changes are possible only at the time of speciation into two daughter lineages), and punctuational (changes occur at the time of speciation but only in one of the daughter lineages). UPGMA clustering and maximum parsimony ("Wagner trees") methods for estimating phylogenies were compared. All trees were based on eight recent OTUs. The three evolutionary context models were found to have the largest influence on accuracy of estimates by both methods. The next most important effect was that of the stemminess × context interaction. Maximum parsimony and UPGMA performed worst under the punctuational models. Under the phyletic model, trees with high stemminess values could be estimated more accurately and balanced trees were slightly easier to estimate than unbalanced ones. Overall, maximum parsimony yielded more accurate trees than UPGMA-but that was expected for these simulations since many more characters than OTUs were used. Our results suggest that the great majority of estimated phylogenetic trees are likely to be quite inaccurate; they underscore the inappropriateness of characterizing current phylogenetic methods as being for reconstruction rather than for estimation.

Re-evaluation of batoid pectoral morphology reveals novel patterns of diversity among major lineages

Batoids (Chondrichthyes: Batoidea) are a diverse group of cartilaginous fishes which comprise a monophyletic sister lineage to all neoselachians or modern sharks. All species in this group possess anteroposteriorly expanded-pectoral fins, giving them a unique disc-like body form. Reliance on pectoral fins for propulsion ranges from minimal (sawfish) to almost complete dependence (skates and rays). A recent study on the diversity of planform pectoral fin shape in batoids compared overall patterns of morphological variation within the group. However, inconsistent pectoral homology prevented the study from accurately representing relationships within and among major batoid taxa. With previous work in mind, we undertook an independent investigation of pectoral form in batoids and evaluated the implications of shape diversity on locomotion and lifestyle, particularly in the skates (Rajoidei) and rays (Myliobatoidei). We used geometric morphometrics with sliding semilandmarks to analyze pectoral fin outlines and also calculate fin aspect ratios (AR), a functional trait linked to locomotion. In agreement with previous work, our results indicated that much of the evolution of batoid pectoral shape has occurred along a morphological axis that is closely related to AR. For species where kinematic data were available, both shape and AR were associated with swimming mode. This work further revealed novel patterns of shape variation among batoids, including strong bimodality of shape in rays, an intermediate location of skate species in the morphospace between benthic/demersal and pelagic rays, and approximately parallel shape trajectories in the benthic/demersal rays and skates. Finally, manipulation of landmarks verified the need for a consistent and accurate definition of homology for the outcome and efficacy of analyses of pectoral form and function in batoids.

Sexual dimorphism in sister species of Leucoraja skate and its relationship to reproductive strategy and life history

Instances of sexual dimorphism occur in a great variety of forms and manifestations. Most skates (Batoidea: Rajoidei) display some level of body shape dimorphism in which the pectoral fins of mature males develop to create a distinct bell-shaped body not found in females. This particular form of dimorphism is present in each of the sister species Leucoraja erinacea and Leucoraja ocellata, but differences between sexes are much greater in the former. In order to understand the nature and potential causes of pectoral dimorphism, we used geometric morphometrics to investigate allometry of fin shape in L. erinacea and L. ocellata and its relationship to the development of reproductive organs, based on previous work on the bonnethead shark, Sphyrna tiburo. We found that allometric trajectories of overall pectoral shape were different in both species of skate, but only L. erinacea varied significantly with respect to endoskeleton development. Male maturation was characterized by a number of sex-specific morphological changes, which appeared concurrently in developmental timing with elongation of cartilage-supported claspers. We suggest that external sexual dimorphism of pectoral fins in skates is a byproduct of skeletal growth needed for clasper development. Further, the magnitude of male shape change appears to be linked to the differential life histories of species. This work reports for the first time that pectoral dimorphism is a persistent feature in rajoid fishes, occurring in varying degrees across several genera. Lastly, our results suggest that pectoral morphology may be useful as a relative indicator of reproductive strategy in some species.

Effects of lowering cerebrospinal fluid pressure on the shape of the peripapillary retina in intracranial hypertension

PURPOSE

To analyze the deformations of the peripapillary retinal pigment epithelium-basement membrane (ppRPE/BM) layer in response to procedures that lower intracranial pressure (ICP). Second, to demonstrate how shape changes may complement the mean retinal nerve fiber layer (RNFL) thickness as a measure of intracranial hypertension (ICH) and papilledema.

METHODS

We used geometric morphometrics on spectral-domain optical coherence tomography images to analyze shape change of the ppRPE/BM layer after several interventions that lower cerebrospinal fluid (CSF) pressure. We also evaluated the effects of pressure-lowering interventions on both the anterior-posterior displacement of ppRPE/BM and the mean RNFL thickness. Forty-one patients with ICH and papilledema were studied before and after lumbar puncture (20), CSF shunt (9), and medical treatment of idiopathic ICH (23). We also compared the shape of 30 normal subjects to 23 patients whose papilledema resolved after medical treatment.

RESULTS

The ppRPE/BM-layer in ICH and papilledema is characterized by an asymmetric anterior deformation that moves posteriorly and becomes more V-shaped after each pressure-lowering intervention. The differences were statistically significant for all three groups. These shape changes also occur in patients with ongoing ICH who have secondary optic atrophy (without papilledema). Posterior displacement at the margin of the ppRPE/BM layer correlated strongly with overall shape changes.

CONCLUSIONS

The subsurface contour of the ppRPE/BM layer is a dynamic property that changes with CSF pressure-lowering interventions. It can supplement the RNFL thickness as an indirect gauge of ICP and is particularly helpful in patients with secondary optic atrophy. Direct measurements of displacement at the basement membrane opening may serve as a more convenient office-based surrogate for shape analysis.

Combining geometric morphometrics, molecular phylogeny, and micropaleontology to assess evolutionary patterns in Mallomonas (Synurophyceae: Heterokontophyta)

Synurophytes, also known as scaled chrysophytes, are ecologically important algae that produce an array of siliceous structures upon which their taxonomy is based. Despite occupying a key position within the photosynthetic heterokonts, the evolutionary history of synurophytes remains poorly constrained. Here, modern and Middle Eocene siliceous scales of the morphotaxon Mallomonas insignis are used as a model to investigate synurophyte evolutionary patterns. Structural details of scale morphology were examined comparatively with scanning electron microscopy and scored for geometric morphometric analyses to assess the stability of shape characters. Although consistent size differences exist (modern scales are larger than Eocene counterparts), the populations cannot be differentiated on the basis of shape or microstructural detail, implying considerable evolutionary stasis in scale morphology. A time-calibrated relaxed molecular clock analysis using a three-gene concatenated data set (27 strains) suggests that the M. insignis lineage predates the available fossil record, having diverged from closest congeneric taxa in the Cretaceous (≥94 Ma). However, the molecular analysis also implies that considerable genetic variability is present within several morphotaxa of Mallomonas, implying that substantial genetic variability has arisen despite the retention of uniform scale morphologies, and resulting in the widespread occurrence of cryptic taxa. Results from the synurophyte lineage are consistent with the notion of protracted ghost ranges (>10 Ma) implied by the molecular phylogenies of other algal groups, together pointing to the paucity of the fossil record of these organisms on these timescales.

Host genes related to paneth cells and xenobiotic metabolism are associated with shifts in human ileum-associated microbial composition

The aim of this study was to integrate human clinical, genotype, mRNA microarray and 16 S rRNA sequence data collected on 84 subjects with ileal Crohn’s disease, ulcerative colitis or control patients without inflammatory bowel diseases in order to interrogate how host-microbial interactions are perturbed in inflammatory bowel diseases (IBD). Ex-vivo ileal mucosal biopsies were collected from the disease unaffected proximal margin of the ileum resected from patients who were undergoing initial intestinal surgery. Both RNA and DNA were extracted from the mucosal biopsy samples. Patients were genotyped for the three major NOD2 variants (Leufs1007, R702W, and G908R) and the ATG16L1T300A variant. Whole human genome mRNA expression profiles were generated using Agilent microarrays. Microbial composition profiles were determined by 454 pyrosequencing of the V3–V5 hypervariable region of the bacterial 16 S rRNA gene. The results of permutation based multivariate analysis of variance and covariance (MANCOVA) support the hypothesis that host mucosal Paneth cell and xenobiotic metabolism genes play an important role in host microbial interactions.

Inflammatory bowel diseases phenotype, C. difficile and NOD2 genotype are associated with shifts in human ileum associated microbial composition

We tested the hypothesis that Crohn's disease (CD)-related genetic polymorphisms involved in host innate immunity are associated with shifts in human ileum-associated microbial composition in a cross-sectional analysis of human ileal samples. Sanger sequencing of the bacterial 16S ribosomal RNA (rRNA) gene and 454 sequencing of 16S rRNA gene hypervariable regions (V1-V3 and V3-V5), were conducted on macroscopically disease-unaffected ileal biopsies collected from 52 ileal CD, 58 ulcerative colitis and 60 control patients without inflammatory bowel diseases (IBD) undergoing initial surgical resection. These subjects also were genotyped for the three major NOD2 risk alleles (Leu1007fs, R708W, G908R) and the ATG16L1 risk allele (T300A). The samples were linked to clinical metadata, including body mass index, smoking status and Clostridia difficile infection. The sequences were classified into seven phyla/subphyla categories using the Naïve Bayesian Classifier of the Ribosome Database Project. Centered log ratio transformation of six predominant categories was included as the dependent variable in the permutation based MANCOVA for the overall composition with stepwise variable selection. Polymerase chain reaction (PCR) assays were conducted to measure the relative frequencies of the Clostridium coccoides - Eubacterium rectales group and the Faecalibacterium prausnitzii spp. Empiric logit transformations of the relative frequencies of these two microbial groups were included in permutation-based ANCOVA. Regardless of sequencing method, IBD phenotype, Clostridia difficile and NOD2 genotype were selected as associated (FDR ≤ 0.05) with shifts in overall microbial composition. IBD phenotype and NOD2 genotype were also selected as associated with shifts in the relative frequency of the C. coccoides--E. rectales group. IBD phenotype, smoking and IBD medications were selected as associated with shifts in the relative frequency of F. prausnitzii spp. These results indicate that the effects of genetic and environmental factors on IBD are mediated at least in part by the enteric microbiota.

Shape analysis of the peripapillary RPE layer in papilledema and ischemic optic neuropathy

PURPOSE

Geometric morphometrics (GM) was used to analyze the shape of the peripapillary retinal pigment epithelium-Bruch's membrane (RPE/BM) layer imaged on the SD-OCT 5-line raster in normal subjects and in patients with papilledema and ischemic optic neuropathy.

METHODS

Three groups of subjects were compared: 30 normals, 20 with anterior ischemic optic neuropathy (AION), and 25 with papilledema and intracranial hypertension. Twenty equidistant semilandmarks were digitized on OCT images of the RPE/BM layer spanning 2500 μm on each side of the neural canal opening (NCO). The data were analyzed using standard GM techniques, including a generalized least-squares Procrustes superimposition, principal component analysis, thin-plate spline (to visualize deformations), and permutation statistical analysis to evaluate differences in shape variables.

RESULTS

The RPE/BM layer in normals and AION have a characteristic V shape pointing away from the vitreous; the RPE/BM layer in papilledema has an inverted U shape, skewed nasally inward toward the vitreous. The differences were statistically significant. There was no significant difference in shapes between normals and AION. Pre- and posttreatment OCTs, in select cases of papilledema, showed that the inverted U-shaped RPE/BM moved posteriorly into a normal V shape as the papilledema resolved with weight loss or shunting.

CONCLUSIONS

The shape difference in papilledema, absent in AION, cannot be explained by disc edema alone. The difference is a consequence of both the translaminar pressure gradient and the material properties of the peripapillary sclera. GM offers a novel way of statistically assessing shape differences of the peripapillary optic nerve head.

Development of sizing structure for fall arrest harness design

Updated harness designs are needed to accommodate diverse populations in the current workforce. This paper determined an improved fall-arrest harness sizing scheme and strap-length configurations for harness design. A 3-D elliptic Fourier analysis (EFA) procedure with 123 coefficients was developed to quantify torso-shape effect on harness fit, based on 3-D data of 108 women and 108 men. The EFA coefficients were then applied to 600 representative body scans from a national database of 2382 participants to establish an improved sizing system. Study outcomes suggested a more upward back D-ring location for women than current unisex designs to accommodate female torso form and mitigate their fit problem. Results also suggested an improved system of three sizes for women and three sizes for men. New harness sizing charts for women and men were proposed accordingly. Using the most current 3-D whole-body digital scanning technology, this study assembled data from a US workforce to establish an improved fall-arrest harness sizing system and strap configurations for men and women. The information is useful for new generation harness designs to reduce the risk of worker injury.

A COMMENT ON PHYLOGENETIC CORRECTION

Statistical methods are now commonly used to take into account the expected lack of independence of observations across different species (due to their phylogenetic relatedness) when computing correlations or regressions among traits. The methods are often interpreted as removing that part of the regression or correlation that is an artifact due to phylogeny and there is an expectation that the corrected regression or correlation coefficients will usually be closer to zero. It is shown here that this is not an accurate way to interpret these methods. The effect of taking phylogeny into account is to reduce the variance of the estimated regression or correlation coefficients. Their means are not because since estimates of regression coefficients are unbiased whether or not the correct phylogeny is taken into account. Estimates of correlations are only slightly biased (and in the opposite direction that many expect).

Bias and error in estimates of mean shape in geometric morphometrics

Sampling experiments were performed to investigate mean square error and bias in estimates of mean shape produced by different geometric morphometric methods. The experiments use the isotropic error model, which assumes equal and independent variation at each landmark. The case of three landmarks in the plane (i.e., triangles) was emphasized because it could be investigated systematically and the results displayed on the printed page. The amount of error in the estimates was displayed as RMSE surfaces over the space of all possible configurations of three landmarks. Patterns of bias were shown as vector fields over this same space. Experiments were also performed using particular combinations of four or more landmarks in both two and three dimensions. It was found that the generalized Procrustes analysis method produced estimates with the least error and no pattern of bias. Averages of Bookstein shape coordinates performed well if the longest edge was used as the baseline. The method of moments (Stoyan, 1990, Model. Biomet. J. 32, 843) used in EDMA (Lele, 1993, Math. Geol. 25, 573) exhibits larger errors. When variation is not small, it also shows a pattern of bias for isosceles triangles with one side much shorter than the other two and for triangles whose vertices are approximately collinear causing them to resemble their own reflections. Similar problems were found for the log-distance method of Rao and Suryawanshi (1996, Proc. Nat. Acad. Sci. 95, 4121). These results and their implications for the application of different geometric morphometric methods are discussed.

Computing the uniform component of shape variation

Any change in shape of a configuration of landmark points in two or three dimensions includes a uniform component, a component that is a wholly linear (affine) transformation. The formulas for estimating this component have been standardized for two-dimensional data but not for three-dimensional data. We suggest estimating the component by way of the complementarity between the uniform component and the space of partial warps. The component can be estimated by regression in either one space or the other: regression on the partial warps, followed by their removal, or regression on a basis for the uniform component itself. Either of the new methods can be used for both two- and three-dimensional landmark data and thus generalize Bookstein's (1996, pages 153-168 in Advances in morphometrics [L. F. Marcus et al., eds.], Plenum, New York) linearized Procrustes formula for estimating the uniform component in two dimensions.

Comparative methods for the analysis of continuous variables: geometric interpretations

This study is concerned with statistical methods used for the analysis of comparative data (in which observations are not expected to be independent because they are sampled across phylogenetically related species). The phylogenetically independent contrasts (PIC), phylogenetic generalized least-squares (PGLS), and phylogenetic autocorrelation (PA) methods are compared. Although the independent contrasts are not orthogonal, they are independent if the data conform to the Brownian motion model of evolution on which they are based. It is shown that uncentered correlations and regressions through the origin using the PIC method are identical to those obtained using PGLS with an intercept included in the model. The PIC method is a special case of PGLS. Corrected standard errors are given for estimates of the ancestral states based on the PGLS approach. The treatment of trees with hard polytomies is discussed and is shown to be an algorithmic rather than a statistical problem. Some of the relationships among the methods are shown graphically using the multivariate space in which variables are represented as vectors with respect to OTUs used as coordinate axes. The maximum-likelihood estimate of the autoregressive parameter, p, has not been computed correctly in previous studies (an appendix with MATLAB code provides a corrected algorithm). The importance of the eigenvalues and eigenvectors of the connection matrix, W, for the distribution of p is discussed. The PA method is shown to have several problems that limit its usefulness in comparative studies. Although the PA method is a generalized least-squares procedure, it cannot be made equivalent to the PGLS method using a phylogenetic model.

Use of two-block partial least-squares to study covariation in shape

The relatively new two-block partial least-squares method for analyzing the covariance between two sets of variables is described and contrasted with the well-known method of canonical correlation analysis. Their statistical properties, type of answers, and visualization techniques are discussed. Examples are given to show its usefulness in comparing two sets of variables--especially when one or both of the sets of variables are shape variables from a geometric morphometric study.

A geometric morphometric assessment of change in midline brain structural shape following a first episode of schizophrenia

BACKGROUND

Previous reports indicate that brain structural abnormalities may be progressive in some patients with schizophrenia. Our study was designed to determine deviations in the shape of midline brain structures at the time of onset of symptoms of schizophrenia and 3-5 years later.

METHODS

Eleven landmarks were located on the midsagittal magnetic resonance imagery brain scans of 55 patients with schizophrenia and 22 nonpsychiatric control individuals. Geometric morphometric methods were used for the extraction of shape variables from landmark coordinates. Permutation tests were used to test the effects of gender, diagnosis, time elapsed since illness onset, and age on brain shape.

RESULTS

The diagnosis-by-time interaction and the effect of gender were significantly different from zero (p<.027 and p <.039, respectively). The effect of time was significant in patients (p <.002), but not in control subjects. Some anatomical abnormalities in mean patient brain morphology seem to be present both at the time of diagnosis and at follow-up. These are similar to anomalies reported by previous geometric morphometrics studies.

CONCLUSIONS

Some previously identified brain abnormalities are detectable at the time of first hospitalization. The rapid change in midline brain morphology in patients with schizophrenia during the subsequent 3-5 years is consistent with either a neurodegenerative disease process or an effect of treatment with psychiatric drugs. There is a sexual dimorphism in brain morphology that might be reduced by schizophrenia.

Ecological character displacement in Plethodon: biomechanical differences found from a geometric morphometric study

Ecological character displacement describes a pattern where morphological differences between sympatric species are enhanced through interspecific competition. Although widely considered a pervasive force in evolutionary ecology, few clear-cut examples have been documented. Here we report a case of ecological character displacement between two salamander species, Plethodon cinereus and Plethodon hoffmani. Morphology was quantified by using linear measurements and landmark-based geometric morphometric methods for specimens from allopatric and sympatric populations from two geographic transects in south-central Pennsylvania, and stomach contents were assayed to quantify food resource use. Morphological variation was also assessed in 13 additional allopatric populations. In both transects, we found significant morphological differentiation between sympatric populations that was associated with a reduction in prey consumption in sympatry and a segregation of prey according to prey size. No trophic morphological or resource use differences were found between allopatric populations, and comparisons of sympatric populations with randomly paired allopatric populations revealed that the observed sympatric morphological differentiation was greater than expected by chance. The major trophic anatomical differences between sympatric populations relates to functional and biomechanical differences in jaw closure: sympatric P. hoffmani have a faster closing jaw, whereas sympatric P. cinereus have a slower, stronger jaw. Because salamanders immobilize prey of different sizes in different ways, and because the observed sympatric biomechanical differences in jaw closure are associated with the differences in prey consumption, the observed character displacement has a functional ecological correlate, and we can link changes in form with changes in function in this apparent example of character displacement.

Statistical power comparisons among alternative morphometric methods

This paper compares the statistical power of various tests that have been proposed to test for equality of shape in two populations. Power surfaces are computed with emphasis on the simplest case of three points in the plane (i.e., landmarks at the vertices of a triangle). Goodall's ([1991] J Roy Stat Soc Serb 53:285-339) F-test was found to have the highest power followed by T(2)-tests using Kendall tangent space coordinates. Power for T(2)-tests using Bookstein shape coordinates was good if the baseline was not the shortest side of the triangle. The Rao and Suryawanshi ([1996] Proc Natl Acad Sci 93:12132-12136 and Rao and Suryawanshi [1998] Proc Natl Acad Sci 95:4121-4125) shape variables had much lower power when triangles were not close to being equilateral. Power surfaces for the EDMA-I T statistic revealed very low power for many shape comparisons including those between very different shapes. Power surface for the EDMA-II Z statistic were also complicated and depended strongly on the choice of baseline used for size scaling. The type I error rate was also often not correct for this method. Results for more than three landmarks are also presented. The implications of the results for practical applications of morphometrics are discussed.

Comparing frontal cranial profiles in archaic and modern homo by morphometric analysis

Archaic and modern human frontal bones are known to be quite distinct externally, by both conventional visual and metric evaluation. Internally this area of the skull has been considerably less well-studied. Here we present results from a comparison of interior, as well as exterior, frontal bone profiles from CT scans of five mid-Pleistocene and Neanderthal crania and 16 modern humans. Analysis was by a new morphometric method, Procrustes analysis of semi-landmarks, that permits the statistical comparison of curves between landmarks. As expected, we found substantial external differences between archaic and modern samples, differences that are mainly confined to the region around the brow ridge. However, in the inner median-sagittal profile, the shape remained remarkably stable over all 21 specimens. This implies that no significant alteration in this region has taken place over a period of a half-million years or more of evolution, even as considerable external change occurred within the hominid clade spanning several species. This confirms that the forms of the inner and outer aspects of the human frontal bone are determined by entirely independent factors, and further indicates unexpected stability in anterior brain morphology over the period during which modern human cognitive capacities emerged. Anat Rec (New Anat): 257:217-224, 1999.

Criteria for the evaluation of alternative environmental monitoring variables: Theory and an application using winter flounder (Pleuronectes americanus) and Dover sole (Microstomus pacificus)

The design of environmental monitoring programs is frequently hampered by a lack of objective, quantitative criteria for evaluating alternative monitoring variables. In this paper we describe two such criteria, which we call samples required - the number of samples required to detect a given change in value - and information imparted - the amount of environmental information revealed by the monitoring variable. We then use these criteria to evaluate fin erosion in winter flounder (Pleuronectes americanus) and Dover sole (Microstomus pacificus) as marine environmental monitoring variables. Two methods for determining the samples required use contaminated and reference areas to estimate the sample statistics of a hypothetical impacted population. The first method is based on the overall difference in the proportions of diseased fish in the reference and hypothetical populations. The second treats the proportion of diseased fish in individual trawls as the variate and determines the samples required based on the mean and variance of the reference and contaminated populations. We use both methods to predict the number of trawls needed to detect an increase of 200% in fin erosion in the reference population. The first method had greater statistical power but assumes spatially homogeneous populations. The second method accounts for environmental patchiness. For Dover sole it predicted 1661 trawls would be needed to detect the 200% increase. An estuarine winter flounder population would require 74 trawls, and an oceanic winter flounder population would require 142.5 trawls. It appears that fin erosion in winter flounder may be a useful indicator of environmental contamination, but several stipulations apply. Migration may inflate the number of diseased fish observed in the reference population, and a more detailed etiology of the disease is required, including an understanding of what contaminants are responsible for manifestation of the disease.