The long-term mechanical integrity of non-reinforced PEEK-OPTIMA polymer for demanding spinal applications: experimental and finite-element analysis

Polyetheretherketone (PEEK) is a novel polymer with potential advantages for its use in demanding orthopaedic applications (e.g. intervertebral cages). However, the influence of a physiological environment on the mechanical stability of PEEK has not been reported. Furthermore, the suitability of the polymer for use in highly stressed spinal implants such as intervertebral cages has not been investigated. Therefore, a combined experimental and analytical study was performed to address these open questions. A quasi-static mechanical compression test was performed to compare the initial mechanical properties of PEEK-OPTIMA polymer in a dry, room-temperature and in an aqueous, 37 degrees C environment (n=10 per group). The creep behaviour of cylindrical PEEK polymer specimens (n=6) was measured in a simulated physiological environment at an applied stress level of 10 MPa for a loading duration of 2000 hours (12 weeks). To compare the biomechanical performance of different intervertebral cage types made from PEEK and titanium under complex loading conditions, a three-dimensional finite element model of a functional spinal unit was created. The elastic modulus of PEEK polymer specimens in a physiological environment was 1.8% lower than that of specimens tested at dry, room temperature conditions (P<0.001). The results from the creep test showed an average creep strain of less than 0.1% after 2000 hours of loading. The finite element analysis demonstrated high strain and stress concentrations at the bone/implant interface, emphasizing the importance of cage geometry for load distribution. The stress and strain maxima in the implants were well below the material strength limits of PEEK. In summary, the experimental results verified the mechanical stability of the PEEK-OPTIMA polymer in a simulated physiological environment, and over extended loading periods. Finite element analysis supported the use of PEEK-OPTIMA for load-bearing intervertebral implants.

Protrusio acetabuli: joint loading with severe pincer impingement and its theoretical implications for surgical therapy

Severe pincer impingement (acetabular protrusio) is an established cause of hip pain and osteoarthritis. The proposed underlying pathomechanism is a dynamic pathological contact of the prominent acetabular rim with the femoral head-neck junction. However, this cannot explain the classically described medial osteoarthritis in these hips. We therefore asked: (1) Does an overload exist in the medial aspect of the protrusio joint? and (2) What is the influence of three contemporary joint-preserving procedures on load distribution in protrusio hips? In vivo force and motion data for walking and standing to sitting were applied to six 3D finite element models (normal, dysplasia, protrusio, acetabular rim trimming, acetabular reorientation, and combined reorientation/rim trimming). Compared with dysplasia, the protrusio joint resulted in opposite patterns of von Mises stress and contact pressure during walking. In protrusio hips, we found an overload at the medial margin of the lunate surface (54% higher than normal). Isolated rim trimming further increased the medial overload (up to 28% higher than protrusio), whereas acetabular reorientation with/without rim trimming reduced stresses by up to 25%. Our results can be used as an adjunct for surgical decision making in the treatment of acetabular protrusio.

Viral load, organ distribution, histopathological lesions, and cytokine mRNA expression in goats infected with a molecular clone of the caprine arthritis encephalitis virus

Caprine arthritis encephalitis virus (CAEV) is a lentivirus of goats that causes persistent infection characterized by the appearance of inflammatory lesions in various organs. To define the sites of persistence, 5 goats were infected with a molecular clone of CAEV, and the viral load was monitored by real-time-PCR and RT-PCR in different sites 8 years after infection. The lymph nodes proved to be an important virus reservoir, with moderate virus replication relative to what is reported for lentiviruses of primates. Mammary gland and milk cells were preferred sites of viral replication. The viral load varied significantly between animals, which points to an important role of the genetic background. We found a clear association between occurrence of histopathological lesions and viral load in specific sites. The mRNA expression analysis of several cytokines did not reveal differences between animals that could explain the considerable individual variations in viral load observed.

Cement distribution, volume, and compliance in vertebroplasty: some answers from an anatomy-based nonlinear finite element study

