The Role

We are looking for a Software Engineer - Computational Biologist with 3+
years' experience in a similar role to design and implement a backend system
to support an online individualised genomics interface. This position is a
great opportunity for an ambitious, self-motivated individual to work in a
demanding, challenging and interesting role.

Position Description:
. Participate in planning, design, and implementation of Equinome back
end systems and technologies.
. Implement interfaces and management tools for back end services.
. Manage, analyse, interpret and visualise large genomics data sets.
. Work closely with scientific team to develop new features and
application enhancements
. Design, develop and manage a genomics research database.

Qualification/Experience:
. Minimum MSc in Computer Science, Genetics, Bioinformatics or in a
related field (A Ph.D qualification would be an advantage).
. Proven 3+ years of experience in similar role.
. Highly proficient in Python, SQL, MySQL.
. Excellent knowledge of mammalian genomics, bioinformatics and
statistical/population genetics.
. Hands-on experience working with large data sets.
. Experience with front-end technologies (HTML/CSS/Javascript) an
advantage.
. Experience in rapid web application development: e.g. Django.
. Knowledge or experience of Unix Scripting and R statistical
programming language would be an advantage.
. Ability to work with minimum supervision to deliver high-quality
code on time.
. Fluency in English and good written and communication skills.
. Meticulous attention to detail.

Applications should be submitted before Friday, 26 September 2014 using the
following link:
http://bit.ly/WgbhxS

Note: Full information and application procedure is available at this link:
http://bit.ly/WgbhxS

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

Eligibility:
M.Sc. in any stream of Life Sciences with minimum 55% marks.

Desirable:
Candidates having experience in Molecular Spectroscopy, Protein Folding and Bioinformatics will be preferred.

Interested candidate may appear in the walk in Interview conducted on September 16, 2014 (Tuesday) 11:00 AM in the Director's Office, Centre for Interdisciplinary Research in Basic Sciences, lamia Millia Islamia, New Delhi-110025.
Candidates are required to bring a set of Xerox copy of their recent CV and qualifying degree (certificate/mark sheet) along with original documents. NoTA/DA will be paid.

For any further information you may e-mail to: mihassan@jmLac.in

Read more at http://jmi.ac.in/upload/advertisement/jobs_cirbs_2014september8.pdf

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

(Deemed University)

Hamdard Nagar, New Delhi – 110 062

R E C R U I T M E N T

(Advertisement No. 5/2014)

Applications on prescribed form are invited for filling up the following teaching positions in the Department of Biotechnology, Faculty of Science in the university. Eligible candidates may apply on or before 30.09.2014.

1. Professor/Associate Professor - One in Pay Band of Rs. 37400-67000+ AGP Rs.10000/9000

2. Assistant Professor - Two in Pay Band of Rs. 15600-39100+ AGP Rs. 6000/-

ASSISTANT PROFESSOR – 02 (including 01 SFS)

Specialization : Bioinformatics

Qualification and Experience :

Ph.D. in Biotechnology or an allied discipline with M.Sc. in Biotechnology/ Biochemistry in the First division or equivalent grade from a recognized University/ Institute.

NET in Life Science or allied discipline in addition to the above qualification.

Experience : At least two years of Post-doctoral teaching and/or research experience in Bioinformatics or relevant field in a UGC recognized Institution of repute or international research institute/ University. Proof of research to be evidenced by publications in SCI-indexed journals of high impact factor as the first or corresponding author.

Note : University may consider exempting candidates from NET, who has been awarded Ph.D. degree from ‘A’ Grade accredited University following the procedure as notified by the UGC in its Regulations of 2009 and adopted by Jamia Hamdard.

For more information: http://www.jamiahamdard.ac.in/PDF/Online%20application%20form%20_Teaching_1.pdf
http://www.jamiahamdard.ac.in/PDF/PBAS.pdf

If you are interested in gene prediction, have a look at GenomeThreader.

http://genometools.org/pub/

Address of the bookmark: http://genometools.org/

Functional and Comparative Genomics

http://jp.senescence.info/
http://www.senescence.info/
http://www.liv.ac.uk/integrative-biology/staff/joao-de-magalhaes/

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

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

Extent of employment: 30 Hours per Week
Duration of employment: 1st of October 2015 to 30th of September 2019
Gross monthly salary and pay grade in terms of collective agreement for university staff (payable 14 times per year): B1, € 1.997,20

Responsibilities
The successful candidate (f/m) will pursue a Ph.D. project related to the interpretation of plant genome and transcriptome sequencing data from next-generation sequencing (NGS) platforms. In particular, the candidate will characterize the unexplored genome of quinoa, a crop plant of long-standing tradition in Latin America. We collaborate with research partners in Austria and abroad, and the candidate’s project will be of central importance in the context of this research network.

Required skills and qualifications
We are looking for a graduate student (f/m) with a Master’s degree in bioinformatics or in a related field, solid programming skills (e.g. developing sequence analysis tools), experience with the analysis of NGS data sets, understanding of lab methods and knowledge of genomics/transcriptomics. The group has successfully performed several projects using NGS technology. We have recently published the reference genome sequence of sugar beet (Dohm et al., Nature, 2014), a crop plant closely related to quinoa (same family, but different genus). Not yet published is a quinoa genome assembly that we have generated, and which will serve as the starting point of the candidate’s project. We are a multidisciplinary team and offer work in a lively and friendly atmosphere, and state-of-the-art computing infrastructure. We are looking forward to expanding our team by a dedicated and strongly motivated person with a distinct interest in the challenges of plant genomics.

Applications can be submitted until: 16th of August 2015

University of Natural Resources and Life Sciences Vienna seeks to increase the number of its female faculty and staff members. Therefore qualified women are strongly encouraged to apply. In case of equal qualification, female candidates will be given preference unless reasons specific to an individual male candidate tilt the balance in his favour.

Please send your job application (incl. letter of motivation, CV, summary of Master’s thesis and contact details for two referees) to Personnel department, University of Natural Resources and Life Sciences, 1190 Vienna, Peter-Jordan-Straße 70; E-Mail: kerstin.buchmueller@boku.ac.at. (Reference code: 59)

We regret that we cannot reimburse applicants travel and lodging expenses incurred as part of the selection and hiring process.

www.boku.ac.at

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