A Data-Analytic Strategy for Protein-Biomarker Discovery: Profiling of High-Dimensional Proteomic Data for Cancer Detection

With recent advances in mass spectrometry techniques, it is now possible to investigate proteins over a wide range of molecular weights in small biological specimens. This advance has generated data-analytic challenges in proteomics, similar to those created by microarray technologies in genetics, namely, discovery of "signature" protein profiles specific to each pathologic state (e.g., normal vs. cancer) or differential profiles between experimental conditions (e.g., treated by a drug of interest vs. untreated) from high-dimensional data. We propose a data analytic strategy for discovering protein biomarkers based on such high-dimensional mass-spectrometry data. A real biomarker-discovery project on prostate cancer is taken as a concrete example throughout the paper: the project aims to identify proteins in serum that distinguish cancer, benign hyperplasia, and normal states of prostate using the Surface Enhanced Laser Desorption/Ionization (SELDI) technology, a recently developed mass spectrometry technique. Our data analytic strategy takes properties of the SELDI mass-spectrometer into account: the SELDI output of a specimen contains about 48,000 (x, y) points where x is the protein mass divided by the number of charges introduced by ionization and y is the protein intensity of the corresponding mass per charge value, x, in that specimen. Given high coefficients of variation and other characteristics of protein intensity measures (y values), we reduce the measures of protein intensities to a set of binary variables that indicate peaks in the y-axis direction in the nearest neighborhoods of each mass per charge point in the x-axis direction. We then account for a shifting (measurement error) problem of the x-axis in SELDI output. After these pre-analysis processing of data, we combine the binary predictors to generate classification rules for cancer, benign hyperplasia, and normal states of prostate. Our approach is to apply the boosting algorithm to select binary predictors and construct a summary classifier. We empirically evaluate sensitivity and specificity of the resulting summary classifiers with a test dataset that is independent from the training dataset used to construct the summary classifiers. The proposed method performed nearly perfectly in distinguishing cancer and benign hyperplasia from normal. In the classification of cancer vs. benign hyperplasia, however, an appreciable proportion of the benign specimens were classified incorrectly as cancer. We discuss practical issues associated with our proposed approach to the analysis of SELDI output and its application in cancer biomarker discovery.

Methodical study of nitrous oxide eddy covariance measurements using quantum cascade laser spectrometery over a Swiss forest

Nitrous oxide fluxes were measured at the Lägeren CarboEurope IP flux site over the multi-species mixed forest dominated by European beech and Norway spruce. Measurements were carried out during a four-week period in October–November 2005 during leaf senescence. Fluxes were measured with a standard ultrasonic anemometer in combination with a quantum cascade laser absorption spectrometer that measured N2O, CO2, and H2O mixing ratios simultaneously at 5 Hz time resolution. To distinguish insignificant fluxes from significant ones it is proposed to use a new approach based on the significance of the correlation coefficient between vertical wind speed and mixing ratio fluctuations. This procedure eliminated roughly 56% of our half-hourly fluxes. Based on the remaining, quality checked N2O fluxes we quantified the mean efflux at 0.8±0.4 μmol m−2 h−1 (mean ± standard error). Most of the contribution to the N2O flux occurred during a 6.5-h period starting 4.5 h before each precipitation event. No relation with precipitation amount could be found. Visibility data representing fog density and duration at the site indicate that wetting of the canopy may have as strong an effect on N2O effluxes as does below-ground microbial activity. It is speculated that above-ground N2O production from the senescing leaves at high moisture (fog, drizzle, onset of precipitation event) may be responsible for part of the measured flux.

Comparing satellite and ground-based observations of paroxysmal degassing events at Etna volcano, Italy

Mount Etna, Italy, is one of the most active volcanoes in the world, and is also regarded as one of the strongest volcanic sources of sulfur dioxide (SO2) emissions to the atmosphere. Since October 2004, an automated ultraviolet (UV) spectrometer network (FLAME) has provided ground-based SO2 measurements with high temporal resolution, providing an opportunity to validate satellite SO2 measurements at Etna. The Ozone Monitoring Instrument (OMI) on the NASA Aura satellite, which makes global daily measurements of trace gases in the atmosphere, was used to compare SO2 amount released by the volcano during paroxysmal lava-fountaining events from 2004 to present. We present the first comparison between SO2 emission rates and SO2 burdens obtained by the OMI transect technique and OMI Normalized Cloud-Mass (NCM) technique and the ground-based FLAME Mini-DOAS measurements. In spite of a good data set from the FLAME network, finding coincident OMI and FLAME measurements proved challenging and only one paroxysmal event provided a good validation for OMI. Another goal of this work was to assess the efficacy of the FLAME network in capturing paroxysmal SO2 emissions from Etna, given that the FLAME network is only operational during daylight hours and some paroxysms occur at night. OMI measurements are advantageous since SO2 emissions from nighttime paroxysms can often be quantified on the following day, providing improved constraints on Etna’s SO2 budget.

