A comparison of three methods of Mendelian randomization when the genetic instrument, the risk factor and the outcome are all binary.

The method of instrumental variable (referred to as Mendelian randomization when the instrument is a genetic variant) has been initially developed to infer on a causal effect of a risk factor on some outcome of interest in a linear model. Adapting this method to nonlinear models, however, is known to be problematic. In this paper, we consider the simple case when the genetic instrument, the risk factor, and the outcome are all binary. We compare via simulations the usual two-stages estimate of a causal odds-ratio and its adjusted version with a recently proposed estimate in the context of a clinical trial with noncompliance. In contrast to the former two, we confirm that the latter is (under some conditions) a valid estimate of a causal odds-ratio defined in the subpopulation of compliers, and we propose its use in the context of Mendelian randomization. By analogy with a clinical trial with noncompliance, compliers are those individuals for whom the presence/absence of the risk factor X is determined by the presence/absence of the genetic variant Z (i.e., for whom we would observe X = Z whatever the alleles randomly received at conception). We also recall and illustrate the huge variability of instrumental variable estimates when the instrument is weak (i.e., with a low percentage of compliers, as is typically the case with genetic instruments for which this proportion is frequently smaller than 10%) where the inter-quartile range of our simulated estimates was up to 18 times higher compared to a conventional (e.g., intention-to-treat) approach. We thus conclude that the need to find stronger instruments is probably as important as the need to develop a methodology allowing to consistently estimate a causal odds-ratio.

Earliest known use of marine resources by neanderthals

Numerous studies along the northern Mediterranean borderland have documented the use of shellfish by Neanderthals but none of these finds are prior to Marine Isotopic Stage 3 (MIS 3). In this paper we present evidence that gathering and consumption of mollusks can now be traced back to the lowest level of the archaeological sequence at Bajondillo Cave (Málaga, Spain), dated during the MIS 6. The paper describes the taxonomical and taphonomical features of the mollusk assemblages from this level Bj19 and briefly touches upon those retrieved in levels Bj18 (MIS 5) and Bj17 (MIS 4), evidencing a continuity of the shellfishing activity that reaches to MIS 3. This evidence is substantiated on 29 datings through radiocarbon, thermoluminescence and U series methods. Obtained dates and paleoenvironmental records from the cave include isotopic, pollen, lithostratigraphic and sedimentological analyses and they are fully coherent with paleoclimate conditions expected for the different stages. We conclude that described use of shellfish resources by Neanderthals (H. neanderthalensis) in Southern Spain started ~150 ka and were almost contemporaneous to Pinnacle Point (South Africa), when shellfishing is first documented in archaic modern humans.

A concerted kinase interplay identifies PPARgamma as a molecular target of ghrelin signaling in macrophages.

The peroxisome proliferator-activator receptor PPARgamma plays an essential role in vascular biology, modulating macrophage function and atherosclerosis progression. Recently, we have described the beneficial effect of combined activation of the ghrelin/GHS-R1a receptor and the scavenger receptor CD36 to induce macrophage cholesterol release through transcriptional activation of PPARgamma. Although the interplay between CD36 and PPARgamma in atherogenesis is well recognized, the contribution of the ghrelin receptor to regulate PPARgamma remains unknown. Here, we demonstrate that ghrelin triggers PPARgamma activation through a concerted signaling cascade involving Erk1/2 and Akt kinases, resulting in enhanced expression of downstream effectors LXRalpha and ABC sterol transporters in human macrophages. These effects were associated with enhanced PPARgamma phosphorylation independently of the inhibitory conserved serine-84. Src tyrosine kinase Fyn was identified as being recruited to GHS-R1a in response to ghrelin, but failure of activated Fyn to enhance PPARgamma Ser-84 specific phosphorylation relied on the concomitant recruitment of docking protein Dok-1, which prevented optimal activation of the Erk1/2 pathway. Also, substitution of Ser-84 preserved the ghrelin-induced PPARgamma activity and responsiveness to Src inhibition, supporting a mechanism independent of Ser-84 in PPARgamma response to ghrelin. Consistent with this, we found that ghrelin promoted the PI3-K/Akt pathway in a Galphaq-dependent manner, resulting in Akt recruitment to PPARgamma, enhanced PPARgamma phosphorylation and activation independently of Ser-84, and increased expression of LXRalpha and ABCA1/G1. Collectively, these results illustrate a complex interplay involving Fyn/Dok-1/Erk and Galphaq/PI3-K/Akt pathways to transduce in a concerted manner responsiveness of PPARgamma to ghrelin in macrophages.

