Anatomically-Driven Soft-Tissue Simulation Strategy for Cranio-Maxillofacial Surgery Using Facial Muscle Template Model

We propose a computationally efficient and biomechanically relevant soft-tissue simulation method for cranio-maxillofacial (CMF) surgery. A template-based facial muscle reconstruction was introduced to minimize the efforts on preparing a patient-specific model. A transversely isotropic mass-tensor model (MTM) was adopted to realize the effect of directional property of facial muscles in reasonable computation time. Additionally, sliding contact around teeth and mucosa was considered for more realistic simulation. Retrospective validation study with postoperative scan of a real patient showed that there were considerable improvements in simulation accuracy by incorporating template-based facial muscle anatomy and sliding contact.

Leaf transmission reduction using moving jaws for dynamic MLC IMRT

The aim of this work is to investigate to what extent it is possible to use the secondary collimator jaws to reduce the transmitted radiation through the multileaf collimator (MLC) during an intensity modulated radiation therapy (IMRT). A method is developed and introduced where the jaws follow the open window of the MLC dynamically (dJAW method). With the aid of three academic cases (Closed MLC, Sliding-gap, and Chair) and two clinical cases (prostate and head and neck) the feasibility of the dJAW method and the influence of this method on the applied dose distributions are investigated. For this purpose the treatment planning system Eclipse and the Research-Toolbox were used as well as measurements within a solid water phantom were performed. The transmitted radiation through the closed MLC leads to an inhomogeneous dose distribution. In this case, the measured dose within a plane perpendicular to the central axis differs up to 40% (referring to the maximum dose within this plane) for 6 and 15 MV. The calculated dose with Eclipse is clearly more homogeneous. For the Sliding-gap case this difference is still up to 9%. Among other things, these differences depend on the depth of the measurement within the solid water phantom and on the application method. In the Chair case, the dose in regions where no dose is desired is locally reduced by up to 50% using the dJAW method instead of the conventional method. The dose inside the chair-shaped region decreased up to 4% if the same number of monitor units (MU) as for the conventional method was applied. The undesired dose in the volume body minus the planning target volume in the clinical cases prostate and head and neck decreased up to 1.8% and 1.5%, while the number of the applied MU increased up to 3.1% and 2.8%, respectively. The new dJAW method has the potential to enhance the optimization of the conventional IMRT to a further step.

Theoretical considerations to the verification of dynamic multileaf collimated fields with a SLIC-type portal image detector

The verification possibilities of dynamically collimated treatment beams with a scanning liquid ionization chamber electronic portal image device (SLIC-EPID) are investigated. The ion concentration in the liquid of a SLIC-EPID and therefore the read-out signal is determined by two parameters of a differential equation describing the creation and recombination of the ions. Due to the form of this equation, the portal image detector describes a nonlinear dynamic system with memory. In this work, the parameters of the differential equation were experimentally determined for the particular chamber in use and for an incident open 6 MV photon beam. The mathematical description of the ion concentration was then used to predict portal images of intensity-modulated photon beams produced by a dynamic delivery technique, the sliding window approach. Due to the nature of the differential equation, a mathematical condition for 'reliable leaf motion verification' in the sliding window technique can be formulated. It is shown that the time constants for both formation and decay of the equilibrium concentration in the chamber is in the order of seconds. In order to guarantee reliable leaf motion verification, these time constants impose a constraint on the rapidity of the image-read out for a given maximum leaf speed. For a leaf speed of 2 cm s(-1), a minimum image acquisition frequency of about 2 Hz is required. Current SLIC-EPID systems are usually too slow since they need about a second to acquire a portal image. However, if the condition is fulfilled, the memory property of the system can be used to reconstruct the leaf motion. It is shown that a simple edge detecting algorithm can be employed to determine the leaf positions. The method is also very robust against image noise.

Reconstruction of the palatal aponeurosis with autogenous fascia lata in secondary radical intravelar veloplasty: a new method

