Beneficial effect of recruitable collaterals: a 10-year follow-up study in patients with stable coronary artery disease undergoing quantitative collateral measurements

BACKGROUND: The prognostic relevance of the collateral circulation is still controversial. The goal of this study was to assess the impact on survival of quantitatively obtained, recruitable coronary collateral flow in patients with stable coronary artery disease during 10 years of follow-up. METHODS AND RESULTS: Eight-hundred forty-five individuals (age, 62+/-11 years), 106 patients without coronary artery disease and 739 patients with chronic stable coronary artery disease, underwent a total of 1053 quantitative, coronary pressure-derived collateral measurements between March 1996 and April 2006. All patients were prospectively included in a collateral flow index (CFI) database containing information on recruitable collateral flow parameters obtained during a 1-minute coronary balloon occlusion. CFI was calculated as follows: CFI = (P(occl) - CVP)/(P(ao) - CVP) where P(occl) is mean coronary occlusive pressure, P(ao) is mean aortic pressure, and CVP is central venous pressure. Patients were divided into groups with poorly developed (CFI < 0.25) or well-grown collateral vessels (CFI > or = 0.25). Follow-up information on the occurrence of all-cause mortality and major adverse cardiac events after study inclusion was collected. Cumulative 10-year survival rates in relation to all-cause deaths and cardiac deaths were 71% and 88%, respectively, in patients with low CFI and 89% and 97% in the group with high CFI (P=0.0395, P=0.0109). Through the use of Cox proportional hazards analysis, the following variables independently predicted elevated cardiac mortality: age, low CFI (as a continuous variable), and current smoking. CONCLUSIONS: A well-functioning coronary collateral circulation saves lives in patients with chronic stable coronary artery disease. Depending on the exact amount of collateral flow recruitable during a brief coronary occlusion, long-term cardiac mortality is reduced to one fourth compared with the situation without collateral supply.

Multilead quantitative electroencephalogram profile and cognitive evoked potentials (P300) in healthy subjects after a single dose of olanzapine

RATIONALE: Olanzapine is an atypical antipsychotic drug with a more favourable safety profile than typical antipsychotics with a hitherto unknown topographic quantitative electroencephalogram (QEEG) profile. OBJECTIVES: We investigated electrical brain activity (QEEG and cognitive event related potentials, ERPs) in healthy subjects who received olanzapine. METHODS: Vigilance-controlled, 19-channel EEG and ERP in an auditory odd-ball paradigm were recorded before and 3 h, 6 h and 9 h after administration of either a single dose of placebo or olanzapine (2.5 mg and 5 mg) in ten healthy subjects. QEEG was analysed by spectral analysis and evaluated in nine frequency bands. For the P300 component in the odd-ball ERP, the amplitude and latency was analysed. Statistical effects were tested using a repeated-measurement analysis of variance. RESULTS: For the interaction between time and treatment, significant effects were observed for theta, alpha-2, beta-2 and beta-4 frequency bands. The amplitude of the activity in the theta band increased most significantly 6 h after the 5-mg administration of olanzapine. A pronounced decrease of the alpha-2 activity especially 9 h after 5 mg olanzapine administration could be observed. In most beta frequency bands, and most significantly in the beta-4 band, a dose-dependent decrease of the activity beginning 6 h after drug administration was demonstrated. Topographic effects could be observed for the beta-2 band (occipital decrease) and a tendency for the alpha-2 band (frontal increase and occipital decrease), both indicating a frontal shift of brain electrical activity. There were no significant changes in P300 amplitude or latency after drug administration. Conclusion: QEEG alterations after olanzapine administration were similar to EEG effects gained by other atypical antipsychotic drugs, such as clozapine. The increase of theta activity is comparable to the frequency distribution observed for thymoleptics or antipsychotics for which treatment-emergent somnolence is commonly observed, whereas the decrease of beta activity observed after olanzapine administration is not characteristic for these drugs. There were no clear signs for an increased cerebral excitability after a single-dose administration of 2.5 mg and 5 mg olanzapine in healthy controls.

A novel quantitative method for analyzing the distributions of nanoparticles between different tissue and intracellular compartments

