GPU accelerated sequence alignment with traceback for GATK HaplotypeCaller

Background

Pairwise sequence alignment is widely used in many biological tools and applications. Existing GPU accelerated implementations mainly focus on calculating optimal alignment score and omit identifying the optimal alignment itself. In GATK HaplotypeCaller (HC), the semi-global pairwise sequence alignment with traceback has so far been difficult to accelerate effectively on GPUs.

Results

We first analyze the characteristics of the semi-global alignment with traceback in GATK HC and then propose a new algorithm that allows for retrieving the optimal alignment efficiently on GPUs. For the first stage, we choose intra-task parallelization model to calculate the position of the optimal alignment score and the backtracking matrix. Moreover, in the first stage, our GPU implementation also records the length of consecutive matches/mismatches in addition to lengths of consecutive insertions and deletions as in the CPU-based implementation. This helps efficiently retrieve the backtracking matrix to obtain the optimal alignment in the second stage.

Conclusions

Experimental results show that our alignment kernel with traceback is up to 80x and 14.14x faster than its CPU counterpart with synthetic datasets and real datasets, respectively. When integrated into GATK HC (alongside a GPU accelerated pair-HMMs forward kernel), the overall acceleration is 2.3x faster than the baseline GATK HC implementation, and 1.34x faster than the GATK HC implementation with the integrated GPU-based pair-HMMs forward algorithm. Although the methods proposed in this paper is to improve the performance of GATK HC, they can also be used in other pairwise alignments and applications.

     

Plant growth regulators improve growth,photosynthesis, mineral nutrient and antioxidant system under cadmium stress in menthol mint ( Mentha arvensis L.)

Menthol mint (Mentha arvensis L.) cultivation is significantly affected by the heavy metals like cadmium (Cd) which also imposes severe health hazards. Two menthol mint cultivars namely Kosi and Kushal were evaluated under Cd stress conditions. Impact of plant growth regulators (PGRs) like salicylic acid (SA), gibberellic acid (GA3) and triacontanol (Tria) on Cd stress tolerance was assessed. Reduced growth, photosynthetic parameters, mineral nutrient concentration, and increased oxidative stress biomarkers like electrolyte leakage, malondialdehyde, and hydrogen peroxide contents were observed under Cd stress. Differential upregulation of proline content and antioxidant activities under Cd stress was observed in both the cultivars. Interestingly, low electrolyte leakage, lipid peroxidation, hydrogen peroxide and Cd concentration in leaves were observed in Kushal compared to Kosi. Among all the PGRs tested, SA proved to be the best in improving Cd-stress tolerance in both the cultivars but Kushal responded better than Kosi.     

Growth of Escherichia coli in iron-enriched medium increases HPI catalase activity

Escherichia coli has two catalases, HPI and HPII. HPI is induced during logarithmic growth in response to low concentrations of hydrogen peroxide. This induction is OxyR-dependent. On the other hand, HPII is not peroxide-inducible but is induced in entry to the stationary phase. We demonstrate here that E. coli displayed higher HPI catalase activity when compared to the cultures that were grown in a normal medium, if grown in a medium supplemented with iron-citrate. Iron supplementation had no effect on HPII catalase. This increase of HPI activity was OxyR-independent and not observed in a Deltafur mutant. The physiological significance of the increase of HPI activity is unclear, but it appears that the katG gene that codes for HPI catalase is among the genes that are regulated by Fur.     

Performance-based selection of likelihood models for phylogeny estimation

Phylogenetic estimation has largely come to rely on explicitly model-based methods. This approach requires that a model be chosen and that that choice be justified. To date, justification has largely been accomplished through use of likelihood-ratio tests (LRTs) to assess the relative fit of a nested series of reversible models. While this approach certainly represents an important advance over arbitrary model selection, the best fit of a series of models may not always provide the most reliable phylogenetic estimates for finite real data sets, where all available models are surely incorrect. Here, we develop a novel approach to model selection, which is based on the Bayesian information criterion, but incorporates relative branch-length error as a performance measure in a decision theory (DT) framework. This DT method includes a penalty for overfitting, is applicable prior to running extensive analyses, and simultaneously compares all models being considered and thus does not rely on a series of pairwise comparisons of models to traverse model space. We evaluate this method by examining four real data sets and by using those data sets to define simulation conditions. In the real data sets, the DT method selects the same or simpler models than conventional LRTs. In order to lend generality to the simulations, codon-based models (with parameters estimated from the real data sets) were used to generate simulated data sets, which are therefore more complex than any of the models we evaluate. On average, the DT method selects models that are simpler than those chosen by conventional LRTs. Nevertheless, these simpler models provide estimates of branch lengths that are more accurate both in terms of relative error and absolute error than those derived using the more complex (yet still wrong) models chosen by conventional LRTs. This method is available in a program called DT-ModSel.     

Transcriptional status of known and novel genes tagged with consensus of 33.15 repeat loci employing minisatellite-associated sequence amplification (MASA) and real-time PCR in water buffalo, Bubalus bubalis

