Obligate intracellular pathogens depend on their host for growth yet must also evade detection by host defenses. Here we investigate host adaptation in two Microsporidia, the specialist Edhazardia aedis and the generalist Vavraia culicis, pathogens of disease vector mosquitoes. Genomic analysis and deep RNA-Seq across infection time courses reveal fundamental differences between these pathogens. E. aedis retains enhanced cell surface modification and signalling capacity, upregulating protein trafficking and secretion dynamically during infection. V. culicis is less dependent on its host for basic metabolites and retains a subset of spliceosomal components, with a transcriptome broadly focused on growth and replication. Transcriptional profiling of mosquito immune responses reveals that response to infection by E. aedis differs dramatically depending on the mode of infection, and that antimicrobial defensins may play a general role in mosquito defense against Microsporidia. This analysis illuminates fundamentally different evolutionary paths and host interplay of specialist and generalist pathogens.

Recent outbreaks of locally transmitted dengue and Zika viruses in Florida have placed more emphasis on integrated vector management plans for Aedes aegypti (L.) and Aedes albopictus Skuse. Adulticiding, primarily with pyrethroids, is often employed for the immediate control of potentially arbovirus-infected mosquitoes during outbreak situations. While pyrethroid resistance is common in Ae. aegypti worldwide and testing is recommended by CDC and WHO, resistance to this class of products has not been widely examined or quantified in Florida. To address this information gap, we performed the first study to quantify both pyrethroid resistance and genetic markers of pyrethroid resistance in Ae. aegypti and Ae. albopictus strains in Florida. Using direct topical application to measure intrinsic toxicity, we examined 21 Ae. aegypti strains from 9 counties and found permethrin resistance (resistance ratio (RR) = 6-61-fold) in all strains when compared to the susceptible ORL1952 control strain. Permethrin resistance in five strains of Ae. albopictus was very low (RR<1.6) even when collected from the same containers producing resistant Ae. aegypti. Characterization of two sodium channel kdr alleles associated with pyrethroid-resistance showed widespread distribution in 62 strains of Ae. aegypti. The 1534 phenylalanine to cysteine (F1534C) single nucleotide polymorphism SNP was fixed or nearly fixed in all strains regardless of RR. We observed much more variation in the 1016 valine to isoleucine (V1016I) allele and observed that an increasing frequency of the homozygous V1016I allele correlates strongly with increased RR (Pearson corr = 0.905). In agreement with previous studies, we observed a very low frequency of three kdr genotypes, IIFF, VIFF, and IIFC. In this study, we provide a statewide examination of pyrethroid resistance, and demonstrate that permethrin resistance and the genetic markers for resistance are widely present in FL Ae. aegypti. Resistance testing should be included in an effective management program. Aedes aegypti (Yellow-fever mosquito) and Aedes albopictus (Asian Tiger mosquito) can vector a variety of arboviruses that cause diseases and are thus a public health concern. Pyrethroid insecticide resistance is common in Ae. aegypti in many locations worldwide and can adversely affect vector control operations. However, the resistance status of these vectors in Florida is largely unreported and recent local transmission of dengue and Zika viruses has made this information critical for effective control operations. In this study, we showed that permethrin resistance and two common SNPs of the voltage gated sodium channel (V1016I and F1534C) previously associated with pyrethroid resistance were widely present in Florida Ae. aegypti strains. We also observed a strong correlation between the dilocus knock down response (kdr) genotype and resistance ratio (RR) as determined by topical application, which suggests, as have others, that kdr frequency may be a useful indicator of resistance in Aedes aegypti.

Puerto Rico (PR) has a long history of vector-borne disease and insecticide-resistant Aedes aegypti (L.). Defining contributing mechanisms behind phenotypic resistance is critical for effective vector control intervention. However, previous studies from PR have each focused on only one mechanism of pyrethroid resistance. This study examines the contribution of P450-mediated enzymatic detoxification and sodium channel target site changes to the overall resistance phenotype of Ae. aegypti collected from San Juan, PR, in 2012. Screening of a panel of toxicants found broad resistance relative to the lab susceptible Orlando (ORL1952) strain. We identified significant resistance to representative Type I, Type II, and nonester pyrethroids, a sodium channel blocker, and a sodium channel blocking inhibitor, all of which interact with the sodium channel. Testing of fipronil, a chloride channel agonist, also showed low but significant levels of resistance. In contrast, the PR and ORL1952 strains were equally susceptible to chlorfenapyr, which has been suggested as an alternative public health insecticide. Molecular characterization of the strain indicated that two common sodium channel mutations were fixed in the population. Topical bioassay with piperonyl butoxide (PBO) indicated cytochrome P450-mediated detoxification accounts for approximately half of the resistance profile. Transcript expression screening of cytochrome P450s and glutathione-S-transferases identified the presence of overexpressed transcripts. This study of Puerto Rican Ae. aegypti with significant contributions from both genetic changes and enzymatic detoxification highlights the necessity of monitoring for resistance but also defining the multiple resistance mechanisms to inform effective mosquito control.

DNA metabarcoding offers new perspectives, especially with regard to the high‐throughput identification and diagnostics of pathogens. Microsporidia are an example of widely distributed, opportunistic and pathogenic microorganisms in which molecular identification is important for both environmental research and clinical diagnostics. We have developed a method for parallel detection of both microsporidian infection and the host species. We designed new primer sets: one specific for the classical Microsporidia (targeting the hypervariable V5 region of small subunit [ssu] rDNA), and a second one targeting a shortened fragment of the COI gene (standard metazoan DNA‐barcode); both markers are well suited for next generation sequencing. Analysis of the ssu rDNA data set representing 607 microsporidian species (120 genera) indicated that the V5 region enables identification of >98% species in the data set (596/607). To test the method, we used microsporidians that infect mosquitoes in natural populations. Using mini‐COI data, all field‐collected mosquitoes were unambiguously assigned to seven species; among them almost 60% of specimens were positive for at least 11 different microsporidian species, including a new microsporidian ssu rDNA sequence (Microsporidium sp. PL01). Phylogenetic analysis showed that this species belongs to one of the two main clades in the Terresporidia. We found a high rate of microsporidian co‐infections (9.4%). The numbers of sequence reads for the operational taxonomic units suggest that the occurrence of Nosema spp. in co‐infections could benefit them; however, this observation should be retested using a more intensive host sampling. Our results show that DNA barcoding is a rapid and cost‐effective method for deciphering sample diversity in greater resolution, including the hidden biodiversity that may be overlooked using classical methodology.