The penetration, translocation, and distribution of ultrafine and nanoparticles in tissues and cells are challenging issues in aerosol research. This article describes a set of novel quantitative microscopic methods for evaluating particle distributions within sectional images of tissues and cells by addressing the following questions: (1) is the observed distribution of particles between spatial compartments random? (2) Which compartments are preferentially targeted by particles? and (3) Does the observed particle distribution shift between different experimental groups? Each of these questions can be addressed by testing an appropriate null hypothesis. The methods all require observed particle distributions to be estimated by counting the number of particles associated with each defined compartment. For studying preferential labeling of compartments, the size of each of the compartments must also be estimated by counting the number of points of a randomly superimposed test grid that hit the different compartments. The latter provides information about the particle distribution that would be expected if the particles were randomly distributed, that is, the expected number of particles. From these data, we can calculate a relative deposition index (RDI) by dividing the observed number of particles by the expected number of particles. The RDI indicates whether the observed number of particles corresponds to that predicted solely by compartment size (for which RDI = 1). Within one group, the observed and expected particle distributions are compared by chi-squared analysis. The total chi-squared value indicates whether an observed distribution is random. If not, the partial chi-squared values help to identify those compartments that are preferential targets of the particles (RDI > 1). Particle distributions between different groups can be compared in a similar way by contingency table analysis. We first describe the preconditions and the way to implement these methods, then provide three worked examples, and finally discuss the advantages, pitfalls, and limitations of this method.

Visualization and quantitative analysis of nanoparticles in the respiratory tract by transmission electron microscopy

ABSTRACT: Nanotechnology in its widest sense seeks to exploit the special biophysical and chemical properties of materials at the nanoscale. While the potential technological, diagnostic or therapeutic applications are promising there is a growing body of evidence that the special technological features of nanoparticulate material are associated with biological effects formerly not attributed to the same materials at a larger particle scale. Therefore, studies that address the potential hazards of nanoparticles on biological systems including human health are required. Due to its large surface area the lung is one of the major sites of interaction with inhaled nanoparticles. One of the great challenges of studying particle-lung interactions is the microscopic visualization of nanoparticles within tissues or single cells both in vivo and in vitro. Once a certain type of nanoparticle can be identified unambiguously using microscopic methods it is desirable to quantify the particle distribution within a cell, an organ or the whole organism. Transmission electron microscopy provides an ideal tool to perform qualitative and quantitative analyses of particle-related structural changes of the respiratory tract, to reveal the localization of nanoparticles within tissues and cells and to investigate the 3D nature of nanoparticle-lung interactions.This article provides information on the applicability, advantages and disadvantages of electron microscopic preparation techniques and several advanced transmission electron microscopic methods including conventional, immuno and energy-filtered electron microscopy as well as electron tomography for the visualization of both model nanoparticles (e.g. polystyrene) and technologically relevant nanoparticles (e.g. titanium dioxide). Furthermore, we highlight possibilities to combine light and electron microscopic techniques in a correlative approach. Finally, we demonstrate a formal quantitative, i.e. stereological approach to analyze the distributions of nanoparticles in tissues and cells.This comprehensive article aims to provide a basis for scientists in nanoparticle research to integrate electron microscopic analyses into their study design and to select the appropriate microscopic strategy.

Magnetic resonance (MR) cholangiography: quantitative and qualitative comparison of 3.0 Tesla with 1.5 Tesla

OBJECTIVES: To determine quantitative and qualitative image quality in patients undergoing magnetic resonance (MR) cholangiography at 3.0 Tesla (T) compared with 1.5 T. MATERIALS AND METHODS: Fifty patients (30 women; mean age, 51 years) underwent MR cholangiography at 1.5 T; another 50 patients (25 women; mean age 51 years) were scanned at 3.0 T. MR sequence protocol consisted of breath-hold single-slice rapid acquisition with relaxation enhancement (RARE) and a respiratory-triggered 3D turbo spin echo (3D TSE) sequence. Maximum intensity projections were generated from the 3D TSE datasets. Contrast-to-noise ratio (CNR) measurements between the common bile duct (CBD), left and right intrahepatic duct (LHD, RHD), and periductal tissue were performed. Three radiologists assessed qualitatively the visibility of the CBD, LHD, and RHD and the overall diagnostic quality. RESULTS: Mean gain in CNR at 3.0 T versus 1.5 T in all 3 locations ranged for the RARE sequence from 7.7% to 38.1% and for the 3D TSE from 0.5% to 26.1% (P > 0.05 for all differences). Qualitative analysis did not reveal any significant difference between the 2 field strengths (P > 0.05). CONCLUSIONS: MR cholangiography at 3.0 T shows a trend toward higher CNR without improving image quality significantly.

Quantitative assessment of soil parameters in Western Tajikistan using a soil spectral library approach

