BOL: Related items

A-allele of SLC24A5 gene is found to be responsible for variation in skin color of South-East Asians and Europeans

Rahul Agarwal — Tue, 12 Nov 2013 21:02:27 -0600

Key finding:

rs1426654 SNP of SLC24A5 gene is decider of skin pigmentation variation in South Asia
rs1426654-A allele is widely spread throughout the Indian subcontinent
Skin pigmentation is also account by the combination of processes like selection and demographic history of populations affected by their language and origin
Sign of positive selection in Europeans, Middle East, Pakistan, Central Asia and North India but not in South India
In European , A-allele is almost reached to fixation

Paper:

http://www.plosgenetics.org/article/info%3Adoi%2F10.1371%2Fjournal.pgen.1003912

Your stress/depression came from ancestor

Rahul Agarwal — Sun, 04 May 2014 18:46:20 -0500

"A study published in Nature Neuroscience finds that stress in early life alters the production of small RNAs, called microRNAs, in the sperm of mice. The mice show depressive behaviours that persist in their progeny."

Source:

http://www.nature.com/news/sperm-rna-carries-marks-of-trauma-1.15049

GeneProf: analysis of high-throughput sequencing experiment

Jit — Wed, 05 Jul 2017 16:47:05 -0500

GeneProf is a web-based, graphical software suite that allows users to analyse data produced using high-throughput sequencing platforms (RNA-seq and ChIP-seq; "Next-Generation Sequencing" or NGS): Next-gen analysis for next-gen data!

Some of GeneProf's highlights include:

Easy-to-use web-based interface:Access your data at any time from any computer with a working internet connection -- no need to install software! (cp. Section 'System Requirements').
Analysis wizards make your life easy:Step-by-step workflows make it easy to analyse high-throughput data within a minimum of hands-on time. (cp. SubConcept 'Analysis Wizards').
Versatile modules:Advanced users and data analysis experts benefit from GeneProf's broad range of analysis modules, which can be combined freely into sophisticated workflows (cp. Concept 'Workflows').
Integrated Analysis:Analysis of ChIP-seq and RNA-seq data in one place, plus support for the integration of other external data (e.g. from microarrays).
Comprehensive Resource:GeneProf provides a comprehensive resource of fully analyzed next-generation sequencing data. Experimental results can be easily accessed and compared and the analysis procedures employed to produce the data are fully transparent (cp. Tutorial 'Examining Public Next-Gen Data..').
Extensibility:Algorithm developers and computer programmers can develop their own modules and extend GeneProf. Existing software can be easily wrapped in the workflow framework (cp. Section 'Module Development: Adding new..') and data from GeneProf may be used externally (cp. Section 'Web API: Retrieving Data from ..').

GeneProf is academic software developed at the Centre for Regenerative Medicine / Institute for Stem Cell Research, University of Edinburgh and has benefited from funding by the Medical Research Council and the EU Framework 7 Project "EuroSyStem".

Address of the bookmark: https://www.geneprof.org/GeneProf/index.jsp

INC-Seq: accurate single molecule reads using nanopore sequencing

Jit — Mon, 27 Nov 2017 10:38:56 -0600

INC-Seq reads enabled accurate species-level classification, identification of species at 0.1 % abundance and robust quantification of relative abundances, providing a cheap and effective approach for pathogen detection and microbiome profiling on the MinION system.

Address of the bookmark: https://github.com/CSB5/INC-Seq

HISAT2: a fast and sensitive alignment program for mapping next-generation sequencing reads

Rahul Nayak — Tue, 08 May 2018 04:27:22 -0500

HISAT2 is a fast and sensitive alignment program for mapping next-generation sequencing reads (both DNA and RNA) to a population of human genomes (as well as to a single reference genome). Based on an extension of BWT for graphs [Sirén et al. 2014], we designed and implemented a graph FM index (GFM), an original approach and its first implementation to the best of our knowledge. In addition to using one global GFM index that represents a population of human genomes, HISAT2 uses a large set of small GFM indexes that collectively cover the whole genome (each index representing a genomic region of 56 Kbp, with 55,000 indexes needed to cover the human population). These small indexes (called local indexes), combined with several alignment strategies, enable rapid and accurate alignment of sequencing reads. This new indexing scheme is called a Hierarchical Graph FM index (HGFM).

