Multi-modal deformable medical image registration

Non-rigid image registration is an essential tool required for overcoming the inherent local anatomical variations that exist between images acquired from different individuals or atlases. Furthermore, certain applications require this type of registration to operate across images acquired from different imaging modalities. One popular local approach for estimating this registration is a block matching procedure utilising the mutual information criterion. However, previous block matching procedures generate a sparse deformation field containing displacement estimates at uniformly spaced locations. This neglects to make use of the evidence that block matching results are dependent on the amount of local information content. This paper presents a solution to this drawback by proposing the use of a Reversible Jump Markov Chain Monte Carlo statistical procedure to optimally select grid points of interest. Three different methods are then compared to propagate the estimated sparse deformation field to the entire image including a thin-plate spline warp, Gaussian convolution, and a hybrid fluid technique. Results show that non-rigid registration can be improved by using the proposed algorithm to optimally select grid points of interest.

Improving fluid registration through white matter segmentation in a twin study design

Robust and automatic non-rigid registration depends on many parameters that have not yet been systematically explored. Here we determined how tissue classification influences non-linear fluid registration of brain MRI. Twin data is ideal for studying this question, as volumetric correlations between corresponding brain regions that are under genetic control should be higher in monozygotic twins (MZ) who share 100% of their genes when compared to dizygotic twins (DZ) who share half their genes on average. When these substructure volumes are quantified using tensor-based morphometry, improved registration can be defined based on which method gives higher MZ twin correlations when compared to DZs, as registration errors tend to deplete these correlations. In a study of 92 subjects, higher effect sizes were found in cumulative distribution functions derived from statistical maps when performing tissue classification before fluid registration, versus fluidly registering the raw images. This gives empirical evidence in favor of pre-segmenting images for tensor-based morphometry.

Fasteners in Composite Plates - Effects of Interfacial Friction

The effects of tangential friction at pin—hole interfaces are appropriately modelled for the analysis of fasteners in large composite (orthotropic) plate loaded along its edges. The pin—hole contact could be of interference, clearance or neat fit. When the plate load is monotonically increased, interference fits give rise to receding contact, whereas clearance fits result in advancing contact. In either case, the changing contact situations lead to non-linear moving boundary value problems. The neat fit comes out as a special case in which the contact and separation regions are invariant with the applied load level and so the problem remains linear. The description of boundary conditions in the presence of tangential friction, will depend on whether the problem is one of advancing or receding contact, advancing contact presenting a special problem. A model is developed for the limiting case of a rigid pin and an ideally rough interface (infinitely large friction coefficient). The non-linearity resulting from the continuously varying proportions of contact and separation at the interface, is handled by an “Inverse Formulation” which was successfully applied earlier by the authors for smooth (zero friction) interfacial conditions. The additional difficulty introduced by advancing contact is handled by adopting a “Marching Solution”. The modelling and the procedure are illustrated in respect of symmetric plate load cases. Numerical results are presented bringing out the effects of interfacial friction and plate orthotropy on load-contact relations and plate stresses.

Non-linear flapping vibrations of rotating blades

The non-linear equations of motion of a rotating blade undergoing extensional and flapwise bending vibration are derived, including non-linearities up to O (ε3). The strain-displacement relationship derived is compared with expressions derived by earlier investigators and the errors and the approximations made in some of those are brought out. The equations of motion are solved under the inextensionality condition to obtain the influence of the amplitude on the fundamental flapwise natural frequency of the rotating blade. It is found that large finite amplitudes have a softening effect on the flapwise frequency and that this influence becomes stronger at higher speeds of rotation.

Influence of environmental and instrumental variables on the non-invasive prediction of Brix in pineapple using near infrared spectroscopy

The Brix content of pineapple fruit can be non-invasively predicted from the second derivative of near infrared reflectance spectra. Correlations obtained using a NIRSystems 6500 spectrophotometer through multiple linear regression and modified partial least squares analyses using a post-dispersive configuration were comparable with that from a pre-dispersive configuration in terms of accuracy (e.g. coefficient of determination, R2, 0.73; standard error of cross validation, SECV, 1.01°Brix). The effective depth of sample assessed was slightly greater using the post-dispersive technique (about 20 mm for pineapple fruit), as expected in relation to the higher incident light intensity, relative to the pre-dispersive configuration. The effect of such environmental variables as temperature, humidity and external light, and instrumental variables such as the number of scans averaged to form a spectrum, were considered with respect to the accuracy and precision of the measurement of absorbance at 876 nm, as a key term in the calibration for Brix, and predicted Brix. The application of post-dispersive near infrared technology to in-line assessment of intact fruit in a packing shed environment is discussed.