Soil degradation is a major problem in the agriculturally dominated country of Tajikistan, which makes it necessary to determine and monitor the state of soils. For this purpose a soil spectral library was established as it enables the determination of soil properties with relatively low costs and effort. A total of 1465 soil samples were collected from three 10x10 km test sites in western Tajikistan. The diffuse reflectance of the samples was measured with a FieldSpec PRO FR from ASD in the spectral range from 380 to 2500 nm in laboratory. 166 samples were finally selected based on their spectral information and analysed on total C and N, organic C, pH, CaCO₃, extractable P, exchangeable Ca, Mg and K, and the fractions clay, silt and sand. Multiple linear regression was used to set up the models. Two third of the chemically analysed samples were used to calibrate the models, one third was used for hold-out validation. Very good prediction accuracy was obtained for total C (R² = 0.76, RMSEP = 4.36 g kg⁻¹), total N (R² = 0.83, RMSEP = 0.30 g kg⁻¹) and organic C (R² = 0.81, RMSEP = 3.30 g kg⁻¹), good accuracy for pH (R² = 0.61, RMSEP = 0.157) and CaCO3(R² = 0.72, RMSEP = 4.63 %). No models could be developed for extractable P, exchangeable Ca, Mg and K, and the fractions clay, silt and sand. It can be concluded that the spectral library approach has a high potential to substitute standard laboratory methods where rapid and inexpensive analysis is required.

Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.

Effect of echo time pair selection on quantitative analysis for adrenal tumor characterization with in-phase and opposed-phase MR imaging: initial experience

PURPOSE: To determine the effect of two pairs of echo times (TEs) for in-phase (IP) and opposed-phase (OP) 3.0-T magnetic resonance (MR) imaging on (a) quantitative analysis prospectively in a phantom study and (b) diagnostic accuracy retrospectively in a clinical study of adrenal tumors, with use of various reference standards in the clinical study. MATERIALS AND METHODS: A fat-saline phantom was used to perform IP and OP 3.0-T MR imaging for various fat fractions. The institutional review board approved this HIPAA-compliant study, with waiver of informed consent. Single-breath-hold IP and OP 3.0-T MR images in 21 patients (14 women, seven men; mean age, 63 years) with 23 adrenal tumors (16 adenomas, six metastases, one adrenocortical carcinoma) were reviewed. The MR protocol involved two acquisition schemes: In scheme A, the first OP echo (approximately 1.5-msec TE) and the second IP echo (approximately 4.9-msec TE) were acquired. In scheme B, the first IP echo (approximately 2.4-msec TE) and the third OP echo (approximately 5.8-msec TE) were acquired. Quantitative analysis was performed, and analysis of variance was used to test for differences between adenomas and nonadenomas. RESULTS: In the phantom study, scheme B did not enable discrimination among voxels that had small amounts of fat. In the clinical study, no overlap in signal intensity (SI) index values between adenomas and nonadenomas was seen (P < .05) with scheme A. However, with scheme B, no overlap in the adrenal gland SI-to-liver SI ratio between adenomas and nonadenomas was seen (P < .05). With scheme B, no overlap in adrenal gland SI index-to-liver SI index ratio between adenomas and nonadenomas was seen (P < .05). CONCLUSION: This initial experience indicates SI index is the most reliable parameter for characterization of adrenal tumors with 3.0-T MR imaging when obtaining OP echo before IP echo. When acquiring IP echo before OP echo, however, nonadenomas can be mistaken as adenomas with use of the SI index value.

Quantitative stress echocardiography in coronary artery disease using contrast-based myocardial blood flow measurements: prospective comparison with coronary angiography

AIM: To test whether quantitative stress echocardiography using contrast-based myocardial blood flow (MBF, ml x min(-1) x g(-1)) measurements can detect coronary artery disease in humans. METHODS: 48 patients eligible for pharmacological stress testing by myocardial contrast echocardiography (MCE) and willing to undergo subsequent coronary angiography were prospectively enrolled in the study. Baseline and adenosine-induced (140 microg x kg(-1) x min(-1)) hyperaemic MBF was analysed according to a three-coronary-artery-territory model. Vascular territories were categorised into three groups with increasing stenosis severity defined as percentage diameter reduction by quantitative coronary angiography. RESULTS: Myocardial blood flow reserve (MBFR)-that is, the ratio of hyperaemic to baseline MBF, was obtained in 128 (89%) territories. Mean (SD) baseline MBF was 1.073 (0.395) ml x min(-1) x g(-1) and did not differ between territories supplied by coronary arteries with mild (<50% stenosis), moderate (50%-74% stenosis) or severe (>or=75% stenosis) disease. Mean (SD) hyperaemic MBF and MBFR were 2.509 (1.078) ml x min(-1) x g(-1) and 2.54 (1.03), respectively, and decreased linearly (r2 = 0.21 and r2 = 0.39) with stenosis severity. ROC analysis revealed that a territorial MBFR <1.94 detected >or=50% stenosis with 89% sensitivity and 92% specificity. CONCLUSION: Quantitative stress testing based on MBF measurements derived from contrast echocardiography is a new method for the non-invasive and reliable assessment of coronary artery disease in humans.

