Sequencing of ZIKV genomes directly from Ae. aegypti and Cx. quinquefasciatus mosquitoes collected during the 2015–16 epidemics in Recife

Zika virus (ZIKV) is a negative sense RNA virus from the Flaviviridae family, which was relatively unknown until the first human epidemic in Micronesia, in 2007. Since then, it spread to French Polynesia and the Americas. Recife, the capital of Pernambuco state and epicenter of the Zika epidemic in Brazil, experienced a large number of microcephaly cases and other congenital abnormalities associated to the ZIKV infection from, 2015 to 16. Evidences suggest that both Aedes aegypti and Culex quinquefasciatus mosquitoes from Recife are capable of replicating and transmitting the virus. Here, we conducted high throughput sequencing of ZIKV genomes directly from Ae. aegypti and Cx. quinquefasciatus mosquitoes collected during the ZIKV epidemics in Recife, in order to investigate the variability and evolution of the virus. We obtained 11 draft ZIKV genomes derived from 5 pools from each Ae. aegypti and Cx. quinquefasciatus species. Genome coverage breadth ranged from 16 to 100% and average depth from 45 to 46,584×. Two of these genomes were obtained from pools of Cx. quinquefasciatus females with no sign of blood in the abdomen. Amino acid substitutions found here were not species-specific. In addition, molecular clock dating estimated that ZIKV draft genomes obtained here were co-circulating in other regions of the country during the epidemics. Overall results highlight that viral mutations and even minor variants can be detected in genomes directly sequenced from mosquito samples and insights about natural viral genomic variability and viral evolution can be useful when designing tools for mosquito control programs. •Zika virus was detected in Aedes aegypti and Culex quinquefasciatus field collected samples;•Eleven complete Zika virus genomes were obtained from these samples;•The detection of the Zika virus in male samples confirms vertical transmission in Ae. aegypti and Cu. quinquefasciatus species.•These genomes belongs to two different Zika lineages that simultaneously cocirculated at the time of the epidemics;•Amino acid changes found were likely not involved in virus adaptation to different vectors species.