Associations with low rates of postpartum glucose screening after gestational diabetes among Indigenous and non-Indigenous Australian women

OBJECTIVES To explore factors associated with postpartum glucose screening among women with Gestational Diabetes Mellitus (GDM). METHODS A retrospective study using linked records from women with GDM who gave birth at Cairns Hospital in Far North Queensland, Australia, from 1 January 2004 to 31 December 2010. RESULTS The rates of postpartum Oral Glucose Tolerance Test (OGTT) screening, while having increased significantly among both Indigenous* and non-Indigenous women from 2004 to 2010 (HR 1.15 per year, 95%CI 1.08-1.22, p<0.0001), remain low, particularly among Indigenous women (10% versus 27%, respectively at six months postpartum). Indigenous women in Cairns had a longer time to postpartum OGTT than Indigenous women in remote areas (HR 0.58, 0.38-0.71, p=0.01). Non-Indigenous women had a longer time to postpartum OGTT if they: were born in Australia (HR 0.76, 0.59-1.00, 0.05); were aged <25 years (HR 0.45, 0.23-0.89, p=0.02); had parity >5 (HR 0.33, 0.12-0.90, p=0.03); smoked (HR 0.48, 0.31-0.76, p=0.001); and did not breastfeed (HR 0.09, 0.01-0.64, p=0.02). CONCLUSIONS Postpartum diabetes screening rates following GDM in Far North Queensland are low, particularly among Indigenous women, with lower rates seen in the regional centre; and among non-Indigenous women with indicators of low socioeconomic status. IMPLICATIONS Strategies are urgently needed to improve postpartum diabetes screening after GDM that reach women most at risk.

Non-linear flapping vibrations of rotating blades

The non-linear equations of motion of a rotating blade undergoing extensional and flapwise bending vibration are derived, including non-linearities up to O (ε3). The strain-displacement relationship derived is compared with expressions derived by earlier investigators and the errors and the approximations made in some of those are brought out. The equations of motion are solved under the inextensionality condition to obtain the influence of the amplitude on the fundamental flapwise natural frequency of the rotating blade. It is found that large finite amplitudes have a softening effect on the flapwise frequency and that this influence becomes stronger at higher speeds of rotation.

The structure of Norway spruce (Picea abies [L.] Karst.) stems in relation to wood properties of sawn timber

An important challenge in forest industry is to get the appropriate raw material out from the forests to the wood processing industry. Growth and stem reconstruction simulators are therefore increasingly integrated in industrial conversion simulators, for linking the properties of wooden products to the three-dimensional structure of stems and their growing conditions. Static simulators predict the wood properties from stem dimensions at the end of a growth simulation period, whereas in dynamic approaches, the structural components, e.g. branches, are incremented along with the growth processes. The dynamic approach can be applied to stem reconstruction by predicting the three-dimensional stem structure from external tree variables (i.e. age, height) as a result of growth to the current state. In this study, a dynamic growth simulator, PipeQual, and a stem reconstruction simulator, RetroSTEM, are adapted to Norway spruce (Picea abies [L.] Karst.) to predict the three-dimensional structure of stems (tapers, branchiness, wood basic density) over time such that both simulators can be integrated in a sawing simulator. The parameterisation of the PipeQual and RetroSTEM simulators for Norway spruce relied on the theoretically based description of tree structure developing in the growth process and following certain conservative structural regularities while allowing for plasticity in the crown development. The crown expressed both regularity and plasticity in its development, as the vertical foliage density peaked regularly at about 5 m from the stem apex, varying below that with tree age and dominance position (Study I). Conservative stem structure was characterized in terms of (1) the pipe ratios between foliage mass and branch and stem cross-sectional areas at crown base, (2) the allometric relationship between foliage mass and crown length, (3) mean branch length relative to crown length and (4) form coefficients in branches and stem (Study II). The pipe ratio between branch and stem cross-sectional area at crown base, and mean branch length relative to the crown length may differ in trees before and after canopy closure, but the variation should be further analysed in stands of different ages and densities with varying site fertilities and climates. The predictions of the PipeQual and RetroSTEM simulators were evaluated by comparing the simulated values to measured ones (Study III, IV). Both simulators predicted stem taper and branch diameter at the individual tree level with a small bias. RetroSTEM predictions of wood density were accurate. For focusing on even more accurate predictions of stem diameters and branchiness along the stem, both simulators should be further improved by revising the following aspects in the simulators: the relationship between foliage and stem sapwood area in the upper stem, the error source in branch sizes, the crown base development and the height growth models in RetroSTEM. In Study V, the RetroSTEM simulator was integrated in the InnoSIM sawing simulator, and according to the pilot simulations, this turned out to be an efficient tool for readily producing stand scale information about stem sizes and structure when approximating the available assortments of wood products.

Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

