Accuracy of femur reconstruction from sparse geometric data using a statistical shape model

Sparse geometric information from limited field-of-view medical images is often used to reconstruct the femur in biomechanical models of the hip and knee. However, the full femur geometry is needed to establish boundary conditions such as muscle attachment sites and joint axes which define the orientation of joint loads. Statistical shape models have been used to estimate the geometry of the full femur from varying amounts of sparse geometric information. However, the effect that different amounts of sparse data have on reconstruction accuracy has not been systematically assessed. In this study, we compared shape model and linear scaling reconstruction of the full femur surface from varying proportions of proximal and distal partial femur geometry in combination with morphometric and landmark data. We quantified reconstruction error in terms of surface-to-surface error as well as deviations in the reconstructed femur’s anatomical coordinate system which is important for biomechanical models. Using a partial proximal femur surface, mean shape model-based reconstruction surface error was 1.8 mm with 0.15° or less anatomic axis error, compared to 19.1 mm and 2.7–5.6° for linear scaling. Similar results were found when using a partial distal surface. However, varying amounts of proximal or distal partial surface data had a negligible effect on reconstruction accuracy. Our results show that given an appropriate set of sparse geometric data, a shape model can reconstruct full femur geometry with far greater accuracy than simple scaling.     

Development and validation of a generic 3D model of the distal femur

The development and validation of a virtual generic 3D model of the distal femur using computer graphical methods is presented. The synthesis of the generic model requires the following steps: acquisition of bony 3D morphology using standard computed tomography (CT) imaging; alignment of 3D models reconstructed from CT images with a common coordinate system; computer graphical sectioning of the models; extraction of bone contours from the image sections; combining and averaging of extracted contours; and 3D reconstruction of the averaged contours.

The generic models reconstructed from the averaged contours of six cadaver femora were validated by comparing their surface geometry on a point to point basis with that of the CT reconstructed reference models. The mean errors ranged from 0.99 to 2.5 mm and were in agreement with the qualitative assessment of the models.     

Investigation of 3D glenohumeral displacements from 3D reconstruction using biplane X-ray images: Accuracy and reproducibility of the technique and preliminary analysis in rotator cuff tear patients

Rotator cuff (RC) tears may be associated with increased glenohumeral instability; however, this instability is difficult to quantify using currently available diagnostic tools. Recently, the three-dimensional (3D) reconstruction and registration method of the scapula and humeral head, based on sequences of low-dose biplane X-ray images, has been proposed for glenohumeral displacement assessment. This research aimed to evaluate the accuracy and reproducibility of this technique and to investigate its potential with a preliminary application comparing RC tear patients and asymptomatic volunteers. Accuracy was assessed using CT scan model registration on biplane X-ray images for five cadaveric shoulder specimens and showed differences ranging from 0.6 to 1.4 mm depending on the direction of interest. Intra- and interobserver reproducibility was assessed through two operators who repeated the reconstruction of five subjects three times, allowing defining 95% confidence interval ranging from ±1.8 to ±3.6 mm. Intraclass correlation coefficient varied between 0.84 and 0.98. Comparison between RC tear patients and asymptomatic volunteers showed differences of glenohumeral displacements, especially in the superoinferior direction when shoulder was abducted at 20° and 45°. This study thus assessed the accuracy of the low-dose 3D biplane X-ray reconstruction technique for glenohumeral displacement assessment and showed potential in biomechanical and clinical research.     

3D虚拟复原精度的差异对头骨测量数值的影响——以Mimics软牛为例

近年来,根据三维软件虚拟复原的头骨来获取测量数据的方法被越来越多地应用在古生物,特别是古人类学的研究中,然而对于三维软件不同精度虚拟复原的头骨,其测量数值是否有差异,研究者并不是很清楚。本文以Mimics软件为例,根据其复原模型简化规则,选择未精简的最佳精度模型作为标准进行配对t检验或非参数检验,通过对43例云南人头骨的顶骨矢状弦长、颅周长、头盖部面积、乳突小房表面积、颅容量、乳突小房体积等六个测量项目的对比和分析,对Mimics软件低、中、高、最佳四种精度3D虚拟复原头骨间的测量差异进行了研究。结果显示：颅周长、头盖部面积、颅容量、乳突小房体积四项的所有简化精度模型的测量数据均与最佳精度模型测量数据的差异具有显著性;而除高精度组外,顶骨矢状弦长及乳突小房表面积的其余精度组测量数据均与最佳精度组差异具有显著性;此外,顶骨矢状弦长、颅周长、头盖部面积、颅容量的简化精度与最佳精度的测量差异占比均小于3%.而乳突小房表面积的低精度与最佳精度测量差异占比可超过50%,乳突小房体积的低精度与最佳精度测量差异占比可超过120%。这一结果提示我们,在测量Mimics复原的三维模型时,体量大差异小的测量项可以在较低精度的复原模型上进行测量;而对头骨内部腔窦这样体量小表面粗糙的结构,复原模型的精度选择及测量数据比较需要格外谨慎。     

