Quantitative Analysis of Kilohertz-Frequency Neurostimulation

Mainstream electrical stimulation therapies, e.g., spinal cord stimulation (SCS) and deep brain stimulation, use pulse trains that are delivered at rates no higher than 200 Hz. In recent years, stimulation of nerve fibers using kilohertz-frequency (KHF) signals has received increased attention due to the potential to penetrate deeper in the tissue and to the ability to block conduction of action potentials. As well, there are a growing number of clinical applications that use KHF waveforms, including transcutaneous electrical stimulation (TES) for overactive bladder and SCS for chronic pain. However, there is a lack of fundamental understanding of the mechanisms of action of KHF stimulation. The goal of this research was to analyze quantitatively KHF neurostimulation.

We implemented a multilayer volume conductor model of TES including dispersion and capacitive effects, and we validated the model with in vitro measurements in a phantom constructed from dispersive materials. We quantified the effects of frequency on the distribution of potentials and fiber excitation. We also quantified the effects of a novel transdermal amplitude modulated signal (TAMS) consisting of a non-zero offset sinusoidal carrier modulated by a square-pulse train. The model revealed that high-frequency signals generated larger potentials at depth than did low frequencies, but this did not translate into lower stimulation thresholds. Both TAMS and conventional rectangular pulses activated more superficial fibers in addition to the deeper, target fibers, and at no frequency did we observe an inversion of the strength-distance relationship. In addition, we performed in vivo experiments and applied direct stimulation to the sciatic nerve of cats and rats. We measured electromyogram and compound action potential activity evoked by pulses, TAMS and modified versions of TAMS in which we varied the amplitude of the carrier. Nerve fiber activation using TAMS showed no difference with respect to activation with conventional pulse for carrier frequencies of 20 kHz and higher, regardless the size of the carrier. Therefore, TAMS with carrier frequencies >20 kHz does not offer any advantage over conventional pulses, even with larger amplitudes of the carrier, and this has implications for design of waveforms for efficient and effective TES.

We developed a double cable model of a dorsal column (DC) fiber to quantify the responses of DC fibers to a novel KHF-SCS signal. We validated the model using in vivo recordings of the strength-duration relationship and the recovery cycle of single DC fibers. We coupled the fiber model to a model of SCS in human and applied the KHF-SCS signal to quantify thresholds for activation and conduction block for different fiber diameters at different locations in the DCs. Activation and block thresholds increased sharply as the fibers were placed deeper in the DCs, and decreased for larger diameter fibers. Activation thresholds were > 5 mA in all cases and up to five times higher than for conventional (~ 50 Hz) SCS. For fibers exhibiting persistent activation, the degree of synchronization of the firing activity to the KHF-SCS signal, as quantified using the vector strength, was low for a broad amplitude range, and the dissimilarity between the activities in pairs of fibers, as quantified using the spike time distance, was high and decreased for more closely positioned fibers. Conduction block thresholds were higher than 30 mA for all fiber diameters at any depth and well above the amplitudes used clinically (0.5 – 5 mA). KHF-SCS appears to activate few, large, superficial fibers, and the activated fibers fire asynchronously to the stimulation signal and to other activated fibers.

The outcomes of this work contribute to the understanding of KHF neurostimulation by establishing the importance of the tissue filtering properties on the distribution of potentials, assessing quantitatively the impact of KHF stimulation on nerve fiber excitation, and developing and validating a detailed model of a DC fiber to characterize the effects of KHF stimulation on DC axons. The results have implications for design of waveforms for efficient and effective nerve fiber stimulation in the peripheral and central nervous system.

Statistical considerations in the kinetic analysis of PET-FDG brain tumour studies

Dynamic positron emission tomography (PET) imaging can be used to track the distribution of injected radio-labelled molecules over time in vivo. This is a powerful technique, which provides researchers and clinicians the opportunity to study the status of healthy and pathological tissue by examining how it processes substances of interest. Widely used tracers include 18F-uorodeoxyglucose, an analog of glucose, which is used as the radiotracer in over ninety percent of PET scans. This radiotracer provides a way of quantifying the distribution of glucose utilisation in vivo. The interpretation of PET time-course data is complicated because the measured signal is a combination of vascular delivery and tissue retention effects. If the arterial time-course is known, the tissue time-course can typically be expressed in terms of a linear convolution between the arterial time-course and the tissue residue function. As the residue represents the amount of tracer remaining in the tissue, this can be thought of as a survival function; these functions been examined in great detail by the statistics community. Kinetic analysis of PET data is concerned with estimation of the residue and associated functionals such as ow, ux and volume of distribution. This thesis presents a Markov chain formulation of blood tissue exchange and explores how this relates to established compartmental forms. A nonparametric approach to the estimation of the residue is examined and the improvement in this model relative to compartmental model is evaluated using simulations and cross-validation techniques. The reference distribution of the test statistics, generated in comparing the models, is also studied. We explore these models further with simulated studies and an FDG-PET dataset from subjects with gliomas, which has previously been analysed with compartmental modelling. We also consider the performance of a recently proposed mixture modelling technique in this study.

