1.Stupid for biologist: Biology is so complex that it will make bioinformatician feel stupid. There are no any universal fixed rules; it can surprise you any time. So be nice to biologists who ask questions and resolve your biological puzzles. Sometime you will have no idea what the hell you were doing either.

2.Puzzling why: Do not hesitate to ask question. Especially. at the beginning of project you will have to ask a lot of questions. Instead of puzzling it out at end check out and clear your doubt even for a single error. It may can leads to wrong conclusion.

3.Running marathon: The most of the biological software’s documentation is always incomplete. In other word they are no more than 95 percent complete. Sometime a single problem can halt your entire project for months. Compilation and running the pipelines in tedious because almost all are interdependent and need proper configuration. I face the same kind of problem with Evolver :( …

4.Folders missing: The pipelines generate lots of data, and we keep them in several folders for future use. But sometime we delete them by mistake and move to recovery…

5.Digging deeper: Digging deeper is fruitful, but some time it can be catastrophic. You may get frustrated or direction less. So keep a biologist with you for rescue …. Sometime an expert computer programmer to handle your server. Remember, the server will always go down when you need it the most.

The most common frustrating  common line: Why do we do this again?

Why Normalize RNA-Seq Data?

Normalization is a crucial step in RNA-Seq analysis for the following reasons:

• Sequencing depth: Different RNA-Seq experiments produce varying numbers of reads, making direct comparisons between samples misleading.

• Gene length: Longer genes inherently generate more reads, irrespective of their actual expression level.

• Bias reduction: Normalization mitigates technical biases, enabling meaningful biological interpretation.

TPM (Transcripts Per Million)

TPM measures the proportion of reads mapped to a transcript, normalized by transcript length and sequencing depth. It is calculated as:

Key Features:

1. Proportionality: TPM values sum to 1,000,000 across all transcripts in a sample, making it easier to compare between samples.

2. Intuitive interpretation: TPM values directly represent the abundance of transcripts in a sample.

3. Preferred for comparisons: TPM facilitates between-sample comparisons better than FPKM.

FPKM (Fragments Per Kilobase Million)

FPKM normalizes read counts by transcript length and sequencing depth, but without enforcing proportionality like TPM. It is defined as:

Key Features:

1. Historical significance: FPKM was one of the first normalization methods used for RNA-Seq.

2. Single-end vs. paired-end: In paired-end sequencing, FPKM becomes RPKM (Reads Per Kilobase Million).

3. Limited utility: FPKM values are not as robust as TPM for cross-sample comparisons due to lack of proportionality.

CPM (Counts Per Million)

CPM normalizes raw read counts by sequencing depth, without considering gene length. It is expressed as:

Key Features:

1. Simplicity: CPM is straightforward and computationally less intensive.

2. Application: Suitable for non-length-dependent analyses, such as comparing total expression levels or differential expression analysis.

3. Gene length agnostic: CPM does not correct for gene length, making it less ideal for measuring expression levels.

When to Use Each Method

• TPM: Best for comparing expression levels between samples, especially when transcript length and sequencing depth vary.

• FPKM: Useful for historical consistency but generally replaced by TPM.

• CPM: Ideal for differential expression analysis when gene length normalization is unnecessary.

Conclusion

Choosing the right normalization method depends on the specific objectives of your RNA-Seq analysis. TPM’s proportionality and robustness make it the preferred choice for most applications, while CPM serves well for differential expression studies. Although FPKM paved the way for RNA-Seq normalization, it has largely been supplanted by TPM in modern workflows. Understanding these methods and their nuances ensures accurate and meaningful interpretations of RNA-Seq data.

References:

1. Li, B., & Dewey, C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics.

2. Trapnell, C., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology.

3. Law, C. W., et al. (2014). voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biology.

More at https://leylab.com/

Summary: The application of protein–protein docking in large-scale interactome analysis is a major challenge in structural bioinformatics and requires huge computing resources. In this work, we present MEGADOCK 4.0, an FFT-based docking software that makes extensive use of recent heterogeneous supercomputers and shows powerful, scalable performance of over 97% strong scaling.

Availability and Implementation: MEGADOCK 4.0 is written in C++ with OpenMPI and NVIDIA CUDA 5.0 (or later) and is freely available to all academic and non-profit users at: http://www.bi.cs.titech.ac.jp/megadock.

Contact: akiyama@cs.titech.ac.jp

Address of the bookmark: http://bioinformatics.oxfordjournals.org/content/early/2014/08/06/bioinformatics.btu532.short

Finally, I find the underlying R code in a file created by GrapheR. For more details read also the package vignette, which is available in English, French and German!

Abstract: Strand NGS includes comprehensive workflows for DNA-Seq, RNA-Seq, Small RNA-Seq, ChIP-Seq, MeDIP-Seq, and Methyl-Seq analysis. Each workflow includes a quality assessment and filter section, followed by a workflow-specific analysis section. The pipeline functionality in Strand NGS allows users to execute a sequence of analysis steps with specific parameters - all without any manual intervention. This simplifies the analysis in large scale sequencing projects where every sample needs to be processed identically.

In this webinar we will discuss the pre-packaged pipelines present in Strand NGS. The packaged pipelines have well-chosen default parameters and are suitable for users analyzing data for the first time in the tool. We will also show how advanced users can customize pipelines and share them with other Strand NGS users. Finally, we will show a brief glimpse of an elaborate pipeline that aligns reads, filters poor-quality matches, computes coverage metrics, identifies variants, checks for sample cross-contamination, and emails quality reports - all from within Strand NGS.

Speaker: Dr. Vamsi Veeramachaneni, Vice President - Bioinformatics, Strand Life Sciences

Details: Session 1: 2:30 PM IST, Session 2 : 10:30 PM IST
Register here: http://www.strand-ngs.com/webinar_registration

More at http://pythonhosted.org/pybedtools/

A powerful toolset for genome arithmetic.http://code.google.com/p/bedtools/

CGHub
https://cghub.ucsc.edu/
Comprehensive data repository; huge data size

EGA
https://www.ebi.ac.uk/ega/
Comprehensive data repository; huge data size

COSMIC
http://cancer.sanger.ac.uk
Largest somatic mutation database; genome sequencing paper curation

CPRG
http://www.broadinstitute.org/software/cprg
Interface for cancer program resources

GDAC
http://gdac.broadinstitute.org/
Data analysis; automatic pipelines; user-friendly reports

SNP500Cancer
http://snp500cancer.nci.nih.gov
Sequence and genotype verification of SNPs

canEvolve
www.canevolve.org/
Comprehensive analysis of tumor profile; Data from 90 studies involving more than 10,000 patients

MethyCancer
http://methycancer.psych.ac.cn
Relationship among DNA methylation, gene expression and cancer

SomamiR
http://compbio.uthsc.edu/SomamiR/
Correlation between somatic mutation and microRNA; genome-wide displaying

cBioPortal
http://www.cbioportal.org/public-portal/
Graphical summaries; gene alteration; processed data; visualization

UCSC Cancer Genomics Browser
https://genome-cancer.soe.ucsc.edu/
Clinical information; gene expression; copy number variation; visualization

CGWB
https://cgwb.nci.nih.gov/
Visualization; gene mutation and variation; automated analysis pipeline

GDSC
http://www.cancerrxgene.org
Drug sensitivity information; drug response information

canSAR
https://cansar.icr.ac.uk/
Multidisciplinary information; drug discovery

NONCODE
http://www.noncode.org/ ncRNAs;
lncRNAs; up-to-date and comprehensive resource

Address of the bookmark: https://sourceforge.net/projects/genoviz/

Essential qualifications for JRF
: M. Sc. in Bioinformatics with experience in Proteomics, genomics and structural biology. Knowledge of programming language, pearl and database – HTML, CSS,php and Java script.
Essential qualifications for Research Associate:
MSc/MTech in Bioinformatics with three years experience or Ph.D in Bioinformatics with experience in proteomics, genomics and structural biology. Knowledge of programming language, perl and database
– HTML, CSS, Java script. NGS sequence assembly and analysis and algorithm designing.
Age limit : 35 years maximum (5 year relaxation for SC/ST and women candidates)
Emoluments:
JRF: 16,000 + 30% HRA
.
Res Assoc: Rs22,000 + 30% HRA
The post is purely temporary in nature and is co-terminus with the project. The appointment would be initially for one year and may be extended further upon satisfactory performance.
Interested candidates
should send the duly filled application forms (format in the following page ) so as to reach on or before 20.9.2014 along with all the relevant documents.

More at http://www.iari.res.in/files/JRF_RA-03092014-20140903-135319.pdf