Question: Find HBV integration sites in human genome.
0
gravatar for Chirag Nepal
4 months ago by
Chirag Nepal2.2k
Copenhagen
Chirag Nepal2.2k wrote:

Hi there,

I want to identify HBV (virus) integration sites in human genome.

I have single end CAGE-seq (~30 nucleotide long) data on HBV patients. I mapped reads to human genome (used Bowtie2 version 2.29), and unmapped reads were mapped to HBV genome.

From the remaining unmapped reads, I want to find reads that partially align to human and partially to HBV. If the data was pair-end, it would have been slightly easier. Can you please suggest how do I systematically (logic how to do it, I can implement it) get this information to find the integration site.

I am thinking of fragmenting (keep same FASTQID) the unmapped reads and map to both human and HBV, and identify which IDs map to both human and HBV.

Any suggestion on how i could do this efficiently. Does using BWA help in this case?

Thank you !!

ADD COMMENTlink modified 4 months ago by Pierre Lindenbaum129k • written 4 months ago by Chirag Nepal2.2k
0
gravatar for Pierre Lindenbaum
4 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum129k wrote:

map your reads with bwa mem against a reference containing both host+virus sequence and then detect the discordant reads: Extracting chimeric reads from mapping

ADD COMMENTlink written 4 months ago by Pierre Lindenbaum129k

Thank you Pierre ! So you suggest to merge genome sequence of human + virus. Then map the unmapped reads using BWA on combined genome assembly. Reads that map partially to human and virus will be flagged as discordant. I think it is a neat idea. What is the flag for discordant reads.

ADD REPLYlink written 4 months ago by Chirag Nepal2.2k

What is the flag for discordant reads.

there is no specific SAM flag, look at the link : Extracting chimeric reads from mapping

ADD REPLYlink written 4 months ago by Pierre Lindenbaum129k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 641 users visited in the last hour