Measurement of the vascular input function in mice for DCE-MRI

DCE-MRI is an important technique in the study of small animal cancer models because its sensitivity to vascular changes opens the possibility of quantitative assessment of early therapeutic response. However, extraction of physiologically descriptive parameters from DCE-MRI data relies upon measurement of the vascular input function (VIF), which represents the contrast agent concentration time course in the blood plasma. This is difficult in small animal models due to artifacts associated with partial volume, inflow enhancement, and the limited temporal resolution achievable with MR imaging. In this work, the development of a suite of techniques for high temporal resolution, artifact resistant measurement of the VIF in mice is described. One obstacle in VIF measurement is inflow enhancement, which decreases the sensitivity of the MR signal to the presence of contrast agent. Because the traditional techniques used to suppress inflow enhancement degrade the achievable spatiotemporal resolution of the pulse sequence, improvements can be achieved by reducing the time required for the suppression. Thus, a novel RF pulse which provides spatial presaturation contemporaneously with the RF excitation was implemented and evaluated. This maximizes the achievable temporal resolution by removing the additional RF and gradient pulses typically required for suppression of inflow enhancement. A second challenge is achieving the temporal resolution required for accurate characterization of the VIF, which exceeds what can be achieved with conventional imaging techniques while maintaining adequate spatial resolution and tumor coverage. Thus, an anatomically constrained reconstruction strategy was developed that allows for sampling of the VIF at extremely high acceleration factors, permitting capture of the initial pass of the contrast agent in mice. Simulation, phantom, and in vivo validation of all components were performed. Finally, the two components were used to perform VIF measurement in the murine heart. An in vivo study of the VIF reproducibility was performed, and an improvement in the measured injection-to-injection variation was observed. This will lead to improvements in the reliability of quantitative DCE-MRI measurements and increase their sensitivity.

High-sensitivity rod photoreceptor input to the blue-yellow color opponent pathway in macaque retina.

Small bistratified cells (SBCs) in the primate retina carry a major blue-yellow opponent signal to the brain. We found that SBCs also carry signals from rod photoreceptors, with the same sign as S cone input. SBCs exhibited robust responses under low scotopic conditions. Physiological and anatomical experiments indicated that this rod input arose from the AII amacrine cell-mediated rod pathway. Rod and cone signals were both present in SBCs at mesopic light levels. These findings have three implications. First, more retinal circuits may multiplex rod and cone signals than were previously thought to, efficiently exploiting the limited number of optic nerve fibers. Second, signals from AII amacrine cells may diverge to most or all of the approximately 20 retinal ganglion cell types in the peripheral primate retina. Third, rod input to SBCs may be the substrate for behavioral biases toward perception of blue at mesopic light levels.

Histamine reduces flash sensitivity of on ganglion cells in the primate retina.

PURPOSE. In Old World primates, the retina receives input from histaminergic neurons in the posterior hypothalamus. They are a subset of the neurons that project throughout the central nervous system and fire maximally during the day. The contribution of these neurons to vision, was examined by applying histamine to a dark-adapted, superfused baboon eye cup preparation while making extracellular recordings from peripheral retinal ganglion cells. METHODS. The stimuli were 5-ms, 560-nm, weak, full-field flashes in the low scotopic range. Ganglion cells with sustained and transient ON responses and two cell types with OFF responses were distinguished; their responses were recorded with a 16-channel microelectrode array. RESULTS. Low micromolar doses of histamine decreased the rate of maintained firing and the light sensitivity of ON ganglion cells. Both sustained and transient ON cells responded similarly to histamine. There were no statistically significant effects of histamine in a more limited study of OFF ganglion cells. The response latencies of ON cells were approximately 5 ms slower, on average, when histamine was present. Histamine also reduced the signal-to-noise ratio of ON cells, particularly in those cells with a histamine-induced increase in maintained activity. CONCLUSIONS. A major action of histamine released from retinopetal axons under dark-adapted conditions, when rod signals dominate the response, is to reduce the sensitivity of ON ganglion cells to light flashes. These findings may relate to reports that humans are less sensitive to light stimuli in the scotopic range during the day, when histamine release in the retina is expected to be at its maximum.

