Evaluation of reference-based two-color methods for measurement of gene expression ratios using spotted cDNA microarrays

Abstract Background Spotted cDNA microarrays generally employ co-hybridization of fluorescently-labeled RNA targets to produce gene expression ratios for subsequent analysis. Direct comparison of two RNA samples in the same microarray provides the highest level of accuracy; however, due to the number of combinatorial pair-wise comparisons, the direct method is impractical for studies including large number of individual samples (e.g., tumor classification studies). For such studies, indirect comparisons using a common reference standard have been the preferred method. Here we evaluated the precision and accuracy of reconstructed ratios from three indirect methods relative to ratios obtained from direct hybridizations, herein considered as the gold-standard. Results We performed hybridizations using a fixed amount of Cy3-labeled reference oligonucleotide (RefOligo) against distinct Cy5-labeled targets from prostate, breast and kidney tumor samples. Reconstructed ratios between all tissue pairs were derived from ratios between each tissue sample and RefOligo. Reconstructed ratios were compared to (i) ratios obtained in parallel from direct pair-wise hybridizations of tissue samples, and to (ii) reconstructed ratios derived from hybridization of each tissue against a reference RNA pool (RefPool). To evaluate the effect of the external references, reconstructed ratios were also calculated directly from intensity values of single-channel (One-Color) measurements derived from tissue sample data collected in the RefOligo experiments. We show that the average coefficient of variation of ratios between intra- and inter-slide replicates derived from RefOligo, RefPool and One-Color were similar and 2 to 4-fold higher than ratios obtained in direct hybridizations. Correlation coefficients calculated for all three tissue comparisons were also similar. In addition, the performance of all indirect methods in terms of their robustness to identify genes deemed as differentially expressed based on direct hybridizations, as well as false-positive and false-negative rates, were found to be comparable. Conclusion RefOligo produces ratios as precise and accurate as ratios reconstructed from a RNA pool, thus representing a reliable alternative in reference-based hybridization experiments. In addition, One-Color measurements alone can reconstruct expression ratios without loss in precision or accuracy. We conclude that both methods are adequate options in large-scale projects where the amount of a common reference RNA pool is usually restrictive.

Peak Last Glacial weathering intensity on the North American continent recorded by the authigenic Hf isotope composition of North Atlantic deep-sea sediments

We have retrieved radiogenic hafnium (Hf) isotope compositions (ɛHf) from authigenic Fe–Mn oxyhydroxides of deep northwest Atlantic sediments deposited over the past 26 ka to investigate the oceanic evidence of changes in dissolved weathering inputs from NE America during the last deglaciation. The extraction of seawater-derived Hf isotopic compositions from Fe–Mn oxyhydroxides is not a standard procedure. Comparisons between the Al/Hf ratios and Hf isotopic compositions of the chemically extracted authigenic phase on the one hand, and those of the corresponding detrital fractions on the other, provide evidence that the composition of past seawater has been reliably obtained for most sampled depths with our leaching procedures. This is endorsed most strongly by data for a sediment core from 4250 m water depth at the deeper Blake Ridge, for which consistent replicates were produced throughout. The Hf isotopic composition of the most recent sample in this core also closely matches that of nearby present day central North Atlantic seawater. Comparison with previously published seawater Nd and Pb isotope compositions obtained on the same cores shows that both Hf and Pb were released incongruently during incipient chemical weathering, but responded differently to the deglacial retreat of the Laurentide Ice Sheet. Hafnium was released more congruently during peak glacial conditions of the Last Glacial Maximum (LGM) and changed to typical incongruent interglacial ɛHf signatures either during or shortly after the LGM. This indicates that some zircon-derived Hf was released to seawater during the LGM. Conversely, there is no clear evidence for an increase in the influence of weathering of Lu-rich mineral phases during deglaciation, possibly since relatively unradiogenic Hf contributions from feldspar weathering were superimposed. While the authigenic Pb isotope signal in the same marine sediment samples traced peak chemical weathering rates on continental North America during the transition to the Holocene a similar incongruent excursion is notably absent in the Hf isotope record. The early change towards more radiogenic ɛHf in relation to the LGM may provide direct evidence for the transition from a cold-based to a warm-based Laurentide Ice Sheet on the Atlantic sector of North America.

