BOL: Related items

HECIL: A Hybrid Error Correction Algorithm for Long Reads with Iterative Learning

Abhimanyu Singh — Tue, 01 Jan 2019 12:01:00 -0600

HECIL—Hybrid Error Correction with Iterative Learning—a hybrid error correction framework that determines a correction policy for erroneous long reads, based on optimal combinations of decision weights obtained from short read alignments.

HECIL’s core algorithm by introducing an iterative learning paradigm that enhances the correction policy at each iteration by incorporating knowledge gathered from previous iterations via data-driven confidence metrics assigned to prior corrections.

Address of the bookmark: https://github.com/NDBL/HECIL

Which of the followings are the best place to study Bioinformatics ?

Reshma Khatun — Sun, 28 Dec 2014 00:20:30 -0600

Bioinformatics is a major growth area and qualified Bioinformaticians are in high demand. An explosion in biological data has resulted from genome projects, next generation sequencing and other 'omics' techniques. Bioinformatics provides the tools to analyse and exploit such data sets.

Can you please suggest me the best place to study bioinformatics ( Grad/PostGrad).

Protocol for De novo Genome Assembly using Illumina Reads

BioStar — Sat, 16 Jan 2021 21:42:11 -0600

In this protocol, we address and describe the de novo assembly method for small to medium-sized genomes.

What is de novo genome assembly?
The method of taking a large number of short DNA sequences and placing them back together to create a reflection of the original chromosomes from which the DNA originated relates to genome assembly. No previous knowledge of the source DNA sequence length, structure or composition is inferred by De novo genome assemblies. The DNA of the target organism is split up into millions of tiny parts and read on a sequencing computer in a genome sequencing experiment. Depending on the sequencing system used, these "reads" range from 20 to 1000 nucleotide base pairs (bp) in length. Usually, length reads of 36 - 150 bp are produced for Illumina style short read sequencing. These reads can be either “single ended” as described above or “paired end.”

Why genome assembly?
In basic research into why and how they live, as well as in applied topics, identifying the DNA sequence of an organism is useful. Awareness of a DNA sequence may be useful in virtually any biological research because of the relevance of DNA to living things. For example, it may be used in medicine to classify, diagnose and eventually improve genetic disorder therapies. Similarly, pathogens study can lead to treatments for infectious diseases.

Raw NGS data
Reads can be saved as a Fasta file as text or in a FastQ file with their attributes. FastQ is the most common read file format since this is what the Illumina sequencing pipeline creates. This will henceforth be the subject of our conversation.

In a nutshell the protocol:
Get the sequence file(s) read from the sequencing machine (s).
Look at the readings - have an idea of what you have and what the standard is like.
If required, raw data cleanup/quality trimming.
Choose an adequate parameter set for assembly.
Assemble the data into scaffolds/contigs.
Examine the assembly performance and determine the efficiency of the assembly.

Read Quality Control:
Check the qualiy with fastQC.
Script
https://bioinformaticsonline.com/snippets/view/42540/install-fastqc-using-conda

Quality trimming/cleanup of read files.
This function trims adapters, barcodes and other contaminants from the reads.
Script
https://bioinformaticsonline.com/snippets/view/42542/trimmomatic-command

Genome Assembly:
The object of this portion of the protocol is to explain the method of assembling the reads trimmed by quality into draft contigs.

spades.py -1 illumina_R1.fastq.gz -2 illumina_R2.fastq.gz --careful --cov-cutoff auto -o result_of_spades_assembly_all_illumina

A significant range of short-read assemblers are available. Everyone with strengths and disadvantages of their own.
Some of the assemblers available include:
Velvet
SOAP-denovo
MIRA
ALLPATHS

Next step is to assess the suitability and what to do with a draft package of contiguous details for the remainder of the study now. Few stuff you can note about the contigs you just created: They're the draft Contigs. Any mis-assemblies can occur.

Mis-assembly checking and assembly metric tools:
QUAST - Quality assessment tool for genome assembly http://bioinf.spbau.ru/quast
Mauve assembly metrics - http://code.google.com/p/ngopt/wiki/How_To_Score_Genome_Assemblies_with_Mauve
InGAP-SV - https://sites.google.com/site/nextgengenomics/ingap and http://ingap.sourceforge.net/
inGAP is also useful for finding structural variants between genomes from read mappings.

