  Batch effects, artificial sources of variation due to experimental design,
are a widespread phenomenon in high throughput data. Therefore, mechanisms for
detection of batch effects are needed requiring comparison of multiple samples.
We apply hierarchical clustering (HC) on DNA k-mer counts of multiple RNA-seq
derived Fastq files. Ideally, HC generated trees reflect experimental treatment
groups and thus may indicate experimental effects, but clustering of
preparation groups indicates the presence of batch effects. In order to provide
a simple applicable tool we implemented sequential analysis of Fastq reads with
low memory usage in an R package (seqTools) available on Bioconductor. DNA
k-mer counts were analysed on 61 Fastq files containing RNA-seq data from two
cell types (dermal fibroblasts and Jurkat cells) sequenced on 8 different
Illumina Flowcells. Results: Pairwise comparison of all Flowcells with
hierarchical clustering revealed strong Flowcell based tree separation in 6 (21
%) and detectable Flowcell based clustering in 17 (60.7 %) of 28 Flowcell
comparisons. In our samples, batch effects were also present in reads mapped to
the human genome. Filtering reads for high quality (Phred >30) did not remove
the batch effects. Conclusions: Hierarchical clustering of DNA k-mer counts
provides a quality criterion and an unspecific diagnostic tool for RNA-seq
experiments.
