Control group between 2 batches have different expression in RNAseq data
2
0
Entering edit mode
@4ba43fb6
Last seen 5 hours ago
Vietnam

Hello, I am working with RNAseq data which was acquired from experiments as following:

  • PBMCs were extracted from blood of healthy donors, then extracting monocytes, pooling monocytes of different healthy donors, then differentiate pooled monocytes into macrophages.
  • After that, Mycobacterium tuberculosis isolates were infected into macrophages, I also included lab strain H37Rv, and no-infection group.
  • Then macrophages were collected RNA at 4h and 24h post-infection.
  • Because the heavy workload, I divided into 2 batches, each batch included 2 samples of no-infection group, 2 samples of H37Rv as control.
  • Macrophages from 2 batches still from the same heathy donors.
  • After collecting all RNA samples, RNA was sequenced at the same time.

When I did PCA, there was big difference between batches, even for no-infection and H37Rv group as picture below

PCA plot for 2 batches

I would like to ask that there are any method for batch effect removal or normalization to combine the data between batches.

I was thinking to normalize Mtb strain to no-infection group in each batch before comparing between isolates.

Thank you in advance.

BatchQC BatchEffect rnaseqGene Normalization • 246 views
ADD COMMENT
2
Entering edit mode
@james-w-macdonald-5106
Last seen 15 hours ago
United States

Since batch is orthogonal to the treatment, you just add a batch factor to your model.

0
Entering edit mode
@61775469
Last seen 5 days ago
United States

ComBat is appropriate for this type of design and BatchQC can be used to evaluate how well ComBat corrects for the batch effect. BatchQC provides a shiny interface where you can normalize, batch correct and view multiple visualization to determine the effectiveness of batch effect correction. Your batch-corrected (and/or normalized) data can then be exported for downstream use.

ADD COMMENT

Login before adding your answer.

Traffic: 721 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6