Bayesian models for the analysis of genetic structure when populations are correlated

Motivation: Population allele frequencies are correlated when populations have a shared history or when they exchange genes. Unfortunately, most models for allele frequency and inference about population structure ignore this correlation. Recent analytical results show that among populations, correlations can be very high, which could affect estimates of population genetic structure. In this study, we propose a mixture beta model to characterize the allele frequency distribution among populations. This formulation incorporates the correlation among populations as well as extending the model to data with different clusters of populations. Results: Using simulated data, we show that in general, the mixture model provides a good approximation of the among-population allele frequency distribution and a good estimate of correlation among populations. Results from fitting the mixture model to a dataset of genotypes at 377 autosomal microsatellite loci from human populations indicate high correlation among populations, which may not be appropriate to neglect. Traditional measures of population structure tend to over-estimate the amount of genetic differentiation when correlation is neglected. Inference is performed in a Bayesian framework.

Molecular interactions and signal transfer between sensory rhodopsin-I and its cognate transducer

SRI is unique among known photoreceptors in that it produces opposite signals depending on the color of light stimuli. Absorption of orange light (587 nm) triggers an attractant response by the cell, whereas absorption of orange light followed by near-UV light (373 run) triggers a repellent response. Using behavioral mutants that exhibit aberrant color-sensing ability, we tested a two-conformation equilibrium model, using FRET and EPR spectroscopy. The essence of the model applied to SRI-HtrI is that the complex exists in a metastable two-conformer equilibrium which is shifted in one direction by orange light absorption (producing an attractant signal) and in the opposite direction by a second UV-violet photon (producing a repellent signal). First, by FRET we found that the E-F cytoplasmic loop of SRI moves toward the RAMP domain of the HtrI transducer during the formation of the orange-light activated signaling state of the complex. This is the first localization of a change in the physical relationship between the receptor and transducer subunits of the complex and provides a structural property of the two proposed conformers that we can monitor. Second, EPR spectra of a spin label probe at this cytoplasmic position showed shifts in the dark in the mutants toward shorter or longer EF loop-RAMP distances, explaining their behavior in terms of their mutations causing pre-stimulus shifts into one or the other conformer. ^ Next, we applied a novel electrophysiological method for monitoring the directionality of proton movement during photoactivation of SRI, to investigate the process of proton transfer in the photoactive site from the chromophore to proton acceptors on both the wildtype and aberrant color-response mutants. We observed an unexpected and critical difference in the two signaling conformations of the SRI-HtrI complex. The finding is that the vectoriality (i.e. movement away or toward the cytoplasm) of the light-induced proton transfer from the chromophore to the protein is opposite in formation of the two conformations. Retinylidene proton transfer is a common critical process in rhodopsins and these results are the first to show differences in vectoriality in a rhodopsin receptor, and to demonstrate functional importance of the direction of proton transfer. ^

Disparities in medical expenditure and utilization among hypertensive men and women in the U.S.: Cross-section and lifetime analysis

