CUDA compatible GPU cards as efficient hardware accelerators for Smith-Waterman sequence alignment

Background Searching for similarities in protein and DNA databases has become a routine procedure in Molecular Biology. The Smith-Waterman algorithm has been available for more than 25 years. It is based on a dynamic programming approach that explores all the possible alignments between two sequences; as a result it returns the optimal local alignment. Unfortunately, the computational cost is very high, requiring a number of operations proportional to the product of the length of two sequences. Furthermore, the exponential growth of protein and DNA databases makes the Smith-Waterman algorithm unrealistic for searching similarities in large sets of sequences. For these reasons heuristic approaches such as those implemented in FASTA and BLAST tend to be preferred, allowing faster execution times at the cost of reduced sensitivity. The main motivation of our work is to exploit the huge computational power of commonly available graphic cards, to develop high performance solutions for sequence alignment. Results In this paper we present what we believe is the fastest solution of the exact Smith-Waterman algorithm running on commodity hardware. It is implemented in the recently released CUDA programming environment by NVidia. CUDA allows direct access to the hardware primitives of the last-generation Graphics Processing Units (GPU) G80. Speeds of more than 3.5 GCUPS (Giga Cell Updates Per Second) are achieved on a workstation running two GeForce 8800 GTX. Exhaustive tests have been done to compare our implementation to SSEARCH and BLAST, running on a 3 GHz Intel Pentium IV processor. Our solution was also compared to a recently published GPU implementation and to a Single Instruction Multiple Data (SIMD) solution. These tests show that our implementation performs from 2 to 30 times faster than any other previous attempt available on commodity hardware. Conclusions The results show that graphic cards are now sufficiently advanced to be used as efficient hardware accelerators for sequence alignment. Their performance is better than any alternative available on commodity hardware platforms. The solution presented in this paper allows large scale alignments to be performed at low cost, using the exact Smith-Waterman algorithm instead of the largely adopted heuristic approaches

What traits are carried on mobile genetic elements, and why?

Although similar to any other organism, prokaryotes can transfer genes vertically from mother cell to daughter cell, they can also exchange certain genes horizontally. Genes can move within and between genomes at fast rates because of mobile genetic elements (MGEs). Although mobile elements are fundamentally self-interested entities, and thus replicate for their own gain, they frequently carry genes beneficial for their hosts and/or the neighbours of their hosts. Many genes that are carried by mobile elements code for traits that are expressed outside of the cell. Such traits are involved in bacterial sociality, such as the production of public goods, which benefit a cell's neighbours, or the production of bacteriocins, which harm a cell's neighbours. In this study we review the patterns that are emerging in the types of genes carried by mobile elements, and discuss the evolutionary and ecological conditions under which mobile elements evolve to carry their peculiar mix of parasitic, beneficial and cooperative genes

PLAST: parallel local alignment search tool for database comparison

Background: Sequence similarity searching is an important and challenging task in molecular biology and next-generation sequencing should further strengthen the need for faster algorithms to process such vast amounts of data. At the same time, the internal architecture of current microprocessors is tending towards more parallelism, leading to the use of chips with two, four and more cores integrated on the same die. The main purpose of this work was to design an effective algorithm to fit with the parallel capabilities of modern microprocessors. Results: A parallel algorithm for comparing large genomic banks and targeting middle-range computers has been developed and implemented in PLAST software. The algorithm exploits two key parallel features of existing and future microprocessors: the SIMD programming model (SSE instruction set) and the multithreading concept (multicore). Compared to multithreaded BLAST software, tests performed on an 8-processor server have shown speedup ranging from 3 to 6 with a similar level of accuracy. Conclusions: A parallel algorithmic approach driven by the knowledge of the internal microprocessor architecture allows significant speedup to be obtained while preserving standard sensitivity for similarity search problems.

Adverse Events and Clinical Correlates in Asian Patients with Atrial Fibrillation and Diabetes Mellitus: A Report from Asia Pacific Heart Rhythm Society Atrial Fibrillation Registry

