The influence of data shape acquisition process and geometric accuracy of the mandible for numerical simulation

Computer-aided technologies have allowed new 3D modelling capabilities and engineering analyses based on experimental and numerical simulation. It has enormous potential for product development, such as biomedical instrumentation and implants. However, due to the complex shapes of anatomical structures, the accuracy of these technologies plays an important key role for adequate and accurate finite element analysis (FEA). The objective of this study was to determine the influence of the geometry variability between two digital models of a human model of the mandible. Two different shape acquisition techniques, CT scan and 3D laser scan, were assessed. A total of 130 points were controlled and the deviations between the measured points of the physical and 3D virtual models were assessed. The results of the FEA study showed a relative difference of 20% for the maximum displacement and 10% for the maximum strain between the two geometries.     

Finite element analysis of peri-implant bone volume affected by stresses around Morse taper implants: effects of implant positioning to the bone crest

Abstract

Objectives: The purpose of the present study was to evaluate the distribution and magnitude of stresses through the bone tissue surrounding Morse taper dental implants at different positioning relative to the bone crest. Materials and Methods: A mandibular bone model was obtained from a computed tomography scan. A three-dimensional (3D) model of Morse taper implant-abutment systems placed at the bone crest (equicrestal) and 2?mm bellow the bone crest (subcrestal) were assessed by finite element analysis (FEA). FEA was carried out on axial and oblique (45°) loading at 150 N relatively to the central axis of the implant. The von Mises stresses were analysed considering magnitude and volume of affected peri-implant bone. Results: On vertical loading, maximum von Mises stresses were recorded at 6-7?MPa for trabecular bone while values ranging from 73 up to 118?MPa were recorded for cortical bone. On oblique loading at the equiquestral or subcrestal positioning, the maximum von Mises stresses ranged from 15 to 21?MPa for trabecular bone while values at 150?MPa were recorded for the cortical bone. On vertical loading, >99.9vol.% cortical bone volume was subjected to a maximum of 2?MPa while von Mises stress values at 15?MPa were recorded for trabecular bone. On oblique loading, >99.9vol.% trabecular bone volume was subjected to maximum stress values at 5?MPa, while von Mises stress values at 35?MPa were recorded for >99.4vol.% cortical bone. Conclusions: Bone volume-based stress analysis revealed that most of the bone volume (>99% by vol) was subjected to significantly lower stress values around Morse taper implants placed at equicrestal or subcrestal positioning. Such analysis is commentary to the ordinary biomechanical assessment of dental implants concerning the stress distribution through peri-implant sites.     

Advanced human carotid plaque progression correlates positively with flow shear stress using follow-up scan data: An in vivo MRI multi-patient 3D FSI study

Although it has been well-accepted that atherosclerosis initiation and early progression correlate negatively with flow wall shear stresses (FSS), increasing evidence suggests mechanisms governing advanced plaque progression are not well understood. Fourteen patients were scanned 2–4 times at 18 month intervals using a histologically validated multi-contrast magnetic resonance imaging (MRI) protocol to acquire carotid plaque progression data. Thirty-two scan pairs (baseline and follow-up scans) were formed with slices matched for model construction and analysis. 3D fluid–structure interaction (FSI) models were constructed and plaque wall stress (PWS) and flow shear stress (FSS) were obtained from all matching lumen data points (400–1000 per plaque; 100 points per matched slice) to quantify correlations with plaque progression measured by vessel wall thickness increase (WTI). Using FSS and PWS data from follow-up scan, 21 out of 32 scan pairs showed a significant positive correlation between WTI and FSS (positive/negative/no significance ratio=21/8/3), and 26 out of 32 scan pairs showed a significant negative correlation between WTI and PWS (positive/negative/no significance ratio=2/26/4). The mean FSS value of lipid core nodes (n=5294) from all 47 plaque models was 63.5 dyn/cm², which was 45% higher than that from all normal vessel nodes (n=27553, p<0.00001). The results from this intensive FSI study indicate that flow shear stress from follow-up scan correlates positively with advanced plaque progression which is different from what has been observed in plaque initiation and early-stage progression. It should be noted that the correlation results do not automatically lead to any causality conclusions.     

