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

Candidate Requirements: A PhD in computer science, computational biology, Bioinformatics, or a related field, along with sufficient programming skills for prototyping. Experience with next generation sequencing data analysis is required. Candidates with MS degree but with relevant work experience can also be considered. The successful candidate must be motivated and capable of working independently as well as in team environment.

Eligible and interested candidates can email your resumes to rohit at strandls dot com

About Strand Life Sciences: Strand was founded in 2000 by computer science and mathematics professors who recognized the need to automate and integrate life science data analysis through an algorithmic and computational approach. Strand’s solutions for life sciences research are robust and easy to use by the most novice user while powerful and configurable for the bioinformatician. Using its award-winning application development platform, AVADIS®, Strand builds innovative products that enable fast and cutting-edge analysis for basic and clinical research, drug discovery and development.

http://www.avadis-ngs.com/careers

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

The NCBI Datasets datasets command line tools are datasets and dataformat .

Use datasets to download biological sequence data across all domains of life from NCBI.

Use dataformat to convert metadata from JSON Lines format to other formats.

Conda download:

https://anaconda.org/conda-forge/ncbi-datasets-cli

Buld Download

 https://www.ncbi.nlm.nih.gov/datasets/builder/?tax_id=29979

Address of the bookmark: https://www.ncbi.nlm.nih.gov/datasets/docs/v1/quickstarts/command-line-tools/

1) Statistical machine learning methods for the analysis of medical
epigenomics data.

2) Sequence analysis algorithms for detection of RNA-DNA interactions.

Applicants should hold a M.Sc . or PhD in Computer Science or related
areas. Experience in the analysis of biological sequences, gene
expression and gene regulation is desirable. The candidate should have
solid programming skills (C, Python and/or R) and acquaintance with
Linux. Experience with high performance computing is a plus. The
working language of the group is English.

The position is based on the German TV-L 13 salary scale, including
all German social benefits (health insurance and pension scheme). The
expected starting date is September 2014. Interested candidates should
send a CV, statement of research interests and the names of three
references to jobs@costalab.org.

More at http://costalab.org/wp/phd-and-postdoc-position-in-computational-biology/

Address of the bookmark: http://www.weizmann.ac.il/Biological_Regulation/IgorUlitsky/PLAR

2. In some perl script I found this
 . . . . . .
 . . . . . .
# It works for me, only God understood how it is working
while (/(<\/[^>]+>)|(<[^>]+>)|(<[^>]+>)$|([^><]+)/go) {
            $startGene=$1;
            $beginChromosome=$2;
    
. . . . . .
 .. . . . . .
}

3. One more interesting message in Perl found …. It will must tickle you bone :)
open(my $fh, "<", "gene.txt")    or kill " Me if you think this is a mistake :$!";

4. From the Perl

  while () {  # "The Mothership Connection is here!"
    print “$_\n”; # Printing the offspring :)

5. Perl message
if ($1) { print “Just found a the error in chromosome !!!, yahoo…”; else { “That is not error, but mutation you moron!”;

6. One genome database curator walk in wine bar asked the bartender:
CREATE TABLE gene IF NOT EXISTS SexOnTheBeach;

Lab page: http://jarvislab.net/

Applications are invited on prescribed format along with self attested copies of the certificates for Faculty Positions in the following areas:
Sciences – Systems Biology, Computer Aided Drug Designing, Statistics, Applied Mathematics, Applied Physics. BioMedical Engineering – BioMechanics, BioMedical Instrumentation.

Last Date : May 10, 2014

Details are available on our website : http://www.iiita.ac.in

http://www.iiita.ac.in/downloads/announcements/uploads/FACULTY_Advertisement_NO-FS-01_2014130.pdf

Install

pip install genomenotebook

How to use

Create a simple genome browser with a search bar. The sequence appears when zooming in.

import genomenotebook as gn

g=gn.GenomeBrowser(genome_path, gff_path, init_pos=10000)
g.show()

Tracks can be added to visualize your favorite genomics data. See Examples for more !!!!

Address of the bookmark: https://dbikard.github.io/genomenotebook/