An introduction to low-level analysis methods of DNA microarray data

This article gives an overview over the methods used in the low--level analysis of gene expression data generated using DNA microarrays. This type of experiment allows to determine relative levels of nucleic acid abundance in a set of tissues or cell populations for thousands of transcripts or loci simultaneously. Careful statistical design and analysis are essential to improve the efficiency and reliability of microarray experiments throughout the data acquisition and analysis process. This includes the design of probes, the experimental design, the image analysis of microarray scanned images, the normalization of fluorescence intensities, the assessment of the quality of microarray data and incorporation of quality information in subsequent analyses, the combination of information across arrays and across sets of experiments, the discovery and recognition of patterns in expression at the single gene and multiple gene levels, and the assessment of significance of these findings, considering the fact that there is a lot of noise and thus random features in the data. For all of these components, access to a flexible and efficient statistical computing environment is an essential aspect.

Electron energy-loss spectroscopy as a tool for elemental analysis in biological speciments.

A transmission electron microscope (TEM) accessory, the energy filter, enables the establishment of a method for elemental microanalysis, the electron energy-loss spectroscopy (EELS). In conventional TEM, unscattered, elastic, and inelastic scattered electrons contribute to image information. Energy-filtering TEM (EFTEM) allows elemental analysis at the ultrastructural level by using selected inelastic scattered electrons. EELS is an excellent method for elemental microanalysis and nanoanalysis with good sensitivity and accuracy. However, it is a complex method whose potential is seldom completely exploited, especially for biological specimens. In addition to spectral analysis, parallel-EELS, we present two different imaging techniques in this chapter, namely electron spectroscopic imaging (ESI) and image-EELS. We aim to introduce these techniques in this chapter with the elemental microanalysis of titanium. Ultrafine, 22-nm titanium dioxide particles are used in an inhalation study in rats to investigate the distribution of nanoparticles in lung tissue.

Quantitative analysis of OCT characteristics in patients with achromatopsia and blue-cone monochromatism

PURPOSE: To quantify optical coherence tomography (OCT) images of the central retina in patients with blue-cone monochromatism (BCM) and achromatopsia (ACH) compared with healthy control individuals. METHODS: The study included 15 patients with ACH, 6 with BCM, and 20 control subjects. Diagnosis of BCM and ACH was established by visual acuity testing, morphologic examination, color vision testing, and Ganzfeld ERG recording. OCT images were acquired with the Stratus OCT 3 (Carl Zeiss Meditec AG, Oberkochen, Germany). Foveal OCT images were analyzed by calculating longitudinal reflectivity profiles (LRPs) from scan lines. Profiles were analyzed quantitatively to determine foveal thickness and distances between reflectivity layers. RESULTS: Patients with ACH and BCM had a mean visual acuity of 20/200 and 20/60, respectively. Color vision testing results were characteristic of the diseases. The LRPs of control subjects yielded four peaks (P1-P4), presumably representing the RPE (P1), the ovoid region of the photoreceptors (P2), the external limiting membrane (ELM) (P3), and the internal limiting membrane (P4). In patients with ACH, P2 was absent, but foveal thickness (P1-P4) did not differ significantly from that in the control subjects (187 +/- 20 vs. 192 +/- 14 microm, respectively). The distance from P1 to P3 did not differ significantly (78 +/- 10 vs. 82 +/- 5 microm) between ACH and controls subjects. In patients with BCM, P3 was lacking, and P2 advanced toward P1 compared with the control subjects (32 +/- 6 vs. 48 +/- 4 microm). Foveal thickness (153 +/- 16 microm) was significantly reduced compared with that in control subjects and patients with ACH. CONCLUSIONS: Quantitative OCT image analysis reveals distinct patterns for controls subjects and patients with ACH and BCM, respectively. Quantitative analysis of OCT imaging can be useful in differentiating retinal diseases affecting photoreceptors. Foveal thickness is similar in both normal subjects and patients with ACH but is decreased in patients with BCM.

NONLINEAR TUBE-FITTING FOR THE ANALYSIS OF ANATOMICAL AND FUNCTIONAL STRUCTURES