Comparison of clinical and economic outcomes between high-intensity care and hospice care for the end-of-life melanoma patients

The intensity of care for patients at the end-of-life is increasing in recent years. Publications have focused on intensity of care for many cancers, but none on melanoma patients. Substantial gaps exist in knowledge about intensive care and its alternative, hospice care, among the advanced melanoma patients at the end of life. End-of-life care may be used in quite different patterns and induce both intended and unintended clinical and economic consequences. We used the Surveillance, Epidemiology, and End Results (SEER)-Medicare linked databases to identify patients aged 65 years or older with metastatic melanoma who died between 2000 and 2007. We evaluated trends and associations between sociodemographic and health services characteristics and the use of hospice care, chemotherapy, surgery, and radiation therapy and costs. Survival, end-of-life costs, and incremental cost-effectiveness ratio were evaluated using propensity score methods. Costs were analyzed from the perspective of Medicare in 2009 dollars. In the first journal Article we found increasing use of surgery for patients with metastatic melanoma from 13% in 2000 to 30% in 2007 (P=0.03 for trend), no significant fluctuation in use of chemotherapy (P=0.43) or radiation therapy (P=0.46). Older patients were less likely to receive radiation therapy or chemotherapy. The use of hospice care increased from 61% in 2000 to 79% in 2007 (P =0.07 for trend). Enrollment in short-term (1-3 days) hospice care use increased, while long-term hospice care (≥ 4 days) remained stable. Patients living in the SEER Northeast and South regions were less likely to undergo surgery. Patients enrolled in long-term hospice care used significantly less chemotherapy, surgery and radiation therapy. In the second journal article, of 611 patients identified for this study, 358 (59%) received no hospice care after their diagnosis, 168 (27%) received 1 to 3 days of hospice care, and 85 (14%) received 4 or more days of hospice care. The median survival time was 181 days for patients with no hospice care, 196 days for patients enrolled in hospice for 1 to 3 days, and 300 days for patients enrolled for 4 or more days (log-rank test, P < 0.001). The estimated hazard ratios (HR) between 4 or more days hospice use and survival were similar within the original cohort Cox proportional hazard model (HR, 0.62; 95% CI, 0.49-0.78, P < 0.0001) and the propensity score-matched model (HR, 0.61; 95% CI, 0.47-0.78, P = 0.0001). Patients with ≥ 4 days of hospice care incurred lower end-of-life costs than the other two groups ($14,298 versus $19,380 for the 1- to 3-days hospice care, and $24,351 for patients with no hospice care; p < 0.0001). In conclusion, Surgery and hospice care use increased over the years of this study while the use of chemotherapy and radiation therapy remained consistent for patients diagnosed with metastatic melanoma. Patients diagnosed with advanced melanoma who enrolled in ≥ 4 days of hospice care experienced longer survival than those who had 1-3 days of hospice or no hospice care, and this longer overall survival was accompanied by lower end-of-life costs.^

Magnesium/Calcium ratios and sea surface temperature estimation for sediments of the Timor Sea

We use a multiproxy approach to monitor changes in the vertical profile of the Indonesian Throughflow as well as monsoonal wind and precipitation patterns in the Timor Sea on glacial-interglacial, precessional, and suborbital timescales. We focus on an interval of extreme climate change and sea level variation: marine isotope (MIS) 6 to MIS 5e. Paleoproductivity fluctuations in the Timor Sea follow a precessional beat related to the intensity of the Australian (NW) monsoon. Paired Mg/Ca and d18O measurements of surface- and thermocline-dwelling planktonic foraminifers (G. ruber and P. obliquiloculata) indicate an increase of >4°C in both surface and thermocline water temperatures during Termination II. Tropical sea surface temperature changed synchronously with ice volume (benthic d18O) during deglaciation, implying a direct coupling of high- and low-latitude climate via atmospheric and/or upper ocean circulation. Substantial cooling and freshening of thermocline waters occurred toward the end of Termination II and during MIS 5e, indicating a change in the vertical profile of the Indonesian Throughflow from surface- to thermocline-dominated flow.