Objectives. To investigate procedural gender equity by assessing predisposing, enabling and need predictors of gender differences in annual medical expenditures and utilization among hypertensive individuals in the U.S. Also, to estimate and compare lifetime medical expenditures among hypertensive men and women in the U.S. ^ Data source. 2001-2004 the Medical Expenditure Panel Survey (MEPS);1986-2000 National Health Interview Survey (NHIS) and National Health Interview Survey linked to mortality in the National Death Index through 2002 (2002 NHIS-NDI). ^ Study design. We estimated total medical expenditure using four equations regression model, specific medical expenditures using two equations regression model and utilization using negative binomial regression model. Procedural equity was assessed by applying the Aday et al. theoretical framework. Expenditures were estimated in 2004 dollars. We estimated hypertension-attributable medical expenditure and utilization among men and women. ^ To estimate lifetime expenditures from ages 20 to 85+, we estimated medical expenditures with cross-sectional data and survival with prospective data. The four equations regression model were used to estimate average annual medical expenditures defined as sum of inpatient stay, emergency room visits, outpatient visits, office based visits, and prescription drugs expenditures. Life tables were used to estimate the distribution of life time medical expenditures for hypertensive men and women at different age and factors such as disease incidence, medical technology and health care cost were assumed to be fixed. Both total and hypertension attributable expenditures among men and women were estimated. ^ Data collection. We used the 2001-2004 MEPS household component and medical condition files; the NHIS person and condition files from 1986-1996 and 1997-2000 sample adult files were used; and the 1986-2000 NHIS that were linked to mortality in the 2002 NHIS-NDI. ^ Principal findings. Hypertensive men had significantly less utilization for most measures after controlling predisposing, enabling and need factors than hypertensive women. Similarly, hypertensive men had less prescription drug (-9.3%), office based (-7.2%) and total medical (-4.5%) expenditures than hypertensive women. However, men had more hypertension-attributable medical expenditures and utilization than women. ^ Expected total lifetime expenditure for average life table individuals at age 20, was $188,300 for hypertensive men and $254,910 for hypertensive women. But the lifetime expenditure that could be attributed to hypertension was $88,033 for men and $40,960 for women. ^ Conclusion. Hypertensive women had more utilization and expenditure for most measures than hypertensive men, possibly indicating procedural inequity. However, relatively higher hypertension-attributable health care of men shows more utilization of resources to treat hypertension related diseases among men than women. Similar results were reported in lifetime analyses.^ Key words: gender, medical expenditures, utilization, hypertension-attributable, lifetime expenditure ^

Network-based genome-wide association study of Crohn's disease

Genome-wide association studies (GWAS) have rapidly become a standard method for disease gene discovery. Many recent GWAS indicate that for most disorders, only a few common variants are implicated and the associated SNPs explain only a small fraction of the genetic risk. The current study incorporated gene network information into gene-based analysis of GWAS data for Crohn's disease (CD). The purpose was to develop statistical models to boost the power of identifying disease-associated genes and gene subnetworks by maximizing the use of existing biological knowledge from multiple sources. The results revealed that Markov random field (MRF) based mixture model incorporating direct neighborhood information from a single gene network is not efficient in identifying CD-related genes based on the GWAS data. The incorporation of solely direct neighborhood information might lead to the low efficiency of these models. Alternative MRF models looking beyond direct neighboring information are necessary to be developed in the future for the purpose of this study.^

Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Network-based analysis of genome-wide scans of natural selection

The genomic era brought by recent advances in the next-generation sequencing technology makes the genome-wide scans of natural selection a reality. Currently, almost all the statistical tests and analytical methods for identifying genes under selection was performed on the individual gene basis. Although these methods have the power of identifying gene subject to strong selection, they have limited power in discovering genes targeted by moderate or weak selection forces, which are crucial for understanding the molecular mechanisms of complex phenotypes and diseases. Recent availability and rapid completeness of many gene network and protein-protein interaction databases accompanying the genomic era open the avenues of exploring the possibility of enhancing the power of discovering genes under natural selection. The aim of the thesis is to explore and develop normal mixture model based methods for leveraging gene network information to enhance the power of natural selection target gene discovery. The results show that the developed statistical method, which combines the posterior log odds of the standard normal mixture model and the Guilt-By-Association score of the gene network in a naïve Bayes framework, has the power to discover moderate/weak selection gene which bridges the genes under strong selection and it helps our understanding the biology under complex diseases and related natural selection phenotypes.^

Los públicos en las redes sociales : nuevas prácticas

