BOL: Related items

Upgrade R 3.0.3

Rahul Nayak — Mon, 10 Mar 2014 11:23:51 -0500

R is a free software programming language and software environment for statistical computing and graphics. The R language is widely used among statisticians and data miners for developing statistical software and data analysis. Polls and surveys of data miners are showing R's popularity has increased substantially in recent years. Recently the new version of R codename “Warm Puppy" have been released.

You can download the latest version from here http://cran.rstudio.com/ . Or, if you are using Windows, you can upgrade to the latest version using the installr package http://cran.r-project.org/web/packages/installr/ . Simply run the following code:

# installing/loading the package:
if(!require(installr)) {
install.packages("installr"); require(installr)} #load / install+load installr

updateR()

I try to keep the installr package updated and useful. If you have any suggestions or remarks on the package, you’re invited to leave a comment below.

If you use the global library system http://www.r-statistics.com/2010/04/changing-your-r-upgrading-strategy-and-the-r-code-to-do-it-on-windows/ , you can run the following in the new version of R:

source("http://www.r-statistics.com/wp-content/uploads/2010/04/upgrading-R-on-windows.r.txt")
New.R.RunMe()

Reference:

http://www.r-statistics.com/2014/03/r-3-0-3-is-released/

A Brief Bioinformatics Tutorial

Jit — Wed, 21 May 2014 12:50:09 -0500

This is about how to use a computer to find what is known about a gene of interest and also how to get new insights about it.

The tutorial is divided in three main parts:

In the Sequence part, you will see how to look efficiently for a particular protein sequence, how to blast it against the database of your choice to find homologues, how to perform a multiple alignment of the homologues you've selected and how to edit this alignment.
The Structure part is about molecular visualization, homology modeling and structural domain prediction.
In the Function part, you will be introduced to you 3 useful servers to investigate the function of a protein. i.e. finding interactors, co-expressed genes, see a phylogenetic profile, easily access papers citing your gene etc ...

During all the three parts, we will use the S. cerevisiae VPS36 protein as an example.

Address of the bookmark: http://www.mrc-lmb.cam.ac.uk/rlw/text/bioinfo_tuto/introduction.html

List of useful machine / ai learning resources !

Rahul Nayak — Tue, 30 Mar 2021 08:56:06 -0500

ML cheatsheet !

https://github.com/remicnrd/ml_cheatsheet

Visual AI / ML

https://setosa.io/ev/

Simple and efficient tools for predictive data analysis

https://scikit-learn.org/stable/

List of GOI approved peer reviewed bioinformatics and computational biology journals

Jit — Tue, 25 Apr 2017 05:03:27 -0500

Unfortunately, we now live in a world where the integrity of peer-reviewed journals is being threatened by the rise of the academic version of fake news – something many call “predatory publishing". Mostly in academic publishing world, "predatory open access publishing" is an exploitative open-access publishing business model that involves charging publication fees to authors without providing the editorial and publishing services associated with legitimate journals (open access or not).

Nearly 20% of the such journals have a flashy impact factor and quick publication time, which are quick give-aways. Interestingly, under contact address, some journal websites do not even provide any address to contact. All of this has led to the emergence of a new and dark market of deceptive publishers that exploit the concept of open access and provide channels for “scientific journal” publication with little or no peer review. For a fee, they will publish almost anything – even if the study was fatally flawed. And these journals provide a forum that can be used as a channel to publish fraudulent “advocacy research.” You can find list of certain such publishers at "Beall's List" http://beallslist.weebly.com/

Keeping all these in mind, Government of India (GOI) decided to approved certain bioinformatics and computational biology journals for your research publication.

Following are the list of GOI validated and peer reviewed bioinformatics and computational biology journals:

NOTE:Each journal details are in following order Tittle\nSource\nSubject.
Point to remember: The list of journals are NOT sorted in any ascending or descending order.

If I missed any other GOI validated bioinformatics journal, then please report me in comment section.

Open Bioinformatics Journal
Scopus
Computer Science; Engineering; Medicine

PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS
WoS
BIOLOGY & BIOCHEMISTRY

Advances and Applications in Bioinformatics and Chemistry
Scopus
Biochemistry, Genetics and Molecular Biology Chemistry; Computer Science

Advances in Bioinformatics
Scopus
Biochemistry, Genetics and Molecular Biology; Computer Science; Engineering

Applied Bioinformatics
Scopus
Agricultural and Biological Sciences; Computer Science

BIOINFORMATICS
WoS & Scopus
COMPUTER SCIENCE