Benthic and planktic foraminiferal stable isotope records combined with Nd isotope ratios of ODP Hole 198-1210B

The Pacific Ocean is the largest water body on Earth, and circulation in the Pacific contributed significantly to climate evolution in the latest Cretaceous, the culmination of a period of long-term cooling. Here, we present new high-resolution late Campanian to Maastrichtian benthic and planktic foraminiferal stable isotope data and a neodymium (Nd) isotope record obtained from sedimentary ferromanganese oxide coatings of Ocean Drilling Program Hole 1210B from the tropical Pacific Ocean (Shatsky Rise). These new records resolve 13 million years in the latest Cretaceous, providing insights into changes in surface and bottom water temperatures and source regions of deep to intermediate waters covering the carbon isotope excursions of the Campanian-Maastrichtian Boundary Event (CMBE) and the Mid-Maastrichtian event (MME). Our new benthic foraminiferal d18O and Nd isotope records together with published Nd isotope data show markedly parallel trends across the studied interval over a broad range of bathyal to abyssal water depths interpreted to reflect changes in the intensity of deep-ocean circulation in the tropical Pacific. In particular, we observe a three-million-year-long period of cooler conditions in the early Maastrichtian (72.5 to 69.5 Ma) when a concomitant change toward less radiogenic seawater Nd isotope signatures probably marks a period of enhanced admixture and northward flow of deep waters with Southern Ocean provenance. We suggest this change to have been triggered by intensified formation and convection of deep waters in the high southern latitudes, a process that weakened during the MME (69.5 to 68.5 Ma). The early Maastrichtian cold interval is closely related to the negative and positive carbon isotope trends of the CMBE and MME. The millions-of-years long duration of these carbon cycle perturbations suggests a tectonic forcing of climatic cooling, possibly related to changes in ocean basin geometry and bathymetry.

(Table 2) Contents of chlorophyll a, organic carbon, nitrogen, and phosphorus and their molar ratios in dried matter of Ahnfeltia tobuchiensis agar containing seaweed from the Izmena Bay, Kunashir Island (South Kurils)

A complex study of influence of various environmental factors on rates of oxygen (M_O2 ), ammonium (M_NH4), and phosphate (M_PO4) metabolism in Ahnfeltia tobuchiensis has been carried out in situ in the Izmena Bay of the Kunashir Island. The following environmental factors have been included into the investigation: photosynthetically active radiation (PAR); ammonium (NH4); phosphate (PO4); and contents of carbon (C), nitrogen (N), phosphorus (P), and chlorophyll a (Chl) in tissue. Population of agar-containing seaweed A. tobuchiensis forms a layer with thickness up to 0.5 m, which occupies about 23.3 km**2; biomass is equal to 125000 tons. Quantitative assessment of organic matter production and nutrient consumption during oxygen metabolism has been carried out for the whole population. It has been shown that daily oxygen metabolism depends on PAR intensity, concentrations of PO4 and NH4 in seawater, and contents of N and P in tissues (r**2=0.78, p<0.001). Average daily NH4 consumption is 0.21 µmol/g of dry weight/hour and depends on NH4 and O2 concentrations in seawater and on ? and Chl a contents in algal tissues (r**2=0.64, p<0.001). Average daily PO4 consumption is 0.01 µmol/g of dry weight/hour and depends on NH4 concentrations in seawater and on P contents in algal tissues (r**2=0.40, p<0.001).