Las redes sociales parecen admitir un modelo de comunicación bidireccional que permitiría la interacción y conversación entre los partícipes de los nuevos procesos de comunicación. Este escenario novedoso ha motivado a indagar en los términos que usualmente se utilizaban para definir a los actores del modelo de comunicación unidireccional que caracterizaba a los medios de comunicación tradicional. Así, los usos de los términos de público y audiencia necesitan de su revisión al plantearse interrogantes sobre las formas precisas de nombrar a los actores en los distintos comportamientos intercambiables que llevan adelante en los nuevos procesos de comunicación que las redes sociales vehiculizan. En primer lugar este artículo examina los aspectos novedosos que las redes sociales admiten, para luego indagar en los conceptos de audiencia y público con el fin de conocer las nuevas acepciones que tienen ambos términos en los procesos comunicacionales actuales.

Isotopic geochemistry of granitoids in the Jinshajiang suture zone, SW China

The Jinshajiang suture zone, located in the eastern part of the Tethyan tectonic domain, is noticeable for a large-scale distribution of Late Jurassic to Triassic granitoids. These granitoids were genetically related to the evolution of the Paleo-Tethys Ocean. The Beiwu, Linong and Lunong granitoids occur in the middle zone of the Jinshajiang Suture Zone, and possess similar geochemical features, indicating they share a common magma source. SIMS zircon U-Pb dating reveals the Beiwu, Linong and Lunong granitic intrusions were emplaced at 233.9±1.4 Ma (2 sigma), 233.1 ±1.4 Ma (2 sigma) and 231.0±1.6 Ma (2 sigma), respectively. All of these granitoids are enriched in abundances of Si (SiO2 =65.2-73.5 wt.%), and large-ion-lithophile-elements (LILEs), but depleted in high-field-strength-elements contents (HFSEs, e.g., Nb, Ta, Ti). In addition, they have low P2O5 contents (0.06-0.11 wt.%), A/CNK values ([molecular Al2O3/(CaO+Na2O+K2O)], mostly<1.1) and 10000Ga/Al ratios (1.7-2.2), consistent with the characteristics of I-type granites. In terms of isotopic compositions, these granitoids have high initial 87Sr/86Sr ratios (0.7078-0.7148), Pb isotopic compositions [(206Pb/204Pb)t=18.213-18.598, (207Pb/204Pb)t=15.637-15.730 and (208Pb/204Pb)t=38.323-38.791], zircon d18O values (7. per mil-9.3 per mil) and negative eNd(t) values (-5.1 to -6.7), suggesting they were predominantly derived from the continental crust. Their Nb/Ta ratios (average value=8.6) are consistent with those of the lower continental crust (LCC). However, variable ?Hf(t) values (-8.6 to +2.8) and the occurrences of mafic microgranular enclaves (MMEs) suggest that mantle-derived melts and lower crustal magmas were involved in the generation of these granitoids. Moreover, the high Pb isotopic ratios and elevated zircon d18O values of these rocks indicate a significant contribution of the upper crustal composition. We propose a model in which the Beiwu, Linong and Lunong granitoids were generated under a late collisional or post-collisional setting. It is possible that this collision was completed before Late Triassic. Decompression induced mantle-derived magmas underplated and provided the heat for the anatexis of the crust. Hybrid melts including mantle-derived and the lower crustal magmas were then generated. The hybrid melts thereafter ascended to a shallow depth and resulted in some degree of sedimentary rocks assimilation. Such three-component mixing magmas source and subsequent fractional crystallization could be responsible for the formation of the Beiwu, Linong and Lunong granitoids.

Lithogeochemical and isotope (Rb-Sr, Sm-Nd, Pb-Pb whole rock and Lu-Hf zircon) composition of samples from the Shackleton Range, East Antarctica

