A comparison of strategies for generating artificial replicates in RNA-seq experiments

Abstract

Due to the overall high costs, technical replicates are usually omitted in RNA-seq experiments, but several methods exist to generate them artificially. Bootstrapping reads from FASTQ-files has recently been used in the context of other NGS analyses and can be used to generate artificial technical replicates. Bootstrapping samples from the columns of the expression matrix has already been used for DNA microarray data and generates a new artificial replicate of the whole experiment. Mixing data of individual samples has been used for data augmentation in machine learning. The aim of this comparison is to evaluate which of these strategies are best suited to study the reproducibility of differential expression and gene-set enrichment analysis in an RNA-seq experiment. To study the approaches under controlled conditions, we performed a new RNA-seq experiment on gene expression changes upon virus infection compared to untreated control samples. In order to compare the approaches for artificial replicates, each of the samples was sequenced twice, i.e. as true technical replicates, and differential expression analysis and GO term enrichment analysis was conducted separately for the two resulting data sets. Although we observed a high correlation between the results from the two replicates, there are still many genes and GO terms that would be selected from one replicate but not from the other. Cluster analyses showed that artificial replicates generated by bootstrapping reads produce it p values and fold changes that are close to those obtained from the true data sets. Results generated from artificial replicates with the approaches of column bootstrap or mixing observations were less similar to the results from the true replicates. Furthermore, the overlap of results among replicates generated by column bootstrap or mixing observations was much stronger than among the true replicates. Artificial technical replicates generated by bootstrapping sequencing reads from FASTQ-files are better suited to study the reproducibility of results from differential expression and GO term enrichment analysis in RNA-seq experiments than column bootstrap or mixing observations. However, FASTQ-bootstrapping is computationally more expensive than the other two approaches. The FASTQ-bootstrapping may be applicable to other applications of high-throughput sequencing.

Transcriptome Profiling after Early Spinal Cord Injury in the Axolotl and Its Comparison with Rodent Animal Models through RNA-Seq Data Analysis

“…If neither was present, the axolotl gene IDs from the current reference genome (Amex60DD) were used as the gene names. To complete three technical replicates from the axolotl data, bootstrap resampling was applied to the matrix read counts in R [22].…”

A reliable and quick method for screening alternative splicing variants for low-abundance genes

“…For example, due to the overall costs of RNA-seq, technical replicates are usually omitted, especially in observational studies that already include a large number of biological replicates. For example, high reproducibility of RNA-seq analyses has been reported, and the use of technical replicates at the same study points could improve the statistical power of an experiment [ 18 ]. To date, RNA-Seq data obtained with short-read sequencers have been selected for transcriptome analysis due to their high fidelity, high coverage, and single-nucleotide resolution.…”