Bioinformatics and Biology Insights
Scopus
Biochemistry, Genetics and Molecular Biology; Computer Science; Mathematics

BMC BIOINFORMATICS
WoS & Scopus
COMPUTER SCIENCE

BRIEFINGS IN BIOINFORMATICS
WoS & Scopus
COMPUTER SCIENCE

Computational systems bioinformatics / Life Sciences Society. Computational Systems Bioinformatics Conference
Scopus
Medicine

Current Bioinformatics
WoS & Scopus
COMPUTER SCIENCE

Current Protocols in Bioinformatics
Scopus
Biochemistry, Genetics and Molecular Biology

JOURNAL OF COMPUTATIONAL INTELLIGENCE IN BIOINFORMATICS
ICI
BIOLOGICAL SCIENCE

Journal of integrative bioinformatics
Scopus
Medicine

Journal of Proteomics and Bioinformatics
Scopus
Biochemistry, Genetics and Molecular Biology; Computer Science

Mathematical Biology and Bioinformatics
Scopus
Engineering; Mathematics

Trends in Bioinfprmatics
Scopus
Computer Science

Eurasip Journal on Bioinformatics and Systems Biology
Scopus
General; Computer Science; Mathematics; Medicine

Evolutionary Bioinformatics
WoS & Scopus
COMPUTER SCIENCE

Genomics, Proteomics and Bioinformatics
Scopus
Biochemistry, Genetics and Molecular Biology;Mathematics

IEEE/ACM Transactions on Computational Biology and Bioinformatics
Scopus
Biochemistry, Genetics and Molecular Biology;Mathematics

IEEE-ACM Transactions on Computational Biology and Bioinformatics
WoS
COMPUTER SCIENCE

International Journal of Bioinformatics Research and Application
Scopus
Biochemistry, Genetics and Molecular Biology; Medicine, Health

International Journal o f Data M ining and Bioinformatics
WoS & Scopus
COMPUTER SCIENCE

IPSJ Transactions on Bioinformatics
Scopus
Biochemistry, Genetics and Molecular Biology;Computer Science

Journal of Bioinformatics and Computational Biology
WoS & Scopus
COMPUTER SCIENCE

Journal of Clinical Bioinformatics
Scopus
Medicine

PLoS Computational Biology
WoS & Scopus
BIOLOGY & BIOCHEMISTRY

Reviews in Computational Chemistry
WoS & Scopus
CHEMISTRY

RSC Theoretical and Computational Chemistry Series
Scopus
Chemistry; Computer Science

Annual Reports in Computational Chemistry
Scopus
Chemistry; Mathematics

Computational and Structural Biotechnology Journal
Scopus
Biochemistry, Genetics and Molecular Biology; Computer Science

Computational and Theoretical Chemistry
WoS & Scopus
CHEMISTRY

COMPUTATIONAL BIOLOGY AND CHEMISTRY
WoS & Scopus
COMPUTER SCIENCE

COMPUTATIONAL CHEMISTRY
WoS
CHEMISTRY

Journal of Theoretical and Computational Chemistry
Scopus
Chemistry; Computer Science

Theoretical and Computational Chemistry
Scopus
Chemistry

Wiley Interdisciplinary Reviews: Computational Molecular Science
Scopus
Biochemistry, Genetics and Molecular Biology;Chemistry; Computer Science; Materials Science; Mathematics

Wiley Interdisciplinary Reviews- Computational Molecular Science
WoS
CHEMISTRY

Interdisciplinary sciences, computational life sciences
Scopus
Medicine

Interdisciplinary Sciences-Computational Life Science
WoS
Biology and Biochemistry

International Journal of Computational Biology and Drug Design
Scopus
Computer Science; Pharmacology, Toxicology and Pharmaceutics

Early Genome Screening: The New Health Horoscope!

LEGE — Thu, 02 Jan 2025 19:44:36 -0600

In an era where precision medicine is reshaping healthcare, genome screening is emerging as the modern equivalent of a health horoscope. It offers insights into our biological "stars," unraveling predispositions to various conditions and empowering individuals with knowledge to navigate their health journeys proactively. But how reliable is this "horoscope," and how does it impact our lives?

Understanding Genome Screening

Genome screening involves analyzing an individual's DNA to identify genetic variations that may influence health and disease susceptibility. This can range from simple single-gene tests to comprehensive whole-genome sequencing. By peering into our genetic blueprint, we can uncover risks for conditions like cancer, diabetes, cardiovascular diseases, and even rare genetic disorders.

