ParaHaplo 2.0: a program package for haplotype-estimation and haplotype-based whole-genome association study using parallel computing

<p>Abstract</p> <p>Background</p> <p>The use of haplotype-based association tests can improve the power of genome-wide association studies. Since the observed genotypes are unordered pairs of alleles, haplotype phase must be inferred. However, estimating haplotype phase is time consuming. When millions of single-nucleotide polymorphisms (SNPs) are analyzed in genome-wide association study, faster methods for haplotype estimation are required.</p> <p>Methods</p> <p>We developed a program package for parallel computation of haplotype estimation. Our program package, ParaHaplo 2.0, is intended for use in workstation clusters using the Intel Message Passing Interface (MPI). We compared the performance of our algorithm to that of the regular permutation test on both Japanese in Tokyo, Japan and Han Chinese in Beijing, China of the HapMap dataset.</p> <p>Results</p> <p>Parallel version of ParaHaplo 2.0 can estimate haplotypes 100 times faster than a non-parallel version of the ParaHaplo.</p> <p>Conclusion</p> <p>ParaHaplo 2.0 is an invaluable tool for conducting haplotype-based genome-wide association studies (GWAS). The need for fast haplotype estimation using parallel computing will become increasingly important as the data sizes of such projects continue to increase. The executable binaries and program sources of ParaHaplo are available at the following address: <url>http://en.sourceforge.jp/projects/parallelgwas/releases/</url></p

Genotype determination for polymorphisms in linkage disequilibrium

<p>Abstract</p> <p>Background</p> <p>Genome-wide association studies with single nucleotide polymorphisms (SNPs) show great promise to identify genetic determinants of complex human traits. In current analyses, genotype calling and imputation of missing genotypes are usually considered as two separated tasks. The genotypes of SNPs are first determined one at a time from allele signal intensities. Then the missing genotypes, i.e., no-calls caused by not perfectly separated signal clouds, are imputed based on the linkage disequilibrium (LD) between multiple SNPs. Although many statistical methods have been developed to improve either genotype calling or imputation of missing genotypes, treating the two steps independently can lead to loss of genetic information.</p> <p>Results</p> <p>We propose a novel genotype calling framework. In this framework, we consider the signal intensities and underlying LD structure of SNPs simultaneously by estimating both cluster parameters and haplotype frequencies. As a result, our new method outperforms some existing algorithms in terms of both call rates and genotyping accuracy. Our studies also suggest that jointly analyzing multiple SNPs in LD provides more accurate estimation of haplotypes than haplotype reconstruction methods that only use called genotypes.</p> <p>Conclusion</p> <p>Our study demonstrates that jointly analyzing signal intensities and LD structure of multiple SNPs is a better way to determine genotypes and estimate LD parameters.</p

Rapid and Accurate Multiple Testing Correction and Power Estimation for Millions of Correlated Markers

With the development of high-throughput sequencing and genotyping technologies, the number of markers collected in genetic association studies is growing rapidly, increasing the importance of methods for correcting for multiple hypothesis testing. The permutation test is widely considered the gold standard for accurate multiple testing correction, but it is often computationally impractical for these large datasets. Recently, several studies proposed efficient alternative approaches to the permutation test based on the multivariate normal distribution (MVN). However, they cannot accurately correct for multiple testing in genome-wide association studies for two reasons. First, these methods require partitioning of the genome into many disjoint blocks and ignore all correlations between markers from different blocks. Second, the true null distribution of the test statistic often fails to follow the asymptotic distribution at the tails of the distribution. We propose an accurate and efficient method for multiple testing correction in genome-wide association studies—SLIDE. Our method accounts for all correlation within a sliding window and corrects for the departure of the true null distribution of the statistic from the asymptotic distribution. In simulations using the Wellcome Trust Case Control Consortium data, the error rate of SLIDE's corrected p-values is more than 20 times smaller than the error rate of the previous MVN-based methods' corrected p-values, while SLIDE is orders of magnitude faster than the permutation test and other competing methods. We also extend the MVN framework to the problem of estimating the statistical power of an association study with correlated markers and propose an efficient and accurate power estimation method SLIP. SLIP and SLIDE are available at http://slide.cs.ucla.edu

