Finite difference time domain (FDTD) method for modeling the effect of switched gradients on the human body in MRI

Numerical modeling of the eddy currents induced in the human body by the pulsed field gradients in MRI presents a difficult computational problem. It requires an efficient and accurate computational method for high spatial resolution analyses with a relatively low input frequency. In this article, a new technique is described which allows the finite difference time domain (FDTD) method to be efficiently applied over a very large frequency range, including low frequencies. This is not the case in conventional FDTD-based methods. A method of implementing streamline gradients in FDTD is presented, as well as comparative analyses which show that the correct source injection in the FDTD simulation plays a crucial rule in obtaining accurate solutions. In particular, making use of the derivative of the input source waveform is shown to provide distinct benefits in accuracy over direct source injection. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and the source injection method has been verified against examples with analytical solutions. Results are presented showing the spatial distribution of gradient-induced electric fields and eddy currents in a complete body model.

Treatment of non-resectable hepatocellular carcinoma with autologous tumor-pulsed dendritic cells

Background: The response of hepatocellular carcinoma (HCC) to therapy is often disappointing and new modalities of treatment are clearly needed. Active immunotherapy based on the injection of autologous dendritic cells (DC) co-cultured ex vivo with tumor antigens has been used in pilot studies in various malignancies such as melanoma and lymphoma with encouraging results. Methods: In the present paper, the preparation and exposure of patient DC to autologous HCC antigens and re-injection in an attempt to elicit antitumor immune responses are described. Results: Therapy was given to two patients, one with hepatitis C and one with hepatitis B, who had large, multiple HCC and for whom no other therapy was available. No significant side-effects were observed. The clinical course was unchanged in one patient, who died a few months later. The other patient, whose initial prognosis was considered poor, is still alive and well more than 3 years later with evidence of slowing of tumor growth based on organ imaging. Conclusions: It is concluded that HCC may be a malignancy worthy of DC trials and sufficient details in the present paper are given for the protocol to be copied or modified. (C) 2002 Blackwell Publishing Asia Pty Ltd.

The pyramidal cell of the sensorimotor cortex of the macaque monkey: Phenotypic variation

Recent studies have revealed striking differences in pyramidal cell structure among cortical regions involved in the processing of different functional modalities. For example, cells involved in visual processing show systematic variation, increasing in morphological complexity with rostral progression from V1 through extrastriate areas. Differences have also been identified between pyramidal cells in somatosensory, motor and prefrontal cortex, but the extent to which the pyramidal cell phenotype may vary between these functionally related cortical regions remains unknown. In the present study we investigated the structure of layer III pyramidal cells in somatosensory and motor areas 3b, 4, 5, 6 and 7b of the macaque monkey. Cells were intracellularly injected in fixed, flat-mounted cortical slices and analysed for morphometric parameters. The size of the basal dendritic arbours, the number of their branches and their spine density were found to vary systematically between areas. Namely, we found a trend for increasing complexity in dendritic arbour structure through areas 3b, 5 and 7b. A similar trend occurred through areas 4 and 6. The differences in arbour structure may determine the number of inputs received by neurons and may thus be an important factor in determining function at the cellular and systems level.

Independent and conjugate eye movements during optokinesis in teleost fish

In response to movements involving a large part of the visual field, the eyes of vertebrates typically show an optokinetic nystagmus, a response in which both eyes are tightly yoked. Using a comparative approach, this study sets out to establish whether fish with independent spontaneous eye movements show independent optokinetic nystagmus in each eye. Two fish with independent spontaneous eye movements, the pipefish Corythoichthyes intestinalis and the sandlance Limnichthyes fasciatus were compared with the butterflyfish Chaetodon rainfordi, which exhibits tightly yoked eye movements. In the butterflyfish a single whole-field stimulus elicits conjugate optokinesis, whereas the sandlance and pipefish show asynchronous optokinetic movements. In a split drum experiment, when both eyes were stimulated in opposite directions with different speeds, both the sandlance and the pipefish compensated independently with each eye. The optokinetic response in the butterflyfish showed some disconjugacy but was generally confused. When one eye was occluded, the seeing eye was capable of driving the occluded eye in both the butterflyfish and the pipefish but not in the sandlance. Monocular occlusion therefore unmasks a link between the two eyes in the pipefish, which is overridden when both eyes receive visual input. The sandlance never showed any correlation between the eyes during optokinesis in all stimulus conditions. This suggests that there are different levels of linkage between the two eyes in the oculomotor system of teleosts, depending on the visual input.

Light and electron microscopic analysis of aquaporin 1-like-immunoreactive amacrine cells in the rat retina