Three distinct, spatially separated crustal terranes have been recognised in the Shackleton Range, East Antarctica: the Southern, Eastern and Northern Terranes. Mafic gneisses from the Southern Terrane provide geochemical evidence for a within-plate, probably back-arc origin of their protoliths. A plume-distal ridge origin in an incipient ocean basin is the favoured interpretation for the emplacement site of these rocks at c. 1850 Ma, which, together with a few ocean island basalts, were subsequently incorporated into an accretionary continental arc/supra-subduction zone tectonic setting. Magmatic underplating resulted in partial melting of the lower crust, which caused high-temperature granulite-facies metamorphism in the Southern Terrane at c. 1710-1680 Ma. Mafic and felsic gneisses there are characterised by isotopically depleted, positive Nd and Hf initials and model ages between 2100 and 2000 Ma. They may be explained as juvenile additions to the crust towards the end of the Palaeoproterozoic. These juvenile rocks occur in a narrow, c. 150 km long E-W trending belt, inferred to trace a suture that is associated with a large Palaeoproterozoic accretionary orogenic system. The Southern Terrane contains many features that are similar to the Australo-Antarctic Mawson Continent and may be its furthermost extension into East Antarctica. The Eastern Terrane is characterised by metagranitoids that formed in a continental volcanic arc setting during a late Mesoproterozoic orogeny at c. 1060 Ma. Subsequently, the rocks experienced high-temperature metamorphism during Pan-African collisional tectonics at 600 Ma. Isotopically depleted zircon grains yielded Hf model ages of 1600-1400 Ma, which are identical to Nd model ages obtained from juvenile metagranitoids. Most likely, these rocks trace the suture related to the amalgamation of the Indo-Antarctic and West Gondwana continental blocks at ~600 Ma. The Eastern Terrane is interpreted as the southernmost extension of the Pan-African Mozambique/Maud Belt in East Antarctica and, based on Hf isotope data, may also represent a link to the Ellsworth-Whitmore Mountains block in West Antarctica and the Namaqua-Natal Province of southern Africa. Geochemical evidence indicates that the majority of the protoliths of the mafic gneisses in the Northern Terrane formed as oceanic island basalts in a within-plate setting. Subsequently the rocks were incorporated into a subduction zone environment and, finally, accreted to a continental margin during Pan-African collisional tectonics. Felsic gneisses there provide evidence for a within-plate and volcanic arc/collisional origin. Emplacement of granitoids occurred at c. 530 Ma and high-temperature, high-pressure metamorphism took place at 510-500 Ma. Enriched Hf and Nd initials and Palaeoproterozoic model ages for most samples indicate that no juvenile material was added to the crust of the Northern Terrane during the Pan-African Orogeny but recycling of older crust or mixing of crustal components of different age must have occurred. Isotopically depleted mafic gneisses, which are spatially associated with eclogite-facies pyroxenites, yielded late Mesoproterozoic Nd model ages. These rocks occur in a narrow, at least 100 km long, E-W trending belt that separates alkaline ocean island metabasalts and within-plate metagranitoids from volcanic arc metabasalts and volcanic arc/syn-collisional metagranitoids in the Northern Terrane. This belt is interpreted to trace the late Neoproterozoic/early Cambrian Pan-African collisional suture between the Australo-Antarctic and the combined Indo-Antarctic/West Gondwana continental blocks that formed during the final amalgamation of Gondwana.

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.

Methane, helium isotope and hydrogen data from the Nibelungen hydrothermal area, Mid Atlantic Ridge