Normalized glandular dose coefficients in mammography,digital breast tomosynthesis and dedicated breast CT

PurposeTo provide mean glandular dose (MGD) estimates via Monte Carlo (MC) simulations as a function of the breast models and scan parameters in mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (BCT).MethodsThe MC code was based on GEANT4 toolkit. The simulated compressed breast was either a cylinder with a semi-circular section or ad hoc shaped for oblique view (MLO). In DBT we studied the influence of breast models and exam parameters on the T-factors (i.e. the conversion factor for the calculation of the MGD in DBT from that for a 0-degree projection), and in BCT we investigated the influence on the MGD estimates of the ion chamber volume used for the air kerma measurements.ResultsIn mammography, a model representative of a breast undergoing an MLO view exam did not produce substantial differences (0.4%) in MGD estimates, when compared to a conventional cranio-caudal (CC) view breast model. The beam half value layer did not present a significant influence on T-factors in DBT (<0.8%), while the skin model presented significant influence on MGD estimates (up to 3.3% at 30 degrees scan angle), increasing for larger scan angles. We derived a correction factor for taking into account the different ion chamber volume used in MGD estimates in BCT.ConclusionsA series of MC code modules for MGD estimates in 2D and 3D breast imaging have been developed in order to take into account the most recent advances in breast models.     

Identifying the Functional Flexion-extension Axis of the Knee: An In-Vivo Kinematics Study

Purpose

This study aimed to calculate the flexion-extension axis (FEA) of the knee through in-vivo knee kinematics data, and then compare it with two major anatomical axes of the femoral condyles: the transepicondylar axis (TEA) defined by connecting the medial sulcus and lateral prominence, and the cylinder axis (CA) defined by connecting the centers of posterior condyles.

Methods

The knee kinematics data of 20 healthy subjects were acquired under weight-bearing condition using bi-planar x-ray imaging and 3D-2D registration techniques. By tracking the vertical coordinate change of all points on the surface of femur during knee flexion, the FEA was determined as the line connecting the points with the least vertical shift in the medial and lateral condyles respectively. Angular deviation and distance among the TEA, CA and FEA were measured.

Results

The TEA-FEA angular deviation was significantly larger than that of the CA-FEA in 3D and transverse plane (3.45° vs. 1.98°, p < 0.001; 2.72° vs. 1.19°, p = 0.002), but not in the coronal plane (1.61° vs. 0.83°, p = 0.076). The TEA-FEA distance was significantly greater than that of the CA-FEA in the medial side (6.7 mm vs. 1.9 mm, p < 0.001), but not in the lateral side (3.2 mm vs. 2.0 mm, p = 0.16).

Conclusion

The CA is closer to the FEA compared with the TEA; it can better serve as an anatomical surrogate for the functional knee axis.     

Multi-patient study for coronary vulnerable plaque model comparisons: 2D/3D and fluid–structure interaction simulations

Several image-based computational models have been used to perform mechanical analysis for atherosclerotic plaque progression and vulnerability investigations. However, differences of computational predictions from those models have not been quantified at multi-patient level. In vivo intravascular ultrasound (IVUS) coronary plaque data were acquired from seven patients. Seven 2D/3D models with/without circumferential shrink, cyclic bending and fluid–structure interactions (FSI) were constructed for the seven patients to perform model comparisons and quantify impact of 2D simplification, circumferential shrink, FSI and cyclic bending plaque wall stress/strain (PWS/PWSn) and flow shear stress (FSS) calculations. PWS/PWSn and FSS averages from seven patients (388 slices for 2D and 3D thin-layer models) were used for comparison. Compared to 2D models with shrink process, 2D models without shrink process overestimated PWS by 17.26%. PWS change at location with greatest curvature change from 3D FSI models with/without cyclic bending varied from 15.07% to 49.52% for the seven patients (average = 30.13%). Mean Max-FSS, Min-FSS and Ave-FSS from the flow-only models under maximum pressure condition were 4.02%, 11.29% and 5.45% higher than those from full FSI models with cycle bending, respectively. Mean PWS and PWSn differences between FSI and structure-only models were only 4.38% and 1.78%. Model differences had noticeable patient variations. FSI and flow-only model differences were greater for minimum FSS predictions, notable since low FSS is known to be related to plaque progression. Structure-only models could provide PWS/PWSn calculations as good approximations to FSI models for simplicity and time savings in calculation.

     

