BOL: Related items

KAD: Assessing genome assemblies using K-mer copies in assemblies and K-mer abundance in Illumina reads

Jit — Fri, 19 Jun 2020 07:34:12 -0500

KAD is designed for evaluating the accuracy of nucleotide base quality of genome assemblies. Briefly, abundance of k-mers are quantified for both sequencing reads and assembly sequences. Comparison of the two values results in a single value per k-mer, K-mer Abundance Difference (KAD), which indicates how well the assembly matches read data for each k-mer.

where, c is the count of a k-mer from reads, m is the mode of counts of read k-mers, and n is the copy of the k-mer in the assembly.

Address of the bookmark: https://github.com/liu3zhenlab/KAD

Bioinformatics Scientist at ICAR Labs

Sun, 21 Dec 2014 23:47:03 -0600

AGRICUL AGRICULTURAL SCIENTISTS RECRUITMENT BOARD TURAL SCIENTISTS RECRUITMENT BOARD
KRISHI ANUSANDHAN BHAVAN-I, PUSA, NEW DELHI-110 012

ADVERTISEMENT NO. 03/2014

PRINCIPAL SCIENTIST

Pay Band: Minimum pay of `43,000 in the PB-4 of `37400-67000/- + RGP of `10,000/-.

Age: The candidates must not have attained the age of 52 years as on 19.01.2015. There shall be no age limit for the Council’s employees.

ICAR-Indian Institute for Agricultural Biotechnology, (IIAB) Ranchi (Jharkhand)

151. Principal Scientist (Bioinformatics) (One post)

SENIOR SCIENTIST/PROGRAMME COORDINATOR

Pay Band: PB-4 of ` 37400-67000/- + RGP of ` 9,000/-.

Age: The candidates must not have attained the age of 47 years as on 19.01.2015. There shall be no age limit for the Council’s employees.

National Institute of Biotic Stress Management, Raipur (Chhattishgarh)

166. Senior Scientist (Bioinformatics) (One post)

IMPORTANT NOTE
I. (i) CLOSING DATE

THE CLOSING DATE FOR RECEIPT OF APPLICATIONS IN AGRICULTURAL SCIENTISTS RECRUITMENT BOARD IS 19.01.2015 (For applications posted from abroad and in the Andaman and Nicobar Islands, Lakshdweep, Minicoy and Amindivi islands, States/ Union Territories in the North-Eastern Region, Ladakh Division of J & K State, Sikkim, Pangi, Sub-division of Chamba, Lahul and Spiti Districts of Himachal Pradesh, the last date for receipt of application will be 02.02.2015). Non receipt of the application by the closing date will result in rejection of the application.

More Info: http://asrb.org.in/administrator/uploads_dir/1418978057english.pdf

Ancient whole genome duplication (WGD) detection tools !

Rahul Nayak — Sun, 07 Mar 2021 00:32:44 -0600

There are two methods for ancient WGD detection, one is collinearity analysis, and the other is based on the Ks distribution map. Among them, Ks is defined as the average number of synonymous substitutions at each synonymous site, and there is also a Ka corresponding to it, which refers to the average number of non-synonymous substitutions at each non-synonymous site.

At present, some people have posted articles about the analysis process of WGD. I searched for the keyword "wgd pipeline" and found the following:

GenoDup: https:// github.com/MaoYafei/GenoDup-Pipeline
https://peerj.com/articles/6303/
WGDdetector: https:// github.com/yongzhiyang2 012/WGDdetector
https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-019-2670-3
wgd: https:// github.com/arzwa/wgd
https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1142-2#Sec1
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-017-0399-x
GeNoGAP https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1142-2
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-017-0399-x
https://github.com/dfguan/purge_dups
https://www.biorxiv.org/content/10.1101/2020.01.24.917997v1

This article introduces the usage of wgd.

Wgd cannot be installed directly with bioconda at present, so it is a little troublesome to install, because it depends on a lot of software. wgd depends on the following software

BLAST
MCL
MUSCLE/MAFFT/PRANK
PAML
PhyML/FastTree
i-ADHoRe

But the good news is that most of the software it depends on can be installed with bioconda

conda create -n wgd python=3.5 blast mcl muscle mafft prank paml fasttree cmake libpng mpi=1.0=mpich
conda activate wgd

Here mpi=1.0=mpich is selected, because i-adhore depends on mpich. If openmpi is installed, an error will appear while loading shared libraries: libmpi_cxx.so.40: cannot open shared object file: No such file or directory

After that, the installation is much simpler

git clone https://github.com/arzwa/wgd.git
cd wgd
pip install .
pip install git+https://github.com/arzwa/wgd.git
For i-ADHoRe, you need to register at http:// bioinformatics.psb.ugent.be /webtools/i-adhore/licensing/Agree to the license to download i-ADHoRe-3.0

Since my miniconda3 installed ~/opt/, the installation path is so~/opt/miniconda3/envs/wgd/

tar -zxvf i-adhore-3.0.01.tar.gz
cd i-adhore-3.0.01
mkdir -p build && cd build
cmake .. -DCMAKE_INSTALL_PREFIX=~/opt/miniconda3/envs/wgd/
make -j 4
make insatall

Take the sugarcane genome Saccharum spontaneum L as an example. The genome is 8-ploid with 32 chromosomes (2n = 4x8 = 32)

Download the tutorial for CDS and GFF annotation files

mkdir -p wgd_tutorial && cd wgd_tutorial
wget http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Sspon.v20190103.cds.fasta.gz
wget http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Sspon.v20190103.gff3.gz
gunzip *.gz

First conda activate wgdstart our analysis environment, and then start the analysis

Step 1 : Use to wgd mclidentify homologous genes in the genome

wgd mcl -n 20 --cds --mcl -s Sspon.v20190103.cds.fasta -o Sspon_cds.out

Step 2 : Use to wgd ksdbuild Ks distribution

wgd ksd --n_threads 80 Sspon_cds.out/Sspon.v20190103.cds.fasta.blast.tsv.mcl Sspon.v20190103.cds.fasta

Step 3 : If the quality of the genome is good, then wgd syncollinearity analysis can be used . It can help us find the collinearity block in the genome and the corresponding anchor point

wgd syn --feature gene --gene_attribute ID \
-ks wgd_ksd/Sspon.v20190103.cds.fasta.ks.tsv \
Sspon.v20190103.gff3 Sspon_cds.out/Sspon.v20190103.cds.fasta.blast.tsv.mcl

For more reading - There are 9 sub-modules in WGD

kde: KDE fitting to the Ks distribution
ksd: Ks distribution construction
mcl: BLASP comparison of All-vs-ALl + MCL classification analysis.
mix: Hybrid modeling of Ks distribution.
pre: preprocess the CDS file
syn: Call I-ADHoRe 3.0 to use GFF files for collinearity analysis
viz: draw histogram and density plot
wf1: Ks standard analysis procedure of the whole genome paranome (paranome), call mcl, ksd and syn
wf2: Ks standard analysis procedure of one-vs-one homologous gene (ortholog), call wcl and kSD

Calling variants in non-diploid systems

Neel — Sat, 26 Jun 2021 15:37:49 -0500

The main challenge associated with non-diploid variant calling is the difficulty in distinguishing between the sequencing noise (abundant in all NGS platforms) and true low frequency variants. Some of the early attempts to do this well have been accomplished on human mitochondrial DNA although the same approaches will work equally good on viral and bacterial genomes (Rebolledo-Jaramillo et al. 2014, Li et al. 2015).

Address of the bookmark: https://training.galaxyproject.org/training-material/topics/variant-analysis/tutorials/non-dip/tutorial.html

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).

MitoZ: a toolkit for animal mitochondrial genome assembly, annotation and visualization

Neel — Tue, 30 Nov 2021 23:23:57 -0600

MitoZ, consisting of independent modules of de novo assembly, findMitoScaf (find Mitochondrial Scaffolds), annotation and visualization, that can generate mitogenome assembly together with annotation and visualization results from HTS raw reads.

https://academic.oup.com/nar/article/47/11/e63/5377471

Address of the bookmark: https://github.com/linzhi2013/MitoZ

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/

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

Short-read assembly using Spades !

Abhimanyu Singh — Mon, 31 Jan 2022 07:18:16 -0600

If we only had Illumina reads, we could also assemble these using the tool Spades.

You can try this here, or try it later on your own data.

Get data

We will use the same Illumina data as we used above:

illumina_R1.fastq.gz: the Illumina forward reads
illumina_R2.fastq.gz: the Illumina reverse reads

Assemble

Run Spades:

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

-1 is input file of forward reads
-2 is input file of reverse reads
--careful minimizes mismatches and short indels
--cov-cutoff auto computes the coverage threshold (rather than the default setting, “off”)
-o is the output directory

Results

Move into the output directory and look at the contigs:

infoseq contigs.fasta

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