Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Should I sequence un-replicated blocks in RNAseq?

    Hi all,

    I am designing a new experiment to sequence transcriptomes from wasps. I have 3 treatment "environments" and 4 replicate populations in each (kept separate throughout), for a total of 12 separate rearing cages.

    I want to compare expression among treatment "environments", but of course the replicate populations should be taken into account. They are separate pools of individuals that have evolved separately. The problem is that my experimental wasps are not very cooperative, and I simply don't have 3 individuals from some of my replicate populations. So, I am trying to decide how to go about sequencing them.

    Can I assume similar variance across populations, and so get a good estimate of expression from non-replicated populations? Is this a terrible idea because they are coming from different genetic backgrounds?

    For some of my replicate populations, I have plenty of individuals. For others, just one. So, I am wondering how useful it will be to sequence those single individuals. As an alternative, I can use my limited resources differently, and just sequence more individuals from fewer replicate populations. This will tell me less about my experiment in general, but more about that population.....

    According to this thread, I can analyze using edgeR if I just have good replication in some groups: http://seqanswers.com/forums/showthread.php?t=46068

    Does anyone have experience with this? In EdgeR or other programs?

    For comparison/context, in previous sequencing from these populations, the wasps were more cooperative and I used the following design to analyze differential expression in DESeq2 (35 libraries in total):


    sel.line <- factor(c(rep("GroupB",11),rep("GroupD",12), rep ("GroupE", 12)))

    block <- factor(c(rep("1", 3), rep("2", 3), rep("3", 3), rep("4", 2),
    rep("1", 3), rep("2", 3), rep("3", 3), rep("4", 3),
    rep("1", 3), rep("2", 3), rep("3", 3), rep("4", 3)))

    ExpDesign3 <- data.frame(
    row.names = colnames(countdata3),
    sel.line = sel.line,
    block = block)

    CDS3 <- DESeqDataSetFromMatrix(countData = countdata3, colData=ExpDesign3, design= ~ block + sel.line)
    "

    In this case, block refers to replicate cages within treatment, and sel.line is the level of "interest" in my study.

    Many thanks,
    Alice

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, Today, 06:37 PM
0 responses
6 views
0 likes
Last Post seqadmin  
Started by seqadmin, Today, 06:07 PM
0 responses
6 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
66 views
0 likes
Last Post seqadmin  
Working...
X