Mesh-based vs. Image-based Statistical Appearance Model of the Human Femur: a Preliminary Comparison Study for the Creation of Finite Element Meshes

Statistical models have been recently introduced in computational orthopaedics to investigate the bone mechanical properties across several populations. A fundamental aspect for the construction of statistical models concerns the establishment of accurate anatomical correspondences among the objects of the training dataset. Various methods have been proposed to solve this problem such as mesh morphing or image registration algorithms. The objective of this study is to compare a mesh-based and an image-based statistical appearance model approaches for the creation of nite element(FE) meshes. A computer tomography (CT) dataset of 157 human left femurs was used for the comparison. For each approach, 30 finite element meshes were generated with the models. The quality of the obtained FE meshes was evaluated in terms of volume, size and shape of the elements. Results showed that the quality of the meshes obtained with the image-based approach was higher than the quality of the mesh-based approach. Future studies are required to evaluate the impact of this finding on the final mechanical simulations.

Statistical model-based segmentation of the proximal femur in digital antero-posterior (AP) pelvic radiographs.

PURPOSE    Segmentation of the proximal femur in digital antero-posterior (AP) pelvic radiographs is required to create a three-dimensional model of the hip joint for use in planning and treatment. However, manually extracting the femoral contour is tedious and prone to subjective bias, while automatic segmentation must accommodate poor image quality, anatomical structure overlap, and femur deformity. A new method was developed for femur segmentation in AP pelvic radiographs. METHODS    Using manual annotations on 100 AP pelvic radiographs, a statistical shape model (SSM) and a statistical appearance model (SAM) of the femur contour were constructed. The SSM and SAM were used to segment new AP pelvic radiographs with a three-stage approach. At initialization, the mean SSM model is coarsely registered to the femur in the AP radiograph through a scaled rigid registration. Mahalanobis distance defined on the SAM is employed as the search criteria for each annotated suggested landmark location. Dynamic programming was used to eliminate ambiguities. After all landmarks are assigned, a regularized non-rigid registration method deforms the current mean shape of SSM to produce a new segmentation of proximal femur. The second and third stages are iteratively executed to convergence. RESULTS    A set of 100 clinical AP pelvic radiographs (not used for training) were evaluated. The mean segmentation error was [Formula: see text], requiring [Formula: see text] s per case when implemented with Matlab. The influence of the initialization on segmentation results was tested by six clinicians, demonstrating no significance difference. CONCLUSIONS    A fast, robust and accurate method for femur segmentation in digital AP pelvic radiographs was developed by combining SSM and SAM with dynamic programming. This method can be extended to segmentation of other bony structures such as the pelvis.

Automatic Scan Planning for Magnetic Resonance Imaging of the Knee Joint

Automatic scan planning for magnetic resonance imaging of the knee aims at defining an oriented bounding box around the knee joint from sparse scout images in order to choose the optimal field of view for the diagnostic images and limit acquisition time. We propose a fast and fully automatic method to perform this task based on the standard clinical scout imaging protocol. The method is based on sequential Chamfer matching of 2D scout feature images with a three-dimensional mean model of femur and tibia. Subsequently, the joint plane separating femur and tibia, which contains both menisci, can be automatically detected using an information-augmented active shape model on the diagnostic images. This can assist the clinicians in quickly defining slices with standardized and reproducible orientation, thus increasing diagnostic accuracy and also comparability of serial examinations. The method has been evaluated on 42 knee MR images. It has the potential to be incorporated into existing systems because it does not change the current acquisition protocol.

Statistical shape analysis via principal factor analysis

Statistical shape analysis techniques commonly employed in the medical imaging community, such as active shape models or active appearance models, rely on principal component analysis (PCA) to decompose shape variability into a reduced set of interpretable components. In this paper we propose principal factor analysis (PFA) as an alternative and complementary tool to PCA providing a decomposition into modes of variation that can be more easily interpretable, while still being a linear efficient technique that performs dimensionality reduction (as opposed to independent component analysis, ICA). The key difference between PFA and PCA is that PFA models covariance between variables, rather than the total variance in the data. The added value of PFA is illustrated on 2D landmark data of corpora callosa outlines. Then, a study of the 3D shape variability of the human left femur is performed. Finally, we report results on vector-valued 3D deformation fields resulting from non-rigid registration of ventricles in MRI of the brain.

Automatic multi-resolution shape modeling of multi-organ structures.

