design: Paired-end 2  101 nt sequencing for juvenile samples and paired-end 2  150 nt sequencing for embryonic samples were conducted by the High Throughput Genomics Core Facility, Biodiversity Research Center, Academia Sinica, Taiwan, and Novogen, United States, respectively