Remote sensing of volcanic plumes using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER).

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) has been used to quantify SO2 emissions from passively degassing volcanoes. This dissertation explores ASTER’s capability to detect SO2 with satellite validation, enhancement techniques and extensive processing of images at a variety of volcanoes. ASTER is compared to the Mini UV Spectrometer (MUSe), a ground based instrument, to determine if reasonable SO2 fluxes can be quantified from a plume emitted from Lascar, Chile. The two sensors were in good agreement with ASTER proving to be a reliable detector of SO2. ASTER illustrated the advantages of imaging a plume in 2D, with better temporal resolution than the MUSe. SO2 plumes in ASTER imagery are not always discernible in the raw TIR data. Principal Component Analysis (PCA) and Decorrelation Stretch (DCS) enhancement techniques were compared to determine how well they highlight a variety of volcanic plumes. DCS produced a consistent output and the composition of the plumes was easy to identify from explosive eruptions. As the plumes became smaller and lower in altitude they became harder to distinguish using DCS. PCA proved to be better at identifying smaller low altitude plumes. ASTER was used to investigate SO2 emissions at Lascar, Chile. Activity at Lascar has been characterized by cyclic behavior and persistent degassing (Matthews et al. 1997). Previous studies at Lascar have primarily focused on changes in thermal infrared anomalies, neglecting gas emissions. Using the SO2 data along with changes in thermal anomalies and visual observations it is evident that Lascar is at the end an eruptive cycle that began in 1993. Declining gas emissions and crater temperatures suggest that the conduit is sealing. ASTER and the Ozone Monitoring Instrument (OMI) were used to determine the annual contribution of SO2 to the troposphere from the Central and South American volcanic arcs between 2000 and 2011. Fluxes of 3.4 Tg/a for Central America and 3.7 Tg/a for South America were calculated. The detection limits of ASTER were explored. The results a proved to be interesting, with plumes from many of the high emitting volcanoes, such as Villarrica, Chile, not being detected by ASTER.

SO2 emissions at Volcan de Colima, 2003-2007

Volcán de Colima has been continuously erupting since the onset of dome growth in 1998. This period of unrest has had 4 prominent periods; 1998-1999, 2003, 2004-2005, and the current dome growth that began in February of 2007. Each of these episodes was marked by lava extrusion forming a dome and lava flows, followed by explosions that destroyed the dome. The Correlation Spectrometer (COSPEC) was used to determine SO2 emission rates on 164 days from May 2003 to February 2007, using both stationary ground based scans and some flight traverses. Scans were separated into the categories of explosive degassing and passive, or background degassing. These scans show variation in the SO2 flow rate from below detection limit (~3 t/d depending on environmental conditions) during background, passive emissions to a peak of 2949 t/d (34 kilograms/second) during an explosion on 9 October, 2004. Both passive and explosive degassing increased when there was lava extrusion in 2004 and with the increased explosive activity in 2005. These two different processes of degassing wax with each other when activity increases and wane together as well, indicating a parallel cyclicity in the volcanic eruption and degassing rates, where the conduit partially seals (pressurizes) between explosions. Colima’s gas and eruptive behavior is compared to similar systems such as Santiaguito and Soufrière Hills, Montserrat. About 2/3 of Colima’s SO2 degassing, amounting to 1.3 x 105 tonnes in 3.74 yrs has come in short lived small (VEI=0-1) vertical explosions that occurred at the rate of 100-3000explosions/ month, and the remaining third has occured in continuous passive degassing. Colima emits sulfur at a rate equivalent to about 0.04 to 0.08 wt % S, similar to other andesitic convergent plate boundary volcanoes. There has been an explosive destruction of the dome in every cycle for that past 5 years, and it is assumed that the current dome which began growth in February, 2007 (just at the end of this study) will be destroyed. Higher emission rates seen in the quiescence of 2006 may have eased the pressure at the time, resulting in the slow effusion of the current dome and lack of explosivity.

Organic Wastewater Chemicals in Silver Bow Creek - Butte to Warm Springs Ponds