Point Distribution Models (PDM) are among the most popular shape description techniques and their usefulness has been demonstrated in a wide variety of medical imaging applications. However, to adequately characterize the underlying modeled population it is essential to have a representative number of training samples, which is not always possible. This problem is especially relevant as the complexity of the modeled structure increases, being the modeling of ensembles of multiple 3D organs one of the most challenging cases. In this paper, we introduce a new GEneralized Multi-resolution PDM (GEM-PDM) in the context of multi-organ analysis able to efficiently characterize the different inter-object relations, as well as the particular locality of each object separately. Importantly, unlike previous approaches, the configuration of the algorithm is automated thanks to a new agglomerative landmark clustering method proposed here, which equally allows us to identify smaller anatomically significant regions within organs. The significant advantage of the GEM-PDM method over two previous approaches (PDM and hierarchical PDM) in terms of shape modeling accuracy and robustness to noise, has been successfully verified for two different databases of sets of multiple organs: six subcortical brain structures, and seven abdominal organs. Finally, we propose the integration of the new shape modeling framework into an active shape-model-based segmentation algorithm. The resulting algorithm, named GEMA, provides a better overall performance than the two classical approaches tested, ASM, and hierarchical ASM, when applied to the segmentation of 3D brain MRI.

PSGL-1 is dispensible for the development of active experimental autoimmune encephalomyelitis in SJL/J mice

The adhesion molecule P-selectin glycoprotein ligand (PSGL)-1 has been suggested to be involved in the immunopathogenesis of multiple sclerosis (MS). However, in C57BL/6 mice PSGL-1 was found to be dispensible for the development of MOG(aa35-55)-induced experimental autoimmune encephalomyelitis (EAE), an animal model for MS. To study, if involvement of PSGL-1 to EAE pathogenesis can be observed in another common mouse model, we backcrossed PSGL-1(-/-) mice for at least 12 generations into the SJL/J background and compared PLP(aa139-151) induced EAE in PSGL-1(-/-) SJL/J mice versus wild-type SJL/J mice. Here, we demonstrate that PSGL-1(-/-) SJL/J mice exhibited EAE pathogenesis indistinguishable from wild-type SJL/J mice. Our present study underscores and emphasizes previous observations that PSGL-1 is dispensible for EAE pathogenesis.

Atypical presentation of a hormonally active adrenocortical tumor in an adolescent leading to delayed diagnosis

Adrenocortical tumors are rare in children and present with variable signs depending on the type of hormone excess. We herein describe the unusual presentation of a child with adrenocortical tumor and introduce the concept of in vitro chemosensitivity testing. CASE REPORT: A 10.5-year-old girl presented with hypertrichosis/hirsutism and weight loss. The weight loss and behavioral problems, associated with halted puberty and growth, led to the initial diagnosis of anorexia nervosa. However, subsequent weight gain but persisting arrest in growth and puberty and the appearance of central fat distribution prompted further evaluation. RESULTS AND FOLLOW-UP: 24h-urine free cortisol was elevated. Morning plasma ACTH was undetectable, while cortisol was elevated and circadian rhythmicity was absent. Thus a hormonally active adrenal cortical tumor (ACT) was suspected. On magnetic resonance imaging (MRI) a unilateral, encapsulated tumor was found which was subsequently removed surgically. Tissue was investigated histologically and for chemosensitivity in primary cell cultures. Although there were some risk factors for malignancy, the tumor was found to be a typical adenoma. Despite this histology, tumor cells survived in culture and were sensitive to cisplatin in combination with gemcitabine or paclitaxel. At surgery, the patient was started on hydrocortisone replacement which was unsuccessfully tapered over 3 months. Full recovery of the hypothalamus-pituitary-adrenal axis occurred only after 3 years. CONCLUSIONS: The diagnosis of a hormonally active adrenocortical tumor is often delayed because of atypical presentation. Cortisol replacement following unilateral tumor excision is mandatory and may be required for months or years. Individualized chemosensitivity studies carried out on primary cultures established from the tumor tissue itself may provide a tool in evaluating the effectiveness of chemotherapeutic drugs in the event that the adrenocortical tumor may prove to be carcinoma.

A protective antiarrhythmic role of ursodeoxycholic acid in an in vitro rat model of the cholestatic fetal heart

Intrahepatic cholestasis of pregnancy may be complicated by fetal arrhythmia, fetal hypoxia, preterm labor, and, in severe cases, intrauterine death. The precise etiology of fetal death is not known. However, taurocholate has been demonstrated to cause arrhythmia and abnormal calcium dynamics in cardiomyocytes. To identify the underlying reason for increased susceptibility of fetal cardiomyocytes to arrhythmia, we studied myofibroblasts (MFBs), which appear during structural remodeling of the adult diseased heart. In vitro, they depolarize rat cardiomyocytes via heterocellular gap junctional coupling. Recently, it has been hypothesized that ventricular MFBs might appear in the developing human heart, triggered by physiological fetal hypoxia. However, their presence in the fetal heart (FH) and their proarrhythmogenic effects have not been systematically characterized. Immunohistochemistry demonstrated that ventricular MFBs transiently appear in the human FH during gestation. We established two in vitro models of the maternal heart (MH) and FH, both exposed to increasing doses of taurocholate. The MH model consisted of confluent strands of rat cardiomyocytes, whereas for the FH model, we added cardiac MFBs on top of cardiomyocytes. Taurocholate in the FH model, but not in the MH model, slowed conduction velocity from 19 to 9 cm/s, induced early after depolarizations, and resulted in sustained re-entrant arrhythmias. These arrhythmic events were prevented by ursodeoxycholic acid, which hyperpolarized MFB membrane potential by modulating potassium conductance. CONCLUSION: These results illustrate that the appearance of MFBs in the FH may contribute to arrhythmias. The above-described mechanism represents a new therapeutic approach for cardiac arrhythmias at the level of MFB.

