BarraCUDA – a Fast Sequence Mapping Software using Graphics Processing Units

High-throughput DNA sequencing (HTS) instruments today are capable of generating millions of sequencing reads in a short period of time, and this represents a serious challenge to current bioinformatics pipeline in processing such an enormous amount of data in a fast and economical fashion. Modern graphics cards are powerful processing units that consist of hundreds of scalar processors in parallel in order to handle the rendering of high-definition graphics in real-time. It is this computational capability that we propose to harness in order to accelerate some of the time-consuming steps in analyzing data generated by the HTS instruments. We have developed BarraCUDA, a novel sequence mapping software that utilizes the parallelism of NVIDIA CUDA graphics cards to map sequencing reads to a particular location on a reference genome. While delivering a similar mapping fidelity as other mainstream programs , BarraCUDA is a magnitude faster in mapping throughput compared to its CPU counterparts. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the mapping throughput. BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the mapping of millions of sequencing reads generated by HTS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology. BarraCUDA is currently available at http://seqbarracuda.sf.net

Image sequence segmentation using mixture-models with uniqueness and spatio-temporal consistency constraints

The use of mixture-model techniques for motion estimation and image sequence segmentation was discussed. The issues such as modeling of occlusion and uncovering, determining the relative depth of the objects in a scene, and estimating the number of objects in a scene were also investigated. The segmentation algorithm was found to be computationally demanding, but the computational requirements were reduced as the motion parameters and segmentation of the frame were initialized. The method provided a stable description, in whichthe addition and removal of objects from the description corresponded to the entry and exit of objects from the scene.

XMolecular Cloning and Sequence Analyses of cDNAs Encoding Seven C2type Lectin2like Protein Subunits from Daboia russellii siamensis

按照Promega 公司的mRNA 提取试剂盒操作手册, 从圆斑蝰蛇( Daboia russellii siamensis ) 的毒腺中提 取mRNA ; 利用RT2PCR 的方法进行体外扩增, 获得C - 型凝集素蛋白的基因, 克隆到pMD182T 载体中。随机挑 选14 个阳性克隆进行核酸测序, 获得7 个编码不同蛇毒C - 型凝集素样蛋白亚基的cDNA , 分别命名为DRS2L1 、 DRS2L2 、DRS2L3 、DRS2L4 、DRS2L5 、DRS2L6 和DRS2L7 。由基因序列推导出的氨基酸序列表明, 克隆到的7 个蛇 毒C - 型凝集素样蛋白的亚基中均有糖识别结构域存在。BLAST 分析显示, 仅有DRS2L1 的蛋白序列与目前已知 的蛇毒C - 型凝集素样蛋白的α亚基相似。序列同源性比较并结合半胱氨酸位点分析, 推测DRS2L1 和DRS2L2 可能分别是圆斑蝰蛇毒Ⅹ因子激活剂的轻链LC2 和LC1 。DRS2L3 和DRS2L4 可能是高分子量的蛇毒C - 型凝集 素样蛋白的β亚基, 而DRS2L5 和DRS2L6 可能是低分子量的蛇毒C - 型凝集素样蛋白的β亚基。DRS2L7 可能是 类似于血小板膜糖蛋白Ib 结合蛋白的β亚基。

BarraCUDA - a fast short read sequence aligner using graphics processing units.

BACKGROUND: With the maturation of next-generation DNA sequencing (NGS) technologies, the throughput of DNA sequencing reads has soared to over 600 gigabases from a single instrument run. General purpose computing on graphics processing units (GPGPU), extracts the computing power from hundreds of parallel stream processors within graphics processing cores and provides a cost-effective and energy efficient alternative to traditional high-performance computing (HPC) clusters. In this article, we describe the implementation of BarraCUDA, a GPGPU sequence alignment software that is based on BWA, to accelerate the alignment of sequencing reads generated by these instruments to a reference DNA sequence. FINDINGS: Using the NVIDIA Compute Unified Device Architecture (CUDA) software development environment, we ported the most computational-intensive alignment component of BWA to GPU to take advantage of the massive parallelism. As a result, BarraCUDA offers a magnitude of performance boost in alignment throughput when compared to a CPU core while delivering the same level of alignment fidelity. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the alignment throughput. CONCLUSIONS: BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the alignment of millions of sequencing reads generated by NGS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology.BarraCUDA is currently available from http://seqbarracuda.sf.net.

Statistical analysis of sequence features of introns with positive transcriptional regulation in yeast genes