Smear slide analysis of sediment core GeoB5836-2

Laminated sediments spanning the last 20,000 years (though not continuously) in the Shaban Deep, a brine-filled basin in the northern Red Sea, were analyzed microscopically and with backscattered electron imagery in order to determine laminae composition with emphasis on the diatomaceous component. Based on this detailed study, we present schematic models to propose paleoflux scenarios for laminae formation at different time-slices. The investigated core (GeoB 5836-2; 26°12.61'N, 35°21.56'E; water depth 1475 m) shows light and dark alternating laminae that are easily distinguishable in the mid-Holocene and at the end of the deglaciation (13-15 ka) period. Light layers are mainly composed of coccoliths, terrigenous material and diatom fragments, while dark layers consist almost exclusively of diatom frustules (monospecific or mixed assemblages). The regularity in the occurrence of coccolith/diatom couplets points to an annual deposition cycle where contrasting seasons and associated plankton blooms are represented (diatoms-fall/winter deposition, coccoliths-summer signal). We propose that, for the past ~15,000 years, the laminations represent two-season annual varves. Strong dissolution of carbonate, with the concomitant loss of the coccolith-rich layer in sediments older than 15 ka, prevents us from presenting a schematic model of annual deposition. However, the diatomaceous component reveals a marked switch in species composition between Last Glacial Maximum (LGM) sediments (dominated by Chaetoceros resting spores) and sediments somewhat younger (18-19 ka; dominated by Rhizosolenia). We propose that different diatom assemblages reflect changing conditions in stratification in the northern Red Sea: Strong stratification conditions, such as during two meltwater pulses at 14.5 and 11.4 ka, are reflected in the sediment by Rhizosolenia layers, while Chaetoceros-dominated assemblages represent deep convection conditions.

(Table T2) Relative signal intensity of intact polar lipids in sediments from ODP Sites 207-1257 and 207-1258

We report results from the analysis of intact polar lipids (IPLs) in sediments from Ocean Drilling Program Sites 1257 and 1258. IPLs, constituting the cell membranes of living organisms, were detected in organic-lean sediments but not in underlying organic-rich black shales. Microbial activity in organic-lean sediments is likely due to sulfate-dependent oxidation of methane whereas difficulties detecting IPLs in black shales are interpreted to result from unfavorable signal-to-noise ratios due to low cell concentrations in combination with extremely high analytical noise created by uncharacterized organic matrix. IPLs found are consistent with a low-diversity community of archaea and bacteria. The concentrations of IPLs are more than one order of magnitude lower than those in Neogene deep subsurface sediments at the Peruvian margin, suggestive of significantly lower cell concentrations in Demerara Rise. This finding is consistent with inferred low rates of subsurface microbial activity.

Raman Analysis of Dilute Aqueous Samples by Localized Evaporation of Submicroliter Droplets on the Tips of Superhydrophobic Copper Wires

Raman analysis of dilute aqueous solutions is normally prevented by their low signal levels. A very general method to increase the concentration to detectable levels is to evaporate droplets of the sample to dryness, creating solid deposits which are then Raman probed. Here, superhydrophobic (SHP) wires with hydrophilic tips have been used as supports for drying droplets, which have the advantage that the residue is automatically deposited at the tip. The SHP wires were readily prepared in minutes using electroless galvanic deposition of Ag onto copper wires followed by modification with a polyfluorothiol (3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,10-heptadecafluoro-1-decanethiol, HDFT). Cutting the coated wires with a scalpel revealed hydrophilic tips which could support droplets whose maximum size was determined by the wire diameter. Typically, 230 μm wires were used to support 0.6 μL droplets. Evaporation of dilute melamine droplets gave solid deposits which could be observed by scanning electron microscopy (SEM) and Raman spectroscopy. The limit of detection for melamine using a two stage evaporation procedure was 1 × 10^-6 mol dm^-3. The physical appearance of dried droplets of sucrose and glucose showed that the samples retained significant amounts of water, even under high vacuum. Nonetheless, the Raman detection limits of sucrose and glucose were 5 × 10^-4 and 2.5 × 10^-3 mol dm^-3, respectively, which is similar to the sensitivity reported for surface-enhanced Raman spectroscopy (SERS) detection of glucose. It was also possible to quantify the two sugars in mixtures at concentrations which were similar to those found in human blood through multivariate analysis.