STUDY DESIGN: The biomechanics of vertebral bodies augmented with real distributions of cement were investigated using nonlinear finite element (FE) analysis. OBJECTIVES: To compare stiffness, strength, and stress transfer of augmented versus nonaugmented osteoporotic vertebral bodies under compressive loading. Specifically, to examine how cement distribution, volume, and compliance affect these biomechanical variables. SUMMARY OF BACKGROUND DATA: Previous FE studies suggested that vertebroplasty might alter vertebral stress transfer, leading to adjacent vertebral failure. However, no FE study so far accounted for real cement distributions and bone damage accumulation. METHODS: Twelve vertebral bodies scanned with high-resolution pQCT and tested in compression were augmented with various volumes of cements and scanned again. Nonaugmented and augmented pQCT datasets were converted to FE models, with bone properties modeled with an elastic, plastic and damage constitutive law that was previously calibrated for the nonaugmented models. The cement-bone composite was modeled with a rule of mixture. The nonaugmented and augmented FE models were subjected to compression and their stiffness, strength, and stress map calculated for different cement compliances. RESULTS: Cement distribution dominated the stiffening and strengthening effects of augmentation. Models with cement connecting either the superior or inferior endplate (S/I fillings) were only up to 2 times stiffer than the nonaugmented models with minimal strengthening, whereas those with cement connecting both endplates (S + I fillings) were 1 to 8 times stiffer and 1 to 12 times stronger. Stress increases above and below the cement, which was higher for the S + I cases and was significantly reduced by increasing cement compliance. CONCLUSION: The developed FE approach, which accounts for real cement distributions and bone damage accumulation, provides a refined insight into the mechanics of augmented vertebral bodies. In particular, augmentation with compliant cement bridging both endplates would reduce stress transfer while providing sufficient strengthening.

Embedding of human vertebral bodies leads to higher ultimate load and altered damage localisation under axial compression

Computer tomography (CT)-based finite element (FE) models of vertebral bodies assess fracture load in vitro better than dual energy X-ray absorptiometry, but boundary conditions affect stress distribution under the endplates that may influence ultimate load and damage localisation under post-yield strains. Therefore, HRpQCT-based homogenised FE models of 12 vertebral bodies were subjected to axial compression with two distinct boundary conditions: embedding in polymethylmethalcrylate (PMMA) and bonding to a healthy intervertebral disc (IVD) with distinct hyperelastic properties for nucleus and annulus. Bone volume fraction and fabric assessed from HRpQCT data were used to determine the elastic, plastic and damage behaviour of bone. Ultimate forces obtained with PMMA were 22% higher than with IVD but correlated highly (R2 = 0.99). At ultimate force, distinct fractions of damage were computed in the endplates (PMMA: 6%, IVD: 70%), cortex and trabecular sub-regions, which confirms previous observations that in contrast to PMMA embedding, failure initiated underneath the nuclei in healthy IVDs. In conclusion, axial loading of vertebral bodies via PMMA embedding versus healthy IVD overestimates ultimate load and leads to distinct damage localisation and failure pattern.

Chronic bovine besnoitiosis: intra-organ parasite distribution, parasite loads and parasite-associated lesions in subclinical cases.

Bovine besnoitiosis caused by Besnoitia besnoiti is a chronic and debilitating disease. The most characteristic clinical signs of chronic besnoitiosis are visible tissue cysts in the scleral conjunctiva and the vagina, thickened skin and a generally poor body condition. However, many seropositive animals remain subclinically infected, and the role that these animals may play in spreading the disease is not known. The aim of the present study was to assess the intra-organ parasite distribution, the parasite load and the parasite-associated lesions in seropositive but subclinically infected animals. These animals were seropositive at the time of several consecutive samplings, had visible tissue cysts in the past and, at time of slaughter, had detectable specific anti-Besnoitia spp. antibody levels, but they did not show evident clinical signs at culling. Thus, histopathological, immunohistochemical and molecular analyses of several samples from the respiratory tract, reproductive tract, other internal organs and skin from six cows were performed. The tissue cysts were located primarily in the upper respiratory tract, i.e., in the rhinarium and larynx/pharynx (four cows), followed by the distal genital tract (vulva/vagina) and the skin of the neck (three and two cows, respectively, out of the four cows with cysts in the respiratory tract). We were unable to detect any parasites in the two remaining cows. Cysts were associated with a significant non-purulent inflammatory infiltrate consisting predominantly of T lymphocytes and activated monocytes/macrophages in two cows. The parasite burden, estimated by quantitative real-time PCR, was very low. It is noteworthy that the only animal that showed a recent increase in the antibody titre had the highest parasite burden and the most conspicuous inflammatory reaction against the cysts. In conclusion, although these cows no longer displayed any visible signs of besnoitiosis, they remained infected. Therefore, cows without visible signs of disease may still be able to transmit the parasite.

