Analytic validation of a gas chromatography-mass spectrometry method for quantification of six amino acids in canine serum samples.

OBJECTIVE To analytically validate a gas concentration of chromatography-mass spectrometry (GC-MS) method for measurement of 6 amino acids in canine serum samples and to assess the stability of each amino acid after sample storage. SAMPLES Surplus serum from 80 canine samples submitted to the Gastrointestinal Laboratory at Texas A&M University and serum samples from 12 healthy dogs. PROCEDURES GC-MS was validated to determine precision, reproducibility, limit of detection, and percentage recovery of known added concentrations of 6 amino acids in surplus serum samples. Amino acid concentrations in serum samples from healthy dogs were measured before (baseline) and after storage in various conditions. RESULTS Intra- and interassay coefficients of variation (10 replicates involving 12 pooled serum samples) were 13.4% and 16.6% for glycine, 9.3% and 12.4% for glutamic acid, 5.1% and 6.3% for methionine, 14.0% and 15.1% for tryptophan, 6.2% and 11.0% for tyrosine, and 7.4% and 12.4% for lysine, respectively. Observed-to-expected concentration ratios in dilutional parallelism tests (6 replicates involving 6 pooled serum samples) were 79.5% to 111.5% for glycine, 80.9% to 123.0% for glutamic acid, 77.8% to 111.0% for methionine, 85.2% to 98.0% for tryptophan, 79.4% to 115.0% for tyrosine, and 79.4% to 110.0% for lysine. No amino acid concentration changed significantly from baseline after serum sample storage at -80°C for ≤ 7 days. CONCLUSIONS AND CLINICAL RELEVANCE GC-MS measurement of concentration of 6 amino acids in canine serum samples yielded precise, accurate, and reproducible results. Sample storage at -80°C for 1 week had no effect on GC-MS results.

Video-based quantification of body movement indicates negative symptoms: a replication

Objective: In schizophrenia, abnormalities in nonverbal behaviors have always been considered as highly relevant. However, due to methodological limitations, nonverbal behavior was rarely quantified objectively. Recent methodological advances now allow a quantification of body movement from ordinary video recordings. We showed that patients’ objectively measured amount of movement in social role-play interactions was closely associated with their symptom profiles (Kupper, Ramseyer, Hoffmann, & Tschacher, Schizophrenia Research 2010). In the present study, a replication of these results in the context of semi-standardized PANSS (Positive and Negative Syndrome Scale) interviews was intended. Methods: 17 patients with schizophrenia were analyzed during the initial 15-min sequence of a videotaped PANSS interview using Motion Energy Analysis (MEA). The amount of patients’ movement was then correlated with their PANSS symptom scores. Results: Sizeable and significant correlations between negative symptoms and reduced movements (r = -.68, p<0.01) and reduced movement speed (r = -.80, p<0.001) were found. Moreover, cognitive symptoms were related to reduced movement speed (r = -.70, p<.01). Conclusion: Negative symptoms were reliably indicated by patients’ nonverbal behavior in psychopathology interviews. Hence, the main result of our earlier study, examining patients’ nonverbal behavior in role play tests, was replicated for the less structured interactions in psychopathological interviews. Results could encourage the use of MEA in a wide range of videotaped social interactions of patients with schizophrenia and other psychiatric disorders.

Development of a rapid column-switching LC-MS/MS method for the quantification of THCCOOH and THCCOOH-glucuronide in whole blood for assessing cannabis consumption frequency

The concentration of 11-nor-9-carboxy-Δ(9)-tetrahydrocannabinol (THCCOOH) in whole blood is used as a parameter for assessing the consumption behavior of cannabis consumers. The blood level of THCCOOH-glucuronide might provide additional information about the frequency of cannabis use. To verify this assumption, a column-switching liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for the rapid and direct quantification of free and glucuronidated THCCOOH in human whole blood was newly developed. The method comprised protein precipitation, followed by injection of the processed sample onto a trapping column and subsequent gradient elution to an analytical column for separation and detection. The total LC run time was 4.5 min. Detection of the analytes was accomplished by electrospray ionization in positive ion mode and selected reaction monitoring using a triple-stage quadrupole mass spectrometer. The method was fully validated by evaluating the following parameters: linearity, lower limit of quantification, accuracy and imprecision, selectivity, extraction efficiency, matrix effect, carry-over, dilution integrity, analyte stability, and re-injection reproducibility. All acceptance criteria were analyzed and the predefined criteria met. Linearity ranged from 5.0 to 500 μg/L for both analytes. The method was successfully applied to whole blood samples from a large collective of cannabis consumers, demonstrating its applicability in the forensic field.