The purpose of this study is to detail and analyze the distribution, concentration, and loads of 5 organic compounds along Silver Bow Creek in Butte, Montana from the Municipal Wastewater treatment plant to the Warm Springs Ponds. The chemicals analyzed include Carbamazepine (pharmaceutical), Miconazole (fungicide) and three antibiotics – Sulfamethoxazole, Thiabendazole, and Ciprofloxacin. This project begins a 2 year study to analyze 6 additional compounds (11 compounds total), to develop an effective method to detail and analyze OWCs using Mass Spectrometer/Liquid chromatography system, and to aid in assessment of aquatic health and ongoing restoration work. The EPA method 1694 was used for analysis

The STEREO/PLASTIC response to solar wind ions (Flight measurements and models)

The Plasma and Supra-Thermal Ion Composition (PLASTIC) instrument is one of four experiment packages on board of the two identical STEREO spacecraft A and B, which were successfully launched from Cape Canaveral on 26 October 2006. During the two years of the nominal STEREO mission, PLASTIC is providing us with the plasma characteristics of protons, alpha particles, and heavy ions. PLASTIC will also provide key diagnostic measurements in the form of the mass and charge state composition of heavy ions. Three measurements (E/qk, time of flight, ESSD) from the pulse height raw data are used to characterize the solar wind ions from the solar wind sector, and part of the suprathermal particles from the wide-angle partition with respect to mass, atomic number and charge state. In this paper, we present a new method for flight data analysis based on simulations of the PLASTIC response to solar wind ions. We present the response of the entrance system / energy analyzer in an analytical form. Based on stopping power theory, we use an analytical expression for the energy loss of the ions when they pass through a thin carbon foil. This allows us to model analytically the response of the time of flight mass spectrometer to solar wind ions. Thus we present a new version of the analytical response of the solid state detectors to solar wind ions. Various important parameters needed for our models were derived, based on calibration data and on the first flight measurements obtained from STEREO-A. We used information from each measured event that is registered in full resolution in the Pulse Height Analysis words and we derived a new algorithm for the analysis of both existing and future data sets of a similar nature which was tested and works well. This algorithm allows us to obtain, for each measured event, the mass, atomic number and charge state in the correct physical units. Finally, an important criterion was developed for filtering our Fe raw flight data set from the pulse height data without discriminating charge states.

Validation of water vapour profiles (version 13) retrieved by the IMK / IAA scientific retrieval processor based on full resolution spectra measured by MIPAS on board Envisat

Vertical profiles of stratospheric water vapour measured by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) with the full resolution mode between September 2002 and March 2004 and retrieved with the IMK/IAA scientific retrieval processor were compared to a number of independent measurements in order to estimate the bias and to validate the existing precision estimates of the MIPAS data. The estimated precision for MIPAS is 5 to 10% in the stratosphere, depending on altitude, latitude, and season. The independent instruments were: the Halogen Occultation Experiment (HALOE), the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS), the Improved Limb Atmospheric Spectrometer-II (ILAS-II), the Polar Ozone and Aerosol Measurement (POAM III) instrument, the Middle Atmospheric Water Vapour Radiometer (MIAWARA), the Michelson Interferometer for Passive Atmospheric Sounding, balloon-borne version (MIPAS-B), the Airborne Microwave Stratospheric Observing System (AMSOS), the Fluorescent Stratospheric Hygrometer for Balloon (FLASH-B), the NOAA frostpoint hygrometer, and the Fast In Situ Hygrometer (FISH). For the in-situ measurements and the ground based, air- and balloon borne remote sensing instruments, the measurements are restricted to central and northern Europe. The comparisons to satellite-borne instruments are predominantly at mid- to high latitudes on both hemispheres. In the stratosphere there is no clear indication of a bias in MIPAS data, because the independent measurements in some cases are drier and in some cases are moister than the MIPAS measurements. Compared to the infrared measurements of MIPAS, measurements in the ultraviolet and visible have a tendency to be high, whereas microwave measurements have a tendency to be low. The results of χ2-based precision validation are somewhat controversial among the comparison estimates. However, for comparison instruments whose error budget also includes errors due to uncertainties in spectrally interfering species and where good coincidences were found, the χ2 values found are in the expected range or even below. This suggests that there is no evidence of systematically underestimated MIPAS random errors.

Light curing through glass ceramics with a second- and a third-generation LED curing unit: effect of curing mode on the degree of conversion of dual-curing resin cements