Genome finishing tools:
Semi-automated gap fillers:
Gap filler - http://www.baseclear.com/landingpages/basetools-a-wide-range-of-bioinformatics-solutions/gapfiller/

IMAGE (V2) - http://sourceforge.net/apps/mediawiki/image2/index.php?title=Main_Page

Genome visualisers and editors:
Artemis - http://www.sanger.ac.uk/resources/software/artemis/
IGV - http://www.broadinstitute.org/igv/

Automated and semi automated annotation tools:
Prokka - https://github.com/tseemann/prokka
RAST - http://www.nmpdr.org/FIG/wiki/view.cgi/FIG/RapidAnnotationServer
JCVI Annotation Service - http://www.jcvi.org/cms/research/projects/annotation-service/

Frequent command use for the analysis are at:

https://bioinformaticsonline.com/blog/view/38765/list-of-tools-frequently-used-while-genome-assembly
https://bioinformaticsonline.com/pages/view/42275/frequent-parameters-for-bioinformatics-tools

Roche has acquired Bina Technologies !!!

Jit — Tue, 30 Dec 2014 09:42:16 -0600

Bina Technologies is a privately held company that provides a big data platform for centralized management and processing of next generation sequencing (NGS) data for the academic and translational research markets. Bina will be integrated into the Roche Sequencing Unit, and will continue to focus on development of their innovative genomic analysis solution.

Roche has acquired Bina Technologies, a privately-owned biotech company based in California. The biotech’s first product was the Bina Box, a platform for secondary genomic analysis, sequence alignment, and variant calling, but since 2012, it has developed other products and platforms.

It is our shared vision with Roche that informatics and data sciences are critical elements of an end-to-end genomics solution. Fast, easy-to-use, scalable, and robust informatics solutions make a big difference in the quality and impact of the work of scientists and researchers. We believe in the future of data-driven, personalized medicine. We are passionate about accelerating that future together with Roche.

Financial details of the deal were not disclosed. For Roche, the move is yet another in a string of acquisitions. Last week (December 18), Roche paid $489 million for antibody maker Dutalys. And earlier this month, Roche bought prenatal testing company Ariosa Diagnostics.

Reference

http://blog.bina.com/news/bina-technologies-acquired-by-roche?&__hssc=109677338.1.1419953400266&__hstc=109677338.b8350f2729889b08f1325906d5236cd3.1419953400266.1419953400266.1419953400266.1&hsCtaTracking=96cac941-9372-4bbf-bacb-3ca6f1ff8cfd|3fce0f18-835b-4086-9345-388880861732

http://www.the-scientist.com/?articles.view/articleNo/41750/title/Roche-Buys-Bioinformatics-Firm/

Perl Special Vars Quick Reference

Abhimanyu Singh — Tue, 07 Feb 2017 05:08:47 -0600