Comparison of 2.5D and 3D Quantification of Femoral Head Coverage in Normal Control Subjects and Patients with Hip Dysplasia.

Hip dysplasia is characterized by insufficient femoral head coverage (FHC). Quantification of FHC is of importance as the underlying goal of the surgery to treat hip dysplasia is to restore a normal acetabular morphology and thereby to improve FHC. Unlike a pure 2D X-ray radiograph-based measurement method or a pure 3D CT-based measurement method, previously we presented a 2.5D method to quantify FHC from a single anteriorposterior (AP) pelvic radiograph. In this study, we first quantified and compared 3D FHC between a normal control group and a patient group using a CT-based measurement method. Taking the CT-based 3D measurements of FHC as the gold standard, we further quantified the bias, precision and correlation between the 2.5D measurements and the 3D measurements on both the control group and the patient group. Based on digitally reconstructed radiographs (DRRs), we investigated the influence of the pelvic tilt on the 2.5D measurements of FHC. The intraclass correlation coefficients (ICCs) for absolute agreement was used to quantify interobserver reliability and intraobserver reproducibility of the 2.5D measurement technique. The Pearson correlation coefficient, r, was used to determine the strength of the linear association between the 2.5D and the 3D measurements. Student's t-test was used to determine whether the differences between different measurements were statistically significant. Our experimental results demonstrated that both the interobserver reliability and the intraobserver reproducibility of the 2.5D measurement technique were very good (ICCs > 0.8). Regression analysis indicated that the correlation was very strong between the 2.5D and the 3D measurements (r = 0.89, p < 0.001). Student's t-test showed that there were no statistically significant differences between the 2.5D and the 3D measurements of FHC on the patient group (p > 0.05). The results of this study provided convincing evidence demonstrating the validity of the 2.5D measurements of FHC from a single AP pelvic radiograph and proved that it could serve as a surrogate for 3D CT-based measurements. Thus it may be possible to use this method to avoid a CT scan for the purpose of estimating 3D FHC in diagnosis and post-operative treatment evaluation of patients with hip dysplasia.

Pulmonary lymphangioleiomyomatosis: analysis of disease manifestation by region-based quantification of lung parenchyma

PURPOSE Lymphangioleiomyomatosis (LAM) is characterized by proliferation of smooth muscle tissue that causes bronchial obstruction and secondary cystic destruction of lung parenchyma. The aim of this study was to evaluate the typical distribution of cystic defects in LAM with quantitative volumetric chest computed tomography (CT). MATERIALS AND METHODS CT examinations of 20 patients with confirmed LAM were evaluated with region-based quantification of lung parenchyma. Additionally, 10 consecutive patients were identified who had recently undergone CT imaging of the lung at our institution, in which no pathologies of the lung were found, to serve as a control group. Each lung was divided into three regions (upper, middle and lower thirds) with identical number of slices. In addition, we defined a "peel" and "core" of the lung comprising the 2 cm subpleural space and the remaining inner lung area. Computerized detection of lung volume and relative emphysema was performed with the PULMO 3D software (v3.42, Fraunhofer MEVIS, Bremen, Germany). This software package enables the quantification of emphysematous lung parenchyma by calculating the pixel index, which is defined as the ratio of lung voxels with a density <-950HU to the total number of voxels in the lung. RESULTS Cystic changes accounted for 0.1-39.1% of the total lung volume in patients with LAM. Disease manifestation in the central lung was significantly higher than in peripheral areas (peel median: 15.1%, core median: 20.5%; p=0.001). Lower thirds of lung parenchyma showed significantly less cystic changes than upper and middle lung areas combined (lower third: median 13.4, upper and middle thirds: median 19.0, p=0.001). CONCLUSION The distribution of cystic lesions in LAM is significantly more pronounced in the central lung compared to peripheral areas. There is a significant predominance of cystic changes in apical and intermediate lung zones compared to the lung bases.

A certified plasmid reference material for the standardisation of BCR-ABL1 mRNA quantification by real-time quantitative PCR.