The process is straightforward: a saliva or blood sample is collected, and advanced sequencing technologies decipher the genetic code. The results provide a personalized health map, guiding lifestyle modifications, preventive measures, or medical interventions.

A Shift from Reactive to Proactive Healthcare

Traditional healthcare often focuses on treating diseases after they manifest. Genome screening flips this model on its head, enabling a shift toward prevention and early intervention. For instance:

Cancer Risk Management: Individuals with BRCA1 or BRCA2 gene mutations can opt for enhanced screening programs or preventive surgeries to mitigate their risk of breast and ovarian cancers.
Cardiovascular Health: Genetic predispositions to conditions like familial hypercholesterolemia can prompt early cholesterol monitoring and lifestyle adjustments.
Rare Diseases: Identifying carriers of genetic disorders can aid in family planning and reduce the incidence of inherited conditions.

The Ethical and Practical Concerns

While genome screening offers incredible promise, it is not without challenges:

Accuracy and Interpretation: Genetic predisposition does not guarantee disease. Misinterpretation of results can lead to unnecessary anxiety or unwarranted medical interventions.
Privacy and Data Security: Genetic data is highly sensitive. Ensuring robust data protection measures is crucial to prevent misuse.
Accessibility and Equity: High costs and limited availability may restrict access to genome screening, exacerbating health disparities.

Balancing Science and Pseudoscience

The comparison of genome screening to horoscopes isn’t entirely unfounded. Both offer predictive insights, but the scientific foundation of genome screening distinguishes it from astrology. Unlike the alignment of celestial bodies, genetic predictions are based on rigorous data and evidence. However, the probabilistic nature of genetic predispositions underscores the importance of interpreting results in conjunction with clinical and lifestyle factors.

The Road Ahead

As genome screening becomes more affordable and integrated into routine healthcare, its potential to transform lives is immense. Policymakers, healthcare providers, and genetic counselors must collaborate to ensure ethical implementation, public awareness, and equitable access.

Imagine a future where your genetic "horoscope" is a trusted guide, not just a prediction. Early genome screening could help chart a healthier path for generations, making it a cornerstone of personalized medicine. After all, our genes might just hold the key to unlocking a future of better health and well-being.

Try R Online

Jitendra Narayan — Tue, 16 Jul 2013 06:15:11 -0500

One of the best R tutorial website, which provide an online interative interface to try and learn R language without any hassle.

Link @ http://tryr.codeschool.com/

R and Bioconductor Tutorial

Jitendra Narayan — Fri, 23 Aug 2013 08:23:59 -0500

This tutorial is intended to introduce users quickly to the basics of R, focusing on a few common tasks that biologists need to perform some basic analysis: load a table, plot some graphs, and perform some basic statistics. More extensive tutorials can be found on the project website and via bioconductor (not covered here).

You can add more tutorial links in comments if found new pages.

Address of the bookmark: http://manuals.bioinformatics.ucr.edu/home/R_BioCondManual

Next generation sequencing in R or bioconductor environment

John Parker — Mon, 02 Jun 2014 18:03:09 -0500

There are many R software and bioconductor packages for NGS data analysis, some of them are as follows

Biostrings

The Biostrings package from Bioconductor provides an advanced environment for efficient sequence management and analysis in R. It contains many speed and memory effective string containers, string matching algorithms, and other utilities, for fast manipulation of large sets of biological sequences. The objects and functions provided by Biostrings form the basis for many other sequence analysis packages. Documentation

IRanges Overview

IRanges provides the low-level infrastructure and containers for handling sets of integer ranges within Bioconductor's BioC-Seq domain. Its classes and methods provide support for many more high-level packages like GenomicRanges, ShortRead, Rsamtools, etc. Documentation

GenomicRanges Overview

The GenomicRanges package serves as the foundation for representing genomic locations within the Bioconductor project. It is built upon the IRanges infrastructure and defines three major data containers - GRanges, GRangesList and GappedAlignments - which are supporting other important BioC-Seq packages including ShortRead, Rsamtools, rtracklayer, GenomicFeatures and BSgenome. Compared to the IRanges container, the GRanges/GRangesList classes are more flexible and extensible to store additional information about sequence ranges, such as chromosome identifiers (sequence space), strand information and annotation data. Documentation

Motif Discovery

cosmo

The cosmo package allows to search a set of unaligned DNA sequences for a shared motif that may function as transcription factor binding site. The algorithm extends the popular motif discovery tool MEME (Bailey and Elkan, 1995) in that it allows the search to be supervised by specifying a set of constraints that the motif to be discovered must satisfy. Documentation

