Extremely low-coverage sequencing and imputation increases power for genome-wide association studies

Genome-wide association studies (GWAS) have proven to be a powerful method to identify common genetic variants contributing to susceptibility to common diseases. Here, we show that extremely low-coverage sequencing (0.1-0.5×) captures almost as much of the common (>5%) and low-frequency (1-5%) variation across the genome as SNP arrays. As an empirical demonstration, we show that genome-wide SNP genotypes can be inferred at a mean r(2) of 0.71 using off-target data (0.24× average coverage) in a whole-exome study of 909 samples. Using both simulated and real exome-sequencing data sets, we show that association statistics obtained using extremely low-coverage sequencing data attain similar P values at known associated variants as data from genotyping arrays, without an excess of false positives. Within the context of reductions in sample preparation and sequencing costs, funds invested in extremely low-coverage sequencing can yield several times the effective sample size of GWAS based on SNP array data and a commensurate increase in statistical power.     

Physiology: does gut hormone PYY3-36 decrease food intake in rodents?

Batterham et al. report that the gut peptide hormone PYY3-36 decreases food intake and body-weight gain in rodents, a discovery that has been heralded as potentially offering a new therapy for obesity. However, we have been unable to replicate their results. Although the reasons for this discrepancy remain undetermined, an effective anti-obesity drug ultimately must produce its effects across a range of situations. The fact that the findings of Batterham et al. cannot easily be replicated calls into question the potential value of an anti-obesity approach that is based on administration of PYY3-36.     

Dynamics of seismogenic volcanic extrusion at Mount St Helens in 2004-05

The 2004-05 eruption of Mount St Helens exhibited sustained, near-equilibrium behaviour characterized by relatively steady extrusion of a solid dacite plug and nearly periodic shallow earthquakes. Here we present a diverse data set to support our hypothesis that these earthquakes resulted from stick-slip motion along the margins of the plug as it was forced incrementally upwards by ascending, solidifying, gas-poor magma. We formalize this hypothesis with a dynamical model that reveals a strong analogy between behaviour of the magma-plug system and that of a variably damped oscillator. Modelled stick-slip oscillations have properties that help constrain the balance of forces governing the earthquakes and eruption, and they imply that magma pressure never deviated much from the steady equilibrium pressure. We infer that the volcano was probably poised in a near-eruptive equilibrium state long before the onset of the 2004-05 eruption.     

Bacterial carbon processing by generalist species in the coastal ocean

The assimilation and mineralization of dissolved organic carbon (DOC) by marine bacterioplankton is a major process in the ocean carbon cycle. However, little information exists on the specific metabolic functions of participating bacteria and on whether individual taxa specialize on particular components of the marine DOC pool. Here we use experimental metagenomics to show that coastal communities are populated by taxa capable of metabolizing a wide variety of organic carbon compounds. Genomic DNA captured from bacterial community subsets metabolizing a single model component of the DOC pool (either dimethylsulphoniopropionate or vanillate) showed substantial overlap in gene composition as well as a diversity of carbon-processing capabilities beyond the selected phenotypes. Our direct measure of niche breadth for bacterial functional assemblages indicates that, in accordance with ecological theory, heterogeneity in the composition and supply of organic carbon to coastal oceans may favour generalist bacteria. In the important interplay between microbial community structure and biogeochemical cycling, coastal heterotrophic communities may be controlled less by transient changes in the carbon reservoir that they process and more by factors such as trophic interactions and physical conditions.     

A voltage-gated proton-selective channel lacking the pore domain

Voltage changes across the cell membrane control the gating of many cation-selective ion channels. Conserved from bacteria to humans, the voltage-gated-ligand superfamily of ion channels are encoded as polypeptide chains of six transmembrane-spanning segments (S1-S6). S1-S4 functions as a self-contained voltage-sensing domain (VSD), in essence a positively charged lever that moves in response to voltage changes. The VSD 'ligand' transmits force via a linker to the S5-S6 pore domain 'receptor', thereby opening or closing the channel. The ascidian VSD protein Ci-VSP gates a phosphatase activity rather than a channel pore, indicating that VSDs function independently of ion channels. Here we describe a mammalian VSD protein (H(V)1) that lacks a discernible pore domain but is sufficient for expression of a voltage-sensitive proton-selective ion channel activity. H(v)1 currents are activated at depolarizing voltages, sensitive to the transmembrane pH gradient, H+-selective, and Zn2+-sensitive. Mutagenesis of H(v)1 identified three arginine residues in S4 that regulate channel gating and two histidine residues that are required for extracellular inhibition of H(v)1 by Zn2+. H(v)1 is expressed in immune tissues and manifests the characteristic properties of native proton conductances (G(vH+)). In phagocytic leukocytes, G(vH+) are required to support the oxidative burst that underlies microbial killing by the innate immune system. The data presented here identify H(v)1 as a long-sought voltage-gated H+ channel and establish H(v)1 as the founding member of a family of mammalian VSD proteins.