`$_`	The default or implicit variable.
`@_`	Subroutine parameters.
`$a` `$b`	sort comparison routine variables.
`@ARGV`	The command-line args.
Regular Expressions
`$`	Regexp parenthetical capture holders.
`$&`	Last successful match (degrades performance).
`${^MATCH}`	Similar to `$&` without performance penalty. Requires /p modifier.
$`	Prematch for last successful match string (degrades performance).
`${^PREMATCH}`	Similar to $` without performance penalty. Requires `/p` modifier.
`$'`	Postmatch for last successful match string (degrades performance).
`${^POSTMATCH}`	Similar to `$'` without performance penalty. Requires `/p` modifier.
`$+`	Last paren match.
`$^N`	Last closed paren match (last submatch).
`@+`	Offsets of ends of successful submatches in scope.
`@-`	Offsets of starts of successful submatches in scope.
`%+`	Like `@+`, but for named submatches.
`%-`	Like `@-`, but for named submatches.
`$^R`	Last regexp (?{code}) result.
`${^RE_DEBUG_FLAGS}`	Current value of regexp debugging flags. See `use re 'debug';`
`${^RE_TRIE_MAXBUF}`	Control memory allocations for RE optimizations for large alternations.
Encoding
`${^ENCODING}`	The object reference to the Encode object, used to convert the source code to Unicode.
`${^OPEN}`	Internal use: \0 separated Input / Output layer information.
`${^UNICODE}`	Read-only Unicode settings.
`${^UTF8CACHE}`	State of the internal UTF-8 offset caching code.
`${^UTF8LOCALE}`	Indicates whether UTF8 locale was detected at startup.
IO and Separators
`$.`	Current line number (or record number) of most recent filehandle.
`$/`	Input record separator.
`$\|`	Output autoflush. 1=autoflush, 0=default. Applies to currently selected handle.
`$,`	Output field separator (lists)
`$\`	Output record separator.
`$"`	Output list separator. (interpolated lists)
`$;`	Subscript separator. (Use a real multidimensional array instead.)
Formats
`$%`	Page number for currently selected output channel.
`$=`	Current page length.
`$-`	Number of lines left on page.
`$~`	Format name.
`$^`	Name of top-of-page format.
`$:`	Format line break characters
`$^L`	Form feed (default "\f").
`$^A`	Format Accumulator
Status Reporting
`$?`	Child error. Status code of most recent system call or pipe.
`$!`	Operating System Error. (What just went 'bang'?)
`%!`	Error number hash
`$^E`	Extended Operating System Error (Extra error explanation).
`$@`	Eval error.
`${^CHILD_ERROR_NATIVE}`	Native status returned by the last pipe close, backtick (`` ) command, successful call to wait() or waitpid(), or from the system() operator.
ID's and Process Information
`$$`	Process ID
`$<`	Real user id of process.
`$>`	Effective user id of process.
`$(`	Real group id of process.
`$)`	Effective group id of process.
`$0`	Program name.
`$^O`	Operating System name.
Perl Status Info
`$]`	Old: Version and patch number of perl interpreter. Deprecated.
`$^C`	Current value of flag associated with -c switch.
`$^D`	Current value of debugging flags
`$^F`	Maximum system file descriptor.
`$^I`	Value of the -i (inplace edit) switch.
`$^M`	Emergency Memory pool.
`$^P`	Internal variable for debugging support.
`$^R`	Last regexp (?{code}) result.
`$^S`	Exceptions being caught. (eval)
`$^T`	Base time of program start.
`$^V`	Perl version.
`$^W`	Status of -w switch
`${^WARNING_BITS}`	Current set of warning checks enabled by `use warnings;`
`$^X`	Perl executable name.
`${^GLOBAL_PHASE}`	Current phase of the Perl interpreter.
`$^H`	Internal use only: Hook into Lexical Scoping.
`%^H`	Internaluse only: Useful to implement scoped pragmas.
`${^TAINT}`	Taint mode read-only flag.
`${^WIN32_SLOPPY_STAT}`	If true on Windows `stat()` won't try to open the file.
Command Line Args
`ARGV`	Filehandle iterates over files from command line (see also `<>`).
`$ARGV`	Name of current file when reading <>
`@ARGV`	List of command line args.
`ARGVOUT`	Output filehandle for -i switch
Miscellaneous
`@F`	Autosplit (-a mode) recipient.
`@INC`	List of library paths.
`%INC`	Keys are filenames, values are paths to modules included via `use, require,` or `do`.
`%ENV`	Hash containing current environment variables
`%SIG`	Signal handlers.
`$[`	Array and substr first element (Deprecated!).

See perlvar for detailed descriptions of each of these (and a few more) special variables.

Type Hinting

Pranjali Yadav — Fri, 09 Jan 2015 22:26:13 -0600

Python creator Guido van Rossum’s proposal for static type-checking annotations is inching closer to reality, and the feature has taken on a new name: type hinting.

Back in August, van Rossum published a proposal on the Python mailing list recommending type-checking annotations as a valuable feature for the next version of Python to improve the performance of editors and IDEs, linter capabilities, standard notation, and refactoring. Van Rossum’s latest proposal, posted late last month, outlined plans to publish a Python Enhancement Proposal (PEP) in early January to put the feature now known as type hinting on track for inclusion in Python 3.5, slated for release this September.

