Clinical versus pre-clinical FE models for vertebral body strength predictions

The finite element analysis is an accepted method to predict vertebral body compressive strength. This study compares measurements obtained from in vitro tests with the ones from two different simulation models: clinical quantitative computer tomography (QCT) based homogenized finite element (hFE) models and pre-clinical high-resolution peripheral QCT-based (HR-pQCT) hFE models. About 37 vertebral body sections were prepared by removing end-plates and posterior elements, scanned with QCT (390/450μm voxel size) as well as HR-pQCT (82μm voxel size), and tested in compression up to failure. Non-linear viscous damage hFE models were created from QCT/HT-pQCT images and compared to experimental results based on stiffness and ultimate load. As expected, the predictability of QCT/HR-pQCT-based hFE models for both apparent stiffness (r2=0.685/0.801r2=0.685/0.801) and strength (r2=0.774/0.924r2=0.774/0.924) increased if a better image resolution was used. An analysis of the damage distribution showed similar damage locations for all cases. In conclusion, HR-pQCT-based hFE models increased the predictability considerably and do not need any tuning of input parameters. In contrast, QCT-based hFE models usually need some tuning but are clinically the only possible choice at the moment.

Improved Localization of Implanted Subdural Electrode Contacts on MRI Using an Elastic Image Fusion Algorithm in Invasive EEG Recording

BACKGROUND: Accurate projection of implanted subdural electrode contacts in presurgical evaluation of pharmacoresistant epilepsy cases by invasive EEG is highly relevant. Linear fusion of CT and MRI images may display the contacts in the wrong position due to brain shift effects. OBJECTIVE: A retrospective study in five patients with pharmacoresistant epilepsy was performed to evaluate whether an elastic image fusion algorithm can provide a more accurate projection of the electrode contacts on the pre-implantation MRI as compared to linear fusion. METHODS: An automated elastic image fusion algorithm (AEF), a guided elastic image fusion algorithm (GEF), and a standard linear fusion algorithm (LF) were used on preoperative MRI and post-implantation CT scans. Vertical correction of virtual contact positions, total virtual contact shift, corrections of midline shift and brain shifts due to pneumencephalus were measured. RESULTS: Both AEF and GEF worked well with all 5 cases. An average midline shift of 1.7mm (SD 1.25) was corrected to 0.4mm (SD 0.8) after AEF and to 0.0mm (SD 0) after GEF. Median virtual distances between contacts and cortical surface were corrected by a significant amount, from 2.3mm after LF to 0.0mm after AEF and GEF (p<.001). Mean total relative corrections of 3.1 mm (SD 1.85) after AEF and 3.0mm (SD 1.77) after GEF were achieved. The tested version of GEF did not achieve a satisfying virtual correction of pneumencephalus. CONCLUSION: The technique provided a clear improvement in fusion of pre- and post-implantation scans, although the accuracy is difficult to evaluate.

Digital Data Collection and Processing in Field and Lab: The Approach of Kinneret Regional Project

The paper showcases the field- and lab-documentation system developed for Kinneret Regional Project, an international archaeological expedition to the Northwestern shore of the Sea of Galilee (Israel) under the auspices of the University of Bern, the University of Helsinki, Leiden University and Wofford College. The core of the data management system is a fully relational, server-based database framework, which also includes time-based and static GIS services, stratigraphic analysis tools and fully indexed document/digital image archives. Data collection in the field is based on mobile, hand-held devices equipped with a custom-tailored stand-alone application. Comprehensive three-dimensional documentation of all finds and findings is achieved by means of total stations and/or high-precision GPS devices. All archaeological information retrieved in the field – including tachymetric data – is synched with the core system on the fly and thus immediately available for further processing in the field lab (within the local network) or for post-excavation analysis at remote institutions (via the WWW). Besides a short demonstration of the main functionalities, the paper also presents some of the key technologies used and illustrates usability aspects of the system’s individual components.

Image registration for triggered and non-triggered DTI of the human kidney: Reduced variability of diffusion parameter estimation

BACKGROUND: To investigate if non-rigid image-registration reduces motion artifacts in triggered and non-triggered diffusion tensor imaging (DTI) of native kidneys. A secondary aim was to determine, if improvements through registration allow for omitting respiratory-triggering. METHODS: Twenty volunteers underwent coronal DTI of the kidneys with nine b-values (10-700 s/mm2 ) at 3 Tesla. Image-registration was performed using a multimodal nonrigid registration algorithm. Data processing yielded the apparent diffusion coefficient (ADC), the contribution of perfusion (FP ), and the fractional anisotropy (FA). For comparison of the data stability, the root mean square error (RMSE) of the fitting and the standard deviations within the regions of interest (SDROI ) were evaluated. RESULTS: RMSEs decreased significantly after registration for triggered and also for non-triggered scans (P < 0.05). SDROI for ADC, FA, and FP were significantly lower after registration in both medulla and cortex of triggered scans (P < 0.01). Similarly the SDROI of FA and FP decreased significantly in non-triggered scans after registration (P < 0.05). RMSEs were significantly lower in triggered than in non-triggered scans, both with and without registration (P < 0.05). CONCLUSION: Respiratory motion correction by registration of individual echo-planar images leads to clearly reduced signal variations in renal DTI for both triggered and particularly non-triggered scans. Secondarily, the results suggest that respiratory-triggering still seems advantageous.J. Magn. Reson. Imaging 2014. (c) 2014 Wiley Periodicals, Inc.

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.