A great deal of experimental studies have shown that many introns of eukaryotic genes function as regulators of transcription. However, comprehensive studies of this problem have not yet been conducted. After checking the transcription frequencies of some Saccharomyces cerevisiae (yeast), genes and their introns, a remarkable phenomenon was discovered that generally the introns of the genes with higher transcription frequencies are longer, and the introns of the genes with lower transcription frequencies are shorter. This suggests that the longer introns of genes with higher transcription frequencies may contain some characteristic sequence structures, which could enhance the transcription of genes. Therefore, two sets of introns of yeast genes were chosen for further study. The transcription frequencies of the first set of genes are higher (>30), and those of the second set of genes are lower (less than or equal to10). Some oligonucleotides are detected by statistically comparative analyses of the occurrence frequencies of oligonucleotides (mainly tetranucleotides and pentanucleotides), whose occurrence frequencies in the first set of introns; are significantly higher than those in the second set of introns, and are also significantly higher than those in the exons flanking the introns of the first set. Some of these extracted oligonucleotides are the same as the regulatory elements of transcription revealed by experimental analyses. Besides, the distributions of these extracted oligonucleotides in the two sets of introns and the exons show that the sequence structures of the first set of introns are favorable for transcription of genes.

Sequence context analysis of 8.2 million single nucleotide polymorphisms in the human genome

We analyzed n-mers (n=3-8) in the local environment of 8,249,446 human SNPs and compared their distribution with that in the genome reference sequences. The results revealed that the short sequences, which contained at least one CpG dinucleotide, occurred

Sequence context analysis in the mouse genome: Single nucleotide polymorphisms and CpG island sequences

A genome-wide view of sequence mutability in mice is still limited, although biologists usually assume the same scenario for mice as for humans. In this study, we examined the sequence context in the local environment of 482,528 mouse single nucleotide po

Modeling of SFR cores with Serpent-DYN3D codes sequence

DYN3D reactor dynamics nodal diffusion code was originally developed for the analysis of Light Water Reactors. In this paper, we demonstrate the feasibility of using DYN3D for modeling of fast spectrum reactors. A homogenized cross sections data library was generated using continuous energy Monte-Carlo code Serpent which provides significant modeling flexibility compared with traditional deterministic lattice transport codes and tolerable execution time. A representative sodium cooled fast reactor core was modeled with the Serpent-DYN3D code sequence and the results were compared with those produced by ERANOS code and with a 3D full core Monte-Carlo solution. Very good agreement between the codes was observed for the core integral parameters and power distribution suggesting that the DYN3D code with cross section library generated using Serpent can be reliably used for the analysis of fast reactors. © 2012 Elsevier Ltd. All rights reserved.

The complete mitochondrial DNA sequence and the phylogenetic position of Achalinus meiguensis (Reptilia: Squamata)

The mitochondrial genome complete sequence of Achalinus meiguensis was reported for the first time in the present study. The complete mitochondrial genome of A. meiguensis is 17239 bp in length and contains 13 protein-coding genes, 22 tRNA, 2 rRNA, and 2 non-coding regions (Control regions). On the basis of comparison with the other complete mitochondrial sequences reported, we explored the characteristic of structure and evolution. For example, duplication control regions independently occurred in the evolutionary history of reptiles; the pseudo-tRNA of snakes occurred in the Caenophidia; snake is shorter than other vertebrates in the length of tRNA because of the truncations of T psi C arm (less than 5 bp) and "DHU" arm. The phylogenic analysis by MP and BI analysis showed that the phylogenetic position of A. meiguensis was placed in Caenophidia as a sister group to other advanced snakes with the exclusion of Acrochordus granulatus which was rooted in the Caenophidia. Therefore we suggested that the subfamily Xenodermatinae, which contains A. meiguensis, should be raised to a family rank or higher rank. At the same time, based on the phylogenic statistic test, the tree of Bayesian was used for estimating the divergence time. The results showed that the divergence time between Henophidia and Caenophidia was 109.50 Mya; 106.18 Mya for divergence between Acrochordus granulatus and the other snakes of the Caenophidia; the divergence time of A. meiguensis was 103 Mya, and Viperidae diverged from the unilateral of Elapidae and Colubridae was 96.06 Mya.

Cloning and sequence analysis of Sox genes in a tetraploid cyprinid fish, Tor douronensis

A PCR survey for Sox genes in a young tetraploid fish Tor douronensis (Teleostei: Cyprinidae) was performed to access the evolutionary fates of important functional genes after genome duplication caused by polyploidization event. Totally 13 Sox genes were obtained in Tor douronensis, which represent SoxB, SoxC and SoxE groups. Phylogenetic analysis of Sox genes in Tor douronensis provided evidence for fish-specific genome duplication, and suggested that Sox19 might be a teleost specific Sox gene member. Sequence analysis revealed most of the nucleotide substitutions between duplicated copies of Sox genes caused by tetraploidization event or their orthologues in other species are silent substitutions. It would appear that the sequences are under purifying selective pressure, strongly suggesting that they represent functional genes and supporting selection against all null allele at either of two duplicated loci of Sox4a, Sox9a and Sox9b. Surprising variations of the intron length and similarities of two duplicated copies of Sox9a and Sox9b, suggest that Tor douronensis might be an allotetraploidy.