Reference

https://quip.com/r69HA9GhGa7J

ARC: pipeline which facilitates iterative, reference guided de novo assemblies

Jit — Thu, 26 Jul 2018 09:20:26 -0500

ARC is a pipeline which facilitates iterative, reference guided de novo assemblies with the intent of:

Reducing time in analysis and increasing accuracy of results by only considering those reads which should assemble together.
Reducing/removing reference bias as compared to mapping based approaches.

The software is designed to work in situations where a whole-genome assembly is not the objective, but rather when the researcher wishes to assemble discreet 'targets' contained within next-generation shotgun sequence data. ARC decomplexifies the traditionally difficult problem of assembly by breaking the reads into small, manageable subsets which can then be assembled quickly and efficiently in parallel. Applications include those in which the researcher wishes to de novo assemble specific content and a set of semi-similar reference targets is available to initialize the assembly process.

https://ibest.github.io/ARC/

Address of the bookmark: https://ibest.github.io/ARC/

Postdoctoral Researcher in Cancer Systems Biology

Mon, 12 Jan 2015 01:44:11 -0600

Postdoctoral Researcher in Cancer Systems Biology
Department of Oncology, Old Road Campus Research Building, Roosevelt Drive, Oxford
Grade 7: £30,434 - £37,394 with a discretionary range to £40,847 p.a.
Applications are invited for a Postdoctoral Researcher in Cancer Systems Biology to join a rapidly developing Bioinformatics Research Core group headed by Dr Anastasia Samsonova. The purpose of the role is to develop and deliver integrative approaches to dissect the complexity of cancer as a genomic disease. The research will focus on development and application of effective strategies for mining and integration of complex human *omics datasets and clinical/phenotypic data in cancer studies.

The role sits at the critical interface between genetics and cancer systems biology, and would suit a candidate who is interested in developing a career at the confluence of Statistics/Data Mining/Machine Learning and Biology. Ideally, you will have experience in development of analytical approaches to high-throughput and multivariate data mining and integration gained through a PhD (or equivalent) in a quantitative subject (eg mathematics, statistics, physics, engineering or computer science).

Experience of statistics and/or machine learning techniques is highly desirable as is evidence of prior experience of developing bioinformatics software and/or analysing complex *omics data sets. You will be able to work as part of a team and independently and deliver results to the required standard and schedule. You should be able to organise and prioritise your own work, as well as have excellent communication skills, both written and oral. The post will involve interactions with collaborators from such diverse fields as applied mathematics, statistics, computer science and medicine.

This is a full-time post, fixed-term until 31 March 2017. For informal enquiries, contact Dr Anastasia Samsonova (bioinformatics@oncology.ox.ac.uk).

All applicants must complete a short application form and upload a CV and supporting statement.

The closing date for applications is 12.00 noon on 26 January 2015.

More at https://www.recruit.ox.ac.uk/pls/hrisliverecruit/erq_jobspec_version_4.display_form

Interactive Market Intelligence

Pranjali Yadav — Mon, 19 Jan 2015 08:20:38 -0600

BioInformatics LLC, a premier research and advisory firm serving the life science industry, has launched groundbreaking, dynamic-data presentation platform, Interactive Market Intelligence— the only cloud-based market research analytics tool for the life science tools industry.

Superior to traditional PDF and PowerPoint reports, Interactive Market Intelligence allows end-users to filter, create and export literally thousands of views of data — all easily obtainable from a set of core metrics that include market, brand, customer and workflow analytics in well-defined segments of the life science market.

The Market for Real-Time PCR is the first in a series of topics to be explored using the Interactive Market Intelligence platform. The primary research analysis is based on a survey of 900+ international scientists performing qPCR in their laboratories.

Key data findings from "The Market for Real-Time PCR": Global market for qPCR in 2015 is estimated to be $3.6B; The average growth in qPCR throughput is expected to be at 9.8% in 2015; 22% of respondents are highly likely to switch primary suppliers of qPCR products; 50% of respondents use pre-designed primer/probe sets.