Approximate Gene Cluster Discovery Problem (AGCDP)
Several genes across different genomes are grouped together as a result of functional dependencies. These set of genes that are kept together are called gene cluster. Identifying the set of clustered genes in the literature is widely studied in the fields of Biology and Computer Science. In the field of computation, biological problems as such are modelled as optimization problems.
A model introduced by Rahman et.al. is the Approximate Gene Cluster Discovery Problem (AGCDP) where the genes are treated as integers and genomes are treated as a string of integers.
The input of AGCDP is the following.
- Set of integer strings which represent the genomes.
- Number of genes expected to be in the cluster
The solution of the problem is a set of genes which minimizes a certain score. Details on the score computation is detailed in the paper by Rahman. They also presented an integer linear programming formulation of the problem.
Randomized Computation from Jhoirene Clemente
Lecture on Randomized Computation for our Computational Complexity class.
This lecture includes a short introduction to randomized algorithms. It includes discussion on probabilistic Turing Machines (PTMs), TMs that make use of random number in its computation, as well as different complexity classes that it recognize. The complexity classes include the Bounded Error Probability Polynomial (BPP), RP, co-RP, and ZPP. Relationship of the said classes where also presented as well as their relationship to known complexity class such as P and NP.
GMAC is a user-level library that implements an Asymmetric Distributed Shared Memory model to be used by CUDA programs. An ADSM model allows CPU code to access data hosted in accelerator (GPU) memory.
GMAC is being developed by the Operating System Group at the Universitat Politecnica de Catalunya and the IMPACT Research Group at theUniversity of Illinois
To install gmac,
1. Download the latest version of gmac to date: gmac-0.0.10.tar.gz
2. Unzip gmac-0.0.10.tar.gz using the following command
gunzip gmac-0.0.10.tar.gz tar -xvf gmac-0.0.10.tar
3. Go inside the directory gmac-0.0.10 and make another directory build to separate the installation files from the source.
cd gmac-0.0.10 mkdir build
4. Go inside the directory build and run configuration file
cd build ../configure
5. Install using make all install command
sudo make all install
I am installing gmac on a 32 bit UNIX machine with gcc 4.3 and CUDA 3.0.
After installing from the Makefile I encountered the following error
/usr/bin/ld: cannot find -lcuda collect2: ld returned 1 exit status make: *** [libgmac.la] Error 1 make: Leaving directory `/home/imacbuntu/Projects/gmac-0.0.10/build/ src/gmac' make: *** [all-recursive] Error 1 make: Leaving directory `/home/imacbuntu/Projects/gmac-0.0.10/build/ src/gmac' make: *** [all] Error 2 make: Leaving directory `/home/imacbuntu/Projects/gmac-0.0.10/build/ src/gmac' make: *** [all-recursive] Error 1 make: Leaving directory `/home/imacbuntu/Projects/gmac-0.0.10/build/ src' make: *** [all] Error 2 make: Leaving directory `/home/imacbuntu/Projects/gmac-0.0.10/build/ src' make: *** [all-recursive] Error 1
The error has something to do with the library of cuda,
when I tried to link the libcuda.so in /usr/lib by performing the following command, installation of gmac became successful.
sudo ln -s /usr/lib/nvidia-current/libcuda.so /usr/lib/libcuda.so Sources: http://code.google.com/p/adsm/
Frontiers in Engineering and Innovation (Part 2)
Every great and…
Every great and deep difficulty bears in itself its own solution. It forces us to change our thinking in order to find it.
Enggineering and Science Forum
Motif Finding on GPUs using Random Projections
CLEMENTE, JHOIRENE. FINDING MOTIFS IN PARALLEL USING RANDOM PROJECTION ON GPUS (Under the direction of HENRY N. ADORNA, Ph.D.) Biological motifs are short patterns that have significant number of occurrences in the set of DNA sequences. These motifs are transcription binding sites that help regulate transcription and therefore gene expression. Detection of these patterns helps in gene function discovery and building regulatory networks. Mutations may occur at random positions of the genome and these patterns are also subject to modifications, making the problem more challenging. A variant called planted (l, d)-motif finding models the detection of these subtle motif patterns in the DNA. However, several algorithms used fail to recover most of the planted (l,d)-motifs. To address this problem, a hybrid algorithm was proposed in the literature which we will refer to as FMURP (Finding Motifs using Random Projection). It uses an initialization method called Projection to avoid being trapped in the local maxima and therefore increases the chance of getting the planted motifs. This algorithm is shown to have a high accuracy on solving motif finding and planted (l,d)-motif finding problem. This research presents a parallel algorithm and implementation of FMURP on Graphics Processing Units(GPUs) using CUDA. It also provides details on the implementation and optimizations done in GPU in order to minimize usage of space. The implementation called CUDA-FMURP was tested on randomly generated (l,d)-motif instances and is shown to have recovered majority of the planted motifs. It is also shown that CUDA-FMURP obtains a maximum speedup of 6.8 using a 512 core GPU with 2.0 compute capability.
Analysis of Gene Expression Data
University of the Philippines Manila (UPM) invited us to talk about our research in Bioinformatics. I will discuss our research on visualization of yeast gene expression data.