CHOLINERGIC MODULATION OF THALAMIC INPUT INTO THE CINGULATE CORTEX

This research demonstrates cholinergic modulation of thalamic input into the limbic cortex. A projection from the mediodorsal thalamus (MD) to the anterior cingulate cortex was defined anatomically and physiologically. Injections of horse-radish peroxidase into the anterior cingulate cortex labels neurons in the lateral, parvocellular, region of MD. Electrical Stimulation of this area produces a complex field potential in the anterior cingulate cortex which was further characterized by current density analysis and single cell recordings.^ The monsynaptic component of the response was identified as a large negative field which is maximal in layer IV of the anterior cingulate cortex. This response shows remarkable tetanic potentiation of frequencies near 7 Hz. During a train of 50 or more stimuli, the response would grow quickly and remain at a fairly stable potentiated level throughout the train.^ Cholinergic modulation of this thalamic response was demonstrated by iontophoretic application of the cholinergic agonist carbachol decreased the effectiveness of the thalamic imput by rapidly attenuation the response during a train of stimuli. The effect was apparently mediated by muscarinic receptors since the effect of carbachol was blocked by atropine but not by hexamethonium.^ To determine the source of the cingulate cortex cholinergic innervation, lesions were made in the anterior and medial thalamus and in the nucleus of the diagonal band of Broca. The effects of these lesions on choline acetyltranferase activity in the cingulate cortex were determined by a micro-radio-enzymatical assay. Only the lesions of the nucleus of the diagonal band significantly decreased the choline acetyltransferase activity in the cingulate cortex regions. Therefore, the diagonal band appears to be a major source of sensory cholinergic innervation and may be involved in gating of sensory information from the thalamus into the limbic cortex. Attempts to modulate the cingulate response to MD stimulation with electrical stimulation of the diagonal band, however were not successful.^

Diversity of neuronal connexins in the mammalian retina

Many neurons in the mammalian retina are electrically coupled by intercellular channels or gap junctions, which are assembled from a family of proteins called connexins. Numerous studies indicate that gap junctions differ in properties such as conductance and tracer permeability. For example, A-type horizontal cell gap junctions are permeable to Lucifer Yellow, but B-type horizontal cell gap junctions are not. This suggests the two cell types express different connexins. My hypothesis is that multiple neuronal connexins are expressed in the mammalian retina in a cell type specific manner. Immunohistochemical techniques and confocal microscopy were used to localize certain connexins within well-defined neuronal circuits. The results of this study can be summarized as follows: AII amacrine cells, which receive direct input from rod bipolar cells, are well-coupled to neighboring AIIs. In addition, AII amacrine cells also form gap junctions with ON cone bipolar cells. This is a complex heterocellular network. In both rabbit and primate retina, connexin36 occurs at dendritic crossings in the AII matrix as well as between AIIs and ON cone bipolar cells. Coupling in the AII network is thought to reduce noise in the rod pathway while AII/bipolar gap junctions are required for the transmission of rod signals to ON ganglion cells. In the outer plexiform layer, connexin36 forms gap junctions between cones and between rods and cones via cone telodendria. Cone to cone coupling is thought to reduce noise and is partly color selective. Rod to cone coupling forms an alternative rod pathway thought to operate at intermediate light intensity. A-type horizontal cells in the rabbit retina are strongly coupled via massive low resistance gap junctions composed from Cx50. Coupling dramatically extends the receptive field of horizontal cells and the modulation of coupling is thought to change the strength of the feedback signal from horizontal cells to cones. Finally, there are other coupled networks, such as B-type horizontal cells and S1/S2 amacrine cells, which do not use either connexin36 or Cx50. These results confirm the hypothesis that multiple neuronal connexins are expressed in the mammalian retina and these connexins are localized to particular retinal circuits. ^

An empirical evaluation of the Random Forests classifier models for variable selection in a large-scale lung cancer case-control study

