Comparative analysis of RNA enrichment methods for preparation of Cryptococcus neoformans RNA sequencing libraries

Ribosomal RNA (rRNA) is the major RNA constituent of cells, therefore most RNA sequencing (RNA-Seq) experiments involve removal of rRNA. This process, called RNA enrichment, is done primarily to reduce cost: without rRNA removal, deeper sequencing would need to be performed to balance the sequencing reads wasted on rRNA. The ideal RNA enrichment method would remove all rRNA without affecting other RNA in the sample. We have tested the performance of three RNA enrichment methods on RNA isolated from Cryptococcus neoformans, a fungal pathogen of humans. We show that the RNase H depletion method unambiguously outperforms the commonly used Poly(A) isolation method: the RNase H method more efficiently depletes rRNA while more accurately recapitulating the expression levels of other RNA observed in an unenriched gold standard. The RNase H depletion method is also superior to the Ribo-Zero depletion method as measured by rRNA depletion efficiency and recapitulation of protein-coding gene expression levels, while the Ribo-Zero depletion method performs moderately better in preserving non-coding RNA (ncRNA). Finally, we have leveraged this dataset to identify novel long non-coding RNA (lncRNA) genes and to accurately map the C. neoformans mitochondrial rRNA genes.