Progressive Image Denoising

Image denoising continues to be an active research topic. Although state-of-the-art denoising methods are numerically impressive and approch theoretical limits, they suffer from visible artifacts.While they produce acceptable results for natural images, human eyes are less forgiving when viewing synthetic images. At the same time, current methods are becoming more complex, making analysis, and implementation difficult. We propose image denoising as a simple physical process, which progressively reduces noise by deterministic annealing. The results of our implementation are numerically and visually excellent. We further demonstrate that our method is particularly suited for synthetic images. Finally, we offer a new perspective on image denoising using robust estimators.

Fully Automatic CT Segmentation for Computer-Assisted Pre-operative Planning of Hip Arthroscopy

Extraction of both pelvic and femoral surface models of a hip joint from CT data for computer-assisted pre-operative planning of hip arthroscopy is addressed. We present a method for a fully automatic image segmentation of a hip joint. Our method works by combining fast random forest (RF) regression based landmark detection, atlas-based segmentation, with articulated statistical shape model (aSSM) based hip joint reconstruction. The two fundamental contributions of our method are: (1) An improved fast Gaussian transform (IFGT) is used within the RF regression framework for a fast and accurate landmark detection, which then allows for a fully automatic initialization of the atlas-based segmentation; and (2) aSSM based fitting is used to preserve hip joint structure and to avoid penetration between the pelvic and femoral models. Validation on 30 hip CT images show that our method achieves high performance in segmenting pelvis, left proximal femur, and right proximal femur surfaces with an average accuracy of 0.59 mm, 0.62 mm, and 0.58 mm, respectively.

A complete-pelvis segmentation framework for image-free total hip arthroplasty (THA): methodology and clinical study

Background Complete-pelvis segmentation in antero-posterior pelvic radiographs is required to create a patient-specific three-dimensional pelvis model for surgical planning and postoperative assessment in image-free navigation of total hip arthroplasty. Methods A fast and robust framework for accurately segmenting the complete pelvis is presented, consisting of two consecutive modules. In the first module, a three-stage method was developed to delineate the left hemipelvis based on statistical appearance and shape models. To handle complex pelvic structures, anatomy-specific information processing techniques were employed. As the input to the second module, the delineated left hemi-pelvis was then reflected about an estimated symmetry line of the radiograph to initialize the right hemi-pelvis segmentation. The right hemi-pelvis was segmented by the same three-stage method, Results Two experiments conducted on respectively 143 and 40 AP radiographs demonstrated a mean segmentation accuracy of 1.61±0.68 mm. A clinical study to investigate the postoperative assessment of acetabular cup orientations based on the proposed framework revealed an average accuracy of 1.2°±0.9° and 1.6°±1.4° for anteversion and inclination, respectively. Delineation of each radiograph costs less than one minute. Conclusions Despite further validation needed, the preliminary results implied the underlying clinical applicability of the proposed framework for image-free THA.

CstF64: cell cycle regulation and functional role in 3' end processing of replication-dependent histone mRNAs.

The 3' end processing of animal replication-dependent histone mRNAs is activated during G1/S-phase transition. The processing site is recognized by stem-loop binding protein and the U7 snRNP, but cleavage additionally requires a heat-labile factor (HLF), composed of cleavage/polyadenylation specificity factor, symplekin, and cleavage stimulation factor 64 (CstF64). Although HLF has been shown to be cell cycle regulated, the mechanism of this regulation is unknown. Here we show that levels of CstF64 increase toward the S phase and its depletion affects histone RNA processing, S-phase progression, and cell proliferation. Moreover, analyses of the interactions between CstF64, symplekin, and the U7 snRNP-associated proteins FLASH and Lsm11 indicate that CstF64 is important for recruiting HLF to histone precursor mRNA (pre-mRNA)-resident proteins. Thus, CstF64 is central to the function of HLF and appears to be at least partly responsible for its cell cycle regulation. Additionally, we show that misprocessed histone transcripts generated upon CstF64 depletion mainly accumulate in the nucleus, where they are targets of the exosome machinery, while a small cytoplasmic fraction is partly associated with polysomes.

Antisense derivatives of U7 and other small nuclear RNAs as tools to modify pre-mRNA splicing patterns