Fast 3D reconstruction of the lower limb using a parametric model and statistical inferences and clinical measurements calculation from biplanar X-rays

In clinical routine, lower limb analysis relies on conventional X-ray (2D view) or computerised tomography (CT) Scan (lying position). However, these methods do not allow 3D analysis in standing position. The aim of this study is to propose a fast and accurate 3D-reconstruction-method based on parametric models and statistical inferences from biplanar X-rays with clinical measurements' (CM) assessment in standing position for a clinical routine use. For the reproducibility study, the 95% CI was under 2.7° for all lower limbs' angular measurements except for tibial torsion, femoral torsion and tibiofemoral rotation ( < 5°). The 95% CI were under 2.5 mm for lower limbs' lengths and 1.5 to 3° for the pelvis' CM. Comparisons between X-rays and CT-scan based 3D shapes in vitro showed mean differences of 1.0 mm (95% CI = 2.4 mm). Comparisons of 2D lower limbs' and 3D pelvis' CM between standing ‘Shifted-Feet’ and ‘Non-Shifted-Feet’ position showed means differences of 0.0 to 1.4°. Significant differences were found only for pelvic obliquity and rotation. The reconstruction time was about 5 min.     

Microcirculation in the murine liver: a computational fluid dynamic model based on 3D reconstruction from in vivo microscopy

The liver is organized in hexagonal functional units – termed lobules – characterized by a rather peculiar blood microcirculation, due to the presence of a tangled network of capillaries – termed sinusoids. A better understanding of the hemodynamics that governs liver microcirculation is relevant to clinical and biological studies aimed at improving our management of liver diseases and transplantation.Herein, we built a CFD model of a 3D sinusoidal network, based on in vivo images of a physiological mouse liver obtained with a 2-photon microscope. The CFD model was developed with Fluent 16.0 (ANSYS Inc., Canonsburg, PA), particular care was taken in imposing the correct boundary conditions representing a physiological state. To account for the remaining branches of the sinusoids, a lumped parameter model was used to prescribe the correct pressure at each outlet. The effect of an adhered cell on local hemodynamics is also investigated for different occlusion degrees.The model here proposed accurately reproduces the fluid dynamics in a portion of the sinusoidal network in mouse liver. Mean velocities and mass flow rates are in agreement with literature values from in vivo measurements. Our approach provides details on local phenomena, hardly described by other computational studies, either focused on the macroscopic hepatic vasculature or based on homogeneous porous medium model.     

Genome-wide variants of Eurasian facial shape differentiation and a prospective model of DNA based face prediction

It is a long-standing question as to which genes define the characteristic facial features among different ethnic groups. In this study, we use Uyghurs, an ancient admixed population to query the genetic bases why Europeans and Han Chinese look different. Facial traits were analyzed based on high-dense 3D facial images; numerous biometric spaces were examined for divergent facial features between European and Han Chinese, ranging from inter-landmark distances to dense shape geometrics. Genome-wide association studies(GWAS) were conducted on a discovery panel of Uyghurs. Six significant loci were identified, four of which, rs1868752, rs118078182, rs60159418 at or near UBASH3B, COL23A1, PCDH7 and rs17868256 were replicated in independent cohorts of Uyghurs or Southern Han Chinese. A prospective model was also developed to predict 3D faces based on top GWAS signals and tested in hypothetic forensic scenarios.     

Statistical shape and appearance models for fast and automated estimation of proximal femur fracture load using 2D finite element models