Adsorption behavior of redox-active suppressor additives: Combined electrochemical and STM studies

The redox chemistry and the related surface phase behavior of Safranine (SAF) and Janus Green B (JGB) have been studied by means of cyclic voltammetry in combination with in situ Scanning Tunneling Microscopy using HOPG (Highly Oriented Pyrolytic Graphite) and single crystalline Cu(1 0 0) as model substrates, both revealing different widths of the accessible potential windows. JGB and SAF serve as prototypical heterocyclic suppressor/leveler additives that are used for the metallization of 3D-TSVs (3D Through Silicon Vias) following a classical "leveling" concept. SAF can be considered as the reductive decomposition product of JGB that is formed at the copper/electrolyte interface upon electroplating. Both additives reveal a pronounced pH-dependent redox-chemistry with redox-transitions lying close to or even beyond the anodic limit of the copper potential window. Affected by these redox-processes are in particular the aromatic cores of those heterocycles that can be (quasi)reversibly reduced by a two electron transfer process within the potential window of copper. Therefore we identify the reduced form of those dyes as the active components for the suppressing/leveling effect in copper plating. STM data clearly shows a dye surface phase behavior that is crucially determined by its potential-dependent redox-chemistry. This will be exemplarily discussed for the SAF dye. On chloride-modified Cu(1 0 0) mono-reduced SAF forms a structurally well-defined monolayer of cationic stacking polymers. However, this coupled anion/cation layer reveals only minor suppressing capabilities with respect to the copper dissolution and deposition processes. Complete reduction of the aromatic heterocycle finally leads to the 3D precipitation of hydrophobic reaction products. 3D clusters of this SAF precipitate are discussed as the active structural motif for the suppressing effect of these dyes. (C) 2011 Elsevier Ltd. All rights reserved.

The S-1/S-2 exciton interaction in 2-pyridone center dot 6-methyl-2-pyridone: Davydov splitting, vibronic coupling, and vibronic quenching

The excitonic splitting between the S-1 and S-2 electronic states of the doubly hydrogen-bonded dimer 2-pyridone center dot 6-methyl-2-pyridone (2PY center dot 6M2PY) is studied in a supersonic jet, applying two-color resonant two-photon ionization (2C-R2PI), UV-UV depletion, and dispersed fluorescence spectroscopies. In contrast to the C-2h symmetric (2-pyridone) 2 homodimer, in which the S-1 <- S-0 transition is symmetry-forbidden but the S-2 <- S-0 transition is allowed, the symmetry-breaking by the additional methyl group in 2PY center dot 6M2PY leads to the appearance of both the S-1 and S-2 origins, which are separated by Delta(exp) = 154 cm(-1). When combined with the separation of the S-1 <- S-0 excitations of 6M2PY and 2PY, which is delta = 102 cm(-1), one obtains an S-1/S-2 exciton coupling matrix element of V-AB, el = 57 cm(-1) in a Frenkel-Davydov exciton model. The vibronic couplings in the S-1/S-2 <- S-0 spectrum of 2PY center dot 6M2PY are treated by the Fulton-Gouterman single-mode model. We consider independent couplings to the intramolecular 6a' vibration and to the intermolecular sigma' stretch, and obtain a semi-quantitative fit to the observed spectrum. The dimensionless excitonic couplings are C(6a') = 0.15 and C(sigma') = 0.05, which places this dimer in the weak-coupling limit. However, the S-1/S-2 state exciton splittings Delta(calc) calculated by the configuration interaction singles method (CIS), time-dependent Hartree-Fock (TD-HF), and approximate second-order coupled-cluster method (CC2) are between 1100 and 1450 cm(-1), or seven to nine times larger than observed. These huge errors result from the neglect of the coupling to the optically active intra-and intermolecular vibrations of the dimer, which lead to vibronic quenching of the purely electronic excitonic splitting. For 2PY center dot 6M2PY the electronic splitting is quenched by a factor of similar to 30 (i.e., the vibronic quenching factor is Gamma(exp) = 0.035), which brings the calculated splittings into close agreement with the experimentally observed value. The 2C-R2PI and fluorescence spectra of the tautomeric species 2-hydroxypyridine center dot 6-methyl-2-pyridone (2HP center dot 6M2PY) are also observed and assigned. (C) 2011 American Institute of Physics.

