Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • aggp11
    replied
    Originally posted by mimi_lupton View Post
    Dear all,
    I have a couple of questions regarding CNV-seq.
    I am having a play around with some DNA sequencing that is from custom capture. I have analysed some BAM files using the program. Some work fine, but some file combinations come up with an error, which is strange as all the data has been produced and analysed in exactly the same way. The error is;

    "Can't use an undefined value as an ARRAY reference at /share/apps/cnv-seq_1.0/bin/cnv-seq.pl line 204, <REF> line 6460211."

    I think this is some kind of Perl script error??

    Also does anyone have any experience of using capture data with CNV-seq? any advice would be greatly appreciated.

    Thanks for you help
    Hello mimi,

    I am trying to work on custom capture data usin cnv-seq. Could you enlighten me on how it worked for you and any suggestions on how to use it would be greatly appreciated.

    Thanks,
    Praful

    Leave a comment:


  • gprakhar
    replied
    Originally posted by cdinz View Post
    I have the same question too. It would be really great if someone can shed some light on it.
    The input type will not vary, it will still be a BAM file aligned to the ref genome.
    The difference comes in the analysis part.
    You will have to take into consideration the Intron-Exon Junctions. As a Intron will be taken as a CNV. Also another problem is that if a CNV breakpiomt lies in the Intron, you can not come to know about it.
    We have tumor Exome data, done by the True-Seq protocol, hence while doing the analysis use the co-ordinates provided by them.

    Leave a comment:


  • DineshCyanam
    replied
    Originally posted by louis7781x View Post
    Thanks for your help,I want to know the input type because the coverage of DNA-seq and RNA-seq is different.So if I use RNA-seq aligned to ref genome and use output file bam to run CNV-seq,Does it apply to analyze copy number varivation?
    I have the same question too. It would be really great if someone can shed some light on it.

    Leave a comment:


  • gprakhar
    replied
    Dear xiechao,

    I am trying to use CNV-seq for 36 length Read size simulated data. I am not detecting any CNV.
    While for the same data with 76bp Read Length, CNV-seq detects CNV.

    Is the read length hard coded in the script?

    Is there a way to use CNV-seq on 36 length read data.

    Thanking you,
    pg

    Leave a comment:


  • mimi_lupton
    replied
    Dear all,
    I have a couple of questions regarding CNV-seq.
    I am having a play around with some DNA sequencing that is from custom capture. I have analysed some BAM files using the program. Some work fine, but some file combinations come up with an error, which is strange as all the data has been produced and analysed in exactly the same way. The error is;

    "Can't use an undefined value as an ARRAY reference at /share/apps/cnv-seq_1.0/bin/cnv-seq.pl line 204, <REF> line 6460211."

    I think this is some kind of Perl script error??

    Also does anyone have any experience of using capture data with CNV-seq? any advice would be greatly appreciated.

    Thanks for you help

    Leave a comment:


  • jflores
    replied
    Hi Tonio,

    You can plot specific coordinates and try to infer from that if your gene of interest has or not gain/loss of copies.

    Script :
    samtools view Patient_test/Parsed_X/chr17.fa/bam/sorted.bam | perl -lane 'print "$F[2]\t$F[3]"' > test.hits
    samtools view Patient_ref/Parsed_X/chr17.fa/bam/sorted.bam | perl -lane 'print "$F[2]\t$F[3]"' > ref.hits

    perl cnv-seq.pl --test test.hits --ref ref.hits --genome chrom17

    data <- read.delim("test.hits-vs-ref.hits.log2-0.6.pvalue-0.001.minw-4.cnv"))
    plot.cnv.chr(data, chromosome=NA, from=NA, to=NA)

    ... For more options about plotting using CNV-seq you can have a look at the cnv.R file

    Regards,

    Rodrigo

    Leave a comment:


  • tonio100680
    replied
    Hi,

    how to view the level of a gene in CNV-seq ?

    Script :
    samtools view Patient_test/Parsed_X/chr17.fa/bam/sorted.bam | perl -lane 'print "$F[2]\t$F[3]"' > test.hits
    samtools view Patient_ref/Parsed_X/chr17.fa/bam/sorted.bam | perl -lane 'print "$F[2]\t$F[3]"' > ref.hits

    perl cnv-seq.pl --test test.hits --ref ref.hits --genome chrom17

    data <- read.delim("test.hits-vs-ref.hits.log2-0.6.pvalue-0.001.minw-4.cnv")
    cnv.print(data)
    cnv.summary(data)
    plot.cnv(data)
    ggsave("sample.pdf")
    [IMG]/home/labo/Sample1.pdf[/IMG]

    Regards

    Leave a comment:


  • louis7781x
    replied
    Originally posted by jflores View Post
    Hi Louis,

    CNV-seq input consists in just two files (reference & test) with just two columns each.
    First column corresponds to the third column in a BAM file, that is the reference sequence name of the alignment (Chr1, Chr2, ..), the second column is the fourth column of a BAM file, that is the corresponding 1-based leftmost mapping position of that read. Input files can look like this:

    1 999
    1 1234
    1 23456
    1 25234

    Full explanation of the input, how to get it, how to run CNV-seq, etc:




    Jose Flores
    Thanks for your help,I want to know the input type because the coverage of DNA-seq and RNA-seq is different.So if I use RNA-seq aligned to ref genome and use output file bam to run CNV-seq,Does it apply to analyze copy number varivation?

    Leave a comment:


  • jflores
    replied
    Originally posted by roman.sergio View Post
    Hello

    Did you already found how to detect the CNV with your data from non control case study?

    I am continuing seraching for a solution.

    Regards,

    Sergio
    Hi Sergio,

    With CNV-seq in mind and for non-control-case studies you can basically choose whatever reference you want. There might be some issues to worry about like what population the reference you choose come from? or what's the coverage .. It will depend on your purpose and the features of your sample to analyze.
    I’ve seen some studies taking the health individual NA10851 from the 1000 genomes data as a control/reference.

    ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/

    Hope can help,
    Jose Flores

    Leave a comment:


  • jflores
    replied
    Originally posted by louis7781x View Post
    Hi,I want to know:CNV-seq input file is DNA-seq or RNA-seq?
    Hi Louis,

    CNV-seq input consists in just two files (reference & test) with just two columns each.
    First column corresponds to the third column in a BAM file, that is the reference sequence name of the alignment (Chr1, Chr2, ..), the second column is the fourth column of a BAM file, that is the corresponding 1-based leftmost mapping position of that read. Input files can look like this:

    1 999
    1 1234
    1 23456
    1 25234

    Full explanation of the input, how to get it, how to run CNV-seq, etc:




    Jose Flores

    Leave a comment:


  • roman.sergio
    replied
    Hello

    Did you already found how to detect the CNV with your data from non control case study?

    I am continuing seraching for a solution.

    Regards,

    Sergio

    Leave a comment:


  • roman.sergio
    replied
    Hi !

    My data are DNA-seq. But there are not from one Case-control study.

    I working on the human Chr 8. I use a reference sequence to simulate my DNA-seq. I just would like to know if some parameters that I modify on my simulation influence the CNV detection on my DNA-seq.

    Some ideas?

    Sergio

    Leave a comment:


  • louis7781x
    replied
    CNV-seq input file is DNA-seq or RNA-seq?

    Originally posted by xiechao View Post
    CNV-seq, a new method to detect copy number variation using high-throughput sequencing.

    more details at BMC Bioinformatics :
    Background DNA copy number variation (CNV) has been recognized as an important source of genetic variation. Array comparative genomic hybridization (aCGH) is commonly used for CNV detection, but the microarray platform has a number of inherent limitations. Results Here, we describe a method to detect copy number variation using shotgun sequencing, CNV-seq. The method is based on a robust statistical model that describes the complete analysis procedure and allows the computation of essential confidence values for detection of CNV. Our results show that the number of reads, not the length of the reads is the key factor determining the resolution of detection. This favors the next-generation sequencing methods that rapidly produce large amount of short reads. Conclusion Simulation of various sequencing methods with coverage between 0.1× to 8× show overall specificity between 91.7 – 99.9%, and sensitivity between 72.2 – 96.5%. We also show the results for assessment of CNV between two individual human genomes.


    The package is at
    http://tiger.dbs.nus.edu.sg/cnv-seq
    Hi,I want to know:CNV-seq input file is DNA-seq or RNA-seq?

    Leave a comment:


  • roman.sergio
    replied
    Hi,

    I am trying to undertand how to build the file for CNV detection with sequence.
    I would like to know wich information is need it on the files.

    Do you have some examples?

    Leave a comment:


  • adamreid
    replied
    Detecting alignment collapse

    Hi,

    I want to check my genome sequence for collapsed, tandemly repeated genes. It seems to me that I could use software for detecting CNVs in that I am looking for significantly increased read depth. CNV-seq requires test and reference sequences, which makes sense to control for sequence-based noise in read depth however I only have my reference genome and one instance of Illumina data. Does anyone have any ideas?

    Adam

    Leave a comment:

Latest Articles

Collapse

  • seqadmin
    Best Practices for Single-Cell Sequencing Analysis
    by seqadmin



    While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
    06-06-2024, 07:15 AM
  • seqadmin
    Latest Developments in Precision Medicine
    by seqadmin



    Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

    Somatic Genomics
    “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
    05-24-2024, 01:16 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 07:49 AM
0 responses
12 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-20-2024, 07:23 AM
0 responses
14 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-17-2024, 06:54 AM
0 responses
16 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-14-2024, 07:24 AM
0 responses
24 views
0 likes
Last Post seqadmin  
Working...
X