cnv detection in cow
1
0
Entering edit mode
12 days ago
hamideh ▴ 20

Hi

I want to get cnv for cow. I used cnv-facet and cnvpytor for this. But if I understand correctly these are designed for humans and I didn't get any response. What do you suggest for me? Or with what changes can I use these for cow?

Thanks for your guidance

cow cnv • 347 views
ADD COMMENT
0
Entering edit mode
12 days ago
dthorbur ★ 2.6k

CNV detection tools like cnv-facet and cnvpytor (which is a python wrapper of CNVnator) have been successfully deployed on many different non-model species. I've used CNVnator on a diploid fish species without any problems. Whilst assumptions about data quality might differ in non-model species, you should be fine with cow which has chromosome level assemblies.

If you're not getting any results I suspect there is an issue with the quality of your data, a mistake in the data processing, or you haven't followed the instructions correctly. Impossible to say without any detail.

Some things that may be useful to know:

  1. What is your input data? (Platform, depth, etc...)
  2. How have you processed/cleaned the data?
  3. What reference genome are you using?
  4. What does the run log say?
  5. If everything looks good until here, then maybe we should look at the commands you used.
ADD COMMENT
0
Entering edit mode

Thank you very much for your answer. My data is wgs and my input file was a bam file. And I run the command and it stopped at length 1034!

cnvpytor -root file.pytor -rd 478m.sorted.bam -T genome.fa

and I received this response.

  • cnvpytor.root - INFO - Reading data for chromosome NKLS02001987.1 with length 1034 was killed
ADD REPLY
0
Entering edit mode

I've not gotten that error before, but I would be wary of a chromosome that is only 1kb. You're not going to get any meaningful results with really small reference contigs. I think you might need some filtering of your reference genome (as well as your input data if you've not already done that).

ADD REPLY
0
Entering edit mode

Can I see the steps to prepare your data for use in cnvpytor? After mapping with bwa, I convert the sam file to bam, sort and index it. Of course, this is the first time I've done cnv. But I've done other analyses. Thank you.

ADD REPLY

Login before adding your answer.

Traffic: 2724 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6