Aquaporin 1 (AQP1; also known as CHIP, a channel-forming integral membrane protein of 28 kDa) is the first protein to be shown to function as a water channel and has been recently shown to be present in the rat retina. We previously showed (Kim et al. [1998] Neurosci Lett 244:52-54) that AQP1-like immunoreactivity is present in a certain population of amacrine cells in the rat retina. This study was conducted to characterize these cells in more detail, With immunocytochemistry using specific antisera against AQP1, whole-mount preparations and 50-mum-thick vibratome sections were examined by light and electron microscopy. These cells were a class of amacrine cells, which had symmetric bistratified dendritic trees ramified in stratum 2 and in the border of strata 3 and 4 of the inner plexiform layer (IPL). Their dendritic field diameters ranged from 90 to 230 mum. Double labeling with antisera against AQP1 and gamma-aminobutyric acid or glycine demonstrated that these AQP1-like-immunoreactive amacrine cells were immunoreactive for glycine. Their most frequent synaptic input was from other amacrine cell processes in both sublaminae a and b of the IPL, followed by a few cone bipolar cells. Their primary targets were other amacrine cells and ganglion cells in both sublaminae a and b of the IPL. In addition, synaptic output Onto bipolar cells was rarely observed in sublamina b of the IPL. Thus, the AQP1 antibody labels a class of glycinergic amacrine cells with small to medium-sized dendritic fields in the rat retina. (C) 2002 Wiley-Liss, Inc.

Diverse synaptic mechanisms generate direction selectivity in the rabbit retina

The synaptic conductance of the On-Off direction-selective ganglion cells was measured during visual stimulation to determine whether the direction selectivity is a property of the circuitry presynaptic to the ganglion cells or is generated by postsynaptic interaction of excitatory and inhibitory inputs. Three synaptic asymmetries were identified that contribute to the generation of direction-selective responses: (1) a presynaptic mechanism producing stronger excitation in the preferred direction, (2) a presynaptic mechanism producing stronger inhibition in the opposite direction, and (3) postsynaptic interaction of excitation with spatially offset inhibition. Although the on- and off-responses showed the same directional tuning, the off-response was generated by all three mechanisms, whereas the on- response was generated primarily by the two presynaptic mechanisms. The results indicate that, within a single neuron, different strategies are used within distinct dendritic arbors to accomplish the same neural computation.

The role of gender in reactions to service failure and recovery

Male and female consumers place different emphasis on elements of the service recovery process. Perceptions were influenced by gender of the service provider and by a match of customer and service provider gender. The study, an experimental design with 712 respondents, found that when service providers, irrespective of gender, display concern and give customers voice and a sizable compensation, both men and women reported more positive attitudes compared with when this was not so. Combinations of high voice with high outcome and high voice with high concern were especially important in positively influencing perceptions of effort, regardless of gender. However, the authors also found that there were significant differences between male and female respondents regarding their perceptions of how service recovery should be handled. Women want their views heard during service recovery attempts and to be allowed to provide input. Men, in contrast, do not view voice as important.

Capillaries within compartments: Microvascular interpretation of dynamic positron emission tomography data

Measurement of exchange of substances between blood and tissue has been a long-lasting challenge to physiologists, and considerable theoretical and experimental accomplishments were achieved before the development of the positron emission tomography (PET). Today, when modeling data from modern PET scanners, little use is made of earlier microvascular research in the compartmental models, which have become the standard model by which the vast majority of dynamic PET data are analysed. However, modern PET scanners provide data with a sufficient temporal resolution and good counting statistics to allow estimation of parameters in models with more physiological realism. We explore the standard compartmental model and find that incorporation of blood flow leads to paradoxes, such as kinetic rate constants being time-dependent, and tracers being cleared from a capillary faster than they can be supplied by blood flow. The inability of the standard model to incorporate blood flow consequently raises a need for models that include more physiology, and we develop microvascular models which remove the inconsistencies. The microvascular models can be regarded as a revision of the input function. Whereas the standard model uses the organ inlet concentration as the concentration throughout the vascular compartment, we consider models that make use of spatial averaging of the concentrations in the capillary volume, which is what the PET scanner actually registers. The microvascular models are developed for both single- and multi-capillary systems and include effects of non-exchanging vessels. They are suitable for analysing dynamic PET data from any capillary bed using either intravascular or diffusible tracers, in terms of physiological parameters which include regional blood flow. (C) 2003 Elsevier Ltd. All rights reserved.

T cell regulation of the immune response to infection in periodontal diseases