Inhibition of no tail (ntl) gene expression in zebrafish by external guide sequence (EGS) technique

External guide sequence (EGS) technique, a branch of ribozyme strategy, can be enticed to cleave the target mRNA by forming a tRNA-like structure. In the present study, no tail gene (ntl), a key gene participating in the formation of normal tail, was used as a target for ribonuclease (RNase) P-mediated gene disruption in zebrafish in vivo. Transient expression of pH1-m3/4 ntl-EGS or pH1-3/4 ntl-EGS produced the full no tail phenotype at long-pec stage in proportion as 24 or 35%, respectively. As is expected that the full-length ntl mRNA of embryos at 50% epiboly stage decreased relative to control when injected the embryos with 3/4 EGS or m3/4 EGS RNA. Interestingly, ntl RNA transcripts, including the cleaved by EGS and the untouched, increased. Taken together, these results indicate that EGS strategy can work in zebrafish in vivo and becomes a potential tool for degradation of targeted mRNAs.

Isolation, characterization and genome sequence of a birnavirus strain from flounder Paralichthys olivaceus in China

A birnavirus strain, Paralichthys olivaceus birnavirus (POBV), was isolated and characterized from cultured flounder in China, and its complete genomic sequence was subsequently determined. The virus could induce cytopathic effects (CPE) in four of seven fish cell lines and was resistant to chloroform, 5-iodo-2'-deoxyuridine, acid and alkaline pH, and heat treatment. Purified virus particles had a typical icosahedral shape, with a diameter of approximately 55-60 nm. The genomic segments A and B of POBV were 3,091 and 2,780 bp in length and shared many of the features of the members of the family Birnaviridae. Segment A contained two partially overlapping ORFs encoding a polyprotein, pVP2-VP4-VP3, and a nonstructural protein, VP5, while segment B had only one ORF encoding for the VP1, a viral RNA-dependent RNA polymerase (RdRp). This is the first report about a birnavirus strain from a new non-salmonid host in China and its complete genome sequence.

Molecular phylogenetics of the family Cyprinidae (Actinopterygii : Cypriniformes) as evidenced by sequence variation in the first intron of S7 ribosomal protein-coding gene: Further evidence from a nuclear gene of the systematic chaos in the family

The family Cyprinidae is the largest freshwater fish group in the world, including over 200 genera and 2100 species. The phylogenetic relationships of major clades within this family are simply poorly understood, largely because of the overwhelming diversity of the group; however, several investigators have advanced different hypotheses of relationships that pre- and post-date the use of shared-derived characters as advocated through phylogenetic systematics. As expected, most previous investigations used morphological characters. Recently, mitochondrial DNA (mtDNA) sequences and combined morphological and mtDNA investigations have been used to explore and advance our understanding of species relationships and test monophyletic groupings. Limitations of these studies include limited taxon sampling and a strict reliance upon maternally inherited mtDNA variation. The present study is the first endeavor to recover the phylogenetic relationships of the 12 previously recognized monophyletic subfamilies within the Cyprinidae using newly sequenced nuclear DNA (nDNA) for over 50 species representing members of the different previously hypothesized subfamily and family groupings within the Cyprinidae and from other cypriniform families as outgroup taxa. Hypothesized phylogenetic relationships are constructed using maximum parsimony and Basyesian analyses of 1042 sites, of which 971 sites were variable and 790 were phylogenetically informative. Using other appropriate cypriniform taxa of the families Catostomidae (Myxocyprinus asiaticus), Gyrinocheilidae (Gyrinocheilus aymonieri), and Balitoridae (Nemacheilus sp. and Beaufortia kweichotvensis) as outgroups, the Cyprinidae is resolved as a monophyletic group. Within the family the genera Raiamas, Barilius, Danio, and Rasbora, representing many of the tropical cyprinids, represent basal members of the family. All other species can be classified into variably supported and resolved monophyletic lineages, depending upon analysis, that are consistent with or correspond to Barbini and Leuciscini. The Barbini includes taxa traditionally aligned with the subfamily Cyprininae sensu previous morphological revisionary studies by Howes (Barbinae, Labeoninae, Cyprininae and Schizothoracinae). The Leuciscini includes six other subfamilies that are mainly divided into three separate lineages. The relationships among genera and subfamilies are discussed as well as the possible origins of major lineages. (c) 2008 Published by Elsevier Inc.