Biochemical and structural properties of Potato virus A VPg

The potato virus A (PVA) genome linked protein (VPg) is a multifunctional protein that takes part in vital infection cycle events such as replication and movement of the virus from cell to cell. VPg is attached to the 5´ end of the genome and is carried in the tip structure of the filamentous virus particle. VPg is also the last protein to be cleaved from the polyprotein. VPg interacts with several viral and host proteins and is phosphorylated at several positions. These features indicate a central role in virus epidemiology and a requirement for an efficient but flexible mechanism for switching between different functions. -- This study examines some of the key VPg functions in more detail. Mutations in the positively charged region from Ala38 to Lys44 affected the NTP binding, uridylylation, and in vitro translation inhibition activities of VPg, whereas in vivo translation inhibition was not affected. Some of the data generated in this study implicated the structural flexibility of the protein in functional activities. VPg lacks a rigid structure, which could allow it to adapt conformationally to different functions as needed. A major finding of this study is that PVA VPg belongs to the class of ´intrinsically disordered proteins´ (IDPs). IDPs are a novel protein class that has helped to explain the observed lack of structure. The existence of IDPs clearly shows that proteins can be functional and adapt a native fold without a rigid structure. Evidence for the intrinsic disorder of VPg was provided by CD spectroscopy, NMR, fluorescence spectroscopy, bioinformatic analysis, and limited proteolytic digestion. The structure of VPg resembles that of a molten globule-type protein and has a hydrophobic core domain. Approximately 50% of the protein is disordered and an α-helical stabilization of these regions has been hypothesized. Surprisingly, VPg structure was stabilized in the presence of anionic lipid vesicles. The stabilization was accompanied by a change in VPg structure and major morphological modifications of the vesicles, including a pronounced increase in the size and appearance of pore or plaque like formations on the vesicle surface. The most likely scenario seems to be an α-helical stabilization of VPg which induces formation of a pore or channel-like structure on the vesicle surface. The size increase is probably due to fusion or swelling of the vesicles. The latter hypothesis is supported by the evident disruption of the vesicles after prolonged incubation with VPg. A model describing the results is presented and discussed in relation to other known properties of the protein.

Stability of large flexible damped spacecraft modeled as elastic continua

Vibrational stability of a large flexible, structurally damped spacecraft subject to large rigid body rotations is analysed modelling the system as an elastic continuum. Using solution of rigid body attitude motion under torque free conditions and modal analysis, the vibrational equations are reduced to ordinary differential equations with time-varying coefficients. Stability analysis is carried out using Floquet theory and Sonin-Polya theorem. The cases of spinning and non-spinning spacecraft idealized as a flexible beam plate undergoing simple structural vibration are analysed in detail. The critical damping required for stabilization is shown to be a function of the spacecraft's inertia ratio and the level of disturbance.

Korteweg-De Vries Equation for Non-Linear Ion-Acoustic Waves in a Plasma with Negative Ions

The theoretical analysis, based on the perturbation technique, of ion-acoustic waves in the vicinity of a Korteweg-de Vries (K-dV) equation derived in a plasma with some negative ions has been made. The investigation shows that the negative ions in plasma with isothermal electrons introduced a critical concentration at which the ion-acoustic wave plays an important role of wave-breaking and forming a precursor while the plasma with non-isothermal electrons has no such singular behaviour of the wave. These two distinct features of ion waves lead to an overall different approach of present study of ion-waves. A distinct feature of non-uniform transition from the nonisothermal case to isothermal case has been shown. Few particular plasma models have been chosen to show the characteristics behaviour of the ion-waves existing in different cases

Long-term results of obstetric brachial plexus surgery

