BOL: Related items

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

Postdoctoral fellowship in Bioinformatics

Sat, 08 Nov 2014 14:41:14 -0600

A two-year post-doctoral position is available in the Biocomputing group of the Sapienza University led by Anna Tramontano to work on either genomics research or structural bioinformatics, focusing on the study of relevant biomedical problems.
The ideal candidate should be motivated and talented, hold a PhD degree, have good programming skills, a grasp of statistical methods and an understanding of biology.
Experience in the development of computational biology methods would be an added value.

Good communication skills and fluency in spoken and written English are required.
Please apply sending a curriculum vitae, the names of at least two referees and a letter of motivation describing past experience and future goals to anna.tramontano@uniroma1.it with subject: “Application for post-doctoral position November 2014 YOUR LAST NAME”

Deadline: No later than November 28th, 2014.
Duration: 2 years

Salary on grant: Commeasured to the experience of the candidate
Contact Person (Referent): Anna Tramontano
Ref. E-Mail: anna.tramontano@uniroma1.it
Group Web Page: http:/www.biocomputing.it

Understanding kmer !

BioStar — Wed, 18 Aug 2021 04:27:51 -0500

What is a k-mer anyway? A k-mer is just a sequence of k characters in a string (or nucleotides in a DNA sequence). Now, it is important to remember that to get all k-mers from a sequence you need to get the first k characters, then move just a single character for the start of the next k-mer and so on. Effectively, this will create sequences that overlap in k-1 positions.

Address of the bookmark: https://bioinfologics.github.io/post/2018/09/17/k-mer-counting-part-i-introduction/

Centre for Systems Biology & Bioinformatics, Panjab University vacancy of Research Fellow

Wed, 12 Nov 2014 06:18:54 -0600

Applications are invited along with complete bio-data and attested copies of certificates of qualifications, experience etc. for the one post of
Research Fellow and one post of Program Assistant under PURSE Grant of the University in Centre for Systems Biology & Bioinformatics, UIEAST, Panjab University, Chandigarh which is tenable till the period of the project

Essential Qualification
For Research Fellow:-
M.Sc. in Systems Biology and Bioinformatics / Life
Sciences with minimum 55% marks.
Preference will be given to NET/GATE/ICMR qualified candidates without fellowship however, candidates who have cleared the Panjab University Ph.D. entrance test in Systems Biology & Bioinformatics will also be eligible.

Applications should be reach on or before 19-11-2014 in the office of the undersigned. Interview will be held on 21-11-2014 in the office of the Coordinator, Centre for Systems Biology & Bioinformatics, South Campus, Block-3, Sector-25, Panjab University, Chandigarh. No TA/DA will be paid.

more at http://jobs.puchd.ac.in/includes/jobs/2014/20141110143634-Advertisement.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

Bioinformatics JRF/RA/SRF position at Institute of Cytology and Preventive Oncology (ICPO)

Wed, 19 Nov 2014 20:16:32 -0600

Institute of Cytology and Preventive Oncology (ICPO) I-7, Sector-39, Noida-201301

Candidates having the below mentioned qualifications may appear for walk in interview at ICPO on 2nd December 2014 between 10.00 AM and 12:00 PM under the below time bound projects under Dr. Subhash M. Agarwal, Scientist C. The post is purely temporary and co-terminus with the project.

Research Assistant (One)
25650/- consolidated
Discovery of EGFR secondary mutant inhibitors using structure based screening approach (ICMR)
Duration: 7 months

Essential: M.Sc./ M.Tech in Bioinformatics or any other related subject with good academic record.

Desirable: Experience in scripting and molecular docking.

Below 30 years

Junior Research Fellow (One)

16,000 + 30% HRA = Rs. 20800/-

Identification of novel inhibitors targeting EGFR using an integrated ligand and structure based approach (DBT)

Duration: 9 months

Essential: M.Sc./ M.Tech in Bioinformatics or any other related subject with good academic record. Candidates with CSIR-UGC / ICMR, NET qualification will be preferred

Desirable: Experience in scripting, QSAR and molecular docking.

Below 28 years

Interested eligible candidates may send their applications with Bio-data by email at (smagarwal@gmail.com) or by post addressed to Dr. Subhash M Agarwal, Scientist C, Institute of Cytology and Preventive Oncology (ICPO) I-7, Sector-39, Noida-201301 so as to reach latest by 1st December, 2014. The candidates may appear for interview at ICPO along with 3 copies of CV, photo and relevant certificates of qualifications in original and reprints of publications at the time of interview. It should be noted that No TA/DA will be paid for the walk in Interview.

Advertisement: www.icpo.org.in/advt-walk-in-interview.docx

maftools

Surabhi Chaudhary — Fri, 17 Dec 2021 03:18:28 -0600

With advances in Cancer Genomics, Mutation Annotation Format (MAF) is being widely accepted and used to store somatic variants detected. The Cancer Genome Atlas Project has sequenced over 30 different cancers with sample size of each cancer type being over 200. Resulting data consisting of somatic variants are stored in the form of Mutation Annotation Format. This package attempts to summarize, analyze, annotate and visualize MAF files in an efficient manner from either TCGA sources or any in-house studies as long as the data is in MAF format.

https://www.bioconductor.org/packages/devel/bioc/vignettes/maftools/inst/doc/maftools.html

Address of the bookmark: https://github.com/PoisonAlien/maftools

Genome Origami

Jitendra Narayan — Fri, 12 Dec 2014 22:48:17 -0600

There are several interesting factoid about our genomes, one of them is their folding. If we stretched out the DNA in a single cell, which is only a few millionths of an inch wide, it would span more than six feet. In other word, the size of six feet DNA fold themself to fit in a few millionths of an inch wide space. These DNA folding is a dynamic process that changes over time (!!). Researchers around the world have been trying to understand how DNA folds itself up so efficiently, and a recent post on the NIH Director’s Blog highlights new research illustrating how the human genome folds inside the cell’s nucleus, as well as how DNA folding affects gene regulation. The research team created this delightful video that demonstrates the principles involved using origami art.

http://bioinformaticsonline.com/videolist/watch/19555/a-3d-map-of-the-human-genome

Researchers have been working to determine how cells regulate gene expression for nearly as long as we’ve known about DNA. How, for example, do nerve cells know to turn off only nerve cell genes and turn off bone cell genes? DNA folding loops are part of the answer. This research team, which published their findings in a paper in Cell http://www.cell.com/cell/abstract/S0092-8674%2814%2901497-4 , found that the number of loops is much lower than expected. There are only 10,000 loops instead of the predicted millions, and they form on/off switches in DNA.

More at http://www.eurekalert.org/pub_releases/2014-12/ru-3mr121114.php

Reference http://www.cell.com/cell/abstract/S0092-8674%2814%2901497-4

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

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