Seqanswers Leaderboard Ad



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

  • Sequenza parameters

    Dear sequenza users,

    I am using Sequenza only recently. After the python script, I ran the 3 functions sequenza.extract(), and sequenza.results() with default parameters. My data is a set of 5 tumor-normal WES pairs.
    Now, I would like to adjust the parameters more precisely. The problem is that I do not clearly understand part of them, even after reading the vignette, the google group and the description of the R functions.

    Sorry for the long post, I do not know a lot sequenza and I am not an expert in segmentation.

    Here are the parameters:
    1) sequenza.extract
    • window=10^6
      "size of windows used when plotting mean and quartile ranges of depth ratios and B-allele frequencies. Smaller windows will take more time to compute"
      Is this parameter useful only for plotting?
      I changed it to 500 and I did not see changes when looking at the genome_view.pdf
    • overlap=1
      "integer specifying the number of overlapping windows"
      If we consider a specific window, it can overlap only 0, 1 or 2 window(s), right?
    • min.type.freq=0.9
      "minimum frequency of aberrant types"
      What does it mean?
    • weighted.mean=TRUE
      "boolean to select if the segments should be calculated using the read depth as weights to calculate depth ratio and B-allele frequency means"
      What does this mean?

    • N.ratio.filter=10
      "Threshold of minimum number of observation of depth ratio in a segment"
      Minimum number of variants in a segment?
    • N.BAF.filter=1
      "threshold of minimum number of observation of B-allele frequency in a segment"
      Minimum number of variants in a segment?
      Why is the default value not the same as for N.ratio.filter?
    • segment.filter=3 10^6
      "threshold segment length (in base pairs) to filter out short segments, that can cause noise when fitting the cellularity and ploidy parameters. The threshold will not affect the allele-specific segmentation"
      Is it the minimum length of a segment?
      What is the usual range? 3 10^6 seems big.
    • ratio.priority=FALSE
      "logical, if TRUE only the depth ratio will be used to determine the copy number state, while the Bf value will be used to determine the number of B-alleles"
      Does this mean that with FALSE, both depth ratio and BAF will be use to determine copy number and only BAF to determine the number of B-alleles?

    My last concern is the gamma parameter, which seems crucial in segmentation.
    Since I have WES data, I chose breaks.method="full" and now I want to determine gamma.pcf and kmin.pcf.
    From what I understood, in cancer data, the range varies often from 15 to 40.
    I started with 40. From 140 to 40, copy number estimations did not change but 2/5 cellularity estimations changed: from 0.46 to 0.97 and from 0.35 to 0.26.
    Do you have a way to determine gamma?
    I will try to use the gamma.plot function of copynumber package.

    I am surprised to get very distant cellularities from what we believe. We are pretty sure to have high purity, most likely >90%, but sequenza returned range between 0.35 to 0.5 with default settings.
    Does someone get close values in its data?

    Sorry for the high number of questions!
    Thank you in advance for your feedback.

  • #2
    It seems that there are not a lot of sequenza users here.
    I will post also the question on the sequenza user google group. If someone is interested, an answer could be there.


    Latest Articles


    • seqadmin
      The Impact of AI in Genomic Medicine
      by seqadmin

      Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
      02-26-2024, 02:07 PM
    • seqadmin
      Multiomics Techniques Advancing Disease Research
      by seqadmin

      New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

      A major leap in the field has
      02-08-2024, 06:33 AM





    Topics Statistics Last Post
    Started by seqadmin, 02-23-2024, 04:11 PM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 02-21-2024, 08:52 AM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 02-20-2024, 08:57 AM
    0 responses
    Last Post seqadmin  
    Started by seqadmin, 02-14-2024, 09:19 AM
    0 responses
    Last Post seqadmin