Velopharyngeal insufficiency in cleft patients with muscular insufficiency detected by nasendoscopy is commonly treated by secondary radical intravelar veloplasty, in which the palatal muscles are reoriented and positioned backwards. The dead space between the retro-displaced musculature and the posterior borders of the palatal bone remains problematic. Postoperatively, the surgically achieved lengthening of the soft palate often diminishes due to scar tissue formation in the dead space, leading to reattachment of the reoriented muscles to the palatal bone and to decreased mobility of the soft palate. To avoid this, the dead space should be restored by a structure imitating the function of the missing palatal aponeurosis. The entire dead space was covered using a double layer of autogenous fascia lata harvested from the lateral thigh, which should allow sufficient and permanent sliding of the retro-positioned musculature. A clinical case of a 9-year-old boy who underwent the operation is reported. Postoperatively, marked functional improvements were observable in speech assessment, nasendoscopy and nasometry. The case reported here suggests that the restoration of the dead space may be beneficial for effective secondary palatal repair. Fascia lata seems to be a suitable graft for this purpose.

Parameter-free extraction of the functional network topology from intracranial EEG recordings: a time-resolved study of graph properties in focal onset seizures

Rationale: Focal onset epileptic seizures are due to abnormal interactions between distributed brain areas. By estimating the cross-correlation matrix of multi-site intra-cerebral EEG recordings (iEEG), one can quantify these interactions. To assess the topology of the underlying functional network, the binary connectivity matrix has to be derived from the cross-correlation matrix by use of a threshold. Classically, a unique threshold is used that constrains the topology [1]. Our method aims to set the threshold in a data-driven way by separating genuine from random cross-correlation. We compare our approach to the fixed threshold method and study the dynamics of the functional topology. Methods: We investigate the iEEG of patients suffering from focal onset seizures who underwent evaluation for the possibility of surgery. The equal-time cross-correlation matrices are evaluated using a sliding time window. We then compare 3 approaches assessing the corresponding binary networks. For each time window: * Our parameter-free method derives from the cross-correlation strength matrix (CCS)[2]. It aims at disentangling genuine from random correlations (due to finite length and varying frequency content of the signals). In practice, a threshold is evaluated for each pair of channels independently, in a data-driven way. * The fixed mean degree (FMD) uses a unique threshold on the whole connectivity matrix so as to ensure a user defined mean degree. * The varying mean degree (VMD) uses the mean degree of the CCS network to set a unique threshold for the entire connectivity matrix. * Finally, the connectivity (c), connectedness (given by k, the number of disconnected sub-networks), mean global and local efficiencies (Eg, El, resp.) are computed from FMD, CCS, VMD, and their corresponding random and lattice networks. Results: Compared to FMD and VMD, CCS networks present: *topologies that are different in terms of c, k, Eg and El. *from the pre-ictal to the ictal and then post-ictal period, topological features time courses that are more stable within a period, and more contrasted from one period to the next. For CCS, pre-ictal connectivity is low, increases to a high level during the seizure, then decreases at offset. k shows a ‘‘U-curve’’ underlining the synchronization of all electrodes during the seizure. Eg and El time courses fluctuate between the corresponding random and lattice networks values in a reproducible manner. Conclusions: The definition of a data-driven threshold provides new insights into the topology of the epileptic functional networks.

Single-Molecule Conductance of Functionalized Oligoynes: Length Dependence and Junction Evolution

We report a combined experimental and theoretical investigation of the length dependence and anchor group dependence of the electrical conductance of a series of oligoyne molecular wires in single-molecule junctions with gold contacts. Experimentally, we focus on the synthesis and properties of diaryloligoynes with n = 1, 2, and 4 triple bonds and the anchor dihydrobenzo[b]thiophene (BT). For comparison, we also explored the aurophilic anchor group cyano (CN), amino (NH2), thiol (SH), and 4-pyridyl (PY). Scanning tunneling microscopy break junction (STM-BJ) and mechanically controllable break junction (MCBJ) techniques are employed to investigate single-molecule conductance characteristics. The BT moiety is superior as compared to traditional anchoring groups investigated so far. BT-terminated oligoynes display a 100% probability of junction formation and possess conductance values which are the highest of the oligoynes studied and, moreover, are higher than other conjugated molecular wires of similar length. Density functional theory (DFT)-based calculations are reported for oligoynes with n = 1−4 triple bonds. Complete conductance traces and conductance distributions are computed for each family of molecules. The sliding of the anchor groups leads to oscillations in both the electrical conductance and the binding energies of the studied molecular wires. In agreement with experimental results, BT-terminated oligoynes are predicted to have a high electrical conductance. The experimental attenuation constants βH range between 1.7 nm−1 (CN) and 3.2 nm−1 (SH) and show the following trend: βH(CN) < βH(NH2) < βH(BT) < βH(PY) ≈ βH(SH). DFT-based calculations yield lower values, which range between 0.4 nm−1 (CN) and 2.2 nm−1 (PY).

