Skip to main content

Exome or Whole Genome Sequencing?

This is one of the most discussed topics especially in Clinical Genomics! Affordability, accuracy, feasibility and of course time consumption - based on these factor mostly, which sequencing technology is more suitable for clinics? Whole Exome Sequencing or Whole Genome Sequencing? (WGS or WES, WGS vs WES) So here's my 2 cents on this discussion! When it comes to DNA sequencing there has always been a raging debate over the choice of Whole Genome Sequencing (WGS) or Whole Exome Sequencing (WXS) for routine use. Whole genome sequencing (WGS), as the name suggests is the process of obtaining the entire genome. In most cases however, this is far from practical and only 95-97% of the genome is covered because it is technically difficult to sequence certain regions of the genome (high GC content, large repeat regions, centromeres, telomeres, etc.) with existing technology. “It’s very fair to say the human genome was never fully sequenced,”  - Craig Venter “The human genome ha...
Post a Comment

Installing Sleuth for RNA-seq analysis

So, I tried installing Sleuth along with Kallisto for RNA-seq data analysis. Getting Kallisto to work was no big deal but Sleuth required a couple of dependencies even before installation.

Following instructions on github and on other online tutorials which seemed pretty straightforward as in any R-Bioconductor based installation, I expected this to go smooth. Much to my surprise, it wasn't a cake walk!

Sleuth requires us to install 'devtools' and 'rhdf5' from within biocLite in R (I think we could google a little about what these packages do).

I had a fresh 64-bit BioLinux-8 installation which is basically a Linux Ubuntu 14.04 based OS, comes bundled with many preconfigured bioinformatics packages including a R version 3.2.0.

Initially running biocLite('rhdf5') didn't throw any error and went all the way with just a couple of minor version warnings but biocLite('devtools') threw numerous errors and just wouldn't work.

After a few minutes of googling, I ended up upgrading the R-version to the latest using the steps here.

But biocLite('devtools') simply wouldn't want to run error-free. Whatever I try, installing individual dependencies manually or try to start R as super user, it simply wouldn't give in.

At one point, there was a warning that said too many errors, terminating installation.
The one big error that kicked off the cascade was

dpkg: error processing r-base-core (--configure): causing build error


I googled again with search terms like "devtools failed to process build dependencies" but without much luck.. and finally I found the answer on this stackoverflow thread.

Basically the best way to go about this was to get out of R on to normal terminal and then update the devtools through gnutls using the usual apt-get command.

Here are the commands to solve the devtools problem.

apt-get -y build-dep libcurl4-gnutls-dev
apt-get -y install libcurl4-gnutls-dev
Once this is done, I went back to R and ran the commands

biocLite('rhdf5')
biocLite('devtools')
biocLite('pachterlab/sleuth')

Worked like magic!! Now I have Sleuth running smoothly! yayy!


_________________________________________________________________________________
Venkatesh Chellappa







Comments

Post a Comment

Popular posts from this blog

Clinical Genomics - what is it?

What is Clinical Genomics? What's the difference between precision medicine and personalised medicine? Genomics is the study of an organism's entire genome (all the DNA for the discussion lets say)! Clinical Genomics is the study of: – the correlation between diseases and human genome – the interaction of medication and human genome – developing disease diagnostic methods from human genome - Finally - coming up with ways to treat that disease based on all this knowledge built with the help of genomics. When we think about genomics in treating a disease like cancer, we consider two aspects. First is the patient’s own genome, which is in every cell in their body. In their genome, they may have mutations that make them more susceptible to a particular type of cancer than someone who doesn’t have that particular type of mutation. For example, we know that individuals with specific inherited mutations in BRCA1 and BRCA2 genes are at a greater risk for bre...
Post a Comment
Read more

Phobias: How do I get rid of this irrational fear of dogs?

Dogs are lovely creatures & are adorable. No doubts about that. But if there is one in close radius, I have an anxiety attack. I do everything to avoid getting near them. Also have dreams of dogs chasing me! How do I get rid of this fear? My partner has the same thing for dogs..Her fear is inversely proportional to the size of the dog - the smaller the dog the more petrified she gets. I am a Veterinarian and to me it is irrational but I would not say it is irrational for her or anyone to have such phobias as each of them have a strong underlying reason embedded in our brain's circuitry.  Here are some things you could do. 1. Visit a dog shelter nearby and try to go near the cages and let the dogs sniff your hand, you would be surprised how those poor doggies will be timid and shy and sometimes even afraid of You. Once you have passed level one, u could perhaps try and stroke their neck or shoulder gently. You will automatically see the tail wagging as a sig...
1 comment
Read more

How to upgrade R version on Ubuntu?

So we all have spent some quality time googling for an easy way to upgrade our R installation to the latest version on Ubuntu! Else why would you be here?! right?! Here are some quick steps to upgrade your R to latest version!  R-installations are tied up to their Ubuntu versions and hence are connected to the Ubuntu Distribution name like Trusty Tahr or Precise Pangolin. The trick here is to send a message to the repository with your latest Ubuntu distribution name, so it will pick up the latest version and update whatever is the new one out!  Just follow these steps, (simply copy-paste these commands in sequence, it will work, fingers crossed!)   1. Get the Ubuntu distribution name matching your current Version of ubuntu  distname=$(lsb_release -c -s)   2. Select the nearest CRAN mirror to you  echo "deb http://mirrors.ebi.ac.uk/CRAN/bin/linux/ubuntu $distname/" | sudo tee -a /etc/apt/sources.list > /dev/null  3. Authenticate do...
Post a Comment
Read more