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

WALK – IN – INTERVIEW @ Agricultural Knowledge Management Unit Indian Agricultural Research Institute, New Delhi-110012

Fri, 12 Dec 2014 21:33:02 -0600

Walk-in-interview for the following temporary positions will be conducted on 20th December 2014 (between 10:00 AM to 01:00 PM) at Agricultural Knowledge Management Unit, A0 block (Ground Floor), LBS Building, Indian Agricultural Research Institute, New Delhi-110012:

1 Dr. A.K.Mishra Coordinator & PI (BTISnet)

Traineeship (two) for one year

Rs. 5000/- (consolidated)

M.Sc. (Bioinformatics) with 60 % marks from a recognized University

20-12-2014 (10:00 AM -11:00 AM)

Studentship (four) for one year

Rs. 2500/- (consolidated)

Final year M.Sc./ M.Tech (Bioinformatics) Students from a recognized University

20-12-2014 (11:00 AM- 1:00 PM)

The positions are purely temporary and co-terminus with the DBT Programme. Eligible candidates are requested to submit the application form in the prescribed format along with original certificates/ documents (Degree, Marks sheets, Work experience, if any) at the time of interview. No TA/DA will be paid. Maximum age limit is 28 years for all positions. Age relaxation of 5 yrs for SC/ST and woman candidates and 3 years for OBC candidates will be given. Canvassing in any form invites disqualification.

Advertisement: http://www.iari.res.in/files/BIC-08122014-20141208-172344.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

Studentship at Nagaland University

Tue, 16 Dec 2014 01:35:00 -0600

Nagaland University
(A Central University Estd. By the Act of Parliament No. 35 of 1989)
Lumami 798 627, Nagaland
DBT Sponsored ‘Bioinformatics Infrastructure Facility’ Centre

Applications in plain paper are invited for the posts of (1) Traineeship (2 Nos.) and (2) Studentship – (2 Nos.) in the DBT funded-Bioinformatics Infrastructure Facility (BIF), Nagaland University, Lumami-798627, Nagaland. Details are given below. Interested candidates may submit the application along with self attested copies of certificates in support of the candidature to Prof. Chitta Ranjan Deb, Coordinator or Dr. L. N. Kakati, Deputy Coordinator, BIF Centre, Nagaland University, Lumami-798627, Nagaland on or before 15th January 2015.

The scanned application with relevant documents may be sent by email attachment to bifnulumami@gmail.com. Shortlisted candidates will be informed by email if called for interview. No TA/DA is admissible for attending the interview.

Traineeship (Two nos.)

Post Graduate degree in any branch of Life Sciences from UGC recognized Universities

Knowledge of computers and bioinformatics

Rs.8000/- p.m. fixed.

6 months

Studentship (Two nos.)

Pursuing Post Graduate degree in any branch of Life Sciences from UGC recognized Universities

Knowledge of computers and bioinformatics

Rs.8000/- p.m. fixed.

6 months

Terms and Conditions:

i) Applicants need to produce all original documents if call for interview.
ii) The posts are purely temporary and the appointment does not confer any entitlement or right over the job and will not be considered as formal service.
iii) No TA & DA will be paid for appearing in the walk-in-interview.
iv) The stipend/salary amount is subject to the sanction of DBT, New Delhi.

Advertisement: http://www.nagauniv.org.in/files/BIF%20Advt.pdf

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

Senior Scientist (Agricultural Bio-informatics) (One post)

Wed, 24 Dec 2014 05:01:02 -0600

Agricultural Scientists Recruitment Board

Qualifications Essential:
Doctoral degree in Bio-informatics/ Biotechnology/Molecular Biology and Biotechnology/Life Sciences/ Computer Sciences with specialization in Bio-informatics including relevant basic sciences with 8 years’ experience in the relevant subject as Scientist/Lecturer or in an equivalent position in the Pay Band- 3 of 15600-39100 with Grade Pay of 5400/6000/7000/8000 having made contribution to research/teachin/extension education as evidenced by published work/innovations and impact.

Doctoral degree in the above subject(s) including relevant basic sciences with minimum 8 years’ experience of high quality post-doctoral research in an institution/organization as evidenced by at least 6 publications in journals with NAAS rating of 7.5 or above.

Desirable:
(i) Specialization and experience: Knowledge of software development and its application in crop bioinformatics, experience in handling ‘omies’data.
(ii) Teaching experience in relevant subject

The Secretary, Agricultural Scientists Recruitment Board, Indian Council of Agricultural Research Krishi Anusandhan Bhavan-I, Pusa New Delhi –110 012, India

Detailed eligibility criteria with how to apply information can be had at:
http://www.indiastudychannel.com/jobs/333467-Indian-Council-of-Agricultural-Research-looking-for-Assistant-Director-General.aspx

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

RStudio

Jitendra Narayan — Sat, 27 Dec 2014 06:50:58 -0600

RStudio IDE is a powerful and productive user interface for R. It’s free and open source, and works great on Windows, Mac, and Linux.

The developers and expert trainers are the authors of several popular R packages, including ggplot2, plyr, lubridate, and others.

More at http://www.rstudio.com/

http://www.rstudio.com/products/RStudio/

Address of the bookmark: http://www.rstudio.com/

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

Smudgeplot: Inference of ploidy and heterozygosity structure using whole genome sequencing data

Neel — Fri, 25 Feb 2022 04:42:09 -0600

This tool extracts heterozygous kmer pairs from kmer count databases and performs gymnastics with them. We are able to disentangle genome structure by comparing the sum of kmer pair coverages (CovA + CovB) to their relative coverage (CovB / (CovA + CovB)). Such an approach also allows us to analyze obscure genomes with duplications, various ploidy levels, etc.

Smudgeplots are computed from raw or even better from trimmed reads and show the haplotype structure using heterozygous kmer pairs. For example:

Address of the bookmark: https://github.com/KamilSJaron/smudgeplot