A geometric approach to study the contact mechanisms in the patellofemoral joint of normal versus patellofemoral pain syndrome subjects

The biomechanics of the patellofemoral (PF) joint is complex in nature, and the aetiology of such manifestations of PF instability as patellofemoral pain syndrome (PFPS) is still unclear. At this point, the particular factors affecting PFPS have not yet been determined. This study has two objectives: (1) The first is to develop an alternative geometric method using a three-dimensional (3D) registration technique and linear mapping to investigate the PF joint contact stress using an indirect measure: the depth of virtual penetration (PD) of the patellar cartilage surface into the femoral cartilage surface. (2) The second is to develop 3D PF joint models using the finite element analysis (FEA) to quantify in vivo cartilage contact stress and to compare the peak contact stress location obtained from the FE models with the location of the maximum PD. Magnetic resonance images of healthy and PFPS subjects at knee flexion angles of 15°, 30° and 45° during isometric loading have been used to develop the geometric models. The results obtained from both approaches demonstrated that the subjects with PFPS show higher PD and contact stresses than the normal subjects. Maximum stress and PD increase with flexion angle, and occur on the lateral side in healthy and on the medial side in PFPS subjects. It has been concluded that the alternative geometric method is reliable in addition to being computationally efficient compared with FEA, and has the potential to assess the mechanics of PFPS with an accuracy similar to the FEA.     

Optimizing cardiac material parameters with a genetic algorithm

Determining the unknown material parameters of intact ventricular myocardium can be challenging due to highly nonlinear material behavior. Previous studies combining a gradient-search optimization procedure with finite element analysis (FEA) were limited to two-dimensional (2D) models or simplified three-dimensional (3D) geometries. Here we present a novel scheme to estimate unknown material parameters for ventricular myocardium by combining a genetic algorithm (GA) with nonlinear finite element analysis. This approach systematically explores the domain of the material parameters. The objective function to minimize was the error between simulated strain data and finite element model strains. The proposed scheme was validated for a 2D problem using a realistic material law for ventricular myocardium. Optimized material parameters were generally within 0.5% of the true values. To demonstrate the robustness of the new scheme, unknown material parameters were also determined for a realistic 3D heart model with an exponential hyperelastic material law. When using strains from two material points, the algorithm converged to parameters within 5% of the true values. We conclude that the proposed scheme is robust when estimating myocardial material parameters in 2D and 3D models.     

Computational simulation study on the anion recognition properties of functionalized tetraphenyl porphyrins

Dedicated bond force constant and bulk modulus of C _n fullerenes (n = 20, 28, 36, 50, 60) are computed using density functional theory (DFT). DFT predicts bond force constants of 611, 648, 675, 686, and 691 N/m, for C₂₀, C₂₈, C₃₆, C₅₀, and C₆₀, respectively, indicating that the bond force constant increases for larger fullerenes. The bulk modulus predicted by DFT increases with decreased fullerene diameter, from 0.874 TPa for C₆₀ to 1.830 TPa for C₂₀. The bond force constants predicted by DFT are then used as an input for finite element analysis (FEA) of the fullerenes, considered as spatial frames in structural models where the bond stiffness is represented by the DFT-computed bond force constant. In agreement with DFT, FEA predicts that smaller fullerenes are stiffer, and underestimates the bulk modulus with respect to DFT. The difference between the FEA and DFT predictions of the bulk modulus decreases as the size of the fullerene increases, from 20.9 % difference for C₂₀ to only 4 % difference for C₆₀. Thus, it is concluded that knowing the appropriate bond force constant, FEA can be used as a plausible approximation to model the elastic behavior of small fullerenes.

     