Ankle joint pressure changes in a pes cavovarus model: supramalleolar valgus osteotomy versus lateralizing calcaneal osteotomy

BACKGROUND: A fixed cavovarus foot deformity can be associated with anteromedial ankle arthrosis due to elevated medial joint contact stresses. Supramalleolar valgus osteotomies (SMOT) and lateralizing calcaneal osteotomies (LCOT) are commonly used to treat symptoms by redistributing joint contact forces. In a cavovarus model, the effects of SMOT and LCOT on the lateralization of the center of force (COF) and reduction of the peak pressure in the ankle joint were compared. METHODS: A previously published cavovarus model with fixed hindfoot varus was simulated in 10 cadaver specimens. Closing wedge supramalleolar valgus osteotomies 3 cm above the ankle joint level (6 and 11 degrees) and lateral sliding calcaneal osteotomies (5 and 10 mm displacement) were analyzed at 300 N axial static load (half body weight). The COF migration and peak pressure decrease in the ankle were recorded using high-resolution TekScan pressure sensors. RESULTS: A significant lateral COF shift was observed for each osteotomy: 2.1 mm for the 6 degrees (P = .014) and 2.3 mm for the 11 degrees SMOT (P = .010). The 5 mm LCOT led to a lateral shift of 2.0 mm (P = .042) and the 10 mm LCOT to a shift of 3.0 mm (P = .006). Comparing the different osteotomies among themselves no significant differences were recorded. No significant anteroposterior COF shift was seen. A significant peak pressure reduction was recorded for each osteotomy: The SMOT led to a reduction of 29% (P = .033) for the 6 degrees and 47% (P = .003) for the 11 degrees osteotomy, and the LCOT to a reduction of 41% (P = .003) for the 5 mm and 49% (P = .002) for the 10 mm osteotomy. Similar to the COF lateralization no significant differences between the osteotomies were seen. CONCLUSION: LCOT and SMOT significantly reduced anteromedial ankle joint contact stresses in this cavovarus model. The unloading effects of both osteotomies were equivalent. More correction did not lead to significantly more lateralization of the COF or more reduction of peak pressure but a trend was seen. CLINICAL RELEVANCE: In patients with fixed cavovarus feet, both SMOT and LCOT provided equally good redistribution of elevated ankle joint contact forces. Increasing the amount of displacement did not seem to equally improve the joint pressures. The site of osteotomy could therefore be chosen on the basis of surgeon's preference, simplicity, or local factors in case of more complex reconstructions.

Active Testing for Face Detection and Localization

We provide a novel search technique which uses a hierarchical model and a mutual information gain heuristic to efficiently prune the search space when localizing faces in images. We show exponential gains in computation over traditional sliding window approaches, while keeping similar performance levels.

A g-factor metric for k-t-GRAPPA- and PEAK-GRAPPA-based parallel imaging.

PURPOSE The aim of this work is to derive a theoretical framework for quantitative noise and temporal fidelity analysis of time-resolved k-space-based parallel imaging methods. THEORY An analytical formalism of noise distribution is derived extending the existing g-factor formulation for nontime-resolved generalized autocalibrating partially parallel acquisition (GRAPPA) to time-resolved k-space-based methods. The noise analysis considers temporal noise correlations and is further accompanied by a temporal filtering analysis. METHODS All methods are derived and presented for k-t-GRAPPA and PEAK-GRAPPA. A sliding window reconstruction and nontime-resolved GRAPPA are taken as a reference. Statistical validation is based on series of pseudoreplica images. The analysis is demonstrated on a short-axis cardiac CINE dataset. RESULTS The superior signal-to-noise performance of time-resolved over nontime-resolved parallel imaging methods at the expense of temporal frequency filtering is analytically confirmed. Further, different temporal frequency filter characteristics of k-t-GRAPPA, PEAK-GRAPPA, and sliding window are revealed. CONCLUSION The proposed analysis of noise behavior and temporal fidelity establishes a theoretical basis for a quantitative evaluation of time-resolved reconstruction methods. Therefore, the presented theory allows for comparison between time-resolved parallel imaging methods and also nontime-resolved methods. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Small quantity but large effect — How minor phases control strain localization in upper mantle shear zones

