Mikko S. Venäläinen

Quantitative Evaluation of the Mechanical Risks Caused by Focal Cartilage Defects in the Knee.

Focal cartilage lesions can proceed to severe osteoarthritis or remain unaltered even for years. A method to identify high risk defects would be of utmost importance to guide clinical decision making and to identify the patients that are at the highest risk for the onset and progression of osteoarthritis. Based on cone beam computed tomography arthrography, we present a novel computational model for evaluating changes in local mechanical responses around cartilage defects. Our model, based on data obtained from a human knee in vivo, demonstrated that the most substantial alterations around the defect, as compared to the intact tissue, were observed in minimum principal (compressive) strains and shear strains. Both strain values experienced up to 3-fold increase, exceeding levels previously associated with chondrocyte apoptosis and failure of collagen crosslinks. Furthermore, defects at the central regions of medial tibial cartilage with direct cartilage-cartilage contact were the most vulnerable to loading. Also locations under the meniscus experienced substantially increased minimum principal strains. We suggest that during knee joint loading particularly minimum principal and shear strains are increased above tissue failure limits around cartilage defects which might lead to osteoarthritis. However, this increase in strains is highly location-specific on the joint surface.

Effect of bone inhomogeneity on tibiofemoral contact mechanics during physiological loading

It is not known how inhomogeneous mechanical properties of bone affect contact mechanics and cartilage response during physiological loading of the knee joint. In this study, a finite element model of a cadaver knee joint was constructed based on quantitative computed tomography (QCT). The mechanical properties of bone were altered and their effect on tibiofemoral contact mechanics and cartilage stresses, strains and pore pressures were evaluated during the first 20% of stance. For this purpose, models with rigid, homogeneous and inhomogeneous bones were created. When bone was modeled to be rigid, the resulting contact pressures were substantially higher in the medial side of the joint, as compared to the non-rigid bones. Similar changes were revealed also in stresses, strains and pore pressures throughout the cartilage depth at the cartilage-cartilage contact area. Furthermore, the mechanical response of medial tibial cartilage was found to be highly dependent on the bone properties. When Young׳s modulus in the model with homogeneous bone was 5GPa, cartilage mechanical response approached to that of the model with inhomogeneous bone. Finally, when the apparent bone mineral densities were decreased globally in the inhomogeneous bone, stresses, strains and pore pressures were decreased at all layers of medial tibial cartilage. Similar changes were observed also in cartilage-cartilage contact area of the lateral compartment but with a lesser extent. These results indicate that during physiological loading Young׳s modulus of bone has a substantial influence on cartilage stresses and strains, especially in the medial compartment.

Importance of Material Properties and Porosity of Bone on Mechanical Response of Articular Cartilage in Human Knee Joint—A Two-Dimensional Finite Element Study

Mechanical behavior of bone is determined by the structure and intrinsic, local material properties of the tissue. However, previously presented knee joint models for evaluation of stresses and strains in joints generally consider bones as rigid bodies or linearly elastic solid materials. The aim of this study was to estimate how different structural and mechanical properties of bone affect the mechanical response of articular cartilage within a knee joint. Based on a cadaver knee joint, a two-dimensional (2D) finite element (FE) model of a knee joint including bone, cartilage, and meniscus geometries was constructed. Six different computational models with varying properties for cortical, trabecular, and subchondral bone were created, while the biphasic fibril-reinforced properties of cartilage and menisci were kept unaltered. The simplest model included rigid bones, while the most complex model included specific mechanical properties for different bone structures and anatomically accurate trabecular structure. Models with different porosities of trabecular bone were also constructed. All models were exposed to axial loading of 1.9 times body weight within 0.2 s (mimicking typical maximum knee joint forces during gait) while free varus-valgus rotation was allowed and all other rotations and translations were fixed. As compared to results obtained with the rigid bone model, stresses, strains, and pore pressures observed in cartilage decreased depending on the implemented properties of trabecular bone. Greatest changes in these parameters (up to -51% in maximum principal stresses) were observed when the lowest modulus for trabecular bone (measured at the structural level) was used. By increasing the trabecular bone porosity, stresses and strains were reduced substantially in the lateral tibial cartilage, while they remained relatively constant in the medial tibial plateau. The present results highlight the importance of long bones, in particular, their mechanical properties and porosity, in altering and redistributing forces transmitted through the knee joint.

Spatial variation of fixed charge density in knee joint cartilage from sodium MRI – Implication on knee joint mechanics under static loading

The effects of fixed charge density (FCD) and cartilage swelling have not been demonstrated on cartilage mechanics on knee joint level before. In this study, we present how the spatial and local variations of FCD affects the mechanical response of the knee joint cartilage during standing (half of the body weight, 13 minutes) using finite element (FE) modeling. The FCD distribution of tibial cartilage of an asymptomatic subject was determined using sodium (23Na) MRI at 7T and implemented into a 3-D FE-model of the knee joint (Subject-specific model, FCD: 0.18±0.08 mEq/ml). Tissue deformation in the Subject-specific model was validated against experimental, in vivo loading of the joint conducted with a MR-compatible compression device. For comparison, models with homogeneous FCD distribution (homogeneous model) and FCD distribution obtained from literature (literature model) were created. Immediately after application of the load (dynamic response), the variations in FCD had minor effects on cartilage stresses and strains. After 13 minutes of standing, the spatial and local variations in FCD had most influence on axial strains. In the superficial tibial cartilage in the Subject-specific model, axial strains were increased up to +13% due to smaller FCD (mean -11%), as compared to the homogeneous model. Compared to the literature model, those were decreased up to -18% due to greater FCD (mean +7%). The findings demonstrate that the spatial and local FCD variations in cartilage modulates strains in knee joint cartilage. Thereby, the results highlight the mechanical importance of site-specific content of proteoglycans in cartilage.