Although T cells have been implicated in the pathogenesis and are considered to be central both in progression and control of the chronic inflammatory periodontal diseases, the precise contribution of T cells to the regulation of tissue destruction has not been fully elucidated. Current dogma suggests that immunity to infection is controlled by distinct T helper 1 (Th1) and T helper 2 (Th2) subsets of T cells classified on the basis of their cytokine profile. Further, a subset of T cells with immunosuppressive function and cytokine profile distinct from Th1 or Th2 has been described and designated as regulatory T cells. Although these regulatory T cells have been considered to maintain self-tolerance resulting in the suppression of auto-immune responses, recent data suggest that these cells may also play a role in preventing infection-induced immunopathology. In this review, the role of functional and regulatory T cells in chronic inflammatory periodontal diseases will be summarized. This should not only provide an insight into the relationship between the immune response to periodontopathic bacteria and disease but should also highlight areas of development for potentially new therapeutic modalities.

Impulse-response function of splanchnic circulation with model-independent constraints: theory and experimental validation

Modeling physiological processes using tracer kinetic methods requires knowledge of the time course of the tracer concentration in blood supplying the organ. For liver studies, however, inaccessibility of the portal vein makes direct measurement of the hepatic dual-input function impossible in humans. We want to develop a method to predict the portal venous time-activity curve from measurements of an arterial time-activity curve. An impulse-response function based on a continuous distribution of washout constants is developed and validated for the gut. Experiments with simultaneous blood sampling in aorta and portal vein were made in 13 anesthetized pigs following inhalation of intravascular [O-15] CO or injections of diffusible 3-O[ C-11] methylglucose (MG). The parameters of the impulse-response function have a physiological interpretation in terms of the distribution of washout constants and are mathematically equivalent to the mean transit time ( T) and standard deviation of transit times. The results include estimates of mean transit times from the aorta to the portal vein in pigs: (T) over bar = 0.35 +/- 0.05 min for CO and 1.7 +/- 0.1 min for MG. The prediction of the portal venous time-activity curve benefits from constraining the regression fits by parameters estimated independently. This is strong evidence for the physiological relevance of the impulse-response function, which includes asymptotically, and thereby justifies kinetically, a useful and simple power law. Similarity between our parameter estimates in pigs and parameter estimates in normal humans suggests that the proposed model can be adapted for use in humans.

Mean anisotropy of homogeneous Gaussian random fields and anisotropic norms of linear translation-invariant operators on multidimensional integer lattices

Sensitivity of output of a linear operator to its input can be quantified in various ways. In Control Theory, the input is usually interpreted as disturbance and the output is to be minimized in some sense. In stochastic worst-case design settings, the disturbance is considered random with imprecisely known probability distribution. The prior set of probability measures can be chosen so as to quantify how far the disturbance deviates from the white-noise hypothesis of Linear Quadratic Gaussian control. Such deviation can be measured by the minimal Kullback-Leibler informational divergence from the Gaussian distributions with zero mean and scalar covariance matrices. The resulting anisotropy functional is defined for finite power random vectors. Originally, anisotropy was introduced for directionally generic random vectors as the relative entropy of the normalized vector with respect to the uniform distribution on the unit sphere. The associated a-anisotropic norm of a matrix is then its maximum root mean square or average energy gain with respect to finite power or directionally generic inputs whose anisotropy is bounded above by a≥0. We give a systematic comparison of the anisotropy functionals and the associated norms. These are considered for unboundedly growing fragments of homogeneous Gaussian random fields on multidimensional integer lattice to yield mean anisotropy. Correspondingly, the anisotropic norms of finite matrices are extended to bounded linear translation invariant operators over such fields.

Predicting cave initiation and propagation

The most widely used method for predicting the onset of continuous caving is Laubscher's caving chart. A detailed examination of this method was undertaken which concluded that it had limitations which may impact on results, particularly when dealing with stronger rock masses that are outside current experience. These limitations relate to inadequate guidelines for adjustment factors to rock mass rating (RMR), concerns about the position on the chart of critical case history data, undocumented changes to the method and an inadequate number of data points to be confident of stability boundaries. A review was undertaken on the application and reliability of a numerical method of assessing cavability. The review highlighted a number of issues, which at this stage, make numerical continuum methods problematic for predicting cavability. This is in particular reference to sensitivity to input parameters that are difficult to determine accurately and mesh dependency. An extended version of the Mathews method for open stope design was developed as an alternative method of predicting the onset of continuous caving. A number of caving case histories were collected and analyzed and a caving boundary delineated statistically on the Mathews stability graph. The definition of the caving boundary was aided by the existence of a large and wide-ranging stability database from non-caving mines. A caving rate model was extrapolated from the extended Mathews stability graph but could only be partially validated due to a lack of reliable data.

Critical comparison of Julius Kruttschnitt Mineral Research Centre (JKMRC) blast fragmentation models