Low viscosity domains such as localized shear zones exert an important control on the geodynamics of the uppermost mantle. Grain size reduction and subsequent strain localization related to a switch from dislocation to diffusion creep is one mechanism to form low viscosity domains. To sustain strain localization, the grain size of mantle minerals needs to be kept small over geological timescales. One way to keep olivine grain sizes small is by pinning of mobile grain boundaries during grain growth by other minerals (second phases). Detailed microstructural studies based on natural samples from three shear zones formed at different geodynamic settings, allowed the derivation of the olivine grain-size dependence on the second-phase content. The polymineralic olivine grain-size evolution with increasing strain is similar in the three shear zones. If the second phases are to pin the mobile olivine grain boundary the phases need to be well mixed before grain growth. We suggest that melt-rock and metamorphic reactions are crucial for the initial phase mixing in mantle rocks. With ongoing deformation and increasing strain, grain boundary sliding combined with mass transfer processes and nucleation of grains promotes phase mixing resulting in fine-grained polymineralic mixtures that deform by diffusion creep. Strain localization due to the presence of volumetrically minor minerals in polymineralic mantle rocks is only important at high strain deformation (ultramylonites) at low temperatures (<~800°C). At smaller strain and stress conditions and/or higher temperatures other parameters like overall energy available to deform a given rock volume, the inheritance of mechanical anisotropies or the presence of water or melts needs to be considered to explain strain localization in the upper mantle.

Granular flow in polymineralic rocks bearing sheet silicates: New evidence from natural examples

Natural deformation in carbonate mylonites bearing sheet silicates occurs via a complex interaction of granular flow and solution transfer processes and involves continuous cycles of dissolution, grain boundary diffusion, nucleation and growth. In this way, new sheet silicates (a) nucleate within voids formed by grain boundary sliding of calcite grains. (b) grow, and (c) rotate towards the shear plane. As a consequence, small mica grains show a wide range of orientations with respect to the shear plane, but moderate to large grains are subparallel both to each other and to the shear plane. Increases of average grain sizes with increasing temperature of sheet silicates in mica-rich layers is more pronounced than in mica-poor layers. In the calcitic matrix however, sheet silicates can only grow via solution-precipitation and mass transfer processes. Therefore, the observed grain size variability indicates drastic differences in mass transfer behavior between the individual layers, which might be related to differences in the fluid flux. Based on these observations, a conceptual model for the microfabric evolution in sheet silicate bearing mylonites is presented. © 2001 Elsevier Science B.V. All rights reserved.

Complete genome of a Puumala virus strain from Central Europe

Puumala virus (PUUV) is one of the predominant hantavirus species in Europe causing mild to moderate cases of haemorrhagic fever with renal syndrome. Parts of Lower Saxony in north-western Germany are endemic for PUUV infections. In this study, the complete PUUV genome sequence of a bank vole-derived tissue sample from the 2007 outbreak was determined by a combined primer-walking and RNA ligation strategy. The S, M and L genome segments were 1,828, 3,680 and 6,550 nucleotides in length, respectively. Sliding-window analyses of the nucleotide sequences of all available complete PUUV genomes indicated a non-homogenous distribution of variability with hypervariable regions located at the 3′-ends of the S and M segments. The overall similarity of the coding genome regions to the other PUUV strains ranged between 80.1 and 84.7 % at the level of the nucleotide sequence and between 89.5 and 98.1 % for the deduced amino acid sequences. In comparison to the phylogenetic trees of the complete coding sequences, trees based on partial segments revealed a general drop in phylogenetic support and a lower resolution. The Astrup strain S and M segment sequences showed the highest similarity to sequences of strains from geographically close sites in the Osnabrück Hills region. In conclusion, a primer-walking-mediated strategy resulted in the determination of the first complete nucleotide sequence of a PUUV strain from Central Europe. Different levels of variability along the genome provide the opportunity to choose regions for analyses according to the particular research question, e.g., large-scale phylogenetics or within-host evolution.

A role for tuned levels of nucleosome remodeler subunit ACF1 during Drosophila oogenesis