Objectives The aim of this study was to measure the degree of conversion (DC) of five dual-curing resin cements after different curing modes with a second- and a third-generation light-emitting diode (LED) curing unit. Additionally, irradiance of both light curing units was measured at increasing distances and through discs of two glass ceramics for computer-aided design/manufacturing (CAD/CAM). Materials and methods Irradiance and spectra of the Elipar FreeLight 2 (Standard Mode (SM)) and of the VALO light curing unit (High Power Mode (HPM) and Xtra Power Mode (XPM)) were measured with a MARC radiometer. Irradiance was measured at increasing distances (control) and through discs (1.5 to 6 mm thickness) of IPS Empress CAD and IPS e.max CAD. DC of Panavia F2.0, RelyX Unicem 2 Automix, SpeedCEM, BisCem, and BeautiCem SA was measured with an attenuated total reflectance–Fourier transform infrared spectrometer when self-cured (negative control) or light cured in SM for 40 s, HPM for 32 s, or XPM for 18 s. Light curing was performed directly (positive control) or through discs of either 1.5- or 3-mm thickness of IPS Empress CAD or IPS e.max CAD. DC was analysed with Kruskal–Wallis tests followed by pairwise Wilcoxon rank sum tests (α = 0.05). Results Maximum irradiances were 1,545 mW/cm2 (SM), 2,179 mW/cm2 (HPM), and 4,156 mW/cm2 (XPM), and all irradiances decreased by >80 % through discs of 1.5 mm, ≥95 % through 3 mm, and up to >99 % through 6 mm. Generally, self-curing resulted in the lowest DC. For some cements, direct light curing did not result in higher DC compared to when light cured through ceramic discs. For other cements, light curing through ceramic discs of 3 mm generally reduced DC. Conclusions Light curing was favourable for dual-curing cements. Some cements were more susceptible to variations in curing mode than others. Clinical relevance When light curing a given cement, the higher irradiances of the third-generation LED curing unit resulted in similar DC compared to the second-generation one, though at shorter light curing times.

Prediction of soil organic carbon for Ethiopian highlands using soil spectroscopy

Soil spectroscopy was applied for predicting soil organic carbon (SOC) in the highlands of Ethiopia. Soil samples were acquired from Ethiopia’s National Soil Testing Centre and direct field sampling. The reflectance of samples was measured using a FieldSpec 3 diffuse reflectance spectrometer. Outliers and sample relation were evaluated using principal component analysis (PCA) and models were developed through partial least square regression (PLSR). For nine watersheds sampled, 20% of the samples were set aside to test prediction and 80% were used to develop calibration models. Depending on the number of samples per watershed, cross validation or independent validation were used.The stability of models was evaluated using coefficient of determination (R2), root mean square error (RMSE), and the ratio performance deviation (RPD). The R2 (%), RMSE (%), and RPD, respectively, for validation were Anjeni (88, 0.44, 3.05), Bale (86, 0.52, 2.7), Basketo (89, 0.57, 3.0), Benishangul (91, 0.30, 3.4), Kersa (82, 0.44, 2.4), Kola tembien (75, 0.44, 1.9),Maybar (84. 0.57, 2.5),Megech (85, 0.15, 2.6), andWondoGenet (86, 0.52, 2.7) indicating that themodels were stable. Models performed better for areas with high SOC values than areas with lower SOC values. Overall, soil spectroscopy performance ranged from very good to good.

ROSINA/DFMS capabilities to measure isotopic ratios in water at comet 67P/Churyumov–Gerasimenko

The likelihood that comets may have delivered part of the water to Earth has been reinforced by the recent observation of the earth-like D/H ratio in Jupiter-family comet 103P/Hartley 2 by Hartogh et al. (2011). Prior to this observation, results from several Oort cloud comets indicated a factor of 2 enrichment of deuterium relative to the abundance at Earth. The European Space Agency’s Rosetta spacecraft will encounter comet 67P/Churyumov-Gerasimenko, another Jupiter-family comet of likely Kuiper belt origin, in 2014 and accompany it from almost aphelion to and past perihelion. Onboard Rosetta is the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) which consists of two mass spectrometers and a pressure sensor [Balsiger et al. 2007]. With its unprecedented mass resolution, for a space-borne instrument, the Double Focusing Mass Spectrometer (DFMS), one of the major subsystems of ROSINA, will be able to obtain unambiguously the ratios of the isotopes in water from in situ measurements in the coma around the comet. We will discuss the performance of this sensor on the basis of measurements of the terrestrial hydrogen and oxygen isotopic ratios performed with the flight spare instrument in the lab. We also show that the instrument on Rosetta is capable of measuring the D/H even in the very low density water background released by the spacecraft. This capability demonstrates that ROSINA should obtain very sensitive measurements of these ratios in the cometary environment. These measurements will allow detection of fractionation as function of the distance from the nucleus as well as fractionation due to mechanisms that are correlated with heliocentric distance.