Mynodbcsv: Lightweight Zero-Config Database Solution for Handling Very Large CSV Files.

Volumes of data used in science and industry are growing rapidly. When researchers face the challenge of analyzing them, their format is often the first obstacle. Lack of standardized ways of exploring different data layouts requires an effort each time to solve the problem from scratch. Possibility to access data in a rich, uniform manner, e.g. using Structured Query Language (SQL) would offer expressiveness and user-friendliness. Comma-separated values (CSV) are one of the most common data storage formats. Despite its simplicity, with growing file size handling it becomes non-trivial. Importing CSVs into existing databases is time-consuming and troublesome, or even impossible if its horizontal dimension reaches thousands of columns. Most databases are optimized for handling large number of rows rather than columns, therefore, performance for datasets with non-typical layouts is often unacceptable. Other challenges include schema creation, updates and repeated data imports. To address the above-mentioned problems, I present a system for accessing very large CSV-based datasets by means of SQL. It's characterized by: "no copy" approach - data stay mostly in the CSV files; "zero configuration" - no need to specify database schema; written in C++, with boost [1], SQLite [2] and Qt [3], doesn't require installation and has very small size; query rewriting, dynamic creation of indices for appropriate columns and static data retrieval directly from CSV files ensure efficient plan execution; effortless support for millions of columns; due to per-value typing, using mixed text/numbers data is easy; very simple network protocol provides efficient interface for MATLAB and reduces implementation time for other languages. The software is available as freeware along with educational videos on its website [4]. It doesn't need any prerequisites to run, as all of the libraries are included in the distribution package. I test it against existing database solutions using a battery of benchmarks and discuss the results.

Autoimmunity-Associated LYP-W620 Does Not Impair Thymic Negative Selection of Autoreactive T Cells.

A C1858T (R620W) variation in the PTPN22 gene encoding the tyrosine phosphatase LYP is a major risk factor for human autoimmunity. LYP is a known negative regulator of signaling through the T cell receptor (TCR), and murine Ptpn22 plays a role in thymic selection. However, the mechanism of action of the R620W variant in autoimmunity remains unclear. One model holds that LYP-W620 is a gain-of-function phosphatase that causes alterations in thymic negative selection and/or thymic output of regulatory T cells (Treg) through inhibition of thymic TCR signaling. To test this model, we generated mice in which the human LYP-W620 variant or its phosphatase-inactive mutant are expressed in developing thymocytes under control of the proximal Lck promoter. We found that LYP-W620 expression results in diminished thymocyte TCR signaling, thus modeling a "gain-of-function" of LYP at the signaling level. However, LYP-W620 transgenic mice display no alterations of thymic negative selection and no anomalies in thymic output of CD4(+)Foxp3(+) Treg were detected in these mice. Lck promoter-directed expression of the human transgene also causes no alteration in thymic repertoire or increase in disease severity in a model of rheumatoid arthritis, which depends on skewed thymic selection of CD4(+) T cells. Our data suggest that a gain-of-function of LYP is unlikely to increase risk of autoimmunity through alterations of thymic selection and that LYP likely acts in the periphery perhaps selectively in regulatory T cells or in another cell type to increase risk of autoimmunity.