Intra-articular calcium phosphate cement: Its fate and impact on joint tissues in a rabbit model

Clinical application of injectable ceramic cement in comminuted fractures revealed penetration of the viscous paste into the joint space. Not much is known on the fate of this cement and its influence on articular tissues. The purpose of this experimental study was to assess these unknown alterations of joint tissues after intra-articular injection of cement in a rabbit knee. Observation periods were from 1 week up to 24 months, with three rabbits per group. Norian SRS cement was injected into one knee joint, the contralateral side receiving the same volume of Ringers' solution. Light microscopic evaluation of histologic sections was performed, investigating the appearance of the cement, inflammatory reactions, and degenerative changes of the articular surface. No signs of pronounced acute or chronic inflammation were visible. The injected cement was mainly found as a single particle, anterior to the cruciate ligaments. It became surrounded by synovial tissues within 4 weeks and showed signs of superficial resorption. In some specimens, bone formation was seen around the cement. Degeneration of the articular surface showed no differences between experimental and control side, and no changes over time became apparent. No major degenerative changes were induced by the injected cement. The prolonged presence of cement still seems to make it advisable to remove radiologically visible amounts from the joint space.

Risks of non-accidental mortality by baseline CD4+ T-cell strata in hepatitis-C-positive and -negative individuals initiating highly active antiretroviral therapy

BACKGROUND: Patients coinfected with hepatitis C virus (HCV) and HIV experience higher mortality rates than patients infected with HIV alone. We designed a study to determine whether risks for later mortality are similar for HCV-positive and HCV-negative individuals when subjects are stratified on the basis of baseline CD4+ T-cell counts. METHODS: Antiretroviral-naive individuals, who initiated highly active antiretroviral therapy (HAART) between 1996 and 2002 were included in the study. HCV-positive and HCV-negative individuals were stratified separately by baseline CD4+ T-cell counts of 50 cell/microl increments. Cox-proportional hazards regression was used to model the effect of these strata with other variables on survival. RESULTS: CD4+ T-cell strata below 200 cells/microl, but not above, imparted an increased relative hazard (RH) of mortality for both HCV-positive and HCV-negative individuals. Among HCV-positive individuals, after adjustment for baseline age, HIV RNA levels, history of injection drug use and adherence to therapy, only CD4+ T-cell strata of <50 cells/microl (RH=4.60; 95% confidence interval [CI] 2.72-7.76) and 50-199 cells/microl (RH=2.49; 95% CI 1.63-3.81) were significantly associated with increased mortality when compared with those initiating therapy at cell counts >500 cells/microl. The same baseline CD4+ T-cell strata were found for HCV-negative individuals. CONCLUSION: In a within-groups analysis, the baseline CD4+ T-cell strata that are associated with increased RHs for mortality are the same for HCV-positive and HCV-negative individuals initiating HAART. However, a between-groups analysis reveals a higher absolute mortality risk for HCV-positive individuals.

A coronary heart disease risk model for predicting the effect of potent antiretroviral therapy in HIV-1 infected men

BACKGROUND: Many HIV-infected patients on highly active antiretroviral therapy (HAART) experience metabolic complications including dyslipidaemia and insulin resistance, which may increase their coronary heart disease (CHD) risk. We developed a prognostic model for CHD tailored to the changes in risk factors observed in patients starting HAART. METHODS: Data from five cohort studies (British Regional Heart Study, Caerphilly and Speedwell Studies, Framingham Offspring Study, Whitehall II) on 13,100 men aged 40-70 and 114,443 years of follow up were used. CHD was defined as myocardial infarction or death from CHD. Model fit was assessed using the Akaike Information Criterion; generalizability across cohorts was examined using internal-external cross-validation. RESULTS: A parametric model based on the Gompertz distribution generalized best. Variables included in the model were systolic blood pressure, total cholesterol, high-density lipoprotein cholesterol, triglyceride, glucose, diabetes mellitus, body mass index and smoking status. Compared with patients not on HAART, the estimated CHD hazard ratio (HR) for patients on HAART was 1.46 (95% CI 1.15-1.86) for moderate and 2.48 (95% CI 1.76-3.51) for severe metabolic complications. CONCLUSIONS: The change in the risk of CHD in HIV-infected men starting HAART can be estimated based on typical changes in risk factors, assuming that HRs estimated using data from non-infected men are applicable to HIV-infected men. Based on this model the risk of CHD is likely to increase, but increases may often be modest, and could be offset by lifestyle changes.