BCRANK

BCRANK is a method that takes a ranked list of genomic regions as input and outputs short DNA sequences that are overrepresented in some part of the list. The algorithm was developed for detecting transcription factor (TF) binding sites in a large number of enriched regions from high-throughput ChIP-chip or ChIP-seq experiments, but it can be applied to any ranked list of DNA sequences. Documentation

rGADEM: Documentation

MotIV: Documentation

ShortRead

The ShortRead package provides input, quality control, filtering, parsing, and manipulation functionality for short read sequences produced by high throughput sequencing technologies. While support is provided for many sequencing technologies, this package is primairly focused on Solexa/Illumina reads. Documentation

Rsamtools

Rsamtools provides functions for parsing and inspecting samtools BAM formatted binary alignment data. SAM/BAM is quickly becoming a universal standard alignment format, and is now supported by a wide variety of alignment tools. Documentation

Samtools Website
BWA (Burrows-Wheeler Alignment) Website

Additional tools for SNP analysis:

snpMatrix

BSgenome

BSgenome provides an object oriented infrastructure for interacting with a Biostring based genome sequence. BSgenome packages exist for many common genomes, and can be created to represent custom genomes. See the "How to forge a BSgenome data package" Vignette for instructions to create a new BSgenome package if a prebuilt package does not exist for your organism. Documentation

rtracklayer

rtracklayer provides an interface for exporting annotation feature data to various genome browsers and file formats (such as GFF). See the Small RNA Profiling exercise for an example of using rtracklayer to visualize alignment coverage. Documentation

biomaRt