High prevalence and lineage diversity of avian malaria in wild populations of great tits (Parus major) and mosquitoes (Culex pipiens).

Avian malaria studies have taken a prominent place in different aspects of evolutionary ecology. Despite a recent interest in the role of vectors within the complex interaction system of the malaria parasite, they have largely been ignored in most epidemiological studies. Epidemiology of the disease is however strongly related to the vector's ecology and behaviour, and there is a need for basic investigations to obtain a better picture of the natural associations between Plasmodium lineages, vector species and bird hosts. The aim of the present study was to identify the mosquito species involved in the transmission of the haemosporidian parasites Plasmodium spp. in two wild populations of breeding great tits (Parus major) in western Switzerland. Additionally, we compared Plasmodium lineages, based on mitochondrial DNA cytochrome b sequences, between the vertebrate and dipteran hosts, and evaluated the prevalence of the parasite in the mosquito populations. Plasmodium spp. were detected in Culex pipiens only, with an overall 6.6% prevalence. Among the six cytochrome b lineages of Plasmodium identified in the mosquitoes, three were also present in great tits. The results provide evidence for the first time that C. pipiens can act as a natural vector of avian malaria in Europe and yield baseline data for future research on the epidemiology of avian malaria in European countries.

Genome-wide association and functional follow-up reveals new loci for kidney function.

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.

Estimating the confidence level of white matter connections obtained with MRI tractography.

BACKGROUND: Since the emergence of diffusion tensor imaging, a lot of work has been done to better understand the properties of diffusion MRI tractography. However, the validation of the reconstructed fiber connections remains problematic in many respects. For example, it is difficult to assess whether a connection is the result of the diffusion coherence contrast itself or the simple result of other uncontrolled parameters like for example: noise, brain geometry and algorithmic characteristics. METHODOLOGY/PRINCIPAL FINDINGS: In this work, we propose a method to estimate the respective contributions of diffusion coherence versus other effects to a tractography result by comparing data sets with and without diffusion coherence contrast. We use this methodology to assign a confidence level to every gray matter to gray matter connection and add this new information directly in the connectivity matrix. CONCLUSIONS/SIGNIFICANCE: Our results demonstrate that whereas we can have a strong confidence in mid- and long-range connections obtained by a tractography experiment, it is difficult to distinguish between short connections traced due to diffusion coherence contrast from those produced by chance due to the other uncontrolled factors of the tractography methodology.

Impacts on the deep-sea ecosystem by a severe coastal storm

Major coastal storms, associated with strong winds, high waves and intensified currents, and occasionally with heavy rains and flash floods, are mostly known because of the serious damage they can cause along the shoreline and the threats they pose to navigation. However, there is a profound lack of knowledge on the deep-sea impacts of severe coastal storms. Concurrent measurements of key parameters along the coast and in the deep-sea are extremely rare. Here we present a unique data set showing how one of the most extreme coastal storms of the last decades lashing the Western Mediterranean Sea rapidly impacted the deep-sea ecosystem. The storm peaked the 26th of December 2008 leading to the remobilization of a shallow-water reservoir of marine organic carbon associated with fine particles and resulting in its redistribution across the deep basin. The storm also initiated the movement of large amounts of coarse shelf sediment, which abraded and buried benthic communities. Our findings demonstrate, first, that severe coastal storms are highly efficient in transporting organic carbon from shallow water to deep water, thus contributing to its sequestration and, second, that natural, intermittent atmospheric drivers sensitive to global climate change have the potential to tremendously impact the largest and least known ecosystem on Earth, the deep-sea ecosystem.

Sea Urchins Predation Facilitates Coral Invasion in a Marine Reserve