We conducted minisatellite-associated sequence amplification (MASA) with an oligo (5' CACCTCTCCACCTGCC 3') based on consensus of 33.15 repeat loci using cDNA from the testis, ovary, spleen, kidney, heart, liver, and lung of water buffalo Bubalus bubalis and uncovered 25 amplicons of six different sizes (1,263, 846/847, 602, 576, 487, and 324 base pairs). These fragments, cloned and sequenced, were found to represent several functional, regulatory, and structural genes. Blast search of all the 25 amplicons showed homologies with 43 transcribing genes across the species. Of these, the 846/847-bp fragment, having homology with the adenylate kinase gene, showed nucleotide changes at six identical places in the ovary and testis. The 1,263; 324; and 487-bp fragments showed homology with the secreted modular calcium binding protein (SMOC-1), leucine-rich repeat neuronal 6A (LRRN6A) mRNA, and human TTTY5 mRNA, respectively. Real-time PCR showed maximum expression of AKL, LRRN6A, and T-cell receptor gamma (TCR-gamma)-like genes in the testis, SMOC-1 in the liver, and the T-cell receptor-like (TCRL) gene in the spleen compared to those used as endogenous control. We construe that these genes have evolved from a common progenitor and conformed to various biological functions during the course of evolution. MASA approach coupled with real-time PCR has potentials to uncover accurate expression of a large number of genes within and across the species circumventing the screening of cDNA library.     

Potential Early Predictors for Outcomes of Experimental Hemorrhagic Shock Induced by Uncontrolled Internal Bleeding in Rats

Uncontrolled hemorrhage, resulting from traumatic injuries, continues to be the leading cause of death in civilian and military environments. Hemorrhagic deaths usually occur within the first 6 hours of admission to hospital; therefore, early prehospital identification of patients who are at risk for developing shock may improve survival. The aims of the current study were: 1. To establish and characterize a unique model of uncontrolled internal hemorrhage induced by massive renal injury (MRI), of different degrees (20-35% unilateral nephrectomy) in rats, 2. To identify early biomarkers those best predict the outcome of severe internal hemorrhage. For this purpose, male Sprague Dawley rats were anesthetized and cannulas were inserted into the trachea and carotid artery. After abdominal laparotomy, the lower pole of the kidney was excised. During 120 minutes, hematocrit, pO₂, pCO₂, base excess, potassium, lactate and glucose were measured from blood samples, and mean arterial pressure (MAP) was measured through arterial tracing. After 120 minutes, blood loss was determined. Statistical prediction models of mortality and amount of blood loss were performed. In this model, the lowest blood loss and mortality rate were observed in the group with 20% nephrectomy. Escalation of the extent of nephrectomy to 25% and 30% significantly increased blood loss and mortality rate. Two phases of hemodynamic and biochemical response to MRI were noticed: the primary phase, occurring during the first 15 minutes after injury, and the secondary phase, beginning 30 minutes after the induction of bleeding. A Significant correlation between early blood loss and mean arterial pressure (MAP) decrements and survival were noted. Our data also indicate that prediction of outcome was attainable in the very early stages of blood loss, over the first 15 minutes after the injury, and that blood loss and MAP were the strongest predictors of mortality.     

Alpha-actinin-4 is selectively required for insulin-induced GLUT4 translocation

Insulin induces GLUT4 translocation to the muscle cell surface. Using differential amino acid labeling and mass spectrometry, we observed insulin-dependent co-precipitation of actinin-4 (ACTN4) with GLUT4 (Foster, L. J., Rudich, A., Talior, I., Patel, N., Huang, X., Furtado, L. M., Bilan, P. J., Mann, M., and Klip, A. (2006) J. Proteome Res. 5, 64-75). ACTN4 links F-actin to membrane proteins, and actin dynamics are essential for GLUT4 translocation. We hypothesized that ACTN4 may contribute to insulin-regulated GLUT4 traffic. In L6 muscle cells insulin, but not platelet-derived growth factor, increased co-precipitation of ACTN4 with GLUT4. Small interfering RNA-mediated ACTN4 knockdown abolished the gain in surface-exposed GLUT4 elicited by insulin but not by platelet-derived growth factor, membrane depolarization, or mitochondrial uncoupling. In contrast, knockdown of alpha-actinin-1 (ACTN1) did not prevent GLUT4 translocation by insulin. GLUT4 colocalized with ACTN4 along the insulin-induced cortical actin mesh and ACTN4 knockdown prevented GLUT4-actin colocalization without impeding actin remodeling or Akt phosphorylation, maintaining GLUT4 in a tight perinuclear location. We propose that ACTN4 contributes to GLUT4 traffic, likely by tethering GLUT4 vesicles to the cortical actin cytoskeleton.     

Serum Levels of Soluble CD26/Dipeptidyl Peptidase-IV in Type 2 Diabetes Mellitus and Its Association with Metabolic Syndrome and Therapy with Antidiabetic Agents in Malaysian Subjects