more at https://ccb.jhu.edu/software/hisat2/index.shtml

Address of the bookmark: https://github.com/infphilo/hisat2

BlasR Mapping single molecule sequencing reads using Basic Local Alignment with Successive Refinement (BLASR): Theory and Application,

Jit — Wed, 23 May 2018 06:54:32 -0500

BLASR (Basic Local Alignment with Successive Refinement) for mapping Single Molecule Sequencing (SMS) reads that are thousands to tens of thousands of bases long with divergence between the read and genome dominated by insertion and deletion error.

Here is how I use the blasr to align PacBio reads to the contigs (target.fasta). The “target.fasta.sa” is the suffix array from “target.fasta” generated by sawriter.

blasr query.fa ./target.fasta -sa ./target.fasta.sa -bestn 40 -maxScore -500 -m 4 -nproc 24 -out target.m4 -maxLCPLength 15

the output format option “-m 4″ generate the alignment coordinate. Not fully documented, but I can explain that to you.

I use a 24 cores / 48G ram server for the alignment. It took about 2 to 3 hours aligning 3G PacBio Reads to 10^6 sequences of short read contigs with a mean 3.5kbp length.

Address of the bookmark: http://bix.ucsd.edu/projects/blasr/

ChopStitch: exon annotation and splice graph construction using transcriptome assembly and whole genome sequencing data

Rahul Nayak — Tue, 03 Jul 2018 04:14:52 -0500

ChopStitch is a new method for finding putative exons and constructing splice graphs using an assembled transcriptome and whole genome shotgun sequencing (WGSS) data. ChopStitch identifies exon-exon boundaries in de novo assembled RNA-seq data with the help of a Bloom filter that represents the k-mer spectrum of WGSS reads. The algorithm also detects base substitutions in transcript sequences corresponding to sequencing or assembly errors, haplotype variations, or putative RNA editing events. The primary output of our tool is a FASTA file containing putative exons. Further, exon edges are interrogated for alternative exon-exon boundaries to detect transcript isoforms, which are reported as splice graphs in dot output format.

Address of the bookmark: https://github.com/bcgsc/ChopStitch

PureCN: copy number calling and SNV classification using targeted short read sequencing

Jit — Thu, 09 Aug 2018 04:09:37 -0500

This package estimates tumor purity, copy number, and loss of heterozygosity (LOH), and classifies single nucleotide variants (SNVs) by somatic status and clonality. PureCN is designed for targeted short read sequencing data, integrates well with standard somatic variant detection and copy number pipelines, and has support for tumor samples without matching normal samples.

Author: Markus Riester [aut, cre], Angad P. Singh [aut]

Maintainer: Markus Riester

Citation (from within R, enter citation("PureCN")):

Riester M, Singh A, Brannon A, Yu K, Campbell C, Chiang D, Morrissey M (2016). “PureCN: Copy number calling and SNV classification using targeted short read sequencing.” Source Code for Biology and Medicine, 11, 13. doi: 10.1186/s13029-016-0060-z.

Address of the bookmark: http://bioconductor.org/packages/release/bioc/html/PureCN.html

Qualimap2: Evaluating next generation sequencing alignment data

Jit — Tue, 11 Sep 2018 04:44:29 -0500

Qualimap 2 is a platform-independent application written in Java and R that provides both a Graphical User Inteface (GUI) and a command-line interface to facilitate the quality control of alignment sequencing data and its derivatives like feature counts.

Supported types of experiments include:

Whole-genome sequencing
Whole-exome sequencing
RNA-seq (speical mode available)
ChIP-seq

Address of the bookmark: http://qualimap.bioinfo.cipf.es/

AMStat: display statistics of large sequence files from next generation sequencing projects

Neel — Fri, 09 Nov 2018 13:34:56 -0600

SAMStat is an efficient C program to quickly display statistics of large sequence files from next generation sequencing projects. When applied to SAM/BAM files all statistics are reported for unmapped, poorly and accurately mapped reads separately. This allows for identification of a variety of problems, such as remaining linker and adaptor sequences, causing poor mapping. Apart from this SAMStat can be used to verify individual processing steps in large analysis pipelines.

Address of the bookmark: http://samstat.sourceforge.net/