Macroalgae is the dominant trophic group on Mediterranean infralittoral rocky bottoms, whereas zooxanthellate corals are extremely rare. However, in recent years, the invasive coral Oculina patagonica appears to be increasing its abundance through unknown means. Here we examine the pattern of variation of this species at a marine reserve between 2002 and 2010 and contribute to the understanding of the mechanisms that allow its current increase. Because indirect interactions between species can play a relevant role in the establishment of species, a parallel assessment of the sea urchin Paracentrotus lividus, the main herbivorous invertebrate in this habitat and thus a key species, was conducted. O. patagonica has shown a 3-fold increase in abundance over the last 8 years and has become the most abundant invertebrate in the shallow waters of the marine reserve, matching some dominant erect macroalgae in abundance. High recruitment played an important role in this increasing coral abundance. The results from this study provide compelling evidence that the increase in sea urchin abundance may be one of the main drivers of the observed increase in coral abundance. Sea urchins overgraze macroalgae and create barren patches in the space-limited macroalgal community that subsequently facilitate coral recruitment. This study indicates that trophic interactions contributed to the success of an invasive coral in the Mediterranean because sea urchins grazing activity indirectly facilitated expansion of the coral. Current coral abundance at the marine reserve has ended the monopolization of algae in rocky infralittoral assemblages, an event that could greatly modify both the underwater seascape and the sources of primary production in the ecosystem.

The structure of Mediterranean rocky reef ecosystems across environmental and human gradients, and conservation implications

Historical exploitation of the Mediterranean Sea and the absence of rigorous baselines makes it difficult to evaluate the current health of the marine ecosystems and the efficacy of conservation actions at the ecosystem level. Here we establish the first current baseline and gradient of ecosystem structure of nearshore rocky reefs at the Mediterranean scale. We conducted underwater surveys in 14 marine protected areas and 18 open access sites across the Mediterranean, and across a 31-fold range of fish biomass (from 3.8 to 118 g m22). Our data showed remarkable variation in the structure of rocky reef ecosystems. Multivariate analysis showed three alternative community states: (1) large fish biomass and reefs dominated by non-canopy algae, (2) lower fish biomass but abundant native algal canopies and suspension feeders, and (3) low fish biomass and extensive barrens, with areas covered by turf algae. Our results suggest that the healthiest shallow rocky reef ecosystems in the Mediterranean have both large fish and algal biomass. Protection level and primary production were the only variables significantly correlated to community biomass structure. Fish biomass was significantly larger in well-enforced no-take marine reserves, but there were no significant differences between multi-use marine protected areas (which allow some fishing) and open access areas at the regional scale. The gradients reported here represent a trajectory of degradation that can be used to assess the health of any similar habitat in the Mediterranean, and to evaluate the efficacy of marine protected areas.

Multiple processes regulate long-term population dynamics of sea urchins on Mediterranean rocky reefs

We annually monitored the abundance and size structure of herbivorous sea urchin populations (Paracentrotus lividus and Arbacia lixula) inside and outside a marine reserve in the Northwestern Mediterranean on two distinct habitats (boulders and vertical walls) over a period of 20 years, with the aim of analyzing changes at different temporal scales in relation to biotic and abiotic drivers. P. lividus exhibited significant variability in density over time on boulder bottoms but not on vertical walls, and temporal trends were not significantly different between the protection levels. Differences in densities were caused primarily by variance in recruitment, which was less pronounced inside the MPA and was correlated with adult density, indicating density-dependent recruitment under high predation pressure, as well as some positive feedback mechanisms that may facilitate higher urchin abundances despite higher predator abundance. Populations within the reserve were less variable in abundance and did not exhibit the hyper-abundances observed outside the reserve, suggesting that predation effects maybe more subtle than simply lowering the numbers of urchins in reserves. A. lixula densities were an order of magnitude lower than P. lividus densities and varied within sites and over time on boulder bottoms but did not differ between protection levels. In December 2008, an exceptionally violent storm reduced sea urchin densities drastically (by 50% to 80%) on boulder substrates, resulting in the lowest values observed over the entire study period, which remained at that level for at least two years (up to the present). Our results also showed great variability in the biological and physical processes acting at different temporal scales. This study highlights the need for appropriate temporal scales for studies to fully understand ecosystem functioning, the concepts of which are fundamental to successful conservation and management.

