Stratified randomization controls better for batch effects in 450K methylation analysis: a cautionary tale

Background: Batch effects in DNA methylation microarray experiments can lead to spurious results if not properly handled during the plating of samples. Methods: Two pilot studies examining the association of DNA methylation patterns across the genome with obesity in Samoan men were investigated for chip- and row-specific batch effects. For each study, the DNA of 46 obese men and 46 lean men were assayed using Illumina's Infinium HumanMethylation450 BeadChip. In the first study (Sample One), samples from obese and lean subjects were examined on separate chips. In the second study (Sample Two), the samples were balanced on the chips by lean/obese status, age group, and census region. We used methylumi, watermelon, and limma R packages, as well as ComBat, to analyze the data. Principal component analysis and linear regression were respectively employed to identify the top principal components and to test for their association with the batches and lean/obese status. To identify differentially methylated positions (DMPs) between obese and lean males at each locus, we used a moderated t-test.Results: Chip effects were effectively removed from Sample Two but not Sample One. In addition, dramatic differences were observed between the two sets of DMP results. After removing'' batch effects with ComBat, Sample One had 94,191 probes differentially methylated at a q-value threshold of 0.05 while Sample Two had zero differentially methylated probes. The disparate results from Sample One and Sample Two likely arise due to the confounding of lean/obese status with chip and row batch effects.Conclusion: Even the best possible statistical adjustments for batch effects may not completely remove them. Proper study design is vital for guarding against spurious findings due to such effects

Structural basis for the interaction of a human small heat shock protein with the 14-3-3 universal signaling regulator

By interacting with hundreds of protein partners, 14-3-3 proteins coordinate vital cellular processes. Phosphorylation of the small heat shock protein, HSPB6, within its intrinsically disordered N-terminal domain activates its interaction with 14-3-3, ultimately triggering smooth muscle relaxation. After analyzing the binding of an HSPB6-derived phosphopeptide to 14-3-3 using isothermal calorimetry and X-ray crystallography, we have determined the crystal structure of the complete assembly consisting of the 14-3-3 dimer and full-length HSPB6 dimer and further characterized this complex in solution using fluorescence spectroscopy, small-angle X-ray scattering, and limited proteolysis. We show that selected intrinsically disordered regions of HSPB6 are transformed into well-defined conformations upon the interaction, whereby an unexpectedly asymmetric structure is formed. This structure provides the first atomic resolution snapshot of a human small HSP in functional state, explains how 14-3-3 proteins sequester their regulatory partners, and can inform the design of small-molecule interaction modifiers to be used as myorelaxants

Adjusting the melting point of a model system via Gibbs-Duhem integration: application to a model of Aluminum

Model interaction potentials for real materials are generally optimized with respect to only those experimental properties that are easily evaluated as mechanical averages (e.g., elastic constants (at T=0 K), static lattice energies and liquid structure). For such potentials, agreement with experiment for the non-mechanical properties, such as the melting point, is not guaranteed and such values can deviate significantly from experiment. We present a method for re-parameterizing any model interaction potential of a real material to adjust its melting temperature to a value that is closer to its experimental melting temperature. This is done without significantly affecting the mechanical properties for which the potential was modeled. This method is an application of Gibbs-Duhem integration [D. Kofke, Mol. Phys.78, 1331 (1993)]. As a test we apply the method to an embedded atom model of aluminum [J. Mei and J.W. Davenport, Phys. Rev. B 46, 21 (1992)] for which the melting temperature for the thermodynamic limit is 826.4 +/- 1.3K - somewhat below the experimental value of 933K. After re-parameterization, the melting temperature of the modified potential is found to be 931.5K +/- 1.5K.Comment: 9 pages, 5 figures, 4 table

Specific Sequences in the N-terminal Domain of Human Small Heat Shock Protein HSPB6 Dictate Preferential Heterooligomerization with the Orthologue HSPB1

Small heat-shock proteins (sHSPs) are a conserved group of molecular chaperones with important roles in cellular proteostasis. Although sHSPs are characterized by their small monomeric weight, they typically assemble into large polydisperse oligomers that vary in both size and shape but are principally composed of dimeric building blocks. These assemblies can include different sHSP orthologues, creating additional complexity that may affect chaperone activity. However, the structural and functional properties of such hetero-oligomers are poorly understood. We became interested in hetero-oligomer formation between human heat-shock protein family B (small) member 1 (HSPB1) and HSPB6, which are both highly expressed in skeletal muscle. When mixed in vitro, these two sHSPs form a polydisperse oligomer array composed solely of heterodimers, suggesting preferential association that is determined at the monomer level. Previously, we have shown that the sHSP N-terminal domains (NTDs), which have a high degree of intrinsic disorder, are essential for the biased formation. Here we employed iterative deletion mapping to elucidate how the NTD of HSPB6 influences its preferential association with HSPB1 and show that this region has multiple roles in this process. First, the highly conserved motif RLFDQXFG is necessary for subunit exchange among oligomers. Second, a site ∼20 residues downstream of this motif determines the size of the resultant hetero-oligomers. Third, a region unique to HSPB6 dictates the preferential formation of heterodimers. In conclusion, the disordered NTD of HSPB6 helps regulate the size and stability of hetero-oligomeric complexes, indicating that terminal sHSP regions define the assembly properties of these proteins

