Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • batch effect in radseq

    Hi all,

    We sent 600 samples of a plant species (no reference genome) for DNA-extraction, library prep and ddRadSeq (paired end, Illumina Hiseq). Firstly, 200 samples were sequenced as pilot, then the others 400 followed the same procedure. The raw reads were then:

    - Demultiplexed and adepter clipped
    - Filtered by restriction enzyme cut site
    - Merged and Clusterized using CD-HIT-EST to form the reference contigs
    - Quality trimmed
    - Aligned against the reference contigs (Bowtie2)
    - Variant discovery and SNP called using Freebayes.

    Then, many filters were applied to remove less informative SNPs or samples:
    - filter out snps with missingness (% of samples with missing data) above 10%
    - filter out samples with missingness (% of snps with missing data) above 10%
    - filter out snps with minor allele frequency (frequency of the less represented allele) below 5%

    This resulted in roughly 500 samples kept and 1500 snps.

    The results show a clear batch effect between the samples sequenced as pilot and those sequenced after. Particularly:

    - The PCA show a clear separation on the first axis (explaining 3% of the total variance) between pilot and rest of the samples.
    - The samples from the pilot show an overall higher level of heterozygosity (12 % in pilot against 9 % rest of samples). This increased observation of heterozygous loci is distributed across many loci.

    I've spent a few weeks trying to identify a technical factor that could explain these differences. The most remarkable technical difference I found between the samples is:

    - I calculated the median read count for each sample, then compared it between pilot and non-pilot samples. The mean of this parameter is higher in the pilot samples and the variance is half in comparison to the other samples. This is shown in the attached figure (AB: pilot samples, CDEF: other samples).

    My questions are:
    - Can a differential sequencing depth cause an increase/decrease of heterozygous call in radseq? Is there a read-processing step that could overcome this problem?
    - Can you see other possible causes of these batch effects?

    thank you in advance

    OS
    Attached Files

  • #2
    The higher read count and higher heterozygosity is certainly suggestive, since more reads will give you a better chance to see and call a second allele, or create an error artifact that looks like a second allele. You might re-do the genotype calls and cap the reads at a nucleotide to 20 and see if the difference is reduced.

    I'd also look at SNPs per read nucleotide to see if there was a difference in quality scores between runs. The paired-end read seems less consistent so check that closely. You might see an increase in SNP at the end of the read, and a greater increase in the pilot, or a spike at a particular position.
    Providing nextRAD genotyping and PacBio sequencing services. http://snpsaurus.com

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    10 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    9 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    67 views
    0 likes
    Last Post seqadmin  
    Working...
    X