We are concerned with the estimation of the exterior surface of tube-shaped anatomical structures. This interest is motivated by two distinct scientific goals, one dealing with the distribution of HIV microbicide in the colon and the other with measuring degradation in white-matter tracts in the brain. Our problem is posed as the estimation of the support of a distribution in three dimensions from a sample from that distribution, possibly measured with error. We propose a novel tube-fitting algorithm to construct such estimators. Further, we conduct a simulation study to aid in the choice of a key parameter of the algorithm, and we test our algorithm with validation study tailored to the motivating data sets. Finally, we apply the tube-fitting algorithm to a colon image produced by single photon emission computed tomography (SPECT)and to a white-matter tract image produced using diffusion tensor `imaging (DTI).

Cardiac image fusion from stand-alone SPECT and CT: clinical experience

Myocardial perfusion imaging with SPECT (SPECT-MPI) and 64-slice CT angiography (CTA) are both established techniques for the noninvasive evaluation of coronary artery disease (CAD). Three-dimensional (3D) SPECT/CT image fusion may offer an incremental diagnostic value by integrating both sets of information. We report our first clinical experiences with fused 3D SPECT/CT in CAD patients. METHODS: Thirty-eight consecutive patients with at least 1 perfusion defect on SPECT-MPI (1-d adenosine stress/rest SPECT with (99m)Tc-tetrofosmin) and 64-slice CTA were included. 3D volume-rendered fused SPECT/CT images were generated and compared with the findings from the side-by-side analysis with regard to coronary lesion interpretation by assigning the perfusion defects to their corresponding coronary lesion. RESULTS: The fused SPECT/CT images added information on pathophysiologic lesion severity in 27 coronary stenoses (22%) of 12 patients (29%) (P<0.001). Among 40 equivocal lesions on side-by-side analysis, the fused interpretation confirmed hemodynamic significance in 14 lesions and excluded functional relevance in 10 lesions. In 3 lesions, assignment of perfusion defect and coronary lesion appeared to be reliable on side-by-side analysis but proved to be inaccurate on fused interpretation. Added diagnostic information by SPECT/CT was more commonly found in patients with stenoses of small vessels (P=0.004) and involvement of diagonal branches (P=0.01). CONCLUSION: In addition to being intuitively convincing, 3D SPECT/CT fusion images in CAD may provide added diagnostic information on the functional relevance of coronary artery lesions.

Abdominal multislice CT for obese patients: effect on image quality and radiation dose in a phantom study

RATIONALE AND OBJECTIVES: To evaluate the effect of a modified abdominal multislice computed tomography (CT) protocol for obese patients on image quality and radiation dose. MATERIALS AND METHODS: An adult female anthropomorphic phantom was used to simulate obese patients by adding one or two 4-cm circumferential layers of fat-equivalent material to the abdominal portion. The phantom was scanned with a subcutaneous fat thickness of 0, 4, and 8 cm using the following parameters (detector configuration/beam pitch/table feed per rotation/gantry rotation time/kV/mA): standard protocol A: 16 x 0.625 mm/1.75/17.5 mm/0.5 seconds/140/380, and modified protocol B: 16 x 1.25 mm/1.375/27.5 mm/1.0 seconds/140/380. Radiation doses to six abdominal organs and the skin, image noise values, and contrast-to-noise ratios (CNRs) were analyzed. Statistical analysis included analysis of variance, Wilcoxon rank sum, and Student's t-test (P < .05). RESULTS: Applying the modified protocol B with one or two fat rings, the image noise decreased significantly (P < .05), and simultaneously, the CNR increased significantly compared with protocol A (P < .05). Organ doses significantly increased, up to 54.7%, comparing modified protocol B with one fat ring to the routine protocol A with no fat rings (P < .05). However, no significant change in organ dose was seen for protocol B with two fat rings compared with protocol A without fat rings (range -2.1% to 8.1%) (P > .05). CONCLUSIONS: Using a modified abdominal multislice CT protocol for obese patients with 8 cm or more of subcutaneous fat, image quality can be substantially improved without a significant increase in radiation dose to the abdominal organs.

Comparison of magnetic resonance imaging of inhaled SF6 with respiratory gas analysis

Magnetic resonance imaging of inhaled fluorinated inert gases ((19)F-MRI) such as sulfur hexafluoride (SF(6)) allows for analysis of ventilated air spaces. In this study, the possibility of using this technique to image lung function was assessed. For this, (19)F-MRI of inhaled SF(6) was compared with respiratory gas analysis, which is a global but reliable measure of alveolar gas fraction. Five anesthetized pigs underwent multiple-breath wash-in procedures with a gas mixture of 70% SF(6) and 30% oxygen. Two-dimensional (19)F-MRI and end-expiratory gas fraction analysis were performed after 4 to 24 inhaled breaths. Signal intensity of (19)F-MRI and end-expiratory SF(6) fraction were evaluated with respect to linear correlation and reproducibility. Time constants were estimated by both MRI and respiratory gas analysis data and compared for agreement. A good linear correlation between signal intensity and end-expiratory gas fraction was found (correlation coefficient 0.99+/-0.01). The data were reproducible (standard error of signal intensity 8% vs. that of gas fraction 5%) and the comparison of time constants yielded a sufficient agreement. According to the good linear correlation and the acceptable reproducibility, we suggest the (19)F-MRI to be a valuable tool for quantification of intrapulmonary SF(6) and hence lung function.

Comprehensive review and application of particle image velocimetry

For a fluid dynamics experimental flow measurement technique, particle image velocimetry (PIV) provides significant advantages over other measurement techniques in its field. In contrast to temperature and pressure based probe measurements or other laser diagnostic techniques including laser Doppler velocimetry (LDV) and phase Doppler particle analysis (PDPA), PIV is unique due to its whole field measurement capability, non-intrusive nature, and ability to collect a vast amount of experimental data in a short time frame providing both quantitative and qualitative insight. These properties make PIV a desirable measurement technique for studies encompassing a broad range of fluid dynamics applications. However, as an optical measurement technique, PIV also requires a substantial technical understanding and application experience to acquire consistent, reliable results. Both a technical understanding of particle image velocimetry and practical application experience are gained by applying a planar PIV system at Michigan Technological University’s Combustion Science Exploration Laboratory (CSEL) and Alternative Fuels Combustion Laboratory (AFCL). Here a PIV system was applied to non-reacting and reacting gaseous environments to make two component planar PIV as well as three component stereographic PIV flow field velocity measurements in conjunction with chemiluminescence imaging in the case of reacting flows. This thesis outlines near surface flow field characteristics in a tumble strip lined channel, three component velocity profiles of non-reacting and reacting swirled flow in a swirl stabilized lean condition premixed/prevaporized-fuel model gas turbine combustor operating on methane at 5-7 kW, and two component planar PIV measurements characterizing the AFCL’s 1.1 liter closed combustion chamber under dual fan driven turbulent mixing flow.

ANALYSIS AND VISUALIZATION OF FLOW FIELDS USING INFORMATION-THEORETIC TECHNIQUES AND GRAPH-BASED REPRESENTATIONS

Three-dimensional flow visualization plays an essential role in many areas of science and engineering, such as aero- and hydro-dynamical systems which dominate various physical and natural phenomena. For popular methods such as the streamline visualization to be effective, they should capture the underlying flow features while facilitating user observation and understanding of the flow field in a clear manner. My research mainly focuses on the analysis and visualization of flow fields using various techniques, e.g. information-theoretic techniques and graph-based representations. Since the streamline visualization is a popular technique in flow field visualization, how to select good streamlines to capture flow patterns and how to pick good viewpoints to observe flow fields become critical. We treat streamline selection and viewpoint selection as symmetric problems and solve them simultaneously using the dual information channel [81]. To the best of my knowledge, this is the first attempt in flow visualization to combine these two selection problems in a unified approach. This work selects streamline in a view-independent manner and the selected streamlines will not change for all viewpoints. My another work [56] uses an information-theoretic approach to evaluate the importance of each streamline under various sample viewpoints and presents a solution for view-dependent streamline selection that guarantees coherent streamline update when the view changes gradually. When projecting 3D streamlines to 2D images for viewing, occlusion and clutter become inevitable. To address this challenge, we design FlowGraph [57, 58], a novel compound graph representation that organizes field line clusters and spatiotemporal regions hierarchically for occlusion-free and controllable visual exploration. We enable observation and exploration of the relationships among field line clusters, spatiotemporal regions and their interconnection in the transformed space. Most viewpoint selection methods only consider the external viewpoints outside of the flow field. This will not convey a clear observation when the flow field is clutter on the boundary side. Therefore, we propose a new way to explore flow fields by selecting several internal viewpoints around the flow features inside of the flow field and then generating a B-Spline curve path traversing these viewpoints to provide users with closeup views of the flow field for detailed observation of hidden or occluded internal flow features [54]. This work is also extended to deal with unsteady flow fields. Besides flow field visualization, some other topics relevant to visualization also attract my attention. In iGraph [31], we leverage a distributed system along with a tiled display wall to provide users with high-resolution visual analytics of big image and text collections in real time. Developing pedagogical visualization tools forms my other research focus. Since most cryptography algorithms use sophisticated mathematics, it is difficult for beginners to understand both what the algorithm does and how the algorithm does that. Therefore, we develop a set of visualization tools to provide users with an intuitive way to learn and understand these algorithms.

BEYOND ROOTS ALONE: NOVEL METHODOLOGIES FOR ANALYZING COMPLEX SOIL AND MINIRHIZOTRON IMAGERY USING IMAGE PROCESSING AND GIS TOOLS

Quantifying belowground dynamics is critical to our understanding of plant and ecosystem function and belowground carbon cycling, yet currently available tools for complex belowground image analyses are insufficient. We introduce novel techniques combining digital image processing tools and geographic information systems (GIS) analysis to permit semi-automated analysis of complex root and soil dynamics. We illustrate methodologies with imagery from microcosms, minirhizotrons, and a rhizotron, in upland and peatland soils. We provide guidelines for correct image capture, a method that automatically stitches together numerous minirhizotron images into one seamless image, and image analysis using image segmentation and classification in SPRING or change analysis in ArcMap. These methods facilitate spatial and temporal root and soil interaction studies, providing a framework to expand a more comprehensive understanding of belowground dynamics.

Antiplatelets versus Anticoagulants for the Treatment of Cervical Artery Dissection: Bayesian Meta-Analysis

OBJECTIVE To compare the effects of antiplatelets and anticoagulants on stroke and death in patients with acute cervical artery dissection. DESIGN Systematic review with Bayesian meta-analysis. DATA SOURCES The reviewers searched MEDLINE and EMBASE from inception to November 2012, checked reference lists, and contacted authors. STUDY SELECTION Studies were eligible if they were randomised, quasi-randomised or observational comparisons of antiplatelets and anticoagulants in patients with cervical artery dissection. DATA EXTRACTION Data were extracted by one reviewer and checked by another. Bayesian techniques were used to appropriately account for studies with scarce event data and imbalances in the size of comparison groups. DATA SYNTHESIS Thirty-seven studies (1991 patients) were included. We found no randomised trial. The primary analysis revealed a large treatment effect in favour of antiplatelets for preventing the primary composite outcome of ischaemic stroke, intracranial haemorrhage or death within the first 3 months after treatment initiation (relative risk 0.32, 95% credibility interval 0.12 to 0.63), while the degree of between-study heterogeneity was moderate (τ(2) = 0.18). In an analysis restricted to studies of higher methodological quality, the possible advantage of antiplatelets over anticoagulants was less obvious than in the main analysis (relative risk 0.73, 95% credibility interval 0.17 to 2.30). CONCLUSION In view of these results and the safety advantages, easier usage and lower cost of antiplatelets, we conclude that antiplatelets should be given precedence over anticoagulants as a first line treatment in patients with cervical artery dissection unless results of an adequately powered randomised trial suggest the opposite.

Estimating Ribeye Area Using a Longitudinal Ultrasound Image of the 12th and 13th Rib Section

A study was performed to determine the feasability of using a measurement of ribeye depth (RED) from a longitudinal ultrasound image to estimate ribeye area (REA). The correlation between RED obtained with ultrasound and REA from a tracing was high for both implanted (r = .49) and non-implanted (r = .45) steers. The mean bias between predicted REA and actual REA was not different from zero. This analysis shows that RED could be an accurate indicator of REA.

Finite element analysis for prediction of bone strength

Article preview View full access options BoneKEy Reports | Review Print Email Share/bookmark Finite element analysis for prediction of bone strength Philippe K Zysset, Enrico Dall'Ara, Peter Varga & Dieter H Pahr Affiliations Corresponding author BoneKEy Reports (2013) 2, Article number: 386 (2013) doi:10.1038/bonekey.2013.120 Received 03 January 2013 Accepted 25 June 2013 Published online 07 August 2013 Article tools Citation Reprints Rights & permissions Abstract Abstract• References• Author information Finite element (FE) analysis has been applied for the past 40 years to simulate the mechanical behavior of bone. Although several validation studies have been performed on specific anatomical sites and load cases, this study aims to review the predictability of human bone strength at the three major osteoporotic fracture sites quantified in recently completed in vitro studies at our former institute. Specifically, the performance of FE analysis based on clinical computer tomography (QCT) is compared with the ones of the current densitometric standards, bone mineral content, bone mineral density (BMD) and areal BMD (aBMD). Clinical fractures were produced in monotonic axial compression of the distal radii, vertebral sections and in side loading of the proximal femora. QCT-based FE models of the three bones were developed to simulate as closely as possible the boundary conditions of each experiment. For all sites, the FE methodology exhibited the lowest errors and the highest correlations in predicting the experimental bone strength. Likely due to the improved CT image resolution, the quality of the FE prediction in the peripheral skeleton using high-resolution peripheral CT was superior to that in the axial skeleton with whole-body QCT. Because of its projective and scalar nature, the performance of aBMD in predicting bone strength depended on loading mode and was significantly inferior to FE in axial compression of radial or vertebral sections but not significantly inferior to FE in side loading of the femur. Considering the cumulated evidence from the published validation studies, it is concluded that FE models provide the most reliable surrogates of bone strength at any of the three fracture sites.

In vitro accuracy evaluation of image-guided robot system for direct cochlear access

HYPOTHESIS A previously developed image-guided robot system can safely drill a tunnel from the lateral mastoid surface, through the facial recess, to the middle ear, as a viable alternative to conventional mastoidectomy for cochlear electrode insertion. BACKGROUND Direct cochlear access (DCA) provides a minimally invasive tunnel from the lateral surface of the mastoid through the facial recess to the middle ear for cochlear electrode insertion. A safe and effective tunnel drilled through the narrow facial recess requires a highly accurate image-guided surgical system. Previous attempts have relied on patient-specific templates and robotic systems to guide drilling tools. In this study, we report on improvements made to an image-guided surgical robot system developed specifically for this purpose and the resulting accuracy achieved in vitro. MATERIALS AND METHODS The proposed image-guided robotic DCA procedure was carried out bilaterally on 4 whole head cadaver specimens. Specimens were implanted with titanium fiducial markers and imaged with cone-beam CT. A preoperative plan was created using a custom software package wherein relevant anatomical structures of the facial recess were segmented, and a drill trajectory targeting the round window was defined. Patient-to-image registration was performed with the custom robot system to reference the preoperative plan, and the DCA tunnel was drilled in 3 stages with progressively longer drill bits. The position of the drilled tunnel was defined as a line fitted to a point cloud of the segmented tunnel using principle component analysis (PCA function in MatLab). The accuracy of the DCA was then assessed by coregistering preoperative and postoperative image data and measuring the deviation of the drilled tunnel from the plan. The final step of electrode insertion was also performed through the DCA tunnel after manual removal of the promontory through the external auditory canal. RESULTS Drilling error was defined as the lateral deviation of the tool in the plane perpendicular to the drill axis (excluding depth error). Errors of 0.08 ± 0.05 mm and 0.15 ± 0.08 mm were measured on the lateral mastoid surface and at the target on the round window, respectively (n =8). Full electrode insertion was possible for 7 cases. In 1 case, the electrode was partially inserted with 1 contact pair external to the cochlea. CONCLUSION The purpose-built robot system was able to perform a safe and reliable DCA for cochlear implantation. The workflow implemented in this study mimics the envisioned clinical procedure showing the feasibility of future clinical implementation.