Blast fragmentation can have a significant impact on the profitability of a mine. An optimum run of mine (ROM) size distribution is required to maximise the performance of downstream processes. If this fragmentation size distribution can be modelled and controlled, the operation will have made a significant advancement towards improving its performance. Blast fragmentation modelling is an important step in Mine to Mill™ optimisation. It allows the estimation of blast fragmentation distributions for a number of different rock mass, blast geometry, and explosive parameters. These distributions can then be modelled in downstream mining and milling processes to determine the optimum blast design. When a blast hole is detonated rock breakage occurs in two different stress regions - compressive and tensile. In the-first region, compressive stress waves form a 'crushed zone' directly adjacent to the blast hole. The second region, termed the 'cracked zone', occurs outside the crush one. The widely used Kuz-Ram model does not recognise these two blast regions. In the Kuz-Ram model the mean fragment size from the blast is approximated and is then used to estimate the remaining size distribution. Experience has shown that this model predicts the coarse end reasonably accurately, but it can significantly underestimate the amount of fines generated. As part of the Australian Mineral Industries Research Association (AMIRA) P483A Mine to Mill™ project, the Two-Component Model (TCM) and Crush Zone Model (CZM), developed by the Julius Kruttschnitt Mineral Research Centre (JKMRC), were compared and evaluated to measured ROM fragmentation distributions. An important criteria for this comparison was the variation of model results from measured ROM in the-fine to intermediate section (1-100 mm) of the fragmentation curve. This region of the distribution is important for Mine to Mill™ optimisation. The comparison of modelled and Split ROM fragmentation distributions has been conducted in harder ores (UCS greater than 80 MPa). Further work involves modelling softer ores. The comparisons will be continued with future site surveys to increase confidence in the comparison of the CZM and TCM to Split results. Stochastic fragmentation modelling will then be conducted to take into account variation of input parameters. A window of possible fragmentation distributions can be compared to those obtained by Split . Following this work, an improved fragmentation model will be developed in response to these findings.

New fluorescence parameters for monitoring photosynthesis in plants - 1. The effect of illumination on the fluorescence parameters of the JIP-test

Chlorophyll fluorescence measurements have a wide range of applications from basic understanding of photosynthesis functioning to plant environmental stress responses and direct assessments of plant health. The measured signal is the fluorescence intensity (expressed in relative units) and the most meaningful data are derived from the time dependent increase in fluorescence intensity achieved upon application of continuous bright light to a previously dark adapted sample. The fluorescence response changes over time and is termed the Kautsky curve or chlorophyll fluorescence transient. Recently, Strasser and Strasser (1995) formulated a group of fluorescence parameters, called the JIP-test, that quantify the stepwise flow of energy through Photosystem II, using input data from the fluorescence transient. The purpose of this study was to establish relationships between the biochemical reactions occurring in PS II and specific JIP-test parameters. This was approached using isolated systems that facilitated the addition of modifying agents, a PS II electron transport inhibitor, an electron acceptor and an uncoupler, whose effects on PS II activity are well documented in the literature. The alteration to PS II activity caused by each of these compounds could then be monitored through the JIP-test parameters and compared and contrasted with the literature. The known alteration in PS II activity of Chenopodium album atrazine resistant and sensitive biotypes was also used to gauge the effectiveness and sensitivity of the JIP-test. The information gained from the in vitro study was successfully applied to an in situ study. This is the first in a series of four papers. It shows that the trapping parameters of the JIP-test were most affected by illumination and that the reduction in trapping had a run-on effect to inhibit electron transport. When irradiance exposure proceeded to photoinhibition, the electron transport probability parameter was greatly reduced and dissipation significantly increased. These results illustrate the advantage of monitoring a number of fluorescence parameters over the use of just one, which is often the case when the F-V/F-M ratio is used.

On the induced electric field gradients in the human body for magnetic stimulation by gradient coils in MRI

Prior theoretical studies indicate that the negative spatial derivative of the electric field induced by magnetic stimulation may he one of the main factors contributing to depolarization of the nerve fiber. This paper studies this parameter for peripheral nerve stimulation (PNS) induced by time.-varying gradient fields during MRI scans. The numerical calculations are based on an efficient, quasi-static, finite-difference scheme and an anatomically realistic human, full-body model. Whole-body cylindrical and planar gradient sets in MRI systems and various input signals have been explored. The spatial distributions of the induced electric field and their gradients are calculated and attempts are made to correlate these areas with reported experimental stimulation data. The induced electrical field pattern is similar for both the planar coils and cylindrical coils. This study provides some insight into the spatial characteristics of the induced field gradients for PNS in MRI, which may be used to further evaluate the sites where magnetic stimulation is likely to occur and to optimize gradient coil design.