The importance of alternative splicing for the diversity of the proteome and the large number of genetic diseases that are due to splicing defects call for methods to modulate alternative splicing decisions. Although splicing can be modulated by antisense oligonucleotides, this approach is confronted with problems of efficient delivery and the need for repeated administrations of large amounts of the oligonucleotides. Therefore we have developed methods allowing us to modulate splicing with the help of modified derivatives of the U7 small nuclear RNA involved in histone RNA 3' end processing. Its nuclear accumulation as a stable ribonucleoprotein particle makes U7 snRNA especially useful for this purpose. In particular, U7 derivatives containing two tandem antisense sequences directed against targets upstream and downstream of an exon can induce the efficient and specific skipping of that exon. U7 expression cassettes have been successfully introduced into a great number of cell lines, primary cells or tissues with the help of lentiviral and adeno-associated viral vectors. Examples of these therapeutic strategies in the fields of β-thalassemia, Duchenne muscular dytrophy and HIV/AIDS are discussed.

Functional importance of conserved nucleotides at the histone RNA 3' processing site.

Histone pre-mRNA 3' processing is controlled by a hairpin element preceding the processing site that interacts with a hairpin-binding protein (HBP) and a downstream spacer element that serves as anchoring site for the U7 snRNP. In addition, the nucleotides following the hairpin and surrounding the processing site (ACCCA'CA) are conserved among vertebrate histone genes. Single to triple nucleotide mutations of this sequence were tested for their ability to be processed in nuclear extract from animal cells. Changing the first four nucleotides had no qualitative and little if any quantitative effects on histone RNA 3' processing in mouse K21 cell extract, where processing of this gene is virtually independent of the HBP. A gel mobility shift assay revealing HBP interactions and a processing assay in HeLa cell extract (where the contribution of HBP to efficient processing is more important) showed that only one of these mutations, predicted to extend the hairpin by one base pair, affected the interaction with HBP. Mutations in the next three nucleotides affected both the cleavage efficiency and the choice of processing sites. Analysis of these novel sites indicated a preference for the nucleotide 5' of the cleavage site in the order A > C > U > G. Moreover, a guanosine in the 3' position inhibited cleavage. The preference for an A is shared with the cleavage/polyadenylation reaction, but the preference order for the other nucleotides is different [Chen F, MacDonald CC, Wilusz J, 1995, Nucleic Acids Res 23:2614-2620].

A synthetic histone pre-mRNA-U7 small nuclear RNA chimera undergoing cis cleavage in the cytoplasm of Xenopus oocytes.

The 3' processing of histone pre-mRNAs is a nuclear event in which the U7 small nuclear ribonucleoprotein (snRNP) participates as an essential trans-acting factor. We have constructed a chimeric histone-U7 RNA that when injected into the cytoplasm of Xenopus laevis oocytes assembles into a snRNP-like particle and becomes cleaved at the correct site(s). RNP assembly is a prerequisite for cleavage, but, since neither the RNA nor the RNP appreciably enter the nucleus, cleavage occurs mostly, if not exclusively, in the cytoplasm. Consistent with this, cleavage also occurs in enucleated oocytes or in oocytes which have been depleted of U7 snRNPs. Thus all necessary components for cleavage must be present in the oocyte cytoplasm. The novel cleavage occurs in cis, involving only a single molecule of chimeric RNA with its associated proteins. This reaction is equally dependent upon base pairing interactions between histone spacer sequences and the 5'-end of the U7 moiety as the natural in trans reaction. These results imply that U7 is the only snRNP required for histone RNA processing. Moreover, the chimeric RNA is expected to be useful for further studies of the cleavage and assembly mechanisms of U7 snRNP.

Biochemical demonstration of complex formation of histone pre-mRNA with U7 small nuclear ribonucleoprotein and hairpin binding factors.

Histone RNA 3' end formation occurs through a specific cleavage reaction that requires, among other things, base-pairing interactions between a conserved spacer element in the pre-mRNA and the minor U7 snRNA present as U7 snRNP. An oligonucleotide complementary to the first 16 nucleotides of U7 RNA can be used to characterize U7 snRNPs from nuclear extracts by native gel electrophoresis. Using similar native gel techniques, we present direct biochemical evidence for a stable association between histone pre-mRNA and U7 snRNPs. Other complexes formed in the nuclear extract are dependent on the 5' cap structure and on the conserved hairpin element of histone pre-mRNA, respectively. However, in contrast to the U7-specific complex, their formation is not required for processing. Comparison of several authentic and mutant histone pre-mRNAs with different spacer sequences demonstrates that the formation and stability of the U7-specific complex closely follows the predicted stability of the potential RNA-RNA hybrid. However, this does not exclude a stabilization of the complex by U7 snRNP structural proteins.

A total variation discontinuous Galerkin approach for image restoration

The focal point of this paper is to propose and analyze a P 0 discontinuous Galerkin (DG) formulation for image denoising. The scheme is based on a total variation approach which has been applied successfully in previous papers on image processing. The main idea of the new scheme is to model the restoration process in terms of a discrete energy minimization problem and to derive a corresponding DG variational formulation. Furthermore, we will prove that the method exhibits a unique solution and that a natural maximum principle holds. In addition, a number of examples illustrate the effectiveness of the method.