BOL: Related items

Perl One liner basics !!

Abhimanyu Singh — Sun, 24 May 2015 09:28:33 -0500

Perl has a ton of command line switches (see perldoc perlrun), but I'm just going to cover the ones you'll commonly need to debug code. The most important switch is -e, for execute (or maybe "engage" :) ). The -e switch takes a quoted string of Perl code and executes it. For example:

$ perl -e 'print "Hello, World!\n"'
Hello, World!

It's important that you use single-quotes to quote the code for -e. This usually means you can't use single-quotes within the one liner code. If you're using Windows cmd.exe or PowerShell, you must use double-quotes instead.

I'm always forgetting what Perl's predefined special variables do, and often test them at the command line with a one liner to see what they contain. For instance do you remember what $^O is?

$ perl -e 'print "$^O\n"'
linux

It's the operating system name. With that cleared up, let's see what else we can do. If you're using a relatively new Perl (5.10.0 or higher) you can use the -E switch instead of -e. This turns on some of Perl's newer features, like say, which prints a string and appends a newline to it. This saves typing and makes the code cleaner:

$ perl -E 'say "$^O"'
linux

Pretty handy! say is a nifty feature that you'll use again and again.

Frequent words problem solution by Perl

Jit — Tue, 09 Jun 2015 23:38:44 -0500

Solved with perl http://rosalind.info/problems/1a/

#Find the most frequent k-mers in a string.
#Given: A DNA string Text and an integer k.
#Return: All most frequent k-mers in Text (in any order).

use strict;
use warnings;

my $string="ACGTTGCATGTCGCATGATGCATGAGAGCT";
my $kmer=4;
my %myHash;
my $max=0;

for (my $aa=0; $aa<=(length($string)-4); $aa++) {
   my $myStr=substr $string, $aa,$kmer;
   #print "$myStr\n";
   my $km=kmerMatch ($string, $myStr, $kmer);
   if ($km > $max) { $max = $km;}
   #print "$km\t$myStr\n";
   $myHash{$myStr}=$km;

}

#Print all key which have matching values
foreach my $name (keys %myHash){
    print "$name " if $myHash{$name} == $max;
}

sub kmerMatch { #Check the exact matching kmers with sliding window
my ($string, $myStr, $kmer)=@_;
my $count=0;
for (my $aa=0; $aa<=(length($string)-4); $aa++) {
   my $myWin=substr $string, $aa,$kmer;
   if ($myWin eq $myStr) {
       #print "$myWin eq $myStr\n";
       $count++;
   }
}
return $count;
}

Ensembl 27

Jit — Tue, 16 Jun 2015 16:10:36 -0500

What is new?

Expansion of Protists and Fungi with hundreds of annotated genomes
Variation data for bread wheat, rice, Aedes aegypti, and Ixodes scapularis
Whole genome alignments for O. longistaminata and T. cacao
Non-coding RNA gene models in Bacteria
New assembly of tomato (version 2.50)
Full support for UCSC Track Hub format for hosting your own data in Ensembl

More at http://www.ensembl.info/blog/2015/06/16/ensembl-genomes-release-27-is-out/

How long have you been a bioinformatics scientist for?

Jit — Tue, 23 Jun 2015 10:55:33 -0500

Most of the researcher have been a scientist whole life, but infact they actually started paying it with at certain time.So, how long have you been in bioinformatics domain now?

OpenCPU

Rahul Nayak — Sun, 05 Jul 2015 18:34:46 -0500

OpenCPU is a system for embedded scientific computing and reproducible research. The OpenCPU server provides a reliable and interoperable HTTP API for data analysis based on R.

The OpenCPU JavaScript client library provides the most seamless integration of R and JavaScript available today.

OpenCPU uses standard R packaging to develop, ship and deploy web applications. Several open source example apps are available from Github.

Installing your own OpenCPU server is super easy and only takes a few minutes.

More at https://www.opencpu.org/

Five key traits to seek out in potential bioinformatics candidates !!!

Jit — Mon, 10 Aug 2015 12:53:50 -0500

Genomics and proteomics data are being collected in bulk, but mostly, traditional biologist don’t know what to do with it. Perhaps this is the reason why (not only this!!! ) computational biologist/bioinformatics scientists are hot commodities in the research world.