Analysis of her1 and her7 Mutants Reveals a Spatio Temporal Separation of the Somite Clock Module

Somitogenesis is controlled by a genetic network consisting of an oscillator (clock) and a gradient (wavefront). The “hairy and Enhancer of Split”- related (her) genes act downstream of the Delta/Notch (D/N) signaling pathway, and are crucial components of the segmentation clock. Due to genome duplication events, the zebrafish genome, possesses two gene copies of the mouse Hes7 homologue: her1 and her7. To better understand the functional consequences of this gene duplication, and to determine possible independent roles for these two genes during segmentation, two zebrafish mutants her1hu2124 and her7hu2526 were analyzed. In the course of embryonic development, her1hu2124 mutants exhibit disruption of the three anterior-most somite borders, whereas her7hu2526 mutants display somite border defects restricted to somites 8 (+/−3) to 17 (+/−3) along the anterior-posterior axis. Analysis of the molecular defects in her1hu2124 mutants reveals a her1 auto regulatory feedback loop during early somitogenesis that is crucial for correct patterning and independent of her7 oscillation. This feedback loop appears to be restricted to early segmentation, as cyclic her1 expression is restored in her1hu2124 embryos at later stages of development. Moreover, only the anterior deltaC expression pattern is disrupted in the presomitic mesoderm of her1hu2124 mutants, while the posterior expression pattern of deltaC remains unaltered. Together, this data indicates the existence of an independent and genetically separable anterior and posterior deltaC clock modules in the presomitic mesdorm (PSM)

Association between high-dose erythropoiesis-stimulating agents, inflammatory biomarkers, and soluble erythropoietin receptors

<p>Abstract</p> <p>Background</p> <p>High-dose erythropoiesis-stimulating agents (ESA) for anemia of chronic kidney disease (CKD) have been associated with adverse clinical outcomes and do not always improve erythropoiesis. We hypothesized that high-dose ESA requirement would be associated with elevated inflammatory biomarkers, decreased adipokines, and increased circulating, endogenous soluble erythropoietin receptors (sEpoR).</p> <p>Methods</p> <p>A cross-sectional cohort of anemic 32 CKD participants receiving ESA were enrolled at a single center and cytokine profiles, adipokines, and sEpoR were compared between participants stratified by ESA dose requirement (usual-dose darbepoetin-α (< 1 μg/kg/week) and high-dose (≥1 μg/kg/week)).</p> <p>Results</p> <p>Baseline characteristics were similar between groups; however, hemoglobin was lower among participants on high-dose (1.4 μg/kg/week) vs usual-dose (0.5 μg/kg/week) ESA.</p> <p>In adjusted analyses, high-dose ESA was associated with an increased odds for elevations in c-reactive protein and interleukin-6 (p < 0.05 for both). There was no correlation between high-dose ESA and adipokines. Higher ESA dose correlated with higher levels of sEpoR (r_s= 0.39, p = 0.03). In adjusted analyses, higher ESA dose (per μcg/kg/week) was associated with a 53% greater odds of sEpoR being above the median (p < 0.05).</p> <p>Conclusion</p> <p>High-dose ESA requirement among anemic CKD participants was associated with elevated inflammatory biomarkers and higher levels of circulating sEpoR, an inhibitor of erythropoiesis. Further research confirming these findings is warranted.</p> <p>Trial registration</p> <p>Clinicaltrials.gov <a href="http://www.clinicaltrials.gov/ct2/show/NCT00526747">NCT00526747</a></p

Hypoxia Disruption of Vertebrate CNS Pathfinding through EphrinB2 Is Rescued by Magnesium