Inferring epidemic contact structure from phylogenetic trees.

Contact structure is believed to have a large impact on epidemic spreading and consequently using networks to model such contact structure continues to gain interest in epidemiology. However, detailed knowledge of the exact contact structure underlying real epidemics is limited. Here we address the question whether the structure of the contact network leaves a detectable genetic fingerprint in the pathogen population. To this end we compare phylogenies generated by disease outbreaks in simulated populations with different types of contact networks. We find that the shape of these phylogenies strongly depends on contact structure. In particular, measures of tree imbalance allow us to quantify to what extent the contact structure underlying an epidemic deviates from a null model contact network and illustrate this in the case of random mixing. Using a phylogeny from the Swiss HIV epidemic, we show that this epidemic has a significantly more unbalanced tree than would be expected from random mixing.

Alterations in Postural Control during the World's Most Challenging Mountain Ultra-Marathon

We investigated postural control (PC) effects of a mountain ultra-marathon (MUM): a 330-km trail run with 24000 m of positive and negative change in elevation. PC was assessed prior to (PRE), during (MID) and after (POST) the MUM in experienced ultra-marathon runners (n = 18; finish time = 126+/-16 h) and in a control group (n = 8) with a similar level of sleep deprivation. Subjects were instructed to stand upright on a posturographic platform over a period of 51.2 seconds using a double-leg stance under two test conditions: eyes open (EO) and eyes closed (EC). Traditional measures of postural stability (center of pressure trajectory analysis) and stabilogram-diffusion analysis (SDA) parameters were analysed. For the SDA, a significantly greater short-term effective diffusion was found at POST compared with PRE in the medio-lateral (ML; Dxs) and antero-posterior (AP) directions (Dys) in runners (p<0.05) The critical time interval (Ctx) in the ML direction was significantly higher at MID (p<0.001) and POST (p<0.05) than at PRE in runners. At MID (p<0.001) and POST (p<0.05), there was a significant difference between the two groups. The critical displacement (Cdx) in the ML was significantly higher at MID and at POST (p<0.001) compared with PRE for runners. A significant difference in Cdx was observed between groups in EO at MID (p<0.05) and POST (p<0.005) in the ML direction and in EC at POST in the ML and AP directions (p<0.05). Our findings revealed significant effects of fatigue on PC in runners, including, a significant increase in Ctx (critical time in ML plan) in EO and EC conditions. Thus, runners take longer to stabilise their body at POST than at MID. It is likely that the mountainous characteristics of MUM (unstable ground, primarily uphill/downhill running, and altitude) increase this fatigue, leading to difficulty in maintaining balance.

Using transcription modules to identify expression clusters perturbed in Williams-Beuren syndrome.

The genetic dissection of the phenotypes associated with Williams-Beuren Syndrome (WBS) is advancing thanks to the study of individuals carrying typical or atypical structural rearrangements, as well as in vitro and animal studies. However, little is known about the global dysregulations caused by the WBS deletion. We profiled the transcriptomes of skin fibroblasts from WBS patients and compared them to matched controls. We identified 868 differentially expressed genes that were significantly enriched in extracellular matrix genes, major histocompatibility complex (MHC) genes, as well as genes in which the products localize to the postsynaptic membrane. We then used public expression datasets from human fibroblasts to establish transcription modules, sets of genes coexpressed in this cell type. We identified those sets in which the average gene expression was altered in WBS samples. Dysregulated modules are often interconnected and share multiple common genes, suggesting that intricate regulatory networks connected by a few central genes are disturbed in WBS. This modular approach increases the power to identify pathways dysregulated in WBS patients, thus providing a testable set of additional candidates for genes and their interactions that modulate the WBS phenotypes.