STATISTICAL METHODS FOR RARE AND COMMON VARIANT ASSOCIATION STUDIES

Complex human diseases are a major challenge for biological research. The goal of my research is to develop effective methods for biostatistics in order to create more opportunities for the prevention and cure of human diseases. This dissertation proposes statistical technologies that have the ability of being adapted to sequencing data in family-based designs, and that account for joint effects as well as gene-gene and gene-environment interactions in the GWA studies. The framework includes statistical methods for rare and common variant association studies. Although next-generation DNA sequencing technologies have made rare variant association studies feasible, the development of powerful statistical methods for rare variant association studies is still underway. Chapter 2 demonstrates two adaptive weighting methods for rare variant association studies based on family data for quantitative traits. The results show that both proposed methods are robust to population stratification, robust to the direction and magnitude of the effects of causal variants, and more powerful than the methods using weights suggested by Madsen and Browning [2009]. In Chapter 3, I extended the previously proposed test for Testing the effect of an Optimally Weighted combination of variants (TOW) [Sha et al., 2012] for unrelated individuals to TOW &ndash F, TOW for Family &ndash based design. Simulation results show that TOW &ndash F can control for population stratification in wide range of population structures including spatially structured populations, is robust to the directions of effect of causal variants, and is relatively robust to percentage of neutral variants. In GWA studies, this dissertation consists of a two &ndash locus joint effect analysis and a two-stage approach accounting for gene &ndash gene and gene &ndash environment interaction. Chapter 4 proposes a novel two &ndash stage approach, which is promising to identify joint effects, especially for monotonic models. The proposed approach outperforms a single &ndash marker method and a regular two &ndash stage analysis based on the two &ndash locus genotypic test. In Chapter 5, I proposed a gene &ndash based two &ndash stage approach to identify gene &ndash gene and gene &ndash environment interactions in GWA studies which can include rare variants. The two &ndash stage approach is applied to the GAW 17 dataset to identify the interaction between KDR gene and smoking status.

A new method to monitor visual field defects caused by photoreceptor degeneration by quantitative optical coherence tomography

PURPOSE: To correlate the dimension of the visual field (VF) tested by Goldman kinetic perimetry with the extent of visibility of the highly reflective layer between inner and outer segments of photoreceptors (IOS) seen in optical coherence tomography (OCT) images in patients with retinitis pigmentosa (RP). METHODS: In a retrospectively designed cross-sectional study, 18 eyes of 18 patients with RP were examined with OCT and Goldmann perimetry using test target I4e and compared with 18 eyes of 18 control subjects. A-scans of raw scan data of Stratus OCT images (Carl Zeiss Meditec, AG, Oberkochen, Germany) were quantitatively analyzed for the presence of the signal generated by the highly reflective layer between the IOS in OCT images. Starting in the fovea, the distance to which this signal was detectable was measured. Visual fields were analyzed by measuring the distance from the center point to isopter I4e. OCT and visual field data were analyzed in a clockwise fashion every 30 degrees , and corresponding measures were correlated. RESULTS: In corresponding alignments, the distance from the center point to isopter I4e and the distance to which the highly reflective signal from the IOS can be detected correlate significantly (r = 0.75, P < 0.0001). The greater the distance in VF, the greater the distance measured in OCT. CONCLUSIONS: The authors hypothesize that the retinal structure from which the highly reflective layer between the IOS emanates is of critical importance for visual and photoreceptor function. Further research is warranted to determine whether this may be useful as an objective marker of progression of retinal degeneration in patients with RP.

Quantitative OCT analysis of idiopathic perifoveal telangiectasia