Can Systemic Interventions Designed to Reduce Reoffending by Youth also Reduce their Victimization?

Previous research indicates considerable overlap between populations of boys who are victimized and boys who victimize others. This study was concerned with whether a systems-focused treatment program designed to address individual and systemic risk factors associated with the perpetration of sexual and violent crimes might also be successful in reducing boys’ victimization by others. Boys adjudicated for sexual offences who received ‘treatment as usual’ (TAU; n = 335) were compared with similarly adjudicated boys who completed the treatment program (n = 200) on their histories of contact with police either as offenders or victims. Despite their higher rates of pre-intervention victimization, the treatment group were victimized less frequently post-intervention than the TAU group. Continued offending was the strongest predictor of victimization post-intervention. These findings suggest that offending and victimization share common risk factors that may be addressed simultaneously within offence-focused treatment

Crop Updates 2010 - Farming Systems

This session covers twenty papers from different authors: Pests and Disease 1. Preserving phosphine for use in Grain Storage Industry, Christopher R Newman, Department of Agriculture and Food Farming Systems Research 2. Demonstrating the benefits of grazing canola in Western Australia, Jonathan England, Stephen Gherardi and Mohammad Amjad, Department of Agriculture and Food 3. Buloke barley yield when pasture-cropped across subtropical perennial pastures, David Ferris, Department of Agriculture and Food, Phil Ward and Roger Lawes, CSIRO 4. Is pasture cropping viable in WA? Grower perceptions and EverCrop initiatives to evaluate, David Ferris, Tim Wiley, Perry Dolling, Department of Agriculture and Food, Philip Barrett-Lennard, Evergreen farming 5. Best-bet management for dual-purpose canola, John Kirkegaard, Susan Sprague, Hugh Dove and Walter Kelman, CSIRO, Canberra, Peter Hamblin, Agritech Research, Young, NSW 6. Pasture in cropping systems – with and without sheep, Brad Nutt and Angelo Loi, Department of Agriculture and Food 7. Can technology substitute for a lupin break? Wayne Parker, Department of Agriculture and Food 8. Canola row spacing with and without long term stubble retention on a sandy clay loam at Merredin, Glen Riethmuller, Department of Agriculture and Food 9. Impact of stubble retention on water balance and crop yield, Phil Ward1, Ken Flower2,3, Neil Cordingley2 and Shayne Micin1, 1CSIRO, Wembley, Western Australia, 2Western Australian No-Till Farmers Association, 3University of Western Australia Analysis and Modelling 10. Using POAMA rainfall forecasts for crop management in South-West WA, Senthold Asseng1, Peter McIntosh2,3, Mike Pook2,3, James Risbey2,3, Guomin Wang3, Oscar Alves3, Ian Foster4, Imma Farre4 and Nirav Khimashia1, 1CSIRO Plant Industry, Perth, 2CSIRO Marine and Atmospheric Research, Hobart, 3Centre for Australian Weather and Climate Research (CAWCR), A partnership between the Australian Bureau of Meteorology and CSIRO, Melbourne, 4Department of Agriculture and Food 11. Adaption to changing climates and variability – results of the Agribusiness Changing Climates regional workshop, Anderson W3, Beard D3, Blake J3, Grieve R1, Lang M3, Lemon J3, McTaggart R3, Gray D3, Price M2 and Stephens D3, 1Roderick Grieve Farm Management Consultants, 2Coffey International P/L, 3Department of Agriculture and Food 12. Farmers’ management of seasonal variability and climate change in WA, DA Beard, DM Gray, P Carmody, Department of Agriculture and Food 13. Is there a value in having a frost forecast for wheat in South-West WA? Imma Farre1, Senthold Asseng2, Ian Foster1 and Doug Abrecht3, 1Department of Agriculture and Food, CSIRO, Floreat, 2CSIRO Plant Industry, Perth 3Department of Agriculture and Food, Centre for Cropping Systems 14. Does buying rainfall pay? Greg Kirk, Planfarm Agricultural Consultants 15. Which region in the WA wheatbelt makes best use of rainfall? Peter Rowe, Bankwest Agribusiness 16. POAMA – the Predictive Ocean-Atmosphere Model for Australia, Guomin Wang and Oscar Alves, Centre for Australian Weather and Climate Research (CAWCR), A partnership between the Australian Bureau of Meteorology and CSIRO, Melbourne 17. Exploring the link between water use efficiency and farm profitability, Cameron Weeks, Planfarm and Peter Tozer, PRT Consulting Precision Agriculture 18. A plethora of paddock information is available – how does it stack up? Derk Bakker, Department of Agriculture and Food 18. Variable rate prescription mapping for lime inputs based on electromagnetic surveying and deep soil testing, Frank D’Emden, Quenten Knight and Luke Marquis, Precision Agronomics, Australia 19. Trial design and analysis using precision agriculture and farmer’s equipment, Roger Lawes, CSIRO Sustainable Ecosystems, Centre for Environment and Life Sciences, Floreat 20. Farmer perspectives of precision agriculture in Western Australia: Issues and the way forward, Dr Roger Mandel, Curtin Universit