Estimating the risk of osteoporotic fractures is an important diagnostic step that needs to be taken before medicinal treatment. Densitometry-based criteria are normally used in clinical practice for this purpose. However, densitometry-based techniques could not explain all low-energy fractures. As patient-specific finite element (FE) models allow for consideration of other parameters (e.g. load conditions) that are known to be associated with fracture, they are considered promising candidates for more accurate fracture risk estimation. Nevertheless, they are often time consuming, expensive, and complex to build and may need the type of expertise that is not normally available in clinical settings. In this study, we report the development of an automated platform for estimating proximal femur fracture loads using patient-specific 2D FE models generated using dual-energy x-ray absorptiometry (DXA) scans. First, a statistical shape and appearance model (SSAM) is built using DXA scans of patients screened for osteoporosis following a low energy fracture. SSAM is then used together with Active Appearance Models (AAM) for automated segmentation of the proximal femur from new unseen DXA scans. The mean point-to-curve error of the automated procedure, i.e. 1.2–1.4 mm, is shown to be only slightly larger than the intra-observer variability of manual segmentation, i.e. 1.0 mm. Moreover, the developed platform automatically meshes the segmented shape, assigns density-based mechanical properties, assigns loads and boundary conditions, submits the 2D FE model for solution, and performs post-processing of the 2D FE simulation data to determine fracture loads. The fracture loads predicted using the manually generated and automatically generated 2D FE models are shown to be very close with a mean difference of around 8.8%. Repeated measures ANOVA showed no significant differences between the fracture loads calculated using FE models manually generated by three independent observers and those calculated using the automatically generated FE models (p>0.05).     

人源HDAC1/2-RbAp46/48核心复合物三维空间构象的电子显微分析

HDAC1、HDAC2和RbAp46、RbAp48是许多重要功能复合物(如NuRD、Sin3等)的核心亚基.这4个亚基在空间上相互作用,形成一个具有去乙酰化酶活性的核心复合物.但该核心复合物的三维空间构象及其对去乙酰化、染色质重塑等功能的可能影响还所知甚少.本研究中,我们包装了含4个亚基的杆状病毒,利用昆虫细胞表达、纯化了HDAC1/2-RbAp46/48核心复合物.在此基础上,利用电子显微镜单颗粒分析方法对该去乙酰化酶核心复合物的三维结构进行了初步解析.结果表明,HDAC1、HDAC2、RbAp46和RbAp48可以形成一个较为稳定均一的复合物,但该复合物中各个亚基并不是以单拷贝、等比例形式存在的.该核心复合物呈现一个非对称的鞍型结构,其背部隆起,大致形成一个三角形,两边分别有一大一小的两翼,两翼中间有个凹槽,直径大约为6 nm,推测为该核心复合物与核小体的结合位置.本研究结果为了解HDAC1/2-RbAp46/48去乙酰化酶复合物各亚基的空间结构组成、与核小体和染色质的可能相互作用以及研究去乙酰化酶活性的作用机理等提供了有益的信息.     

Accuracy of computer-aided geometric 3D reconstruction based on histological serial microgrinding preparation

Purpose: For our research on computer-optimised and automated cochlear implant surgery, we pursue a model-based approach to overcome the limitations of currently available clinical imaging modalities. A serial cross section preparation procedure has been developed and evaluated concerning accuracy to serve for modelling of a digital anatomic atlas to make delicate soft tissue structures available for pre-operative planning.

Methods: A special grinding tool was developed allowing the setting of a specific amount of abrasion as equidistant slice thickness was considered a crucial step. Additionally, each actual abrasion was accurately measured and used during three-dimensional reconstruction of the serial cross-sectional images obtained via digital photo documentation after each microgrinding step. A well-known reference object was prepared using this procedure and evaluated in terms of accuracy.

Results: Reconstruction of the whole sample was achieved with an error less than 0.4%, and the edge lengths in the direction of abrasion could be reconstructed with an average error of 0.6 ± 0.3 mm; both prove the realisation of equidistant abrasion. Using artificial registration fiducials and a custom-made algorithm for image alignment, parallelism and rectangularity could be preserved with average errors less than 0.4° ± 0.3°.

Conclusion: We present a systematic, practicable and reliable method for the geometrically accurate reconstruction of anatomical structures, which is especially suitable for the middle and inner ear anatomy including soft tissue structures. For the first time, the quality of such a reconstruction process has been quantified and successfully proven for its usability.     

Image based weighted center of proximity versus directly measured knee contact location during simulated gait

