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
    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
  • 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

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-14-2024, 06:13 AM
0 responses
32 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-08-2024, 08:03 AM
0 responses
71 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-07-2024, 08:13 AM
0 responses
79 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-06-2024, 09:51 AM
0 responses
68 views
0 likes
Last Post seqadmin  
Working...
X