Aims. To evaluate the adverse events (and its clinical correlates) in a large prospective cohort of Asian patients with atrial fibrillation (AF) and diabetes mellitus (DM). Material and Methods. We recruited patients with atrial fibrillation (AF) from the Asia-Pacific Heart Rhythm Society (APHRS) AF Registry and included those for whom the diabetic mellitus (DM) status was known. We used Cox-regression analysis to assess the 1-year risk of all-cause death, thromboembolic events, acute coronary syndrome, heart failure and major bleeding. Results. Of 4058 patients (mean age 68.5 ± 11.8 years; 34.4% females) considered for this analysis, 999 (24.6%) had DM (age 71 ± 11 years, 36.4% females). Patients with DM had higher mean CHA2DS2-VASc (2.3 ± 1.6 vs. 4.0 ± 1.5, p < 0.001) and HAS-BLED (1.3 ± 1.0 vs. 1.7 ± 1.1, p < 0.001) risk scores and were less treated with rhythm control strategies compared to patients without DM (18.7% vs. 22.0%). After 1-year of follow-up, patients with DM had higher incidence of all-cause death (4.9% vs. 2.3%, p < 0.001), cardiovascular death (1.3% vs. 0.4%, p = 0.003), and major bleeding (1.8% vs. 0.9%, p = 0.002) compared to those without DM. On Cox regression analysis, adjusted for age, sex, heart failure, coronary and peripheral artery diseases and previous thromboembolic event, DM was independently associated with a higher risk of all-cause death (HR 1.48, 95% CI 1.00–2.19), cardiovascular death (HR 2.33, 95% CI 1.01–5.40), and major bleeding (HR 1.91, 95% 1.01–3.60). On interaction analysis, the impact of DM in determining the risk of all-cause death was greater in young than in older patients (p int = 0.010). Conclusions. Given the high rates of adverse outcomes in these Asian AF patients with DM, efforts to optimize the management approach of these high-risk patients in a holistic or integrated care approach are needed

Feasibility of Tomotherapy-Based Image-Guided Radiotherapy to Reduce Aspiration Risk in Patients with Non-Laryngeal and Non-Pharyngeal Head and Neck Cancer

PURPOSE: The study aims to assess the feasibility of Tomotherapy-based image-guided radiotherapy (IGRT) to reduce the aspiration risk in patients with non-laryngeal and non-hypopharyngeal cancer. A retrospective review of 48 patients undergoing radiation for non-laryngeal and non-hypopharyngeal head and neck cancers was conducted. All patients had a modified barium swallow (MBS) prior to treatment, which was repeated one month following radiotherapy. Mean middle and inferior pharyngeal dose was recorded and correlated with the MBS results to determine aspiration risk. RESULTS: Mean pharyngeal dose was 23.2 Gy for the whole group. Two patients (4.2%) developed trace aspiration following radiotherapy which resolved with swallowing therapy. At a median follow-up of 19 months (1-48 months), all patients were able to resume normal oral feeding without aspiration. CONCLUSION AND CLINICAL RELEVANCE: IGRT may reduce the aspiration risk by decreasing the mean pharyngeal dose in the presence of large cervical lymph nodes. Further prospective studies with IGRT should be performed in patients with non-laryngeal and non-hypopharyngeal head and neck cancers to verify this hypothesis

Propagation of Respiratory Aerosols by the Vuvuzela

Vuvuzelas, the plastic blowing horns used by sports fans, recently achieved international recognition during the FIFA World Cup soccer tournament in South Africa. We hypothesised that vuvuzelas might facilitate the generation and dissemination of respiratory aerosols. To investigate the quantity and size of aerosols emitted when the instrument is played, eight healthy volunteers were asked to blow a vuvuzela. For each individual the concentration of particles in expelled air was measured using a six channel laser particle counter and the duration of blowing and velocity of air leaving the vuvuzela were recorded. To allow comparison with other activities undertaken at sports events each individual was also asked to shout and the measurements were repeated while using a paper cone to confine the exhaled air. Triplicate measurements were taken for each individual. The mean peak particle counts were 658×103 per litre for the vuvuzela and 3.7×103 per litre for shouting, representing a mean log10 difference of 2.20 (95% CI: 2.03,2.36; p<0.001). The majority (>97%) of particles captured from either the vuvuzela or shouting were between 0.5 and 5 microns in diameter. Mean peak airflows recorded for the vuvuzela and shouting were 6.1 and 1.8 litres per second respectively. We conclude that plastic blowing horns (vuvuzelas) have the capacity to propel extremely large numbers of aerosols into the atmosphere of a size able to penetrate the lower lung. Some respiratory pathogens are spread via contaminated aerosols emitted by infected persons. Further investigation is required to assess the potential of the vuvuzela to contribute to the transmission of aerosol borne diseases. We recommend, as a precautionary measure, that people with respiratory infections should be advised not to blow their vuvuzela in enclosed spaces and where there is a risk of infecting others