BackgroundA soluble form of CD26/dipeptidyl peptidase-IV (sCD26/DPP-IV) induces DPP-IV enzymatic activity that degrades incretin. We investigated fasting serum levels of sCD26/DPP-IV and active glucagon-like peptide-1 (GLP-1) in Malaysian patients with type 2 diabetes mellitus (T2DM) with and without metabolic syndrome (MetS), as well as the associations between sCD26/DPP-IV levels, MetS, and antidiabetic therapy.MethodsWe assessed sCD26/DPP-IV levels, active GLP-1 levels, body mass index (BMI), glucose, insulin, A1c, glucose homeostasis indices, and lipid profiles in 549 Malaysian subjects (including 257 T2DM patients with MetS, 57 T2DM patients without MetS, 71 non-diabetics with MetS, and 164 control subjects without diabetes or metabolic syndrome).ResultsFasting serum levels of sCD26/DPP-IV were significantly higher in T2DM patients with and without MetS than in normal subjects. Likewise, sCD26/DPP-IV levels were significantly higher in patients with T2DM and MetS than in non-diabetic patients with MetS. However, active GLP-1 levels were significantly lower in T2DM patients both with and without MetS than in normal subjects. In T2DM subjects, sCD26/DPP-IV levels were associated with significantly higher A1c levels, but were significantly lower in patients using monotherapy with metformin. In addition, no significant differences in sCD26/DPP-IV levels were found between diabetic subjects with and without MetS. Furthermore, sCD26/DPP-IV levels were negatively correlated with active GLP-1 levels in T2DM patients both with and without MetS. In normal subjects, sCD26/DPP-IV levels were associated with increased BMI, cholesterol, and LDL-cholesterol (LDL-c) levels.ConclusionSerum sCD26/DPP-IV levels increased in T2DM subjects with and without MetS. Active GLP-1 levels decreased in T2DM patients both with and without MetS. In addition, sCD26/DPP-IV levels were associated with Alc levels and negatively correlated with active GLP-1 levels. Moreover, metformin monotherapy was associated with reduced sCD26/DPP-IV levels. In normal subjects, sCD26/DPP-IV levels were associated with increased BMI, cholesterol, and LDL-c.     

In vitro acclimatization of date palm (Phoenix dactylifera L.) plantlets: A quantitative comparison of epicuticular leaf wax as a function of polyethylene glycol treatment

Wax deposits on leaf surfaces ofin vitro-grown plantlets,in vitro plantlets treated with polyethylene glycol and greenhouse-grown seedlings from five cultivars of date palm (Phoenix dactylifera L.) were extracted and quantified. Significant variations among treatments and cultivars were obtained. Greenhouse-grown plants had the greatest wax deposits followed by the acclimatized plantlets.In vitro plantlets had an average of 15% of the wax of greenhouse plants. Cultivar and plant age differences had a significant effect on the quantity of wax deposits. Greenhouse seedlings of Majhool, Deglet Noor and Khadraoui (cultivars grown under irrigation) had less wax accumulation than Zahidi and Sayer, dryland cultivars.The increase in wax deposition as a result of polyethylene glycol treatment, explains in part, the decreased water loss observed in acclimatized plantlets when transferredex vitro. Abbreviations EW epicuticular wax - NAA naphthalene acetic acid - PEG polyethylene glycol     

Metabolic flux analysis of<Emphasis Type="Italic"> pykF</Emphasis> gene knockout<Emphasis Type="Italic"> Escherichia coli</Emphasis> based on <Superscript>13</Superscript>C-labeling experiments together with measurements of enzyme activities and intracellular metabolite concentrations

Metabolic flux analysis based on ¹³C-labeling experiments followed by the measurement of intracellular isotope distribution using both 2D NMR and GC-MS was carried out to investigate the effect of pyruvate kinase (pyk) gene knockout on the metabolism of Escherichia coli in continuous culture. In addition, the activities of 16 enzymes, and the concentrations of 5 intracellular metabolites, were measured as a function of time in batch culture as well as continuous culture. It was found that flux through phosphoenol pyruvate carboxylase and malic enzyme were up-regulated in the pykF^– mutant as compared with the wild type, and acetate formation was significantly reduced in the mutant. In addition, flux through the phosphofructose kinase pathway was reduced and that through the oxidative pentose phosphate (PP) pathway increased in the mutant. This was evidenced by the corresponding enzyme activities, and the increase in the concentrations of phosphoenol pyruvate, glucose-6-phosphate and 6-phosphogluconate, etc. It was also found for continuous cultivation that the enzyme activities of the oxidative PP and Entner-Doudoroff pathways increased as the dilution rate increased for the pykF^– mutant. To clarify the metabolism quantitatively, it was found to be quite important to integrate the information on intracellular metabolic flux distribution, enzyme activities and intracellular metabolite concentrations.