Diagnostic Efficacy of 18F-FDG PET/MRI in Peripheral Nerve Injury Models

The aim of this study was to evaluate the diagnostic efficacy of ¹⁸F-FDG PET/MRI in two different peripheral neuropathic pain models using the injured rat sciatic nerves. Twelve rats, with operation on left sciatic nerves, were evenly divided into three groups: sham surgery (control group), crushing injury and chronic constriction injury (CCI) (experimental groups). The nerve damage was assessed at 3 weeks postoperatively using following methods: paw withdrawal threshold values (RevWT), maximum standardized uptake values on PET/MRI images (SUVR), and counting the number of myelinated axons in proximal and distal sites of nerve injury (MAxR). The results were quantified and statistically analyzed. Compared to the control group, the crushing injury demonstrated significant differences in RevWT (p < 0.0001) and SUVR (p = 0.027) and the CCI group demonstrated significant differences in RevWT (p < 0.0001), SUVR (p = 0.001) and MAxR (p = 0.048). There were no significant differences between the two experimental groups for all assessments. Correlation analysis demonstrated that RevWT and SUVR assessments were highly correlated (r = -− 0.710, p = 0.010), and SUVR and MAxR were highly correlated (r = 0.611, p = 0.035). However, there was no significant correlation between RevWT and MAxR. The PET scan may be a valuable imaging modality to enable noninvasive, objective diagnosis of neuropathic pain caused by peripheral nerve injury. Also, MRI fused with PET may help clarify the anatomic location of soft tissue structures, including the peripheral nerves.

     

Assessing the transferability of hydraulic habitat models for juvenile Atlantic salmon

Hydraulic habitat models are logistically and technically challenging and expensive to produce. They are therefore frequently transferred between rivers, often with different environmental conditions, without validation. Although studies have recognised problems with model transfer, few have assessed the consequences for model predictions. This study investigated the local (within sub-catchment) transfer of hydraulic habitat models developed for Atlantic salmon (Salmo salar) fry. Two adjacent reaches were chosen for the study, each containing pool, riffle, glide and run habitats. Detailed topographic surveys were used to develop hydraulic models for the study reaches. Substrate and cover were characterised using transects. Seasonal fish habitat use was characterised by electrofishing. Uniform saturation stocking was employed to ensure standard starting densities across all habitats. Generalised additive models were fitted to fry abundance data. Depth, velocity, dominant substrate and cover were used as predictor variables. The distribution of depth, velocity, dominant substrate and cover differed between the two reaches, but was consistent within reaches across seasons. Substrate exhibited the greatest inter-reach difference. Velocity and depth were the most important individual predictors of fry abundance, with the highest densities observed at moderate velocities (∼0.4 m s⁻¹) and low depths (∼0.1 m) across all seasons. When models were transferred locally between reaches, those that were adjusted for inter-reach differences in mean fry abundance (i.e. those predicting changes in relative abundance) performed better than those transferred without adjustment. Complex models that included substrate and cover typically explained some of the variation in abundance in the transferred reach but performed less well than models containing hydraulic parameters alone. This potentially reflected inter-reach differences in the distribution of substrate and cover. This study suggests that (1) uniform stocking is useful for examining habitat-abundance relationships free from the influence of patchy egg deposition; (2) habitat models should be developed at sites offering maximum environmental complexity at a local level; (3) scientists and managers should avoid transferring models between locations with different environmental characteristics, especially in the absence of model validation; (4) complex models should be avoided, and those containing hydraulic variables alone should be considered, if predictions of habitat quality are to be made at new sites.     

