Haplotype-resolved or phased sequence assembly provides a complete picture of genomes and complex genetic variations. However, current phased assembly algorithms either fail to generate chromosome-scale phasing or require pedigree information, which limits their application. We present a method that leverages long accurate reads and long-range conformation data for single individuals to generate chromosome-scale phased assembly within a day. Applied to three public human genomes, PGP1, HG002, and NA12878, our method produced haplotype-resolved assemblies with contig NG50 up to 25 Mb and phased ∼99.5% of heterozygous sites to 98–99% accuracy, outperforming trio-based approach in terms of both contiguity and phasing completeness. We demonstrate the importance of chromosome-scale phased assemblies to discover structural variants, including thousands of new transposon insertions, and of highly polymorphic and medically important regions such as HLA and KIR. Our improved method will enable high-quality precision medicine and facilitate new studies of individual haplotype variation and population diversity.
See our preprint here: https://doi.org/10.1101/810341.
DipAsm requires docker as DeepVariant uses it. Users need to make sure docker is installed and the docker service is started.
mkdir -p dipasm
cd dipasm
git clone https://github.com/shilpagarg/DipAsm.git
cd DipAsm/docker
docker build -t dipasm .
cd ../../..
docker run -it --rm -v $PWD/dipasm/DipAsm:/wd/dipasm/DipAsm/ -e HOSTWD=$PWD/dipasm/DipAsm -v /var/run/docker.sock:/var/run/docker.sock dipasm:latest /bin/bashThe docker run -it will start an interactive docker container session. You will be in
the virtual container envrionment which have the preinstall DipAsm and testing data.
You can run the test for DipAsm within the docker container environment by:
cd /wd/dipasm/DipAsm
bash test.sh | bash
ls test_output/out/assemble/test-H?.fasta # final assemblyOnce you enter the virtual Docker environment with the docker command line
shown above, go into the /wd/dipasm/DipAsm directory and run python pipeline.py. You will see:
Usage: pipeline.py [-h] --hic-path PATH --pb-path PATH --sample NAME [--female]
--prefix STR
Optional arguments:
-h, --help show this help message and exit
--hic-path PATH Use Hi-C data from this path. Should be named by *1.fastq
and *2.fastq.
--pb-path PATH Use PacBioCCS data from this path. All fastq will be
used.
--sample NAME Sample name to put for Read Group of BAM and Sample of
VCF.
--prefix STR Prefix name for the experiment, for example "refBased",
"ragooBased".
Example:
python pipeline.py --hic-path data/hic --pb-path data/pacbiocss --sample PGP1 --prefix asm
This pipeline produces phased assemblies as sample_output/prefix/assemble/sample-H?.fasta.
DipAsm depends on Peregrine, 3d-dna, minimap2, DeepVariant, whatshap and hapcut2.
React
Typescript
Django
Laravel
Facebook
Microsoft
Google
Alibaba
D3
Tencent