The biomaRt package, provides an interface to a growing collection of databases implementing the BioMart software suite (http:// www.biomart.org). The package enables online retrieval of large amounts of data in a uniform way without the need to know the underlying database schemas. This data is retrieved automatically via the Internet, so it's recommended that you cache the data locally, or check versions if your code will be adversely affected by updates to these data. Documentation

ChIP-Seq Analysis Packages

Bioconductor provides various packages for analyzing and visualizing ChIP-Seq data. Only a small selection of these packages is introduced here. Additional useful introductions to this topic are: BioC ChIP-seq Case Study and BioC ChIP-Seq.

chipseq

The chipseq package combines a variety of HT-Seq packages to a pipeline for ChIP-Seq data analysis. Documentation

BayesPeak

BayesPeak is a peak calling package for identifying DNA binding sites of proteins in ChIP-Seq experiments. Its algorithm uses hidden Markov models (HMM) and Bayesian statistical methods. The following sample code introduces the identification of peaks with the BayesPeak package as well as the incorporation of read coverage information obtained by the chipseq package. Documentation [ Publication ]

PICS

The PICS package applies probabilistic inference to aligned-read ChIP-Seq data in order to identify regions bound by transcription factors. PICS identifies enriched regions by modeling local concentrations of directional reads, and uses DNA fragment length prior information to discriminate closely adjacent binding events via a Bayesian hierarchical t-mixture model. The following sample code uses the test data set from the above BayesPeak package in order to compare the results from both methods by identifying their consensus peak set. Documentation [ Publication ]

ChIPpeakAnno

The ChIPpeakAnno package provides. batch annotation of the peaks identified from either ChIP-seq or ChIP-chip experiments. It includes functions to retrieve the sequences around peaks, obtain enriched Gene Ontology (GO) terms, find the nearest gene, exon, miRNA or custom features such as most conserved elements and other transcription factor binding sites supplied by users. The package leverages the biomaRt, IRanges, Biostrings, BSgenome, GO.db, multtest and stat packages. Documentation

Additional ChIP-Seq Packages

DiffBind: Documentation

MOSAICS: Documentation

iSeq: Documentation

ChIPseqR: Documentation

ChiPsim: Documentation

CSAR: Documentation

ChIP-Seq Pipeline: PICS, rGADEM and MotIV (developer web site)

SPP: ChIP-seq processing pipeline

SPP Tutorial

MACS

SIPeS

RNA-Seq Analysis

Counting Reads that Overlap with Annotation Ranges

The GenomicRanges package provides support for importing into R short read alignment data in BAM format (via Rsamtools) and associating them with genomic feature ranges, such as exons or genes. This way one can quantify the number of reads aligning to annotated genomic regions. The package defines general purpose containers for storing genomic intervals as well as more specialized containers for storing alignments against a reference genome. The two main functions for read counting provided by this infrastructure are countOverlaps and summarizeOverlaps. For their proper usage, it is important to read the corresponding PDF manual. Documentation

Differential Gene Expression Analysis with DESeq

The DESeq package contains functions to call differentially expressed genes (DEGs) in count tables based on a model using the negative binomial distribution. It expects as input a data frame with the raw read counts per region/gene of interest (rows) for each test sample (columns). Such a count table can be imported into R or generated from BAM alignment files using the countOverlaps function as introduced above. Documentation

Differential Gene Expression Analysis with edgeR

The edgeR package uses empirical Bayes estimation and exact tests based on the negative binomial distribution to call differentially expressed genes (DEGs) in count data.

Documentation

A variety of additional R packages are available for normalizing RNA-Seq read count data and identifying differentially expressed genes (DEG):

easyRNASeq (simplifies read counting per genome feature)

DEXSeq (Inference of differential exon usage); parathyroidSE explains how to generate exon read counts in R

DEGseq

baySeq (also see: segmentSeq)

Genominator (Bullard et al. 2010)

Detection of Alternative Splice Junctions

Another utility of RNA-Seq experiments is the analysis of splice junctions. The following software suggestions provide this utility:

ERANGE
TopHat

SpliceMap

SplitSeek

DNA-Methylation Data Analysis

methylPipe
bsseq
BiSeq
Much more under BiocViews

HT-Seq Data Visualization

ggbio: ggplot2 extension for genomics data (online manual) Gviz: Plotting data and annotation information along genomic coordinates HilbertVis: Hilbert genome plots

GenomeGraphs: Plotting genomic information from Ensembl

TileQC: Flow Cell Quality Visualization

rtracklayer: R interface to genome browsers

genoPlotR: Plotting maps of genes and genomes

Genominator: Tools for storing, accessing, analyzing and visualizing genomic data.

To install all packages

source("http://bioconductor.org/biocLite.R")
biocLite()
biocLite(c("ShortRead", "Biostrings", "IRanges", "BSgenome", "rtracklayer", "biomaRt", "chipseq", "ChIPpeakAnno", "Rsamtools", "BayesPeak", "PICS", "GenomicRanges", "DESeq", "edgeR", "leeBamViews", "GenomicFeatures", "BSgenome.Celegans.UCSC.ce2"))

translate2R

Rahul Nayak — Fri, 21 Nov 2014 01:16:06 -0600

After their presentation at the international “user!” conference, data analysis specialist eoda starts the public alpha testing of translate2R. With the start of alpha testing the innovative migration solution by the company hailing from Kassel discards the working title “translateR” and takes on the final product brand name “translate2R”. translate2R is a service for the automated translation of SPSS® syntax to R code, therefore supporting data analysts with a quick and low-risk migration to R.

The manual translation of many, frequently rather complex SPSS scripts often presents itself as a tedious and error-prone task, and represents a rather large obstacle for many analysts and companies to migrate to a modern, open source data management and analysis tool like R. With translate2R this hurdle will be diminished substantially.

Find at https://service.eoda.de/translater/?lang=en

A guide for complete R beginners !

Archana Malhotra — Fri, 20 Feb 2015 23:36:46 -0600

R-language: http://www.r-project.org

BioConductor: http://www.bioconductor.org

Advantages of R

Free!
Powerful, many libraries have been created to perform application specific tasks. e.g. analysis of microarray experiments and Next-Gen sequencing (bioconductor: including Bioseq group).
Presentation quality graphics
- Save as a png, pdf or svg
History
- What you do can be saved for the next time you use R.
- Ability to turn it into an automated script to perform again and again on different data

Disadvantages

Lack of a comprehensive graphical user interface, but two do exist: However some do exist: R commander: http://socserv.mcmaster.ca/jfox/Misc/Rcmdr/ and Limma-gui (microarrays) : http://bioinf.wehi.edu.au/limmaGUI/

Preparation

(Optional) Download and save the tutorial data set from
- http://bioinformatics.knowledgeblog.org/wp-content/uploads/bioinf/kerr/data.tsv
- Start R (type R on a Linux or Mac terminal, or find the starting link from PC)

Getting More Help

Project Home page
- http://www.r-project.org/
- Check out the ‘introduction to R’, which is a much more in depth guide .
- Also R has a built-in help system (see later)

Working directory

This is the directory used to store your data and results. It is useful if it is also the directory where your input data is stored.

Mac/Linux: this is the directory where you typed in R
PC: Change using the change working directory option