In fact, there are huge demands for expert biological data analyst. It’s a fairly new (not exactly) hot area, these bioinformatician are invaluable because they know and understand the significance of biological data for your research and how you can use it for better understanding of biological problems.

The bioinformatics can discover biological patterns and stories in genomic and proteomics data. They can develop the pipeline needed to properly collect, store and analyse it.

Once your research group is ready to make a larger investment and hire a bioinformatician to gain a competitive edge, there are several key traits to seek out in potential candidates. The best bioinformatician are:

1. Highly Skilled - programming skills, experience with the biological software and tools.

The biological data won’t illuminate much if the scientist analysing it doesn’t possess practical programming skills, experience with the biological software and tools and a thorough understanding of basic biological stuff. A solid background in mathematics and statistics is also an indispensable trait.

2. Insight - Real vision, robust understanding and deep insight.

In order to hire the best bioinformatics and computational biologist scientist for your needs, it is always recommended and mostly practiced by the recruiters, to ask each contender to write and develop a sample script/presentation based on a specific set of data you provide. Then, explore the approaches used to deal with data provided and pick up those candidates who convey real vision, robust understanding and deep insight.

3. Energetic – Curiosity to explore

Mostly natural curiosity and enthusiasm for solving big biological problems coupled with an ability to transform data into a scientific stories may place one candidate above the rest. In addition to achieve that, the bioinformatician should be agile enough to quickly modify their methods to suit changes within a particular research.

4. Researcher – Publications

Look for someone who has a keen sense and understanding of concern biological problems. You can judge it by looking at previously published papers and data. It is always recommended to have a look at GitHub and other repository for codes written by her/him.

5. Impressive communicator - Insight that can’t be expressed is worthless.

Good bioinformatics scientists are able to uncover biological patterns and are willing to explain those patterns in clear and helpful ways through thoughtful and open communication. In other words, they should must have good scientific writing skills. A computational biologis/bioinformatician should know how to present the data and tell a scientific story through numbers/images.

Cancer research database

Neel — Tue, 01 Sep 2015 17:36:31 -0500

Researchers in Andhra Pradesh have developed a database to identify genes that are common in tumours to provide their colleagues with easy access to insights into the genetic alterations in cancer.

The database, hosted at the Sri Venkateswara University (SVU) in Tirupati, will integrate information on cancer genes and markers with experimental data.

The Cancer Gene Markers Database (CGMD) is meant to help scientists better understand tumour genes and markers at a molecular level by combining data with literature on treatment regimen and recent advances in cancer therapy.

The database is free to access, and already includes 309 genes and 206 markers that correspond to 40 different human cancers. Accompanying literature comes from databases such as the United States’ National Center for Biotechnology Information and the Kyoto Encyclopedia of Genes and Genomes. It also includes experimental data from PubMed.

In a paper published last month in Nature Scientific Reports, the researchers from SVU’s department of animal biotechnology, describes the need for a database for different genes and markers along with their molecular characteristics and pathway associations.

MAKER

Jitendra Narayan — Sun, 07 Feb 2016 15:59:24 -0600

MAKER is a portable and easily configurable genome annotation pipeline.Its purpose is to allow smaller eukaryotic and prokaryotic genome projects to independently annotate their genomes and to create genome databases. MAKER identifies repeats, aligns ESTs and proteins to a genome, produces ab-initio gene predictions and automatically synthesizes these data into gene annotations having evidence-based quality values.

More at http://www.yandell-lab.org/software/maker.html

Address of the bookmark: http://www.yandell-lab.org/software/maker.html

liftover

Jitendra Narayan — Mon, 08 Feb 2016 15:45:03 -0600

Convenient conversions between genome assemblie. The liftover package makes it easy to remap genomic coordinates to a different genome assembly.

More at https://github.com/aaronwolen/liftover

https://www.bioconductor.org/help/workflows/liftOver/

Address of the bookmark: https://github.com/aaronwolen/liftover

svelter

Jitendra Prajapati — Mon, 29 Feb 2016 17:33:15 -0600

This software is designed to identify both simple and complex rearrangements from paired-end sequencing data. Users could ran it easily by just alling SVelter.py with proper parameters. It's also possible to ran it on multiple cores by calling different sub-functions separately.

More at https://github.com/mills-lab/svelter/

Address of the bookmark: https://github.com/mills-lab/svelter/