Remote energetic neutral atom imaging of electric potential over a lunar magnetic anomaly

The formation of electric potential over lunar magnetized regions is essential for understanding fundamental lunar science, for understanding the lunar environment, and for planning human exploration on the Moon. A large positive electric potential was predicted and detected from single point measurements. Here, we demonstrate a remote imaging technique of electric potential mapping at the lunar surface, making use of a new concept involving hydrogen neutral atoms derived from solar wind. We apply the technique to a lunar magnetized region using an existing dataset of the neutral atom energy spectrometer SARA/CENA on Chandrayaan-1. Electrostatic potential larger than +135 V inside the Gerasimovic anomaly is confirmed. This structure is found spreading all over the magnetized region. The widely spread electric potential can influence the local plasma and dust environment near the magnetic anomaly.

A Texture-based land cover classification for the delineation of a shifting cultivation landscape in the Lao PDR using landscape metrics

The delineation of shifting cultivation landscapes using remote sensing in mountainous regions is challenging. On the one hand, there are difficulties related to the distinction of forest and fallow forest classes as occurring in a shifting cultivation landscape in mountainous regions. On the other hand, the dynamic nature of the shifting cultivation system poses problems to the delineation of landscapes where shifting cultivation occurs. We present a two-step approach based on an object-oriented classification of Advanced Land Observing Satellite, Advanced Visible and Near-Infrared Spectrometer (ALOS AVNIR) and Panchromatic Remote-sensing Instrument for Stereo Mapping (ALOS PRISM) data and landscape metrics. When including texture measures in the object-oriented classification, the accuracy of forest and fallow forest classes could be increased substantially. Based on such a classification, landscape metrics in the form of land cover class ratios enabled the identification of crop-fallow rotation characteristics of the shifting cultivation land use practice. By classifying and combining these landscape metrics, shifting cultivation landscapes could be delineated using a single land cover dataset.

High precision determination of the terrestrial 40K abundance

Recent improvements in the precision of mass spectrometric measurements have reduced the uncertainty of K-Ar and 39Ar-40Ar ages measured on geological materials. Now the major sources of uncertainty are the uncertainties on the 40K decay constant and the absolute abundance of 40K. In order to improve on this situation we determined the abundance of the 40K isotope in terrestrial standards. A ThermoFischer Triton+ thermal ionization mass spectrometer was used for K isotope ratio measurements of the NIST K standard reference materials SRM 918b and SRM 985. Ion beams were measured in Faraday cups with amplifiers equipped with 1E10, 1E11 and 1E12 Ω resistors. Three measurement protocols were used: (A) dynamic measurement with in-run fractionation correction by normalization to the IUPAC recommended isotope ratio 41K/39K = 0.0721677; (B) total evaporation; (C) a modified total evaporation with interblock baseline measurements. Different measurement protocols were combined with different loading procedures. The best results were obtained by loading samples on single tantalum filaments with 0.1M H3PO4. The total ion yields (ionization + transmission) were tested for the evaporation procedures (B) and (C) and ranged up to 48 %. The resulting best estimate for the 40K/39K ratio is 0.000 125 116 ± 57 (2σ), corresponding to 40K/K = (1.1668 ± 8; 2σ) x 10-4.

Improving a gas ion source for C-14 AMS

For more than 4 years, gaseous samples of 1-50 mu g carbon have been routinely measured with the gas ion source of the small AMS (Accelerator Mass Spectrometer) facility MICADAS (Mini CArbon DAting System) at ETH Zurich. The applied measurement technique offers a simple and fast way of C-14 measurements without the need of sample graphitization. A major drawback of gaseous C-14 measurements, however, is the relatively low negative ion current, which results in longer measurement times and lower precision compared to graphitized samples. In December 2009, a new, improved Cs sputter ion source was installed at MICADAS and we began to optimize conditions for the measurement of gaseous samples. C-12(-) currents from the new ion source were improved from initially 3 to 12-15 mu A for routine measurements and the negative ion yield was increased by a factor of 2, reaching 8 on average during routine operation. Moreover, the new measurement settings enable a doubled CO2 flow, thus substantially reducing measurement times. The achieved performance allows closing the sample size gap between gaseous and solid samples and makes the gas ion source a promising tool for dating with a measurement precision of 5 parts per thousand on samples as small as 50 mu g carbon. (C) 2012 Elsevier B.V. All rights reserved.