Serial quantification of BCR-ABL1 mRNA is an important therapeutic indicator in chronic myeloid leukaemia, but there is a substantial variation in results reported by different laboratories. To improve comparability, an internationally accepted plasmid certified reference material (CRM) was developed according to ISO Guide 34:2009. Fragments of BCR-ABL1 (e14a2 mRNA fusion), BCR and GUSB transcripts were amplified and cloned into pUC18 to yield plasmid pIRMM0099. Six different linearised plasmid solutions were produced with the following copy number concentrations, assigned by digital PCR, and expanded uncertainties: 1.08±0.13 × 10(6), 1.08±0.11 × 10(5), 1.03±0.10 × 10(4), 1.02±0.09 × 10(3), 1.04±0.10 × 10(2) and 10.0±1.5 copies/μl. The certification of the material for the number of specific DNA fragments per plasmid, copy number concentration of the plasmid solutions and the assessment of inter-unit heterogeneity and stability were performed according to ISO Guide 35:2006. Two suitability studies performed by 63 BCR-ABL1 testing laboratories demonstrated that this set of 6 plasmid CRMs can help to standardise a number of measured transcripts of e14a2 BCR-ABL1 and three control genes (ABL1, BCR and GUSB). The set of six plasmid CRMs is distributed worldwide by the Institute for Reference Materials and Measurements (Belgium) and its authorised distributors (https://ec.europa.eu/jrc/en/reference-materials/catalogue/; CRM code ERM-AD623a-f).

Quantification of biotic responses to rapid climatic changes around the Younger Dryas — a synthesis

To assess the presence or absence of lags in biotic responses to rapid climatic changes, we: (1) assume that the δ18O in biogenically precipitated carbonates record global or hemispheric climatic change at the beginning and at the end of the Younger Dryas without any lag at our two study sites of Gerzensee and Leysin, Switzerland; (2) derive a time scale by correlating the δ18O record from these two sites with the δ18O record of the GRIP ice core; (3) measure δ18O records in ostracods and molluscs to check the record in the bulk samples and to detect possible hydrological changes; (4) analyse at Gerzensee and Leysin as well as at two additional sites (that lack carbonates and hence a δ18O record) pollen, plant macrofossils, chironomids, beetles and other insects, and Cladocera; (5) estimate our sampling resolution using the GRIP time scale for the isotope stratigraphies and the biostratigraphies; and (6) summarise the major patterns of compositional change in the biostratigraphies by principal component analysis or correspondence analysis. We conclude that, at the major climatic shifts at the beginning and end of the Younger Dryas, hardly any biotic lags occur (within the sampling resolution of 8–30 years) and that upland vegetation responded as fast as aquatic invertebrates. We suggest that the minor climatic changes associated with the Gerzensee and Preboreal oscillations were weakly recorded in the biostratigraphies at the lowland site, but were more distinct at higher altitudes. Individualistic responses of plant and animal species to climatic change may reflect processes in individuals (e.g. productivity and phenology), in populations (e.g. population dynamics), in spatial distributions (e.g. migrations), and in ecosystems (e.g. trophic state). We suggest that biotic responses may be telescoped together into relatively short periods (50 to 150 years), perhaps disrupting functional interactions among species and thus destabilising ecosystems.

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Renal Cryoablation: Investigation of Periprocedural Visualization Tools and Treatment Response Quantification

Cryoablation for small renal tumors has demonstrated sufficient clinical efficacy over the past decade as a non-surgical nephron-sparing approach for treating renal masses for patients who are not surgical candidates. Minimally invasive percutaneous cryoablations have been performed with image guidance from CT, ultrasound, and MRI. During the MRI-guided cryoablation procedure, the interventional radiologist visually compares the iceball size on monitoring images with respect to the original tumor on separate planning images. The comparisons made during the monitoring step are time consuming, inefficient and sometimes lack the precision needed for decision making, requiring the radiologist to make further changes later in the procedure. This study sought to mitigate uncertainty in these visual comparisons by quantifying tissue response to cryoablation and providing visualization of the response during the procedure. Based on retrospective analysis of MR-guided cryoablation patient data, registration and segmentation algorithms were investigated and implemented for periprocedural visualization to deliver iceball position/size with respect to planning images registered within 3.3mm with at least 70% overlap and a quantitative logit model was developed to relate perfusion deficit in renal parenchyma visualized in verification images as a result of iceball size visualized in monitoring images. Through retrospective study of 20 patient cases, the relationship between likelihood of perfusion loss in renal parenchyma and distance within iceball was quantified and iteratively fit to a logit curve. Using the parameters from the logit fit, the margin for 95% perfusion loss likelihood was found to be 4.28 mm within the iceball. The observed margin corresponds well with the clinically accepted margin of 3-5mm within the iceball. In order to display the iceball position and perfusion loss likelihood to the radiologist, algorithms were implemented to create a fast segmentation and registration module which executed in under 2 minutes, within the clinically-relevant 3 minute monitoring period. Using 16 patient cases, the average Hausdorff distance was reduced from 10.1mm to 3.21 mm with average DSC increased from 46.6% to 82.6% before and after registration.

Solar radiation over and under sea ice and photographic quantification of Ice algal aggregates during the POLARSTERN cruise ARK-XXVII/3 (IceArc) in summer 2012

The ice cover of the Arctic Ocean has been changing dramatically in the last decades and the consequences for the sea-ice associated ecosystem remain difficult to assess. Algal aggregates underneath sea ice have been described sporadically but the frequency and distribution of their occurrence is not well quantified. We used upward looking images obtained by a remotely operated vehicle (ROV) to derive estimates of ice algal aggregate biomass and to investigate their spatial distribution. During the IceArc expedition (ARK-XXVII/3) of RV Polarstern in late summer 2012, different types of algal aggregates were observed floating underneath various ice types in the Central Arctic basins. Our results show that the floe scale distribution of algal aggregates in late summer is very patchy and determined by the topography of the ice underside, with aggregates collecting in dome shaped structures and at the edges of pressure ridges. The buoyancy of the aggregates was also evident from analysis of the aggregate size distribution. Different approaches used to estimate aggregate biomass yield a wide range of results. This highlights that special care must be taken when upscaling observations and comparing results from surveys conducted using different methods or on different spatial scales.

Chemical characterization and quantification of the brown algal storage compound laminarin

In search of a meaningful stress indicator for Fucus vesiculosus we found that the often used quantitative determination procedures for the polysaccharide laminarin (beta-1,3-glucan) result in different kind of problems, uncertainties and limitations. This chemical long-term storage form of carbon enables perennial brown algae in seasonally fluctuating ecosystems to uncouple growth from photosynthesis. Because of this high ecological relevance a reliable and precise method for determination and quantification of laminarin is needed. Therefore, a simple, cold water extraction method coupled to a new quantitative liquid chromatography-mass spectrometrical method (LC-MS) was developed. Laminarin was determined in nine out of twelve brown algal species, and its expected typical molar mass distribution of 2000-7000 Da was confirmed. Furthermore, laminarin consisted of a complex mixture of different chemical forms, since fifteen chemical laminarin species with distinct molecular weights were measured in nine species of brown algae. Laminarin concentrations in the algal tissues ranged from 0.03 to 0.86% dry weight (DW). The direct chemical characterization and quantification of laminarin by LC-MS represents a powerful method to verify the biochemical and ecological importance of laminarin for brown algae. Single individuals of Laminaria hyperborea, L. digitata, Saccharina latissima, F. serratus, F. vesiculosus, F. spiralis, Himanthalia elongata, Cystoseira tamariscifolia, Pelvetia canaliculata, Ascophyllum nodosum, Halidrys siliquosa and Dictyota dichotoma were collected in fall (18.11.2013) during spring low tide from the shore of Finavarra, Co. Clare, west coast of Ireland (53° 09' 25'' N, 09° 06' 58'' W). After sampling, the different algae were immediately transported to the lab, lyophilized and sent to the University of Rostock. Laminarin was extracted with cold ultrapure water from the algal samples. Before extraction they were ground to < 1 mm grain size with an analytical mill (Ika MF 10 Basic). The algal material (approx. 1.5 g DW) was extracted in ultrapure water (8 mL) on a shaker (250 rpm) for 5 h. After the addition of surplus ultrapure water (4 mL) and shaking manually, 1 mL of the sample was filter centrifuged (45 µm) at 14,000 rpm (Hettich Mikro 22 R). The slightly viscous supernatant was free of suspended material and converted into a microvial (300 µL) for further analysis. The extracts were analyzed using liquid chromatography-mass spectrometry (LC-MS) analysis (LTQ Velos Pro ion trap spectrometer with Accela HPLC, Thermo Scientific). Laminarin species were separated on a KinetexTM column (2.6 µm C18, 150 x 3 mm). The mobile phase was 90 % ultrapure water and 10 % acetonitrile, run isocratically at a flow rate of 0.2 mL min-1. MS was working in ESI negative ion mode in a mass range of 100 - 4000 amu. Glucose contents were determined after extraction using high-performance liquid chromatography (HPLC). Extracted samples were analyzed in an HPLC (SmartLine, Knauer GmbH) equipped with a SUPELCOGELTM Ca column (30 x 7,8 mm without preColumn) and RI-detector (S2300 PDA S2800). Water was used as eluent at a flow rate of 0.8 mL min-1 at 75 °C. Glucose was quantified by comparison of the retention time and peak area with standard solutions using ChromGate software. Mannitol was extracted from three subsamples of 10-20 mg powdered alga material (L. hyperborea, L. digitata, S. latissima, F. serratus, F. vesiculosus, F. spiralis, H. elongata, P. canaliculata, A. nodosum, H. siliquosa) and quantified, following the HPLC method described by Karsten et al. (1991). For analyzing carbon and nitrogen contents, dried algal material was ground to powder and three subsamples of 2 mg from each alga thalli were loaded and packed into tin cartridges (6×6×12 mm). The packages were combusted at 950 °C and the absolute contents of C and N were automatically quantified in an elemental analyzer (Elementar Vario EL III, Germany) using acetanilide as standard according to Verardo et al. (1990).

Results of bulk sediment X-ray diffraction analysis and quantification of mineral phases based on the RockJock and on the QUAX quantitative analysis

We have performed quantitative X-ray diffraction (qXRD) analysis of 157 grab or core-top samples from the western Nordic Seas between (WNS) ~57°-75°N and 5° to 45° W. The RockJock Vs6 analysis includes non-clay (20) and clay (10) mineral species in the <2 mm size fraction that sum to 100 weight %. The data matrix was reduced to 9 and 6 variables respectively by excluding minerals with low weight% and by grouping into larger groups, such as the alkali and plagioclase feldspars. Because of its potential dual origins calcite was placed outside of the sum. We initially hypothesized that a combination of regional bedrock outcrops and transport associated with drift-ice, meltwater plumes, and bottom currents would result in 6 clusters defined by "similar" mineral compositions. The hypothesis was tested by use of a fuzzy k-mean clustering algorithm and key minerals were identified by step-wise Discriminant Function Analysis. Key minerals in defining the clusters include quartz, pyroxene, muscovite, and amphibole. With 5 clusters, 87.5% of the observations are correctly classified. The geographic distributions of the five k-mean clusters compares reasonably well with the original hypothesis. The close spatial relationship between bedrock geology and discrete cluster membership stresses the importance of this variable at both the WNS-scale and at a more local scale in NE Greenland.

Mineral quantification of a Jordan Cretaceous outcrop (Section GM3, Ghawr Al Mazar, Ghor al Mazrar)

Cyclic fluctuations in global sea level during epochs of warm greenhouse climate have remained enigmatic, because absence or subordinate presence of polar ice during these periods precludes an explanation by glacio-eustatic forcing. An alternative concept suggests that the water-bearing potential of groundwater aquifers is equal to that of ice caps and that changes in the dynamic balance of aquifer charge versus discharge, as a function of the temperature-related intensity of the hydrological cycle, may have driven eustasy during warm climates. However, this idea has long been neglected for two reasons: 1) the large storage potential of subsurface aquifers was confused with the much smaller capacity of rivers and lakes and 2) empirical data were missing that document past variations in the hydrological cycle in relation to eustasy. In the present study we present the first empirical evidence for changes in precipitation, continental weathering intensity and evaporation that correlate with astronomically (long obliquity) forced sea-level cycles during the warmest period of the Cretaceous (Cenomanian-Turonian). We compare sequence-stratigraphic data with changes in the terrigenous mineral assemblage in a low-latitude marine sedimentary sequence from the equatorial humid belt at the South-Tethyan margin (Levant carbonate platform, Jordan), thereby avoiding uncertainties from land-ocean correlations. Our data indicate covariance between cycles in weathering and sea level: predominantly chemical weathering under wet climate conditions is reflected by dominance of weathering products (clays) in deposits that represent sea-level fall (aquifer charge > discharge). Conversely, preservation of weathering-sensitive minerals (feldspars, epidote and pyroxenes) in transgressive sediments reflects decreased continental weathering due to dryer climate (aquifer discharge > charge). Based on our results we suggest that aquifer-eustasy represents a viable alternative to glacio-eustasy as a driver of cyclic 3rd-order sea-level fluctuations during the middle Cretaceous greenhouse climate, and it may have been a pervasive process throughout Earth history.

Radiocarbon measurements and loose-sediment quantification of seafloor sediment samples off northern Mauritania

Heterozoan carbonates are typical for extratropical sedimentary systems. However, under mesotrophic to eutrophic conditions, heterozoan carbonates also form in tropical settings. Nevertheless, such heterozoan tropical sedimentary systems are rare in the modern world and therefore are only poorly understood to date. Here a carbonate depositional system is presented where nutrient-rich upwelling waters push onto a wide shelf. These waters warm up in the shelf, giving rise to the production and deposition of tropical heterozoan facies. The carbonate facies on this shelf are characterized by a mixture of tropical and cosmopolitan biogenic sedimentary grains. Study of facies and taxonomy are the key for identifying and characterizing tropical heterozoan carbonates and for distinguishing them from their coolwater counterparts, in particular in the past where the oceanography cannot be determined directly.