Brief communication: comparing loading scenarios in lower first molar supporting bone structure using 3D finite element analysis

Finite element analysis (FEA) is a widespread technique to evaluate the stress/strain distributions in teeth or dental supporting tissues. However, in most studies occlusal forces are usually simplified using a single vector (i.e., point load) either parallel to the long tooth axis or oblique to this axis. In this pilot study we show how lower first molar occlusal information can be used to investigate the stress distribution with 3D FEA in the supporting bone structure. The LM₁ and the LP²‐LM¹ of a dried modern human skull were scanned by μCT in maximum intercuspation contact. A kinematic analysis of the surface contacts between LM₁ and LP²‐LM¹ during the power stroke was carried out in the occlusal fingerprint analyzer (OFA) software to visualize contact areas during maximum intercuspation contact. This information was used for setting the occlusal molar loading to evaluate the stress distribution in the supporting bone structure using FEA. The output was compared to that obtained when a point force parallel to the long axis of the tooth was loaded in the occlusal basin. For the point load case, our results indicate that the buccal and lingual cortical plates do not experience notable stresses. However, when the occlusal contact areas are considered, the disto‐lingual superior third of the mandible experiences high tensile stresses, while the medio‐lingual cortical bone is subjected to high compressive stresses. Developing a more realistic loading scenario leads to better models to understand the relationship between masticatory function and mandibular shape and structures. Am J Phys Anthropol, 2012. © 2011 Wiley Periodicals, Inc.     

A customized fixation plate with novel structure designed by topological optimization for mandibular angle fracture based on finite element analysis

Background

The purpose of this study was to design a customized fixation plate for mandibular angle fracture using topological optimization based on the biomechanical properties of the two conventional fixation systems, and compare the results of stress, strain and displacement distributions calculated by finite element analysis (FEA).

Methods

A three-dimensional (3D) virtual mandible was reconstructed from CT images with a mimic angle fracture and a 1 mm gap between two bone segments, and then a FEA model, including volume mesh with inhomogeneous bone material properties, three loading conditions and constraints (muscles and condyles), was created to design a customized plate using topological optimization method, then the shape of the plate was referenced from the stress concentrated area on an initial part created from thickened bone surface for optimal calculation, and then the plate was formulated as “V” pattern according to dimensions of standard mini-plate finally. To compare the biomechanical behavior of the “V” plate and other conventional mini-plates for angle fracture fixation, two conventional fixation systems were used: type A, one standard mini-plate, and type B, two standard mini-plates, and the stress, strain and displacement distributions within the three fixation systems were compared and discussed.

Results

The stress, strain and displacement distributions to the angle fractured mandible with three different fixation modalities were collected, respectively, and the maximum stress for each model emerged at the mandibular ramus or screw holes. Under the same loading conditions, the maximum stress on the customized fixation system decreased 74.3, 75.6 and 70.6% compared to type A, and 34.9, 34.1, and 39.6% compared to type B. All maximum von Mises stresses of mandible were well below the allowable stress of human bone, as well as maximum principal strain. And the displacement diagram of bony segments indicated the effect of treatment with different fixation systems.

Conclusions

The customized fixation system with topological optimized structure has good biomechanical behavior for mandibular angle fracture because the stress, strain and displacement within the plate could be reduced significantly comparing to conventional “one mini-plate” or “two mini-plates” systems. The design methodology for customized fixation system could be used for other fractures in mandible or other bones to acquire better mechanical behavior of the system and improve stable environment for bone healing. And together with SLM, the customized plate with optimal structure could be designed and fabricated rapidly to satisfy the urgent time requirements for treatment.

     

Influence of Head Motion on the Accuracy of 3D Reconstruction with Cone-Beam CT: Landmark Identification Errors in Maxillofacial Surface Model