Background: Brachial plexus birth palsy (BPBP) most often occurs as a result of foetal-maternal disproportion. The C5 and C6 nerve roots of the brachial plexus are most frequently affected. In contrast, roots from the C7 to Th1 that result in total injury together with C5 and C6 injury, are affected in fewer than half of the patients. BPBP was first described by Smellie in 1764. Erb published his classical description of the injury in 1874 and his name became linked with the paralysis that is associated with upper root injury. Since then, early results of brachial plexus surgery have been reasonably well documented. However, from a clinical point of view not all primary results are maintained and there is also a need for later follow-up results. In addition most of the studies that are published emanate from highly specialized clinics and no nation wide epidemiological reports are available. One of the plexus injuries is the avulsion type, in which the nerve root or roots are ruptured at the neural cord. It has been speculated whether this might cause injury to the whole neural system or whether shoulder asymmetry and upper limb inequality results in postural deformities of the spine. Alternatively, avulsion could manifest as other signs and symptoms of the whole musculoskeletal system. In addition, there is no available information covering activities of daily living after obstetric brachial plexus surgery. Patients and methods: This was a population-based cross-sectional study on all patients who had undergone brachial plexus surgery with at least 5 years of follow-up. An incidence of 3.05/1000 for BPBP was obtained from the registers for this study period. A total of 1706 BPBP patients needing hospital treatment out of 1 717 057 newborns were registered in Finland between 1971 and 1997 inclusive. Of these BPBP patients, 124 (7.3%) underwent brachial plexus surgery at a mean age of 2.8 months (range: 0.4―13.2 months). Surgery was most often performed by direct neuroraphy after neuroma resection (53%). Depending on the phase of the study, 105 to 112 patients (85-90%) participated in a clinical and radiological follow-up assessment. The mean follow up time exceeded 13 years (range: 5.0―31.5 years). Functional status of the upper extremity was evaluated using Mallet, Gilbert and Raimondi scales. Isometric strength of the upper limb, sensation of the hand and stereognosis were evaluated for both the affected and unaffected sides then the differences and their ratios were calculated and recorded. In addition to the upper extremity, assessment of the spine and lower extremities were performed. Activities of daily living (ADL), participation in normal physical activities, and the use of physiotherapy and occupational therapy were recorded in a questionnaire. Results: The unaffected limb functioned as the dominant hand in all, except four patients. The mean length of the affected upper limb was 6 cm (range: 1-13.5 cm) shorter in 106 (95%) patients. Shoulder function was recorded as a mean Mallet score of 3 (range: 2―4) which was moderate. Both elbow function and hand function were good. The mean Gilbert elbow scale value was 3 (range: -1―5) and the mean Raimondi hand scale was 4 (range:1―5). One-third of the patients experienced pain in the affected limb including all those patients (n=9) who had clavicular non-union resulting from surgery. A total of 61 patients (57%) had an active shoulder external rotation of less than 0° and an active elbow extension deficiency was noted in 82 patients (77%) giving a mean of 26° (range: 5°―80°). In all, expect two patients, shoulder external rotation strength at a mean ratio 35% (range: 0―83%) and in all patients elbow flexion strength at a mean ratio of 41% (range: 0―79%) were impaired compared to the unaffected side. According to radiographs, incongruence of the glenohumeral joint was noted in 15 (16%) patients, whereas incongruence of the radiohumeral joint was found in 20 (21%) patients. Fine sensation was normal for 34/49 (69%) patients with C5-6 injury, for 15/31 (48%) with C5-7 and for only 8/25 (32%) of patients with total injury. Loss of protective sensation or absent sensation was noted in some palmar areas of the hand for 12/105 patients (11%). Normal stereognosis was recorded for 88/105 patients (84%). No significant inequalities in leg length were found and the incidence of structural scoliosis (1.7%) did not differ from that of the reference population. Nearly half of the patients (43%) had asynchronous motion of the upper limbs during gait, which was associated with impaired upper limb function. Data obtained from the completed questionnaires indicated that two thirds (63%) of the patients were satisfied with the functional outcome of the affected hand although one third of all patients needed help with ADL. Only a few patients were unable to participate in physical activities such as: bicycling, cross-country skiing or swimming. However, 71% of the patients reported problems related to the affected upper limb, such as muscle weakness and/or joint stiffness during the aforementioned activities. Incongruity of the radiohumeral joints, extent of the injury, avulsion type injury, age less than three months of age at the time of plexus surgery and inexperience of the surgeon was related to poor results as determined by multivariate analyses. Conclusions: Most of the patients had persistent sequelae, especially of shoulder function. Almost all measurements for the total injury group were poorer compared with those of the C5-6 type injury group. Most of the patients had asymmetry of the shoulder region and a shorter affected upper limb, which is a probable reason for having an abnormal gait. However, BPBP did not have an effect on normal growth of the lower extremities or the spine. Although, participation in physical activities was similar to that of the normal population, two-thirds of the patients reported problems. One-third of the patients needed help with ADL. During the period covered by this study, 7.3% BPBP of patients that needed hospital treatment had a brachial plexus operation, which amounts to fewer than 10 operations per year in Finland. It seems that better results of obstetric plexus surgery and more careful follow-up including opportunities for late reconstructive procedures will be expected, if the treatment is solely concentrated on by a few specialised teams.

