Paired-end reads were filtered and trimmed using fastp (version 0.21.0 with default parameters) (1). Reads with quality inferior to 20 were removed from the data sets
Data sets were aligned against the Arabidopsis thaliana genome (TAIR10) using STAR (version 2.7.2a) allowing two mismatches (--outFilterMismatchNmax 2).
Using featureCounts (version 1.6.4), reads were count the paired-end fragments overlapping with TEs, reads mapping on several TE weighted.
We then computed TPM using the read counts obtained from featureCounts for TEs.
TEs were considered to be upregulated if:
mutant lines showed a 2-fold up-regulation as compared to Col in both biological replicates,
and had a value of TPM ≥ 5.
Heatmap was made with R using pheatmap
S. Chen, Y. Zhou, Y. Chen, J. Gu, fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884-i890 (2018).
K. Panda et al., Full-length autonomous transposable elements are preferentially targeted by expression-dependent forms of RNA-directed DNA methylation. Genome Biol 17, 170 (2016).
Readers should cite both the Bio-protocol preprint and the original research article where this protocol was used:
Davarinejad, H, Huang, Y, Couture, J and Jacob, Y(2026). Transposon analysis from RNA sequencing. Bio-protocol Preprint. bio-protocol.org/prep3008.
Davarinejad, H., Huang, Y., Mermaz, B., LeBlanc, C., Poulet, A., Thomson, G., Joly, V., Muñoz, M., Arvanitis-Vigneault, A., Valsakumar, D., Villarino, G., Ross, A., Rotstein, B. H., Alarcon, E. I., Brunzelle, J. S., Voigt, P., Dong, J., Couture, J. and Jacob, Y.(2022). The histone H3.1 variant regulates TONSOKU-mediated DNA repair during replication. Science 375(6586). DOI: 10.1126/science.abm5320
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this
article to respond.
Post a Question
0 Q&A
Spinning
This protocol preprint was submitted via the "Request
a Protocol" track.