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

Please visit our site at www.arraygen.com

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

JRF (Maximum Tenure‐ Five Years: 2yrs as JRF and 3yrs as SRF)
A first class Master’s Degree in Bioinformatics (likely 2 posts)
Around Rs 16,000/ Plus 30% HRA (as per rules of funding agency)

Applications are invited from candidates possessing the above qualifications. The upper age limit is as on the last date for receipt of application. (5 years relaxation to SC/ST candidates, 3 years to OBC candidates, and other entitled categories as per Govt rules). Actual No. of vacancies may vary.

Application form can be download from the website www.drdo.gov.in and E Mailed to inmashrd@gmail.com.
Last date to apply by email is 1700 hrs on 15 Oct 2014
Incomplete applications are liable to be rejected.
Confirmation will be sent to short-listed candidates through email only
Antecedents of selected candidates will be verified.
Written Test will be conducted from 0930-1030 hrs. Latecomers will not be considered.
Candidates will be required to produce certificates/testimonials in original at the time of interview.
It may please be noted that offer of Fellowship does not confer on fellows any right for absorption in DRDO.
Candidates should carry photocopy of Application form sent by email with them.
No TA/DA will be paid for attending interview & on joining.
Last date to apply by email is 1700 hrs on 15 Oct 2014

More at http://drdo.gov.in/drdo/English/jrf29092014.pdf
http://drdo.gov.in/drdo/English/index.jsp?pg=inmas29092014.jsp

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

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

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

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

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

Current projects in the Morota lab include developing kernel-based whole-genome prediction and kernel-based genome-wide association models, polygenic modeling of binary traits, reexamining the results from quantitative genetics analysis in light of functional annotation, and extending kernel methods (such as GBLUP and RKHS) specifically tailored for diverse types of emerging omics data.

In addition, candidates will be expected to leverage opportunities to interact with faculty in animal genetics and biometrics at the UNL in the areas of bioinformatics, breeding, functional genomics, quantitative genetics, and molecular genetics.

Candidates should have a B.S. or M.S. degree in quantitative disciplines with strong background and interest in statistical computing.
The starting date is Fall 2015.
For more information about research in the Morota lab at the UNL, visit: http://www.morotalab.org

A letter of interest in the position, C.V., and contact information for
three references should be emailed to Gota Morota at .
Review of applications will begin immediately, and continue until the
positions are filled. Informal inquiries are also welcome.

Also, please see: http://animalscience.unl.edu/anscprospectivegraduatestudents