Workflow assessing the effect of gait alterations on stresses in the medial tibial cartilage - combined musculoskeletal modelling and finite element analysis

Knee osteoarthritis (KOA) is most common in the medial tibial compartment. We present a novel method to study the effect of gait modifications and lateral wedge insoles (LWIs) on the stresses in the medial tibial cartilage by combining musculoskeletal (MS) modelling with finite element (FE) analysis. Subject’s gait was recorded in a gait laboratory, walking normally, with 5° and 10° LWIs, toes inward (‘Toe in’), and toes outward (‘Toe out wide’). A full lower extremity MRI and a detailed knee MRI were taken. Bones and most soft tissues were segmented from images, and the generic bone architecture of the MS model was morphed into the segmented bones. The output forces from the MS model were then used as an input in the FE model of the subject’s knee. During stance, LWIs failed to reduce medial peak pressures apart from Insole 10° during the second peak. Toe in reduced peak pressures by −11% during the first peak but increased them by 12% during the second. Toe out wide reduced peak pressures by −15% during the first and increased them by 7% during the second. The results show that the work flow can assess the effect of interventions on an individual level. In the future, this method can be applied to patients with KOA.

Prediction of Adult Dyslipidemia Using Genetic and Childhood Clinical Risk FactorsCLINICAL PERSPECTIVE: The Cardiovascular Risk in Young Finns Study

Background— Dyslipidemia is a major modifiable risk factor for cardiovascular disease. We examined whether the addition of novel single-nucleotide polymorphisms for blood lipid levels enhances the prediction of adult dyslipidemia in comparison to childhood lipid measures. Methods and Results— Two thousand four hundred and twenty-two participants of the Cardiovascular Risk in Young Finns Study who had participated in 2 surveys held during childhood (in 1980 when aged 3–18 years and in 1986) and at least once in a follow-up study in adulthood (2001, 2007, and 2011) were included. We examined whether inclusion of a lipid-specific weighted genetic risk score based on 58 single-nucleotide polymorphisms for low-density lipoprotein cholesterol, 71 single-nucleotide polymorphisms for high-density lipoprotein cholesterol, and 40 single-nucleotide polymorphisms for triglycerides improved the prediction of adult dyslipidemia compared with clinical childhood risk factors. Adjusting for age, sex, body mass index, physical activity, and smoking in childhood, childhood lipid levels, and weighted genetic risk scores were associated with an increased risk of adult dyslipidemia for all lipids. Risk assessment based on 2 childhood lipid measures and the lipid-specific weighted genetic risk scores improved the accuracy of predicting adult dyslipidemia compared with the approach using only childhood lipid measures for low-density lipoprotein cholesterol (area under the receiver-operating characteristic curve 0.806 versus 0.811; P=0.01) and triglycerides (area under the receiver-operating characteristic curve 0.740 versus area under the receiver-operating characteristic curve 0.758; P<0.01). The overall net reclassification improvement and integrated discrimination improvement were significant for all outcomes. Conclusions— The inclusion of weighted genetic risk scores to lipid-screening programs in childhood could modestly improve the identification of those at highest risk of dyslipidemia in adulthood.

Serum Proteomic Profiling to Identify Biomarkers of Premature Carotid Atherosclerosis

To evaluate the presence of serum protein biomarkers associated with the early phases of formation of carotid atherosclerotic plaques, label-free quantitative proteomics analyses were made for serum samples collected as part of The Cardiovascular Risk in Young Finns Study. Samples from subjects who had an asymptomatic carotid artery plaque detected by ultrasound examination (N = 43, Age = 30–45 years) were compared with plaque free controls (N = 43) (matched for age, sex, body weight and systolic blood pressure). Seven proteins (p < 0.05) that have been previously linked with atherosclerotic phenotypes were differentially abundant. Fibulin 1 proteoform C (FBLN1C), Beta-ala-his-dipeptidase (CNDP1), Cadherin-13 (CDH13), Gelsolin (GSN) and 72 kDa type IV collagenase (MMP2) were less abundant in cases, whereas Apolipoproteins C-III (APOC3) and apolipoprotein E (APOE) were more abundant. Using machine learning analysis, a biomarker panel of FBLN1C, APOE and CDH13 was identified, which classified cases from controls with an area under receiver-operating characteristic curve (AUROC) value of 0.79. Furthermore, using selected reaction monitoring mass spectrometry (SRM-MS) the decreased abundance of FBLN1C was verified. In relation to previous associations of FBLN1C with atherosclerotic lesions, the observation could reflect its involvement in the initiation of the plaque formation, or represent a particular risk phenotype.