Harderian gland transcriptome
each gland was conducted for the total RNA extraction with Qiagens RNeasy Mini kit (Venlo, The Netherlands). The extraction followed the protocol for purification of total RNA from animal tissues provided by manufacturer. Then, a total RNA extracted from each species were used for library preparation on Macrogen, Inc. (Seoul, South Korea) using the TruSeq stranded total RNA library prep with Ribo-Zero Globin kit and finally sequenced in one lane on the Illumina HiSeq 2000 (Illumina, Inc., San Diego, CA, USA) with 100-base-pair(bp) paired-end reads (PE) (4-5 Gb).