PurposeThe purpose of this study was to investigate the influence of head motion on the accuracy of three-dimensional (3D) reconstruction with cone-beam computed tomography (CBCT) scan.ResultsRendered surface models with head motion were similar to the control model in appearance; however, the landmark identification errors showed larger values in models with head motion than in the control. In particular, the Porion in the horizontal rotation models presented statistically significant differences (P < .05). Statistically significant difference in the errors between the right and left side landmark was present in the left side rotation which was opposite direction to the scanner rotation (P < .05).ConclusionsPatient movement during CBCT scan might cause landmark identification errors on the 3D surface model in relation to the direction of the scanner rotation. Clinicians should take this into consideration to prevent patient movement during CBCT scan, particularly horizontal movement.     

Dosimetric characterization of 3D printed bolus at different infill percentage for external photon beam radiotherapy

Background and purpose3D printing is rapidly evolving and further assessment of materials and technique is required for clinical applications. We evaluated 3D printed boluses with acrylonitrile butadiene styrene (ABS) and polylactide (PLA) at different infill percentage.Material and methodsA low-cost 3D printer was used. The influence of the air inclusion within the 3D printed boluses was assessed thoroughly both with treatment planning system (TPS) and with physical measurements. For each bolus, two treatment plans were calculated with Monte Carlo algorithm, considering the computed tomography (CT) scan of the 3D printed bolus or modelling the 3D printed bolus as a virtual bolus structure with a homogeneous density. Depth dose measurements were performed with Gafchromic films.ResultsHigh infill percentage corresponds to high density and high homogeneity within bolus material. The approximation of the bolus in the TPS as a homogeneous material is satisfying for infill percentages greater than 20%. Measurements performed with PLA boluses are more comparable to the TPS calculated profiles. For boluses printed at 40% and 60% infill, the discrepancies between calculated and measured dose distribution are within 5%.Conclusions3D printing technology allows modulating the shift of the build-up region by tuning the infill percentage of the 3D printed bolus in order to improve superficial target coverage.     

Fatigue creep damage at the cement–bone interface: An experimental and a micro-mechanical finite element study

The goal of this study was to quantify the micromechanics of the cement–bone interface under tensile fatigue loading using finite element analysis (FEA) and to understand the underlying mechanisms that play a role in the fatigue behavior of this interface. Laboratory cement–bone specimens were subjected to a tensile fatigue load, while local displacements and crack growth on the specimen's surface were monitored. FEA models were created from these specimens based upon micro-computed tomography data. To accurately model interfacial gaps at the interface between the bone and cement, a custom-written erosion algorithm was applied to the bone model. A fatigue load was simulated in the FEA models while monitoring the local displacements and crack propagation. The results showed the FEA models were able to capture the general experimental creep damage behavior and creep stages of the interface. Consistent with the experiments, the majority of the deformation took place at the contact interface. Additionally, the FEA models predicted fatigue crack patterns similar to experimental findings. Experimental surface cracks correlated moderately with FEA surface cracks (r²=0.43), but did not correlate with the simulated crack volume fraction (r²=0.06). Although there was no relationship between experimental surface cracks and experimental creep damage displacement (r²=0.07), there was a strong relationship between the FEA crack volume fraction and the FEA creep damage displacement (r²=0.76). This study shows the additional value of FEA of the cement–bone interface relative to experimental studies and can therefore be used to optimize its mechanical properties.     

Algorithm and validation of a computer method for quantifying attachment locus of glenohumeral ligament in vivo

The aim of this work is to validate an algorithm that quantifies the locus of glenohumeral ligament (GHL) attachments on glenohumeral joint (GHJ) bones.

A computed tomography scan of a GHJ was segmented to reconstruct the humerus, scapula, anatomical neck (AN) and glenoid rim (GR) into 3D meshes of interconnecting nodal vectors. These were applied to construct a ‘clock face’ coordinate system in which 3 o'clock points anteriorly.

Based on the assigned clock face coordinate frame and the fitted plane, the error between the fitted plane and the actual bony node was quantified through manual data extraction. This was tested on 50 specimens.