The brain’s default state interacts with working memory processes in a complex and load-dependent manner

Recently, many studies about a network active during rest and deactivated during tasks emerged in the literature: the default mode network (DMN). Spatial and temporal DMN features are important markers for psychiatric diseases. Another prominent indicator of cognitive functioning, yielding information about the mental condition in health and disease, is working memory (WM) processing. In EEG studies, frontal-midline theta power has been shown to increase with load during WM retention in healthy subjects. From these findings, the conclusion can be drawn that an increase in resting state DMN activity may go along with an increase in theta power in high-load WM conditions. We followed this hypothesis in a study on 17 healthy subjects performing a visual Sternberg WM task. The DMN was obtained by a BOLD-ICA approach and its dynamics represented by the percent-strength during pre-stimulus periods. DMN dynamics were temporally correlated with EEG theta spectral power from retention intervals. This so-called covariance mapping yielded the spatial distribution of the theta EEG fluctuations associated with the dynamics of the DMN. In line with previous findings, theta power was increased at frontal-midline electrodes in high- versus low-load conditions during early WM retention. However, load-dependent correlations of DMN with theta power resulted in primarily positive correlations in low-load conditions, while during high-load conditions negative correlations of DMN activity and theta power were observed at frontal-midline electrodes. This DMN-dependent load effect reached significance during later retention. Our results show a complex and load-dependent interaction of pre-stimulus DMN activity and theta power during retention, varying over the course of the retention period. Since both, WM performance and DMN activity, are markers of mental health, our results could be important for further investigations of psychiatric populations.

The brain’s pre-stimulus default state interacts with working memory in a load-dependent manner

Recently, multiple studies showed that spatial and temporal features of a task-negative default mode network (DMN) (Greicius et al., 2003) are important markers for psychiatric diseases (Balsters et al., 2013). Another prominent indicator of cognitive functioning, yielding information about the mental condition in health and disease, is working memory (WM) processing. In EEG and MEG studies, frontal-midline theta power has been shown to increase with load during WM retention in healthy subjects (Brookes et al., 2011). Negative correlations between DMN activity and theta amplitude have been found during resting state (Jann et al., 2010) as well as during WM (Michels et al., 2010). Likewise, WM training resulted in higher resting state theta power as well as increased small-worldness of the resting brain (Langer et al., 2013). Further, increased fMRI connectivity between nodes of the DMN correlated with better WM performance (Hampson et al., 2006). Hence, the brain’s default state might influence it’s functioning during task. We therefore hypothesized correlations between pre-stimulus DMN activity and EEG-theta power during WM maintenance, depending on the WM load. 17 healthy subjects performed a Sternberg WM task while being measured simultaneously with EEG and fMRI. Data was recorded within a multicenter-study: 12 subjects were measured in Zurich with a 64-channels MR-compatible system (Brain Products) in a 3T Philips scanner, 5 subjects with a 96-channel MR-compatible system (Brain Products) in a 3T Siemens Scanner in Bern. The DMN components was obtained by a group BOLD-ICA approach over the full task duration (figure 1). The subject-wise dynamics were obtained by back-reconstructed onto each subject’s fMRI data and normalized to percent signal change values. The single trial pre-stimulus-DMN activation was then temporally correlated with the single trial EEG-theta (3-8 Hz) spectral power during retention intervals. This so-called covariance mapping (Jann et al., 2010) yielded the spatial distribution of the theta EEG fluctuations during retention associated with the dynamics of the pre-stimulus DMN. In line with previous findings, theta power was increased at frontal-midline electrodes in high- versus low-load conditions during early WM retention (figure 2). However, correlations of DMN with theta power resulted in primarily positive correlations in low-load conditions, while during high-load conditions negative correlations of DMN activity and theta power were observed at frontal-midline electrodes. This DMN-dependent load effect reached significance in the middle of the retention period (TANOVA, p<0.05) (figure 3). Our results show a complex and load-dependent interaction of pre-stimulus DMN activity and theta power during retention, varying over time. While at a more global, load-independent view pre-stimulus DMN activity correlated positively with theta power during retention, the correlation was inversed during certain time windows in high-load trials, meaning that in trials with enhanced pre-stimulus DMN activity theta power decreases during retention. Since both WM performance and DMN activity are markers of mental health our results could be important for further investigations of psychiatric populations.