Random Forests™ is reported to be one of the most accurate classification algorithms in complex data analysis. It shows excellent performance even when most predictors are noisy and the number of variables is much larger than the number of observations. In this thesis Random Forests was applied to a large-scale lung cancer case-control study. A novel way of automatically selecting prognostic factors was proposed. Also, synthetic positive control was used to validate Random Forests method. Throughout this study we showed that Random Forests can deal with large number of weak input variables without overfitting. It can account for non-additive interactions between these input variables. Random Forests can also be used for variable selection without being adversely affected by collinearities. ^ Random Forests can deal with the large-scale data sets without rigorous data preprocessing. It has robust variable importance ranking measure. Proposed is a novel variable selection method in context of Random Forests that uses the data noise level as the cut-off value to determine the subset of the important predictors. This new approach enhanced the ability of the Random Forests algorithm to automatically identify important predictors for complex data. The cut-off value can also be adjusted based on the results of the synthetic positive control experiments. ^ When the data set had high variables to observations ratio, Random Forests complemented the established logistic regression. This study suggested that Random Forests is recommended for such high dimensionality data. One can use Random Forests to select the important variables and then use logistic regression or Random Forests itself to estimate the effect size of the predictors and to classify new observations. ^ We also found that the mean decrease of accuracy is a more reliable variable ranking measurement than mean decrease of Gini. ^

Intensive care noise and mean arterial blood pressure in ELBW neonates

Objective. Loud noises in neonatal intensive care units (NICUs) may impede growth and development for extremely low birthweight (ELBW, < 1000 grams) newborns. The objective of this study was to measure the association between NICU sound levels and ELBW neonates' arterial blood pressure to determine whether these newborns experience noise-induced stress. ^ Methods. Noise and arterial blood pressure recordings were collected for 9 ELBW neonates during the first week of life. Sound levels were measured inside the incubator, and each subject's arterial blood pressures were simultaneously recorded for 15 minutes (at 1 sec intervals). Time series cross-correlation functions were calculated for NICU noise and mean arterial blood pressure (MABP) recordings for each subject. The grand mean noise-MABP cross-correlation was calculated for all subjects and for lower and higher birthweight groups for comparison. ^ Results. The grand mean noise-MABP cross-correlation for all subjects was mostly negative (through 300 sec lag time) and nearly reached significance at the 95% level at 111 sec lag (mean r = -0.062). Lower birthweight newborns (454-709 g) experienced significant decreases in blood pressure with increasing NICU noise after 145 sec lag (peak r = -0.074). Higher birthweight newborns had an immediate negative correlation with NICU sound levels (at 3 sec lag, r = -0.071), but arterial blood pressures increased to a positive correlation with noise levels at 197 sec lag (r = 0.075). ^ Conclusions. ELBW newborns' arterial blood pressure was influenced by NICU noise levels during the first week of life. Lower birthweight newborns may have experienced an orienting reflex to NICU sounds. Higher birthweight newborns experienced an immediate orienting reflex to increasing sound levels, but arterial blood pressure increased approximately 3 minutes after increases in noise levels. Increases in arterial blood pressure following increased NICU sound levels may result from a stress response to noise. ^

TIME SERIES ANALYSIS AS INPUT FOR PREDICTIVE MODELING: PREDICTING CARDIAC ARREST IN A PEDIATRIC INTENSIVE CARE UNIT