The Chromatin Accessibility Complex (CHRAC) consists of the ATPase ISWI, the large ACF1 subunit and a pair of small histone-like proteins, CHRAC-14/16. CHRAC is a prototypical nucleosome sliding factor that mobilizes nucleosomes to improve the regularity and integrity of the chromatin fiber. This may facilitate the formation of repressive chromatin. Expression of the signature subunit ACF1 is restricted during embryonic development, but remains high in primordial germ cells. Therefore, we explored roles for ACF1 during Drosophila oogenesis. ACF1 is expressed in somatic and germline cells, with notable enrichment in germline stem cells and oocytes. The asymmetrical localization of ACF1 to these cells depends on the transport of the Acf1 mRNA by the Bicaudal-D/Egalitarian complex. Loss of ACF1 function in the novel Acf1(7) allele leads to defective egg chambers and their elimination through apoptosis. In addition, we find a variety of unusual 16-cell cyst packaging phenotypes in the previously known Acf1(1) allele, with a striking prevalence of egg chambers with two functional oocytes at opposite poles. Surprisingly, we found that the Acf1(1) deletion - despite disruption of the Acf1 reading frame - expresses low levels of a PHD-bromodomain module from the C-terminus of ACF1 that becomes enriched in oocytes. Expression of this module from the Acf1 genomic locus leads to packaging defects in the absence of functional ACF1, suggesting competitive interactions with unknown target molecules. Remarkably, a two-fold overexpression of CHRAC (ACF1 and CHRAC-16) leads to increased apoptosis and packaging defects. Evidently, finely tuned CHRAC levels are required for proper oogenesis.

Properties of granular analogue model materials: A community wide survey

We report the material properties of 26 granular analogue materials used in 14 analogue modelling laboratories. We determined physical characteristics such as bulk density, grain size distribution, and grain shape, and performed ring shear tests to determine friction angles and cohesion, and uniaxial compression tests to evaluate the compaction behaviour. Mean grain size of the materials varied between c. 100 and 400 μm. Analysis of grain shape factors shows that the four different classes of granular materials (14 quartz sands, 5 dyed quartz sands, 4 heavy mineral sands and 3 size fractions of glass beads) can be broadly divided into two groups consisting of 12 angular and 14 rounded materials. Grain shape has an influence on friction angles, with most angular materials having higher internal friction angles (between c. 35° and 40°) than rounded materials, whereas well-rounded glass beads have the lowest internal friction angles (between c. 25° and 30°). We interpret this as an effect of intergranular sliding versus rolling. Most angular materials have also higher basal friction angles (tested for a specific foil) than more rounded materials, suggesting that angular grains scratch and wear the foil. Most materials have an internal cohesion in the order of 20–100 Pa except for well-rounded glass beads, which show a trend towards a quasi-cohesionless (C < 20 Pa) Coulomb-type material. The uniaxial confined compression tests reveal that rounded grains generally show less compaction than angular grains. We interpret this to be related to the initial packing density after sifting, which is higher for rounded grains than for angular grains. Ring-shear test data show that angular grains undergo a longer strain-hardening phase than more rounded materials. This might explain why analogue models consisting of angular grains accommodate deformation in a more distributed manner prior to strain localisation than models consisting of rounded grains.

Deformation mechanisms in second-phase affected microstructures and their energy balance

Based on the relationship Zener parameter (Z=second-phase size/second-phase volume fraction) vs. calcite grain size (dg), second-phase controlled aggregates and microstructures that are weakly affected by second-phases are discriminated. The latter are characterized by large but constant grain sizes, high calcite grain boundary fractions and crystallographic preferred orientations (CPO), while calcite grain size and calcite grain boundary fraction decrease continuously and CPO weakens with decreasing Z in second-phase controlled microstructures. These observations suggest that second-phase controlled microstructures predominantly deform via granular flow because pinning of calcite grain boundaries reduces the efficiency of dynamic recrystallization favoring mass transfer processes and grain boundary sliding. In contrast, the balance of grain size reduction and growth by dynamic recrystallization maintains a steady state grain size in microstructures that are only weakly affected by second-phases promoting a predominance of dislocation creep. With increasing temperature, the relationship between Z and dg persists but the calcite grain size increases continuously. Based on microstructures, the energy of each modifying process is calculated and its relative contribution is compared with energies of the competing processes (surface energy, dragging energy, dynamic recrystallization energy). The steady state microstructures result from a temperature-dependent energy minimization procedure of the system.