The mechanisms of hypoxic injury to the developing human brain are poorly understood, despite being a major cause of chronic neurodevelopmental impairments. Recent work in the invertebrate Caenorhabditis elegans has shown that hypoxia causes discrete axon pathfinding errors in certain interneurons and motorneurons. However, it is unknown whether developmental hypoxia would have similar effects in a vertebrate nervous system. We have found that developmental hypoxic injury disrupts pathfinding of forebrain neurons in zebrafish (Danio rerio), leading to errors in which commissural axons fail to cross the midline. The pathfinding defects result from activation of the hypoxia-inducible transcription factor (hif1) pathway and are mimicked by chemical inducers of the hif1 pathway or by expression of constitutively active hif1α. Further, we found that blocking transcriptional activation by hif1α helped prevent the guidance defects. We identified ephrinB2a as a target of hif1 pathway activation, showed that knock-down of ephrinB2a rescued the guidance errors, and showed that the receptor ephA4a is expressed in a pattern complementary to the misrouting axons. By targeting a constitutively active form of ephrinB2a to specific neurons, we found that ephrinB2a mediates the pathfinding errors via a reverse-signaling mechanism. Finally, magnesium sulfate, used to improve neurodevelopmental outcomes in preterm births, protects against pathfinding errors by preventing upregulation of ephrinB2a. These results demonstrate that evolutionarily conserved genetic pathways regulate connectivity changes in the CNS in response to hypoxia, and they support a potential neuroprotective role for magnesium

Animal abuse and intimate partner violence: researching the link and its significance in Ireland - a veterinary perspective

Research on domestic violence has established a substantial association between intimate partner abuse and the abuse of children within the home. It is only recently however, that researchers have demonstrated the correlation between non-accidental injury in animals, and abuse of women by their intimate male partners. A growing body of evidence suggests that animal abuse can be an early indicator for other forms of violent behaviour. This research includes the responses of a sample of 23 women using refuge services in the Republic of Ireland. It investigates the connection between domestic violence and animal abuse, and ascertains if there is sufficient support service for animals and people relevant to domestic abuse. In the survey population, 57% of women reported witnessing one or more forms of abuse, or threats of abuse, of their pets. Five of which were reported to have resulted in the death of the pet. Eighty seven per cent of women felt a facility to accommodate pets would have made their decision to leave the family home easier. Four women disclosed that lack of such a service and concern for the welfare of their companion animals caused them to remain in their abusive relationships for longer than they felt appropriate. Nine families placed pets in the care of family or friends, one woman is unaware of the fate of her pet, while the pets of six families remained with the abusive male after his partner entered a refuge. The majority of women felt unable to talk to anyone about their fears for their pets' welfare. Many felt that there is no service which can provide temporary accommodation for womens' pets while they are in refuge. The results obtained support those found elsewhere in larger studies in the USA and UK, and demonstrate an association of animal abuse in households where there is reported domestic violence

Expression Analysis of PAC1-R and PACAP Genes in Zebrafish Embryos

This study describes the expression of the pituitary adenylate cyclase-activating polypeptide (PACAP1 and PACAP2) and PAC1 receptor genes (PAC1a-R and PAC1b-R) in the brain of zebrafish (Danio rerio) during development. In situ hybridization of the 24- and 48-hpf embryos revealed that PACAP genes were expressed in the telencephalon, the diencephalon, the rhombencephalon, and the neurons in the dorsal part of the spinal cord. PACAP2 mRNA appears to be the most abundant form during brain development. The two PAC1-R subtypes showed a similar expression pattern: mRNAs were detected in the forebrain, the thalamus, and the rhombencephalon. However, in the tectum, only PAC1b-R gene was detected. These results suggest that, in fish, PACAP may play a role in brain development

Purification and In Situ Immobilization of Papain with Aqueous Two-Phase System

Papain was purified from spray-dried Carica papaya latex using aqueous two-phase system (ATPS). Then it was recovered from PEG phase by in situ immobilization or preparing cross-linked enzyme aggregates (CLEAs). The Plackett-Burman design and the central composite design (CCD) together with the response surface methodology (RSM) were used to optimize the APTS processes. The highly purified papain (96–100%) was achieved under the optimized conditions: 40% (w/w) 15 mg/ml enzyme solution, 14.33–17.65% (w/w) PEG 6000, 14.27–14.42% (w/w) NaH2PO4/K2HPO4 and pH 5.77–6.30 at 20°C. An in situ enzyme immobilization approach, carried out by directly dispersing aminated supports and chitosan beads into the PEG phase, was investigated to recover papain, in which a high immobilization yield (>90%) and activity recovery (>40%) was obtained. Moreover, CLEAs were successfully used in recovering papain from PEG phase with a hydrolytic activity hundreds times higher than the carrier-bound immobilized papain