Hydrothermal emission of mantle helium appears to be directly related to magma production rate, but other processes can generate methane and hydrogen on mid-ocean ridges. In an on-going effort to characterize these processes in the South Atlantic, the flux and distribution of these gases were investigated in the vicinity of a powerful black smoker recently discovered at 8°17.9' S, 13°30.4' W. The vent lies on the shoulder of an oblique offset in the Mid-Atlantic Ridge and discharges high concentrations of methane and hydrogen. Measurements during expeditions in 2004 and 2006 show that the ratio of CH4 to 3He in the neutrally buoyant plume is quite high, 4 x 10**8. The CTD stations were accompanied by velocity measurements with lowered acoustic Doppler current profilers (LADCP), and from these data we estimate the methane transport to have been 0.5 mol/sec in a WSW-trending plume that seems to develop during the ebb tidal phase. This transport is an order of magnitude greater than the source of CH4 calculated from its concentration in the vent fluid and the rise height of the plume. From this range of methane fluxes, the source of 3He is estimated to be between 0.14 and 1.2 nmol/sec. In either case, the 3He source is significantly lower than expected from the spreading rate of the Mid-Atlantic Ridge. From the inventory of methane in the rift valley adjacent to the vent, it appears that the average specific rate of oxidation is 2.6 to 23/yr, corresponding to a turnover time between 140 and 16 days. Vertical profiles of methane in the surrounding region often exhibited Gaussian-like distributions, and the variances appear to increase with distance from the vent. Using a Gaussian plume model, we obtained a range of vertical eddy diffusivities between 0.009 and 0.08 m2m2/sec. These high values may be due to tidally driven internal waves across the promontory on which the vent is located.

Volcanic rocks from the submarine Vityaz Ridge: their age, chemical and isotopic compositions

Results of geological studies at the submarine Vityaz Ridge carried out during cruises 37 and 41 of R/V Akademik Lavrent'ev in 2005 and 2006 are reported. The studied area is located at an near-island trench of the slope in the central part of the Kuril Island arc. Morphologically it consists of two parts: an inner volcanic arc represented by the Great Kuril Range and an outer arc corresponding to the submarine Vityaz Ridge. Diverse rocks composing the basement and the sedimentary cover of the ridge were recovered by dredging. Based on K-Ar dating and geochemistry, volcanics were divided into Paleocene, Eocene, late Oligocene, and Pliocene-Pleistocene complexes. Each of the complexes reflects a tectonomagmatic stage in the ridge evolution. Geochemical and isotope data on the volcanics indicate contribution of ancient crustal material in the magma source and, correspondingly, formation of this structure on the continental basement. Two-stage model ages (TDM2) vary in a wide range from zero values in mafic rocks to 0.77 Ga in felsic varieties, pointing to presence of Precambrian protolith in the source of the felsic rocks of the Vityaz Ridge. The Pliocene-Pleistocene volcanics are classed with tholeiitic, calc-alkaline, and subalkaline series, which differ in alkali contents and REE fractionation. Values of (La/Sm)_n and (La/Yb)_n ratios vary from 0.74 and 0.84 in the tholeiitic varieties to 1.19 and 1.44 in the calc-alkaline and 2.32 and 3.73 in the subalkaline rocks. All three varieties occur within the same volcanic edifices and formed during differentiation of magmatic melts that were channeled along fault zones from the mantle source slightly enriched in crustal component.

Ba, Mn, and Mg record derived from Hoeglundina elegans of sediment core SLA-9

The evolution of productivity, redox conditions, temperature, and ventilation during the deposition of an Aegean sapropel (S1) is independently constrained using bulk sediment composition and high-resolution single specimen benthic foraminiferal trace metal and stable isotope data. The occurrence of benthic foraminifer, Hoeglundina elegans (H. elegans), through a shallow water (260 m) sapropel, permits for the first time a comparison between dissolved and particulate concentrations of Ba and Mn and the construction of a Mg/Ca-based temperature record through sapropel S1. The simultaneous increase in sedimentary Ba and incorporated Ba in foraminiferal test carbonate, (Ba/Ca)H. elegans, points to a close coupling between Ba cycling and export productivity. During sapropel deposition, sedimentary Mn content ((Mn/Al)sed) is reduced, corresponding to enhanced Mn2+ mobilization from sedimentary Mn oxides under suboxic conditions. The consequently elevated dissolved Mn2+ concentrations are reflected in enhanced (Mn/Ca)H. elegans levels. The magnitude and duration of the sapropel interruption and other short-term cooling events are constrained using Mg/Ca thermometry. Based on integrating productivity and ventilation records with the temperature record, we propose a two-mode hysteresis model for sapropel formation.