A tagging SNP in INSIG2 is associated with obesity-related phenotypes among Samoans

<p>Abstract</p> <p>Background</p> <p>A genome wide association study found significant association of a sequence variant, rs7566605, in the insulin-induced gene 2 (<it>INSIG2</it>) with obesity. However, the association remained inconclusive in follow-up studies. We tested for association of four tagging SNPs (tagSNPs) including this variant with body mass index (BMI) and abdominal circumference (ABDCIR) in the Samoans of the Western Pacific, a population with high levels of obesity.</p> <p>Methods</p> <p>We studied 907 adult Samoan participants from a longitudinal study of adiposity and cardiovascular disease risk in two polities, American Samoa and Samoa. Four tagSNPs were identified from the Chinese HapMap database based on pairwise <it>r</it>^{<it>2 </it>}of ≥0.8 and minor allele frequency of ≥0.05. Genotyping was performed using the TaqMan assay. Tests of association with BMI and ABDCIR were performed under the additive model.</p> <p>Results</p> <p>We did not find association of rs7566605 with either BMI or ABDCIR in any group of the Samoans. However, the most distally located tagSNPs in Intron 3 of the gene, rs9308762, showed significant association with both BMI (p-value 0.024) and ABDCIR (p-value 0.009) in the combined sample and with BMI (p-value 0.038) in the sample from Samoa.</p> <p>Conclusion</p> <p>Although rs7566605 was not significantly associated with obesity in our study population, we can not rule out the involvement of <it>INSIG2 </it>in obesity related traits as we found significant association of another tagSNP in <it>INSIG2 </it>with both BMI and ABDCIR. This study suggests the importance of comprehensive assessment of sequence variants within a gene in association studies.</p

Wolbachia infections that reduce immature insect survival: Predicted impacts on population replacement

<p>Abstract</p> <p>Background</p> <p>The evolutionary success of <it>Wolbachia </it>bacteria, infections of which are widespread in invertebrates, is largely attributed to an ability to manipulate host reproduction without imposing substantial fitness costs. Here, we describe a stage-structured model with deterministic immature lifestages and a stochastic adult female lifestage. Simulations were conducted to better understand <it>Wolbachia </it>invasions into uninfected host populations. The model includes conventional <it>Wolbachia </it>parameters (the level of cytoplasmic incompatibility, maternal inheritance, the relative fecundity of infected females, and the initial <it>Wolbachia </it>infection frequency) and a new parameter termed relative larval viability (<it>RLV</it>), which is the survival of infected larvae relative to uninfected larvae.</p> <p>Results</p> <p>The results predict the <it>RLV </it>parameter to be the most important determinant for <it>Wolbachia </it>invasion and establishment. Specifically, the fitness of infected immature hosts must be close to equal to that of uninfected hosts before population replacement can occur. Furthermore, minute decreases in <it>RLV </it>inhibit the invasion of <it>Wolbachia </it>despite high levels of cytoplasmic incompatibility, maternal inheritance, and low adult fitness costs.</p> <p>Conclusions</p> <p>The model described here takes a novel approach to understanding the spread of <it>Wolbachia </it>through a population with explicit dynamics. By combining a stochastic female adult lifestage and deterministic immature/adult male lifestages, the model predicts that even those <it>Wolbachia </it>infections that cause minor decreases in immature survival are unlikely to invade and spread within the host population. The results are discussed in relation to recent theoretical and empirical studies of natural population replacement events and proposed applied research, which would use <it>Wolbachia </it>as a tool to manipulate insect populations.</p