Sedimentology and Lead and Strontium isotope ratios of Miocene sediments from the Kerguelen Plateau

Characterization of sediment from Ocean Drilling Program Site 745, representing the East Kerguelen Ridge sediment drift, addresses important issues surrounding the timing of Miocene to present East Antarctic ice sheet stability and oceanic environmental change. Our results show three periods of greatly enhanced accumulation of Antarctic-derived sediment, at 6.4-5.9 Ma, 4.9-4.4 Ma and 1.1-0.8 Ma, potentially indicative of warmer, less stable ice sheets at these times. Conversely, the accumulation of Antarctic-derived material is comparatively less during the middle of the Pliocene warm epoch (4.8-3.2 Ma). The deep flow forming the Kerguelen drift was stronger during the latest Miocene and earliest Pliocene and has decreased in intensity continuously since then.

Paleomagnetic, stable oxygen isotope ratios, and U/Th activities of ODP Site 172-1062

Geomagnetic excursions are recognized as intrinsic features of the Earth's magnetic field. High-resolution records of field behaviour, captured in marine sedimentary cores, present an opportunity to determine the temporal and geometric character of the field during geomagnetic excursions and provide constraints on the mechanisms producing field variability. We present here the highest resolution record yet published of the Blake geomagnetic excursion (~125 ka) measured in three cores from Ocean Drilling Program (ODP) Site 1062 on the Blake-Bahama Outer Ridge. The Blake excursion has a controversial structure and timing but these cores have a sufficiently high sedimentation rate (~10cm/ka) to allow detailed reconstruction of the field behaviour at this site during the excursion. Palaeomagnetic measurements of the cores reveal rapid transitions (<500 yr) between the contemporary stable normal polarity and a completely reversed state of long duration which spans a stratigraphic interval of 0.7 m. We determine the duration of the reversed state during the Blake excursion using oxygen isotope stratigraphy, combined with 230Th excess measurements to assess variations in the sedimentation rates through the sections of interest. This provides an age and duration for the Blake excursion with greater accuracy and with constrained uncertainty. We date the directional excursion as falling between 129 and 122 ka with a duration for the deviation of 6.5±1.3 kyr. The long duration of this interval and the fully reversed field suggest the existence of a pseudo-stable, reversed dipole field component during the excursion and challenge the idea that excursions are always of short duration.

Separating ITCZ- and ENSO-related rainfall changes in the Galápagos over the last 3 kyr using D/H ratios of multiple lipid biomarkers

We present a 3000-yr rainfall reconstruction from the Galápagos Islands that is based on paired biomarker records from the sediment of El Junco Lake. Located in the eastern equatorial Pacific, the climate of the Galápagos Islands is governed by movements of the Intertropical Convergence Zone (ITCZ) and the El Niño-Southern Oscillation (ENSO). We use a novel method for reconstructing past ENSO- and ITCZ-related rainfall changes through analysis of molecular and isotopic biomarker records representing several types of plants and algae that grow under differing climatic conditions. We propose that ?D values of dinosterol, a sterol produced by dinoflagellates, record changes in mean rainfall in El Junco Lake, while dD values of C34 botryococcene, a hydrocarbon unique to the green alga Botryococcus braunii, record changes in rainfall associated with moderate-to-strong El Niño events. We use these proxies to infer changes in mean rainfall and El Niño-related rainfall over the past 3000 yr. During periods in which the inferred change in El Niño-related rainfall opposed the change in mean rainfall, we infer changes in the amount of ITCZ-related rainfall. Simulations with an idealized isotope hydrology model of El Junco Lake help illustrate the interpretation of these proxy reconstructions. Opposing changes in El Niño- and ITCZ-related rainfall appear to account for several of the largest inferred hydrologic changes in El Junco Lake. We propose that these reconstructions can be used to infer changes in frequency and/or intensity of El Niño events and changes in the position of the ITCZ in the eastern equatorial Pacific over the past 3000 yr. Comparison with El Junco Lake sediment grain size records indicates general agreement of inferred rainfall changes over the late Holocene.