Functional analysis of the ATM-p53-p21 pathway in the LRF CLL4 trial: blockade at the level of p21 is associated with short response duration.

PURPOSE: This study sought to establish whether functional analysis of the ATM-p53-p21 pathway adds to the information provided by currently available prognostic factors in patients with chronic lymphocytic leukemia (CLL) requiring frontline chemotherapy. EXPERIMENTAL DESIGN: Cryopreserved blood mononuclear cells from 278 patients entering the LRF CLL4 trial comparing chlorambucil, fludarabine, and fludarabine plus cyclophosphamide were analyzed for ATM-p53-p21 pathway defects using an ex vivo functional assay that uses ionizing radiation to activate ATM and flow cytometry to measure upregulation of p53 and p21 proteins. Clinical endpoints were compared between groups of patients defined by their pathway status. RESULTS: ATM-p53-p21 pathway defects of four different types (A, B, C, and D) were identified in 194 of 278 (70%) samples. The type A defect (high constitutive p53 expression combined with impaired p21 upregulation) and the type C defect (impaired p21 upregulation despite an intact p53 response) were each associated with short progression-free survival. The type A defect was associated with chemoresistance, whereas the type C defect was associated with early relapse. As expected, the type A defect was strongly associated with TP53 deletion/mutation. In contrast, the type C defect was not associated with any of the other prognostic factors examined, including TP53/ATM deletion, TP53 mutation, and IGHV mutational status. Detection of the type C defect added to the prognostic information provided by TP53/ATM deletion, TP53 mutation, and IGHV status. CONCLUSION: Our findings implicate blockade of the ATM-p53-p21 pathway at the level of p21 as a hitherto unrecognized determinant of early disease recurrence following successful cytoreduction.

Damping Torque Analysis of Power Systems with DFIGs for Wind Power Generation

A grid-connected DFIG for wind power generation can affect power system small-signal angular stability in two ways: by changing the system load flow condition and dynamically interacting with synchronous generators (SGs). This paper presents the application of conventional method of damping torque analysis (DTA) to examine the effect of DFIG’s dynamic interactions with SGs on the small-signal angular stability. It shows that the effect is due to the dynamic variation of power exchange between the DFIG and power system and can be estimated approximately by the DTA. Consequently, if the DFIG is modelled as a constant power source when the effect of zero dynamic interactions is assumed, the impact of change of load flow brought about by the DFIG can be determined. Thus the total effect of DFIG can be estimated from the result of DTA added on that of constant power source model. Applications of the DTA method proposed in the paper are discussed. An example of multi-machine power systems with grid-connected DFIGs are presented to demonstrate and validate the DTA method proposed and conclusions obtained in the paper.

I/Q imbalance in two-way AF relaying: Performance analysis and detection mode switch

This paper studies the impact of in-phase and quadrature-phase imbalance (IQI) in two-way amplify-and-forward (AF) relaying systems. In particular, the effective signal-to-interference-plus-noise ratio (SINR) is derived for each source node, considering four different linear detection schemes, namely, uncompensated (Uncomp) scheme, maximal-ratio-combining (MRC), zero-forcing (ZF) and minimum mean-square error (MMSE) based schemes. For each proposed scheme, the outage probability (OP) is investigated over independent, non-identically distributed Nakagami-m fading channels, and exact closed-form expressions are derived for the first three schemes. Based on the closed-form OP expressions, an adaptive detection mode switching scheme is designed for minimizing the OP of both sources. An important observation is that, regardless of the channel conditions and transmit powers, the ZF-based scheme should always be selected if the target SINR is larger than 3 (4.77dB), while the MRC-based scheme should be avoided if the target SINR is larger than 0.38 (-4.20dB).

Energy efficiency analysis of rank-1 Ricean fading MIMO channels

This paper studies the energy efficiency (EE) of a point-to-point rank-1 Ricean fading multiple-input-multiple-output (MIMO) channel. In particular, a tight lower bound and an asymptotic approximation for the EE of the considered MIMO system are presented, under the assumption that the channel is unknown at the transmitter and perfectly known at the receiver. Moreover, the effects of different system parameters, namely, transmit power, spectral efficiency (SE), and number of transmit and receive antennas, on the EE are analytically investigated. An important observation is that, in the high signal-to-noise ratio regime and with the other system parameters fixed, the optimal transmit power that maximizes the EE increases as the Ricean-K factor increases. On the contrary, the optimal SE and the optimal number of transmit antennas decrease as K increases.

