<p>The simplest thing to do is to trim the 150 bp fasts so that they are 100 long.  I don't think there is an easy way to correct for the fact that the 150 bp long reads will have a higher unambiguous rate of alignment and gene assignment than the 100 bp long reads.</p>
<p>If you have a mix of all the experimental conditions across all the lengths, you can include length as an element in the design, so that DESeq2 may be able to model the difference between the two batches.</p>

CC BY-SA 4.0
2021-03-08T21:53:11.747Z
2021-03-08T21:53:11.747Z
swbarnes2
<p>I've done this by doing a RPkM-like transformation which I call VSTPk. Transformations that adjust for transcript length will correct a little bit for the noise associated with mapping short reads to transcripts of varying lengths:</p>
<p><a href="https://bioinformatics.stackexchange.com/a/15218/73">https://bioinformatics.stackexchange.com/a/15218/73</a></p>
<p>Another, probably better, alternative is to use Salmon to map reads to transcripts, which applies its own mapping-based correction to read counts. Salmon produces a statistic that relates to the mappable length of a transcript, and normalises counts based on how many molecules the model predicts came from each transcript isoform.</p>
<p>Unfortunately all these normalisation techniques require good gene models, which are difficult to obtain for non-model organisms.</p>

CC BY-SA 4.0
2021-03-09T03:11:52.820Z
2021-03-09T03:11:52.820Z
gringer
<p>I have 20 samples out of which 14 are 100 bp in length and 6 are 150 bp. Is there a way to normalize the read length for cross-sample differential expression comparison? What would be the best way to approach this?</p>
<p>Note: The samples are NOT from reference/model species.</p>

CC BY-SA 4.0
2021-03-09T03:12:50.017Z
Are these paired-end reads (i.e. 2x100bp; 2x150bp), or single-end reads? That will make quite a big difference in terms of the mappability of the reads to transcripts.
gringer
<rna-seq><deseq2><normalization><edger><salmon>
Bulk RNA-Seq Read Length Normalization across different samples
# Bulk RNA-Seq Read Length Normalization across different samples

(Asked by: username_0 on 2021-03-08)

I have 20 samples out of which 14 are 100 bp in length and 6 are 150 bp. Is there a way to normalize the read length for cross-sample differential expression comparison? What would be the best way to approach this?

Note: The samples are NOT from reference/model species.

**Tags:** rna-seq, deseq2, normalization, edger, salmon
> **Comments:**
> **username_1**: Are these paired-end reads (i.e. 2x100bp; 2x150bp), or single-end reads? That will make quite a big difference in terms of the mappability of the reads to transcripts.  

---

## 2 answers

### Answer by username_2 on 2021-03-08. Score: 1

The simplest thing to do is to trim the 150 bp fasts so that they are 100 long. I don't think there is an easy way to correct for the fact that the 150 bp long reads will have a higher unambiguous rate of alignment and gene assignment than the 100 bp long reads.

If you have a mix of all the experimental conditions across all the lengths, you can include length as an element in the design, so that DESeq2 may be able to model the difference between the two batches.


---

### Answer by username_1 on 2021-03-09. Score: 0

I've done this by doing a RPkM-like transformation which I call VSTPk. Transformations that adjust for transcript length will correct a little bit for the noise associated with mapping short reads to transcripts of varying lengths:

<https://bioinformatics.stackexchange.com/a/15218/73>

Another, probably better, alternative is to use Salmon to map reads to transcripts, which applies its own mapping-based correction to read counts. Salmon produces a statistic that relates to the mappable length of a transcript, and normalises counts based on how many molecules the model predicts came from each transcript isoform.

Unfortunately all these normalisation techniques require good gene models, which are difficult to obtain for non-model organisms.


---