Delta 66Zn ratios and age of carbonates from ODP Site 138-849 (Table 1)

The carbonate fraction of sediment core ODP 849, leg 138, located in the eastern equatorial Pacific, mostly consisting of coccoliths, was separated and analyzed for its Zn isotopic composition. The overall variation in Zn isotopic composition, as determined by multiple-collector, magnetic-sector, inductively coupled plasma mass spectrometry, was found to be on the order of 1? (expressed in delta66Zn, where deltaxZn=[(xZn/64Zn)sample/(xZn/64Zn)standard -1]*10**3 and x=66, 67 or 68) over the last 175 ka. The analytical precision was 0.04 per mil and the overall reproducibility was usually better than 0.07 per mil. The Zn isotopic composition signal exhibits several marked peaks and a high-frequency variability. A periodogram of the delta66Zn signal showed two periodicities of 35.2 and 21.2 ka. We suggest that the latter is caused by the precession of the Earth's axis of rotation. The periodogram exhibits a minimum at 41.1 ka, thus showing that the Zn isotopic composition is independent of the obliquity in the eastern equatorial Pacific. The range of delta66Zn values observed for the carbonate fraction of ODP 849 overlaps with the range observed for Fe-Mn nodules in the world's oceans, which suggests that seawater/carbonate Zn isotope fractionation is weak. We therefore assume that most of the Zn isotope variability is a result of the selective entrainment of the light isotopes by organic matter in the surface ocean. The ODP 849 delta66Zn record seems to follow the changes in the insolation cycles. Changes in the late summer/fall equatorial insolation modulate the intensity of the equatorial upwelling, hence the mixing between deep and surface waters. We propose that during decreased summer/fall equatorial insolation, when a steep thermocline can develop (El Niño-like conditions), the surface waters cannot be replenished by deep waters and become depleted in the lighter Zn isotopes by biological activity, thus resulting in the progressive increase of the delta66Zn values of the carbonate shells presumably in equilibrium with surface seawater.

Pore water isotopic ratios of CO2 and SO4, and sediment CH4 content of ODP Leg 164 sites

A unique set of geochemical pore-water data, characterizing the sulfate reduction and uppermost methanogenic zones, has been collected at the Blake Ridge (offshore southeastern North America) from Ocean Drilling Program (ODP) Leg 164 cores and piston cores. The d13C values of dissolved CO2 (sum CO2) are as 13C-depleted as -37.7 per mil PDB (Site 995) at the sulfate-methane interface, reflecting a substantial contribution of isotopically light carbon from methane. Although the geochemical system is complex and difficult to fully quantify, we use two methods to constrain and illustrate the intensity of anaerobic methane oxidation in Blake Ridge sediments. An estimate using a two-component mixing model suggests that ~24% of the carbon residing in the sum CO2 pool is derived from biogenic methane. Independent diagenetic modeling of a methane concentration profile (Site 995) indicates that peak methane oxidation rates approach 0.005 µmol/cm**3/yr, and that anaerobic methane oxidation is responsible for consuming ~35% of the total sulfate flux into the sediments. Thus, anaerobic methane oxidation is a significant biogeochemical sink for sulfate, and must affect interstitial sulfate concentrations and sulfate gradients. Such high proportions of sulfate depletion because of anaerobic methane oxidation are largely undocumented in continental rise sediments with overlying oxic bottom waters. We infer that the additional amount of sulfate depleted through anaerobic methane oxidation, fueled by methane flux from below, causes steeper sulfate gradients above methane-rich sediments. Similar pore water chemistries should occur at other methane-rich, continental-rise settings associated with gas hydrates.

Stable isotope ratios of planktonic foraminifera from the West Equatorial Pacific

