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

AIRVF: a filtering toolbox for precise variant calling in Ion Torrent sequencing

BioStar — Fri, 22 Dec 2017 00:31:06 -0600

AIRVF that works on flowgram, raw and mapped reads and called variants to reduce artifact-driven false variant calls. Tests on sequencing data of standard reference material showed up to ∼98% reduction of false variants when combined to conventional public pipelines and ∼48% to the in-house commercial solution, with a minimal loss of sensitivity

The program with a detailed manual is available at https://sourceforge.net/projects/airvf/

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

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/

gSearch: a fast and flexible general search tool for whole-genome sequencing

Jit — Mon, 06 Aug 2018 17:19:15 -0500

gSearch compares sequence variants in the Genome Variation Format (GVF) or Variant Call Format (VCF) with a pre-compiled annotation or with variants in other genomes. Its search algorithms are subsequently optimized and implemented in a multi-threaded manner.

Address of the bookmark: http://ml.ssu.ac.kr/gSearch/index.html

SimLoRD: A read simulator for third generation sequencing reads

Aaryan Lokwani — Wed, 22 Aug 2018 10:40:27 -0500

SimLoRD is a read simulator for third generation sequencing reads and is currently focused on the Pacific Biosciences SMRT error model.

Reads are simulated from both strands of a provided or randomly generated reference sequence.

The reference can be read from a FASTA file or randomly generated with a given GC content. It can consist of several chromosomes, whose structure is respected when drawing reads. (Simulation of genome rearrangements may be incorporated at a later stage.)
The read lengths can be determined in four ways: drawing from a log-normal distribution (typical for genomic DNA), sampling from an existing FASTQ file (typical for RNA), sampling from a a text file with integers (RNA), or using a fixed length
Quality values and number of passes depend on fragment length.
Provided subread error probabilities are modified according to number of passes
Outputs reads in FASTQ format and alignments in SAM format

Address of the bookmark: https://bitbucket.org/genomeinformatics/simlord/

sim3C: Read-pair simulation of 3C-based sequencing methodologies (HiC, Meta3C, DNase-HiC)

Jit — Tue, 13 Nov 2018 07:25:38 -0600

Required python modules

biopython
intervaltree
numpy
scipy
tqdm
PyYAML

Address of the bookmark: https://github.com/cerebis/sim3C

Katuali is a flexible consensus pipeline implemented in Snakemake to basecall, assemble, and polish Oxford Nanopore Technologies' sequencing data

Jit — Tue, 22 Jan 2019 06:26:55 -0600

Run a pipeline processing fast5s to a consensus in a single command.
Recommended fixed "standard" and "fast" pipelines.
Interchange basecaller, assembler, and consensus components of the pipelines simply by changing the target filepath.
Seemless distribution of tasks over local or distributed compute.
Highly configurable.
Open source (Mozilla Public License 2.0).

Documentation can be found at https://nanoporetech.github.io/katuali/.

Address of the bookmark: https://github.com/nanoporetech/katuali

DeepVariant : an analysis pipeline that uses a deep neural network to call genetic variants from next-generation DNA sequencing data.

Jit — Sat, 25 Jan 2020 13:28:09 -0600

DeepVariant is an analysis pipeline that uses a deep neural network to call genetic variants from next-generation DNA sequencing data.

DeepVariant is an analysis pipeline that uses a deep neural network to call genetic variants from next-generation DNA sequencing data. DeepVariant relies on Nucleus, a library of Python and C++ code for reading and writing data in common genomics file formats (like SAM and VCF) designed for painless integration with the TensorFlow machine learning framework.

https://ai.googleblog.com/2017/12/deepvariant-highly-accurate-genomes.html

https://www.biorxiv.org/content/10.1101/092890v6

Address of the bookmark: https://github.com/google/deepvariant