Total RNA quality was assessed by the HAIB using an Agilent 2100 Bioanalyzer and an Invitrogen Qubit to ensure quality for the libraries for each time point. A total of seven libraries were prepared, one for each time point of 2, 12, 24, 36, 48, 60, and 72 h post-eclosion. The Illumina RNA Sample Prep Kits for mRNA Seq was used to prepare the libraries, which were run using the paired-end 50 nt read module (HAIB), which also conducted the base calling, barcode parsing, and removal of low quality reads. Further cleaning of adapter was performed using Trimmomatic 10 (link). The Cx. quinquefasciatus genome from Vectorbase 11 was used to map the cleaned reads using Tophat 12 (link) and gene expression estimations were performed using Cufflinks, and differential gene expression was tested using Cuffdiff as time series data 13 . After analysis, only genes with expression values ≥1, as measured in number of fragments mapped for per kilo (every thousand) bases of gene length for every million fragments sequenced (FPKM), were retained for expression comparisons 14 . All data have been submitted to the Gene Expression Omnibus at NCBI as accession #GSE51327.