Measuring Global Credibility with Application to Local Sequence Alignment

Computational biology is replete with high-dimensional (high-D) discrete prediction and inference problems, including sequence alignment, RNA structure prediction, phylogenetic inference, motif finding, prediction of pathways, and model selection problems in statistical genetics. Even though prediction and inference in these settings are uncertain, little attention has been focused on the development of global measures of uncertainty. Regardless of the procedure employed to produce a prediction, when a procedure delivers a single answer, that answer is a point estimate selected from the solution ensemble, the set of all possible solutions. For high-D discrete space, these ensembles are immense, and thus there is considerable uncertainty. We recommend the use of Bayesian credibility limits to describe this uncertainty, where a (1−α)%, 0≤α≤1, credibility limit is the minimum Hamming distance radius of a hyper-sphere containing (1−α)% of the posterior distribution. Because sequence alignment is arguably the most extensively used procedure in computational biology, we employ it here to make these general concepts more concrete. The maximum similarity estimator (i.e., the alignment that maximizes the likelihood) and the centroid estimator (i.e., the alignment that minimizes the mean Hamming distance from the posterior weighted ensemble of alignments) are used to demonstrate the application of Bayesian credibility limits to alignment estimators. Application of Bayesian credibility limits to the alignment of 20 human/rodent orthologous sequence pairs and 125 orthologous sequence pairs from six Shewanella species shows that credibility limits of the alignments of promoter sequences of these species vary widely, and that centroid alignments dependably have tighter credibility limits than traditional maximum similarity alignments

Poor adherence to antibiotic prescribing guidelines in acute otitis media—obstacles, implications, and possible solutions

Many countries now have guidelines on the clinical management of acute otitis media. In almost all, the public health goal of containing acquired resistance in bacteria through reduced antibiotic prescribing is the main aim and basis for recommendations. Despite some partial short-term successes, clinical activity databases and opinion surveys suggest that such restrictive guidelines are not followed closely, so this aim is not achieved. Radical new solutions are needed to tackle irrationalities in healthcare systems which set the short-term physician–patient relationship against long-term public health. Resolving this opposition will require comprehensive policy appraisal and co-ordinated actions at many levels, not just dissemination of evidence and promotion of guidelines. The inappropriate clinical rationales that underpin non-compliance with guidelines can be questioned by evidence, but also need specific developments promoting alternative solutions, within a framework of whole-system thinking. Promising developments would be (a) physician training modules on age-appropriate analgesia and on detection plus referral of rare complications like mastoiditis, and (b) vaccination against the most common and serious bacterial pathogens

NGF and proNGF Regulate Functionally Distinct mRNAs in PC12 Cells: An Early Gene Expression Profiling

The biological activities of NGF and of its precursor proNGF are quite distinct, due to different receptor binding profiles, but little is known about how proNGF regulates gene expression. Whether proNGF is a purely pro-apoptotic molecule and/or simply a “less potent NGF” is still a matter of debate. We performed experiments to address this question, by verifying whether a proNGF specific transcriptional signature, distinct from that of NGF, could be identified. To this aim, we studied gene expression regulation by proNGF and NGF in PC12 cells incubated for 1 and 4 hours with recombinant NGF and proNGF, in its wild-type or in a furin-cleavage resistant form. mRNA expression profiles were analyzed by whole genome microarrays at early time points, in order to identify specific profiles of NGF and proNGF. Clear differences between the mRNA profiles modulated by the three neurotrophin forms were identified. NGF and proNGF modulate remarkably distinct mRNA expression patterns, with the gene expression profile regulated by NGF being significantly more complex than that by proNGF, both in terms of the total number of differentially expressed mRNAs and of the gene families involved. Moreover, while the total number of genes modulated by NGF increases dramatically with time, that by proNGFs is unchanged or reduced. We identified a subset of regulated genes that could be ascribed to a “pure proNGF” signalling, distinct from the “pure NGF” one. We also conclude that the composition of mixed NGF and proNGF samples, when the two proteins coexist, influences the profile of gene expression. Based on this comparison of the gene expression profiles regulated by NGF and its proNGF precursor, we conclude that the two proteins activate largely distinct transcriptional programs and that the ratio of NGF to proNGF in vivo can profoundly influence the pattern of regulated mRNAs