Risk charts to guide targeted HIV-1 viral load monitoring of ART: development and validation in patients from resource-limited settings.

BACKGROUND HIV-1 RNA viral load (VL) testing is recommended to monitor antiretroviral therapy (ART) but not available in many resource-limited settings. We developed and validated CD4-based risk charts to guide targeted VL testing. METHODS We modeled the probability of virologic failure up to 5 years of ART based on current and baseline CD4 counts, developed decision rules for targeted VL testing of 10%, 20% or 40% of patients in seven cohorts of patients starting ART in South Africa, and plotted cut-offs for VL testing on colour-coded risk charts. We assessed the accuracy of risk chart-guided VL testing to detect virologic failure in validation cohorts from South Africa, Zambia and the Asia-Pacific. FINDINGS 31,450 adult patients were included in the derivation and 25,294 patients in the validation cohorts. Positive predictive values increased with the percentage of patients tested: from 79% (10% tested) to 98% (40% tested) in the South African, from 64% to 93% in the Zambian and from 73% to 96% in the Asia-Pacific cohorts. Corresponding increases in sensitivity were from 35% to 68% in South Africa, from 55% to 82% in Zambia and from 37% to 71% in Asia-Pacific. The area under the receiver-operating curve increased from 0.75 to 0.91 in South Africa, from 0.76 to 0.91 in Zambia and from 0.77 to 0.92 in Asia Pacific. INTERPRETATION CD4-based risk charts with optimal cut-offs for targeted VL testing may be useful to monitor ART in settings where VL capacity is limited.

Implementation and Operational Research: Risk Charts to Guide Targeted HIV-1 Viral Load Monitoring of ART: Development and Validation in Patients From Resource-Limited Settings.

BACKGROUND HIV-1 RNA viral load (VL) testing is recommended to monitor antiretroviral therapy (ART) but not available in many resource-limited settings. We developed and validated CD4-based risk charts to guide targeted VL testing. METHODS We modeled the probability of virologic failure up to 5 years of ART based on current and baseline CD4 counts, developed decision rules for targeted VL testing of 10%, 20%, or 40% of patients in 7 cohorts of patients starting ART in South Africa, and plotted cutoffs for VL testing on colour-coded risk charts. We assessed the accuracy of risk chart-guided VL testing to detect virologic failure in validation cohorts from South Africa, Zambia, and the Asia-Pacific. RESULTS In total, 31,450 adult patients were included in the derivation and 25,294 patients in the validation cohorts. Positive predictive values increased with the percentage of patients tested: from 79% (10% tested) to 98% (40% tested) in the South African cohort, from 64% to 93% in the Zambian cohort, and from 73% to 96% in the Asia-Pacific cohort. Corresponding increases in sensitivity were from 35% to 68% in South Africa, from 55% to 82% in Zambia, and from 37% to 71% in Asia-Pacific. The area under the receiver operating curve increased from 0.75 to 0.91 in South Africa, from 0.76 to 0.91 in Zambia, and from 0.77 to 0.92 in Asia-Pacific. CONCLUSIONS CD4-based risk charts with optimal cutoffs for targeted VL testing maybe useful to monitor ART in settings where VL capacity is limited.

Distance from Sub-Saharan Africa Predicts Mutational Load in Diverse Human Genomes

The Out-of-Africa (OOA) dispersal ∼50,000 y ago is characterized by a series of founder events as modern humans expanded into multiple continents. Population genetics theory predicts an increase of mutational load in populations undergoing serial founder effects during range expansions. To test this hypothesis, we have sequenced full genomes and high-coverage exomes from seven geographically divergent human populations from Namibia, Congo, Algeria, Pakistan, Cambodia, Siberia, and Mexico. We find that individual genomes vary modestly in the overall number of predicted deleterious alleles. We show via spatially explicit simulations that the observed distribution of deleterious allele frequencies is consistent with the OOA dispersal, particularly under a model where deleterious mutations are recessive. We conclude that there is a strong signal of purifying selection at conserved genomic positions within Africa, but that many predicted deleterious mutations have evolved as if they were neutral during the expansion out of Africa. Under a model where selection is inversely related to dominance, we show that OOA populations are likely to have a higher mutation load due to increased allele frequencies of nearly neutral variants that are recessive or partially recessive.