PURPOSE: To identify and quantitate specific changes in optical coherence tomography (OCT) images of patients with type 2 idiopathic perifoveal telangiectasia (IPT). METHODS: In a prospectively designed, observational, case-control study, 28 eyes of 14 consecutive patients with IPT were examined with OCT and compared with eyes of 14 unaffected control subjects. Light reflectivity profiles of raw scan data of OCT images were quantitatively analyzed for differences in distance between different retinal reflectivity layers and their respective reflectivities. Maculae were examined in four separate regions: (1) central fovea, (2) nasal perifovea, (3) temporal perifovea, and (4) outside the fovea. RESULTS: Retinal thinning, shortening of the photoreceptor outer segments and loss of reflectivity of the photoreceptor ellipsoid region were found in the central foveal region as well as the nasal and temporal perifoveal regions in eyes with IPT. In addition, increased reflectivity of the outer nuclear layer was found in a sharply demarcated area of the inferotemporal perifoveal region in all affected eyes. Retinal tissue located more than 2000 mum away from the foveola was indistinguishable from that in normal eyes. CONCLUSIONS: Quantitative OCT analysis shows unique and specific changes in the photoreceptors of the central macula in IPT which can be detected from first clinical presentation. These changes may be of use as an additional diagnostic tool. Correlation of the findings in the outer nuclear layer with histologic studies may help identify the nature of the reflectivity increase and define more clearly the type of damage sustained by the photoreceptors in this condition.

Quantitative isolation of mouse and human intestinal intraepithelial lymphocytes by elutriation centrifugation

Intestinal intraepithelial lymphocytes (IEL) are specialized subsets of T cells with distinct functional capacities. While some IEL subsets are circulating, others such as CD8alphaalpha TCRalphabeta IEL are believed to represent non-circulating resident T cell subsets [Sim, G.K., Intraepithelial lymphocytes and the immune system. Adv. Immunol., 1995. 58: 297-343.]. Current methods to obtain enriched preparations of intraepithelial lymphocytes are mostly based on Percoll density gradient or magnetic bead-based technologies [Lundqvist, C., et al., Isolation of functionally active intraepithelial lymphocytes and enterocytes from human small and large intestine. J. Immunol. Methods, 1992. 152(2): 253-263.]. However, these techniques are hampered by a generally low yield of isolated cells, and potential artifacts due to the interference of the isolation procedure with subsequent functional assays, in particular, when antibodies against cell surface markers are required. Here we describe a new method for obtaining relatively pure populations of intestinal IEL (55-75%) at a high yield (>85%) by elutriation centrifugation. This technique is equally suited for the isolation and enrichment of intraepithelial lymphocytes of both mouse and human origin. Time requirements for fractionating cell suspensions by elutriation centrifugation are comparable to Percoll-, or MACS-based isolation procedures. Hence, the substantially higher yield and the consistent robust enrichment for intraepithelial lymphocytes, together with the gentle treatment of the cells during elutriation that does not interfere with subsequent functional assays, are important aspects that are in favor of using this elegant technology to obtain unmanipulated, unbiased populations of intestinal intraepithelial lymphocytes, and, if desired, also of pure epithelial cells.

Simultaneous quantitative detection of relevant biomarkers in breast cancer by quantitative real-time PCR

The assessment of ERa, PgR and HER2 status is routinely performed today to determine the endocrine responsiveness of breast cancer samples. Such determination is usually accomplished by means of immunohistochemistry and in case of HER2 amplification by means of fluorescent in situ hybridization (FISH). The analysis of these markers can be improved by simultaneous measurements using quantitative real-time PCR (Qrt-PCR). In this study we compared Qrt-PCR results for the assessment of mRNA levels of ERa, PgR, and the members of the human epidermal growth factor receptor family, HER1, HER2, HER3 and HER4. The results were obtained in two independent laboratories using two different methods, SYBR Green I and TaqMan probes, and different primers. By linear regression we demonstrated a good concordance for all six markers. The quantitative mRNA expression levels of ERa, PgR and HER2 also strongly correlated with the respective quantitative protein expression levels prospectively detected by EIA in both laboratories. In addition, HER2 mRNA expression levels correlated well with gene amplification detected by FISH in the same biopsies. Our results indicate that both Qrt-PCR methods were robust and sensitive tools for routine diagnostics and consistent with standard methodologies. The developed simultaneous assessment of several biomarkers is fast and labor effective and allows optimization of the clinical decision-making process in breast cancer tissue and/or core biopsies.

Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on 210Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for 226Ra, 210Pb and 137Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of 210Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the 137Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the 137Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments