design: rRNA depletion was prepared using three QIAGEN probes kits: Fast Select rRNA Plant, Worm, and 5S/16S/23S. NEBNext Ultra II RNA Library Preparation Kit from Illumina was used following manufacturer's recommendations. Sequencing libraries were multiplexed and clustered onto flow cells and then sequenced using 2x150bp paired end configuration. Image analysis and base calling were conducted with Illumina Control software. Raw sequences data was converted into .fastq files and de-multiplexed using Illuina bcl2fastq 2.17 software with one mis-match allowed for index sequence identification.