Differential dissolution affects the isotopic composition of different species of planktonic foraminifera in different ways. In the two species studied here in cores from Ontong Java Plateau, the less resistant species, Globigerinoides sacculifer, is more readily affected at a shallower depth than the more resistant species, Pulleniatina obliquiloculata (2.9 versus 3.4 km), but shows a smaller and less predictable response to partial dissolution (0.2 to 0.3 per mil versus 0.6 to 0.7 per mil). Comparison of isotopic values from the last glacial period with those from the late Holocene indicates that the apparent dissolution effect is considerably reduced during the last glacial, presumably due to reduced dissolution intensity during glacial time. A change in the level of the lysocline of about 400 m is suggested. In the published isotope records from Pacific cores V28-238 and V28-239, the dissolution-generated difference in delta18O (noted previously by Shackleton and Opdyke [1976]) is seen to describe a mid-Brunhes dissolution maximum, between 300 and 500 thousand years ago. This mid-Brunhes dissolution excursion is well known from the Pacific and the Indian Ocean.

AMS 14C dating points, isotopic compositions, Mg/Ca ratios and sea surface temperatures of sediment cores GeoB10029-4 and GeoB10038-4

Quantifying the spatial and temporal sea surface temperature (SST) and salinity changes of the Indo-Pacific Warm Pool is essential to understand the role of this region in connection with abrupt climate changes particularly during the last deglaciation. In this study we reconstruct SST and seawater d18O of the tropical eastern Indian Ocean for the past 40,000 years from two sediment cores (GeoB 10029-4, 1°30'S, 100°08'E, and GeoB 10038-4, 5°56'S, 103°15'E) retrieved offshore Sumatra. Our results show that annual mean SSTs increased about 2-3 °C at 19,000 years ago and exhibited southern hemisphere-like timing and pattern during the last deglaciation. Our SST records together with other Mg/Ca-based SST reconstructions around Indonesia do not track the monsoon variation since the last glacial period, as recorded by terrestrial monsoon archives. However, the spatial SST heterogeneity might be a result of changing monsoon intensity that shifts either the annual mean SSTs or the seasonality of G. ruber towards the warmer or the cooler season at different locations. Seawater d18O reconstructions north of the equator suggest fresher surface conditions during the last glacial and track the northern high-latitude climate change during the last deglaciation. In contrast, seawater ?18O records south of the equator do not show a significant difference between the last glacial period and the Holocene, and lack Bølling-Allerød and Younger Dryas periods suggestive of additional controls on annual mean surface hydrology in this part of the Indo-Pacific Warm Pool.

Carbon/oxygen ratios, clay composition and gas-hydrate saturation of ODP Leg 164 sites

The analyses of downhole log data from Ocean Drilling Program (ODP) boreholes on the Blake Ridge at Sites 994, 995, and 997 indicate that the Schlumberger geochemical logging tool (GLT) may yield useful gas hydrate reservoir data. In neutron spectroscopy downhole logging, each element has a characteristic gamma ray that is emitted from a given neutron-element interaction. Specific elements can be identified by their characteristic gamma-ray signature, with the intensity of emission related to the atomic elemental concentration. By combining elemental yields from neutron spectroscopy logs, reservoir parameters including porosities, lithologies, formation fluid salinities, and hydrocarbon saturations (including gas hydrate) can be calculated. Carbon and oxygen elemental data from the GLT was used to determine gas hydrate saturations at all three sites (Sites 994, 995, and 997) drilled on the Blake Ridge during Leg 164. Detailed analyses of the carbon and oxygen content of various sediments and formation fluids were used to construct specialized carbon/oxygen ratio (COR) fan charts for a series of hypothetical gas hydrate accumulations. For more complex geologic systems, a modified version of the standard three-component COR hydrocarbon saturation equation was developed and used to calculate gas hydrate saturations on the Blake Ridge. The COR-calculated gas hydrate saturations (ranging from about 2% to 14% bulk volume gas hydrate) from the Blake Ridge compare favorably to the gas hydrate saturations derived from electrical resistivity log measurements.