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

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 download from the CRAN mir
Post a Comment
Read more

How do I get into Veterinary School?

What advice would you give to any aspiring candidate who is seeking to become a veterinarian? As we have seen in the recent past, Boromir tends to be pessimistic about most of the things we could do and say. But do not let Boromir discourage you. Here are some tips for you to prepare yourself to start a career as a Vet. Future Vet Website is an excellent place to start if you are looking for organised information. They advise you briefly about the necessary steps while you are in High School and as a College undergraduate. They also have lots of other interesting materials to read. Otherwise here are my simple tips -  Firstly ask yourself whether you are an animal lover - there's no point in being scared of animals, hate animals or indifferent about animals if you want to be their guardian angel for the rest of your life! Secondly be sure you want to be a vet for the right reason - to care, nurture, protect and save! and not for that Maserati you always wanted to
Post a Comment
Read more

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