Expansion load: recessive mutations and the role of standing genetic variation

Expanding populations incur a mutation burden – the so-called expansion load. Previous studies of expansion load have focused on codominant mutations. An important consequence of this assumption is that expansion load stems exclusively from the accumulation of new mutations occurring in individuals living at the wave front. Using individual-based simulations, we study here the dynamics of standing genetic variation at the front of expansions, and its consequences on mean fitness if mutations are recessive. We find that deleterious genetic diversity is quickly lost at the front of the expansion, but the loss of deleterious mutations at some loci is compensated by an increase of their frequencies at other loci. The frequency of deleterious homozygotes therefore increases along the expansion axis, whereas the average number of deleterious mutations per individual remains nearly constant across the species range. This reveals two important differences to codominant models: (i) mean fitness at the front of the expansion drops much faster if mutations are recessive, and (ii) mutation load can increase during the expansion even if the total number of deleterious mutations per individual remains constant. We use our model to make predictions about the shape of the site frequency spectrum at the front of range expansion, and about correlations between heterozygosity and fitness in different parts of the species range. Importantly, these predictions provide opportunities to empirically validate our theoretical results. We discuss our findings in the light of recent results on the distribution of deleterious genetic variation across human populations and link them to empirical results on the correlation of heterozygosity and fitness found in many natural range expansions.

Size-frequency distribution of boulders ≥7 m on comet 67P/Churyumov-Gerasimenko

Aims. We derive for the first time the size-frequency distribution of boulders on a comet, 67P/Churyumov-Gerasimenko (67P), computed from the images taken by the Rosetta/OSIRIS imaging system. We highlight the possible physical processes that lead to these boulder size distributions. Methods. We used images acquired by the OSIRIS Narrow Angle Camera, NAC, on 5 and 6 August 2014. The scale of these images (2.44−2.03 m/px) is such that boulders ≥7 m can be identified and manually extracted from the datasets with the software ArcGIS. We derived both global and localized size-frequency distributions. The three-pixel sampling detection, coupled with the favorable shadowing of the surface (observation phase angle ranging from 48° to 53°), enables unequivocally detecting boulders scattered all over the illuminated side of 67P. Results. We identify 3546 boulders larger than 7 m on the imaged surface (36.4 km2), with a global number density of nearly 100/km2 and a cumulative size-frequency distribution represented by a power-law with index of −3.6 +0.2/−0.3. The two lobes of 67P appear to have slightly different distributions, with an index of −3.5 +0.2/−0.3 for the main lobe (body) and −4.0 +0.3/−0.2 for the small lobe (head). The steeper distribution of the small lobe might be due to a more pervasive fracturing. The difference of the distribution for the connecting region (neck) is much more significant, with an index value of −2.2 +0.2/−0.2. We propose that the boulder field located in the neck area is the result of blocks falling from the contiguous Hathor cliff. The lower slope of the size-frequency distribution we see today in the neck area might be due to the concurrent processes acting on the smallest boulders, such as i) disintegration or fragmentation and vanishing through sublimation; ii) uplifting by gas drag and consequent redistribution; and iii) burial beneath a debris blanket. We also derived the cumulative size-frequency distribution per km2 of localized areas on 67P. By comparing the cumulative size-frequency distributions of similar geomorphological settings, we derived similar power-law index values. This suggests that despite the selected locations on different and often opposite sides of the comet, similar sublimation or activity processes, pit formation or collapses, as well as thermal stresses or fracturing events occurred on multiple areas of the comet, shaping its surface into the appearance we see today.

Benford's Law and Fraud Detection. Facts and Legends

Is Benford's law a good instrument to detect fraud in reports of statistical and scientific data? For a valid test the probability of "false positives" and "false negatives" has to be low. However, it is very doubtful whether the Benford distribution is an appropriate tool to discriminate between manipulated and non-manipulated estimates. Further research should focus more on the validity of the test and test results should be interpreted more carefully.