I/Q imbalance in two-way AF relaying: Power allocation and performance analysis

We investigate the performance of dual-hop two-way amplify-and-forward (AF) relaying in the presence of inphase and quadrature-phase imbalance (IQI) at the relay node. In particular, the effective signal-to-interference-plus-noise ratio (SINR) at both sources is derived. These SINRs are used to design an instantaneous power allocation scheme, which maximizes the minimum SINR of the two sources under a total transmit power constraint. The solution to this optimization problem is analytically determined and used to evaluate the outage probability (OP) of the considered two-way AF relaying system. Both analytical and numerical results show that IQI can create fundamental performance limits on two-way relaying, which cannot be avoided by simply improving the channel conditions.

Evaluation of Two Ultrasonic Systems for Analysis of Porosity in Ceramic

The Ultrasound Laboratory of the Nuclear Engineering Institute (LABUS / IEN) has developed an ultrasonic technique to measure porosity in nuclear fuel pellets (UO2). By difficulties related to the handling of UO2 pellets, Alumina (Al2O3) pellets have been used in preliminary tests, until a methodology for tests with pellets of UO2 could be defined. In a previous work, in which a contact ultrasonic technique was used, good results were obtained to measure the porosity of Alumina pellets. In the current studies, it was found that the frequency spectrum of an ultrasonic pulse is very sensitive to the porosity of the medium in which it propagates. In order to define the most appropriate experimental apparatus for using immersion technique in future tests, two ultrasonic systems, available in LABUS, which permit to work with the ultrasonic pulse in the frequency domain were evaluated . One system was the Explorer II (Matec INSTRUMENTS) and the other the ultrasonic pulse generator Epoch 4 Plus (Panametrics) coupled with an oscilloscope TDS 3032B (Tektronix). For this evaluation, several frequency spectra were obtained with the two equipment, by the passage of the ultrasonic wave in the same pellet of Alumina. This procedure was performed on four different days, on each day 12 ultrasonic signals were acquired, one signal every 10 minutes, with each apparatus. The results were compared and analyzed as regard the repeatability of the frequency spectra obtained.

Analysis of Ultrasonic Techniques for the Characterization of Microfiltration Polymeric Membranes

The use of polymeric membranes is extremely important in several industries such as nuclear, biotechnology, chemical and pharmaceutical. In the nuclear area, for instance, systems based on membrane separation technologies are currently being used in the treatment of radioactive liquid effluent, and new technologies using membranes are being developed at a great rate. The knowledge of the physical characteristics of these membranes, such as, pore size and the pore size distribution, is very important to the membranes separation processes. Only after these characteristics are known is it possible to determine the type and to choose a particular membrane for a specific application. In this work, two ultrasonic non destructive techniques were used to determine the porosity of membranes: pulse echo and transmission. A 25 MHz immersion transducer was used. Ultrasonic signals were acquired, for both techniques, after the ultrasonic waves passed through a microfiltration polymeric membrane of pore size of 0.45 μm and thickness of 180 μm. After the emitted ultrasonic signal crossed the membrane, the received signal brought several information on the influence of the membrane porosity in the standard signal of the ultrasonic wave. The ultrasonic signals were acquired in the time domain and changed to the frequency domain by application of the Fourier Fast Transform (FFT), thus generating the material frequency spectrum. For the pulse echo technique, the ultrasonic spectrum frequency changed after the ultrasonic wave crossed the membrane. With the transmission technique there was only a displacement of the ultrasonic signal at the time domain.

Flexible strategies for association analysis with genomic phenotypes

Thesis (Ph.D.)--University of Washington, 2016-08

On independent component analysis applied to unmixing hyperspectral data

One of the most challenging task underlying many hyperspectral imagery applications is the spectral unmixing, which decomposes a mixed pixel into a collection of reectance spectra, called endmember signatures, and their corresponding fractional abundances. Independent Component Analysis (ICA) have recently been proposed as a tool to unmix hyperspectral data. The basic goal of ICA is to nd a linear transformation to recover independent sources (abundance fractions) given only sensor observations that are unknown linear mixtures of the unobserved independent sources. In hyperspectral imagery the sum of abundance fractions associated to each pixel is constant due to physical constraints in the data acquisition process. Thus, sources cannot be independent. This paper address hyperspectral data source dependence and its impact on ICA performance. The study consider simulated and real data. In simulated scenarios hyperspectral observations are described by a generative model that takes into account the degradation mechanisms normally found in hyperspectral applications. We conclude that ICA does not unmix correctly all sources. This conclusion is based on the a study of the mutual information. Nevertheless, some sources might be well separated mainly if the number of sources is large and the signal-to-noise ratio (SNR) is high.