Quantifying morphological variation is critical for conducting anatomical research. Three‐dimensional geometric morphometric (3D GM) landmark analyses quantify shape using homologous Cartesian coordinates (landmarks). Setting up a high‐density landmark set and placing it on all specimens, however, can be a time‐consuming task. Weighted spherical harmonics (SPHARM) provides an alternative method for analyzing the shape of such objects. Here we compare sliding semilandmark and SPHARM analyses of the calcaneus of
Recent methodological advances have led to a rapid expansion of evolutionary studies employing three‐dimensional landmark‐based geometric morphometrics (GM). GM methods generally enable researchers to capture and compare complex shape phenotypes, and to quantify their relationship to environmental gradients. However, some recent studies have shown that the common, inexpensive, and relatively rapid two‐dimensional GM methods can distort important information and produce misleading results because they cannot capture variation in the depth (
- NSF-PAR ID:
- 10455725
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Ecology and Evolution
- Volume:
- 10
- Issue:
- 23
- ISSN:
- 2045-7758
- Page Range / eLocation ID:
- p. 13297-13311
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract Gorilla gorilla gorilla (n = 20),Pan troglodytes troglodytes (n = 20), andHomo sapiens (n = 20) to determine whether the SPHARM and sliding semilandmark analyses capture comparable levels of shape variation. We also compare both the sliding semilandmark and SPHARM analyses to a novel combination of the two methods, here termed SPHARM–sliding. In SPHARM–sliding, the vertices of the surface models produced from the SPHARM analysis (that are the same in number and relative location) are used as the starting landmark positions for a sliding semilandmark analysis. Calcaneal shape variation quantified by all three analyses was summarized using separate principal components analyses. Results were compared using the root mean square (RMS) and maximum distance between surface models of species averages scaled (up) to centroid size created from each analysis. The average RMS was 0.23 mm between sliding semilandmark and SPHARM average surface models, 0.19 mm between SPHARM and SPHARM sliding average surface models, and 0.22 mm between sliding semilandmark and SPHARM sliding average surface models. Although results indicate that all three analyses are comparable methods for 3D shape analysis, there are advantages and disadvantages to each. While the SPHARM analysis is less time‐intensive, it is unable to capture the same level of detail around the sharp edges of articular facets on average surface models as the sliding semilandmark analysis. The SPHARM analysis also does not allow for individual articular facets to be analyzed in isolation. SPHARM–sliding, however, captures the same level of detail as the sliding semilandmark analysis, and (as in the sliding semilandmark analysis) allows for the evaluation of individual portions of bone. SPHARM is a comparable method to a 3D GM analysis for small, irregularly shaped bones, such as the calcaneus, and SPHARM–sliding allows for an expedited set up process for a sliding semilandmark analysis. -
null (Ed.)Abstract Phytoplankton biomass is routinely estimated using relationships between cell volume and carbon (C) and nitrogen (N) content that have been defined using diverse plankton that span orders of magnitude in size. Notably, volume has traditionally been estimated with geometric approximations of cell shape using cell dimensions from planar two-dimensional (2D) images, which requires assumptions about the third, depth dimension. Given advances in image processing, we examined how cell volumes determined from three-dimensional (3D), confocal images affected established relationships between phytoplankton cell volume and C and N content. Additionally, we determined that growth conditions could result in 30–40% variation in cellular N and C. 3D phytoplankton cell volume measurements were on average 15% greater than the geometric approximations from 2D images. Volume method variation was minimal compared to both intraspecific variation in volumes (~30%) and the 50-fold variation in elemental density among species. Consequently, C:vol and N:vol relationships were unaltered by volume measurement method and growth environment. Recent advances in instrumentation, including those for at sea and autonomous applications can be used to estimate plankton biomass directly. Going forward, we recommend instrumentation that permits species identification alongside size and shape characteristics for plankton biomass estimates.more » « less
-
Abstract Whether it is swimming, walking, eating, or jumping, motions are a fundamental way in which organisms interact with their environment. Understanding how morphology contributes to motion is a primary focus of kinematic research and is necessary for gaining insights into the evolution of functional systems. However, an element that is largely missing from traditional analyses of motion is the spatial context in which they occur. We explore an application of geometric morphometrics (GM) for analyzing and comparing motions to evaluate the outputs of biomechanical linkage models. We focus on a common model for oral jaw mechanics of perciform fishes, the fourbar linkage, using GM to summarize motion as a trajectory of shape change. Two traits derived from trajectories capture the total kinesis generated by a linkage (trajectory length) and the kinematic asynchrony (KA) of its mobile components (trajectory nonlinearity). Oral jaw fourbar data from two subfamilies of Malagasy cichlids were used to generate form–function landscapes, describing broad features of kinematic diversity. Our results suggest that kinesis and KA have complex relationships with fourbar morphology, each displaying a pattern in which different shapes possess equivalent kinematic trait values, known as many-to-one mapping of form-to-function. Additionally, we highlight the observation that KA captures temporal differences in the activation of motion components, a feature of kinesis that has long been appreciated but was difficult to measure. The methods used here to study fourbar linkages can also be applied to more complex biomechanical models and broadly to motions of live organisms. We suggest that they provide a suitable alternative to traditional approaches for evaluating linkage function and kinematics.more » « less
-
more » « less
Abstract This study evaluates the methods of identifying the height
zi of the top of the convective boundary layer (CBL) during winter (December and January) over the Great Lakes and nearby land areas using observations taken by the University of Wyoming King Air research aircraft during the Lake-Induced Convection Experiment (1997/98) and Ontario Winter Lake-effect Systems (2013/14) field campaigns. Since CBLs facilitate vertical mixing near the surface, the most direct measurement ofzi is that above which the vertical velocity turbulent fluctuations are weak or absent. Thus, we usezi from the turbulence method as the “reference value” to whichzi from other methods, based on bulk Richardson number (Rib ), liquid water content, and vertical gradients of potential temperature, relative humidity, and water vapor mixing ratio, are compared. The potential temperature gradient method using a threshold value of 0.015 K m−1for soundings over land and 0.011 K m−1for soundings over lake provided the estimates ofzi that are most consistent with the turbulence method. The Rib threshold-based method, commonly used in numerical simulation studies, underestimatedzi . Analyzing the methods’ performance on the averaging windowz avgwe recommend usingz avg= 20 or 50 m forzi estimations for lake-effect boundary layers. The present dataset consists of both cloudy and cloud-free boundary layers, some having decoupled boundary layers above the inversion top. Because cases of decoupled boundary layers appear to be formed by nearby synoptic storms, we recommend use of the more general term, elevated mixed layers.Significance Statement The depth
zi of the convective atmospheric boundary layer (CBL) strongly influences precipitation rates during lake-effect snowstorms (LES). However, variouszi approximation methods produce significantly different results. This study utilizes extensive concurrently collected observations by project aircraft during two LES field studies [Lake-Induced Convection Experiment (Lake-ICE) and OWLeS] to assess howzi from common estimation methods compare with “reference”zi derived from turbulent fluctuations, a direct measure of CBL mixing. For soundings taken both over land and lake; with cloudy or cloud-free conditions, potential temperature gradient (PTG) methods provided the best agreement with the referencezi . A method commonly employed in numerical simulations performed relatively poorly. Interestingly, the PTG method worked equally well for “coupled” and elevated decoupled CBLs, commonly associated with nearby cyclones. -
Abstract Modern computational and imaging methods are revolutionizing the fields of comparative morphology, biomechanics, and ecomorphology. In particular, imaging tools such as X-ray micro computed tomography (µCT) and diffusible iodine-based contrast enhanced CT allow observing and measuring small and/or otherwise inaccessible anatomical structures, and creating highly accurate three-dimensional (3D) renditions that can be used in biomechanical modeling and tests of functional or evolutionary hypotheses. But, do the larger datasets generated through 3D digitization always confer greater power to uncover functional or evolutionary patterns, when compared with more traditional methodologies? And, if so, why? Here, we contrast the advantages and challenges of using data generated via (3D) CT methods versus more traditional (2D) approaches in the study of skull macroevolution and feeding functional morphology in bats. First, we test for the effect of dimensionality and landmark number on inferences of adaptive shifts during cranial evolution by contrasting results from 3D versus 2D geometric morphometric datasets of bat crania. We find sharp differences between results generated from the 3D versus some of the 2D datasets (xy, yz, ventral, and frontal), which appear to be primarily driven by the loss of critical dimensions of morphological variation rather than number of landmarks. Second, we examine differences in accuracy and precision among 2D and 3D predictive models of bite force by comparing three skull lever models that differ in the sources of skull and muscle anatomical data. We find that a 3D model that relies on skull µCT scans and muscle data partly derived from diceCT is slightly more accurate than models based on skull photographs or skull µCT and muscle data fully derived from dissections. However, the benefit of using the diceCT-informed model is modest given the effort it currently takes to virtually dissect muscles from CT scans. By contrasting traditional and modern tools, we illustrate when and why 3D datasets may be preferable over 2D data, and vice versa, and how different methodologies can complement each other in comparative analyses of morphological function and evolution.more » « less