The first manuscript, entitled "Time-Series Analysis as Input for Clinical Predictive Modeling: Modeling Cardiac Arrest in a Pediatric ICU" lays out the theoretical background for the project. There are several core concepts presented in this paper. First, traditional multivariate models (where each variable is represented by only one value) provide single point-in-time snapshots of patient status: they are incapable of characterizing deterioration. Since deterioration is consistently identified as a precursor to cardiac arrests, we maintain that the traditional multivariate paradigm is insufficient for predicting arrests. We identify time series analysis as a method capable of characterizing deterioration in an objective, mathematical fashion, and describe how to build a general foundation for predictive modeling using time series analysis results as latent variables. Building a solid foundation for any given modeling task involves addressing a number of issues during the design phase. These include selecting the proper candidate features on which to base the model, and selecting the most appropriate tool to measure them. We also identified several unique design issues that are introduced when time series data elements are added to the set of candidate features. One such issue is in defining the duration and resolution of time series elements required to sufficiently characterize the time series phenomena being considered as candidate features for the predictive model. Once the duration and resolution are established, there must also be explicit mathematical or statistical operations that produce the time series analysis result to be used as a latent candidate feature. In synthesizing the comprehensive framework for building a predictive model based on time series data elements, we identified at least four classes of data that can be used in the model design. The first two classes are shared with traditional multivariate models: multivariate data and clinical latent features. Multivariate data is represented by the standard one value per variable paradigm and is widely employed in a host of clinical models and tools. These are often represented by a number present in a given cell of a table. Clinical latent features derived, rather than directly measured, data elements that more accurately represent a particular clinical phenomenon than any of the directly measured data elements in isolation. The second two classes are unique to the time series data elements. The first of these is the raw data elements. These are represented by multiple values per variable, and constitute the measured observations that are typically available to end users when they review time series data. These are often represented as dots on a graph. The final class of data results from performing time series analysis. This class of data represents the fundamental concept on which our hypothesis is based. The specific statistical or mathematical operations are up to the modeler to determine, but we generally recommend that a variety of analyses be performed in order to maximize the likelihood that a representation of the time series data elements is produced that is able to distinguish between two or more classes of outcomes. The second manuscript, entitled "Building Clinical Prediction Models Using Time Series Data: Modeling Cardiac Arrest in a Pediatric ICU" provides a detailed description, start to finish, of the methods required to prepare the data, build, and validate a predictive model that uses the time series data elements determined in the first paper. One of the fundamental tenets of the second paper is that manual implementations of time series based models are unfeasible due to the relatively large number of data elements and the complexity of preprocessing that must occur before data can be presented to the model. Each of the seventeen steps is analyzed from the perspective of how it may be automated, when necessary. We identify the general objectives and available strategies of each of the steps, and we present our rationale for choosing a specific strategy for each step in the case of predicting cardiac arrest in a pediatric intensive care unit. Another issue brought to light by the second paper is that the individual steps required to use time series data for predictive modeling are more numerous and more complex than those used for modeling with traditional multivariate data. Even after complexities attributable to the design phase (addressed in our first paper) have been accounted for, the management and manipulation of the time series elements (the preprocessing steps in particular) are issues that are not present in a traditional multivariate modeling paradigm. In our methods, we present the issues that arise from the time series data elements: defining a reference time; imputing and reducing time series data in order to conform to a predefined structure that was specified during the design phase; and normalizing variable families rather than individual variable instances. The final manuscript, entitled: "Using Time-Series Analysis to Predict Cardiac Arrest in a Pediatric Intensive Care Unit" presents the results that were obtained by applying the theoretical construct and its associated methods (detailed in the first two papers) to the case of cardiac arrest prediction in a pediatric intensive care unit. Our results showed that utilizing the trend analysis from the time series data elements reduced the number of classification errors by 73%. The area under the Receiver Operating Characteristic curve increased from a baseline of 87% to 98% by including the trend analysis. In addition to the performance measures, we were also able to demonstrate that adding raw time series data elements without their associated trend analyses improved classification accuracy as compared to the baseline multivariate model, but diminished classification accuracy as compared to when just the trend analysis features were added (ie, without adding the raw time series data elements). We believe this phenomenon was largely attributable to overfitting, which is known to increase as the ratio of candidate features to class examples rises. Furthermore, although we employed several feature reduction strategies to counteract the overfitting problem, they failed to improve the performance beyond that which was achieved by exclusion of the raw time series elements. Finally, our data demonstrated that pulse oximetry and systolic blood pressure readings tend to start diminishing about 10-20 minutes before an arrest, whereas heart rates tend to diminish rapidly less than 5 minutes before an arrest.