Solution structure of O-glycosylated C-terminal leucine zipper domain of human salivary mucin (MUC7)

Solution structures of a 23 residue glycopeptide II (KIS* RFLLYMKNLLNRIIDDMVEQ, where * denotes the glycan Gal-beta-(1-3)-alpha-GalNAc) and its deglycosylated counterpart I derived from the C-terminal leucine zipper domain of low molecular weight human salivary mucin (MUC7) were studied using CD, NMR spectroscopy and molecular modeling. The peptide I was synthesized using the Fmoc chemistry following the conventional procedure and the glycopeptide II was synthesized incorporating the O-glycosylated building block (N alpha-Fmoc-Ser-[Ac-4,-beta-D-Gal-(1,3)-Ac(2)alpha-D-GalN(3)]-OPfp) at the appropriate position in stepwise assembly of peptide chain. Solution structures of these glycosylated and nonglycosylated peptides were studied in water and in the presence of 50% of an organic cosolvent, trifluoroethanol (TFE) using circular dichroism (CD), and in 50% TFE using two-dimensional proton nuclear magnetic resonance (2D H-1 NMR) spectroscopy. CD spectra in aqueous medium indicate that the apopeptide I adapts, mostly, a beta-sheet conformation whereas the glycopeptide II assumes helical structure. This transition in the secondary structure, upon glycosylation, demonstrates that the carbohydrate moiety exerts significant effect on the peptide backbone conformation. However, in 50% TFE both the peptides show pronounced helical structure. Sequential and medium range NOEs, C alpha H chemical shift perturbations, (3)J(NH:C alpha H) couplings and deuterium exchange rates of the amide proton resonances in water containing 50% TFE indicate that the peptide I adapts alpha-helical structure from Ile2-Val21 and the glycopeptide II adapts alpha-helical structure from Ser3-Glu22. The observation of continuous stretch of helix in both the peptides as observed by both NMR and CD spectroscopy strongly suggests that the C-terminal domain of MUC7 with heptad repeats of leucines or methionine residues may be stabilized by dimeric leucine zipper motif. The results reported herein may be invaluable in understanding the aggregation (or dimerization) of MUC7 glycoprotein which would eventually have implications in determining its structure-function relationship.

Low-Complexity Near-ML Decoding of Large Non-Orthogonal STBCs using Reactive Tabu Search

Non-orthogonal space-time block codes (STBC) with large dimensions are attractive because they can simultaneously achieve both high spectral efficiencies (same spectral efficiency as in V-BLAST for a given number of transmit antennas) as well as full transmit diversity. Decoding of non-orthogonal STBCs with large dimensions has been a challenge. In this paper, we present a reactive tabu search (RTS) based algorithm for decoding non-orthogonal STBCs from cyclic division algebras (CDA) having largedimensions. Under i.i.d fading and perfect channel state information at the receiver (CSIR), our simulation results show that RTS based decoding of 12 X 12 STBC from CDA and 4-QAM with 288 real dimensions achieves i) 10(-3) uncoded BER at an SNR of just 0.5 dB away from SISO AWGN performance, and ii) a coded BER performance close to within about 5 dB of the theoretical MIMO capacity, using rate-3/4 turbo code at a spectral efficiency of 18 bps/Hz. RTS is shown to achieve near SISO AWGN performance with less number of dimensions than with LAS algorithm (which we reported recently) at some extra complexity than LAS. We also report good BER performance of RTS when i.i.d fading and perfect CSIR assumptions are relaxed by considering a spatially correlated MIMO channel model, and by using a training based iterative RTS decoding/channel estimation scheme.