To understand the mechanical consequences of knee injury requires a detailed analysis of the effect of that injury on joint contact mechanics during activities of daily living. Three-dimensional (3D) knee joint geometric models have been combined with knee joint kinematics to dynamically estimate the location of joint contact during physiological activities—using a weighted center of proximity (WCoP) method. However, the relationship between the estimated WCoP and the actual location of contact has not been defined. The objective of this study was to assess the relationship between knee joint contact location as estimated using the image-based WCoP method, and a directly measured weighted center of contact (WCoC) method during simulated walking. To achieve this goal, we created knee specific models of six human cadaveric knees from magnetic resonance imaging. All knees were then subjected to physiological loads on a knee simulator intended to mimic gait. Knee joint motion was captured using a motion capture system. Knee joint contact stresses were synchronously recorded using a thin electronic sensor throughout gait, and used to compute WCoC for the medial and lateral plateaus of each knee. WCoP was calculated by combining knee kinematics with the MRI-based knee specific model. Both metrics were compared throughout gait using linear regression. The anteroposterior (AP) location of WCoP was significantly correlated with that of WCoC on both tibial plateaus in all specimens (p<0.01, 95% confidence interval of Pearson?s coefficient r>0), but the correlation was not significant in the mediolateral (ML) direction for 4/6 knees (p>0.05). Our study demonstrates that while the location of joint contact obtained from 3D knee joint contact model, using the WCoP method, is significantly correlated with the location of actual contact stresses in the AP direction, that relationship is less certain in the ML direction.     

A strategy for determining the orientations of refractory particles for reconstruction from cryo-electron micrographs with particular reference to round, smooth-surfaced, icosahedral viruses

Cryo-electron microscopy and three-dimensional image reconstruction are powerful tools for analyzing icosahedral virus capsids at resolutions that now extend below 1 nm. However, the validity of such density maps depends critically on correct identification of the viewing geometry of each particle in the data set. In some cases-for example, round capsids with low surface relief-it is difficult to identify orientations by conventional application of the two most widely used approaches-"common lines" and model-based iterative refinement. We describe here a strategy for determining the orientations of such refractory specimens. The key step is to determine reliable orientations for a base set of particles. For each particle, a list of candidate orientations is generated by common lines: correct orientations are then identified by computing a single-particle reconstruction for each candidate and then systematically matching their reprojections with the original images by visual criteria and cross-correlation analysis. This base set yields a first-generation reconstruction that is fed into the model-based procedure. This strategy has led to the structural determination of two viruses that, in our hands, resisted solution by other means.     

A topology optimization based model of bone adaptation

A novel topology optimization model based on homogenization methods was developed for predicting bone density distribution and anisotropy, assuming the bone structure to be a self-optimizing biological material which maximizes its own structural stiffness. The feasibility and efficiency of this method were tested on a 2D model for a proximal femur under single and multiple loading conditions. The main aim was to compute homogenized optimal designs using an optimal laminated microstructure. The computational results showed that high bone density levels are distributed along the diaphysis and form arching struts within the femoral head. The pattern of bone density distribution and the anisotropic bone behavior predicted by the model in the multiple load case were both in good agreement with the structural architecture and bone density distribution occurring in natural femora. This approach provides a novel means of understanding the remodeling processes involved in fracture repair and the treatment of bone diseases.     

Microdosimetry in 3D realistic mitochondria phantoms: Geant4 Monte Carlo tracking of 250 keV photons in phantoms reconstructed from microscopic images

Mitochondria are considered to be sensitive radiation targets since they control processes vital to the cell’s functioning. These organelles are starting to get attention and some studies are investigating the radiation dose inside them. In previous studies, mitochondria are represented as simple ellipsoids inside the cell not taking into consideration the complexity of their shape. In this study, realistic phantoms are built based on deconvolved widefield fluorescent microscopic images of the mitochondrial networks of fibroblast cells. The phantoms are imported into Geant4 as tessellated volumes taking into account the geometrical complexity of these organelles. Irradiation with 250 keV photons is performed and the lineal energy is calculated. The lineal energy distributions inside the produced phantoms are compared with those calculated inside simple volumes, a sphere and an ellipsoid, where the effect of the shape and volume is clearly seen on lineal energies.     

The application of 3D printing patient specific instrumentation model in total knee arthroplasty