Mean algorithm quantification errors for GHL attachments were 4.8 (SD 2.2 mm) and 4.5 mm (1.7 mm) for the humerus and glenoid, respectively. Further studies would apply this to investigate GHL length changes during function and may suggest how these structures should be handled during surgical repairs.     

Comparison between two low profile attachments for implant mandibular overdentures

ObjectiveIn this research it was aimed to evaluate stress distribution on the implants supporting a complete overdenture in addition to compare between two different types of low-profile attachments for implant-retained mandibular overdenture with two techniques (with/without using connecting bar).Materials and methodsTwo 3D finite element models were constructed simulating supported lower complete overdenture with two implants and with two implants and bar. Where, models components were modeled in 3D on commercial general purpose CAD/CAM software. Four runs were carried out, two runs on each model, as linear static analysis.ResultsUsing bar is generally preferred for mucosa and cortical bone, while its effect can be considered as negligible on overdenture. On the other hand, it slightly increases the stresses on spongy bone. Using bar ensures the same level of energy transfer to the spongy bone and increases its maximum Von Mises stresses by about 50%. In addition, increase in maximum Von Mises stress was noticed by about 1% on cortical bone.ConclusionUsing bar is not recommended for patients with flat ridge.     

Creation of accurate 3D models of harbor porpoises (Phocoena phocoena) using 3D photogrammetry

Creating accurate 3D models of marine mammals is valuable for assessment of body condition, computational fluids dynamics models of locomotion, and for education. However, the methods for creating 3D models are not well-developed. We used photography and video to create 3D photogrammetry models of harbor porpoises (Phocoena phocoena). We accessed one live adult female (155.5 cm total length), and two dead animals, one juvenile (110 cm total length) and one calf (88 cm total length). We accessed the two dead individuals through a stranding network in Germany, and the live individual through the Fjord and Baelt research center in Denmark. For all porpoises, we used still photographs from hand-held cameras, drone video, and synchronized GoPro videos to create 3D photogrammetric models. We used Blender software, and other 3D reconstruction software, to recreate the 3D body meshes, and confirmed the accuracy of each of the 3D body meshes by comparing digital measures on the 3D models to original measures taken on the specimens. We also provide a colored, animated version of the live harbor porpoise for educational purposes. These open-access 3D models can be used to develop methods to study body morphometrics and condition, and to study bioenergetics and locomotion costs.     

Ethnic-specific differences in vitamin d status is associated with adiposity

Background

Low circulating 25 hydroxyvitamin D [25(OH)D] concentrations are common in obesity (BMI ≥30 kg/m²) and a negative relationship with body fat distribution has recently been reported. Ethnic-specific differences in body fat distribution have been described with South Asians are reported to have greater visceral adipose tissue (VAT), which could influence circulating 25(OH)D concentrations. The objective of this study is to investigate the relationship between plasma 25(OH)D, adiposity, and body fat distribution in Europeans and South Asians.

Methods/Principal Findings

187 Europeans and 192 South Asians were assessed for demographics, anthropometrics, and plasma 25(OH)D concentrations. Subcutaneous adipose tissue (SAT) and VAT were quantified by CT scan, and percent body fat by DEXA. Data were assessed by general linear models. South Asians had lower (P<0.001) plasma 25(OH)D concentrations and higher VAT (P = 0.04) than Europeans. Plasma 25(OH)D concentrations were negatively (P<0.05) associated with BMI, waist circumference, percent body fat, total adipose tissue, VAT, and SAT in unadjusted models and negatively (P<0.05) associated with VAT, SAT, and percent body fat after adjusting for BMI, ethnicity, age, and season of blood collection in males and females. When percent body fat, VAT, and SAT were included in the same model, only VAT remained negatively (P<0.05) associated with plasma 25(OH)D concentrations. Ethnicity remained significant in all models (P<0.001).

Conclusion

Compared to other adipose tissue compartments, VAT may have a distinct role in determining plasma 25(OH)D concentrations, which may account for the lower levels in South Asians.