The application of 3D printing patient specific instrumentation model in total knee arthroplasty was explored to improve the operative accuracy and safety of artificial total knee arthroplasty. In this study, a total of 52 patients who need knee replacement were selected as the study objects, and 52 patients were divided into experimental group and control group. First, the femoral mechanical-anatomical angle (FMAA), lateral femoral angle (LFA), hip-knee-ankle angle (HKA), femorotibial angle (FTA) of research objects in both groups were measured. Then, the blood loss during the operations, drainage volume after operations, total blood loss, hidden blood loss, and hemoglobin decrease of the experiment group and the control group were measured and calculated. Finally, the postoperative outcomes of patients who underwent total knee arthroplasty were evaluated. The results showed that before the operations, in the PSI group, the femoral mechanical-anatomical angle (FMAA) was (6.9 ± 2.4)°, the lateral femoral angle (LFA) was (82.4 ± 1.6)°, the hip-knee-ankle angle (HKA) was (166.4 ± 1.4)°, and the femorotibial angle (FTA) was (179.5 ± 7.3)°. In the CON group, the FMAA was (5.8 ± 2.4)°, the LFA was (81.3 ± 2.1)°, the HKA was (169.5 ± 1.9)°, and the FTA was (185.4 ± 5.4)°. The differences in these data between the two groups were not statistically significant (P > 0.05). After the operations, in the PSI group, the total blood loss, the hidden blood loss, and the hemoglobin (Hb) decrease were respectively (420.2 ± 210.5), (240.5 ± 234.5), and (1.7 ± 0.9); in the CON group, the total blood loss, the hidden blood loss, and the Hb decrease were respectively (782.1 ± 340.4), (450.9 ± 352.6), and (2.9 ± 1.0). These data of both groups were statistically significant (P < 0.05). Therefore, it can be seen that the 3D printing patient specific instrumentation model can effectively simulate the lower limb coronal force line and was highly consistent of the preoperative software simulation plan. In addition, the random interviews of patients who underwent total knee arthroplasty showed that the knees of patients had recovered well. The application of 3D printing patient specific instrumentation model in artificial total knee arthroplasty can effectively improve the operative accuracy and safety, and the clinical therapeutic effects were significant.     

Predicted structure model of Bungarotoxin from Bungarus fasciatus snake

Snake venoms are cocktails comprising combinations of different proteins, peptides, enzymes and toxins. Snake toxins have diverse characteristics having different molecular configuration, structure and mode of action. Many toxins derived from snake venom have distinct pharmacological activities. Venom from Bungarus fasciatus (commonly known as banded krait) is a species of elapid snake found on the South East Asia and Indian sub-continent, mainly contains neurotoxins. Beta bungartotoxin is the major fraction of Bungarus venom and particularly act pre-synaptically by obstructing neurotransmitter release. This toxin in other snake species functionally forms a heterodimer containing two different subunits (A and B). Dimerization of these two chains is a pre-requisite for the proper functionality of this protein. However, B. fasciatus bungartotoxin contains only B chain and their structural orientation in yet to be resolved. Therefore, it is of interest to describe the predicted structure model of the toxin for functional insights. In this work we analyzed the neurotoxic nature, their alignments, secondary and three dimensional structures, functions, active sites and stability with the help of different bioinformatical tools. A comprehensive analysis of the predicted model provides approaching to the functional interpretation of its molecular action.     

Paurodontella parapitica n. sp. (Nematoda: Hexatylina,Sphaerularioidea) from Kermanshah Province,Western Iran

Paurodontella parapitica n. sp., collected from the rhizosphere of an apple tree in Kermanshah province, western Iran, is described. The new species is characterized by a body length of 505 to 723 µm (females) and 480 to 600 µm (males), lip region continuous by depression; 7 to 8 μm broad, 3 to 4 µm high, stylet length 7 to 9 µm or 1 to 1.3 times the lip region diameter, short postuterine sac of 4 to 6 μm long, lateral fields with five to six incisures; outer incisures crenated and inner incisures weakly crenated, excretory pore situated 90 to 100 µm from anterior end; functional males common in the population, with spicules 24 to 26 μm long. Tail of both sexes similar, almost straight and elongate-conoid. The new species resembles in morphology and morphometrics to four known species of the genus, namely P. apitica, P. minuta, P. myceliophaga, and P. sohailai. The results of phylogenetic analyses based on sequences of D2/D3 expansion region of 28S rRNA gene revealed this genus is polyphyletic in four different clades in Tylenchid.