A differential transcriptional profile by Culex quinquefasciatus larvae resistant to Lysinibacillus sphaericus IAB59 highlights genes and pathways associated with the resistance phenotype

Background: The study of the mechanisms by which larvae of the Culex quinquefasciatus mosquito survive exposure to the entomopathogen Lysinibacillus sphaericus has benefited substantially from the generation of laboratory-selected colonies resistant to this bacterium. One such colony, RIAB59, was selected after regular long-term exposure of larvae to the L. sphaericus IAB59 strain. This strain is characterized by its ability to produce the well known Binary (Bin) toxin, and the recently characterized Cry48Aa/Cry49Aa toxin, able to kill Bin-resistant larvae. Resistance to Bin is associated with the depletion of its receptor, Cqm1 α-glucosidase, from the larvae midgut. This study aimed to identify novel molecules and pathways associated with survival of the RIAB59 larvae and the resistance phenotype.

Methods: A transcriptomic approach and bioinformatic tools were used to compare the profiles derived from the midguts of larvae resistant and susceptible to L. sphaericus IAB59.

Results: The RNA-seq profiles identified 1355 differentially expressed genes (DEGs), with 673 down- and 682 upregulated transcripts. One of the most downregulated DEGs was cqm1, which validates the approach. Other strongly downregulated mRNAs encode the enzyme pantetheinase, apolipoprotein D, lipases, heat-shock proteins and a number of lesser known and hypothetical polypeptides. Among the upregulated DEGs, the top most encodes a peroxisomal enzyme involved in lipid metabolism, while others encode enzymes associated with juvenile hormone synthesis, ion channels, DNA binding proteins and defense polypeptides. Further analyses confirmed a strong downregulation of several enzymes involved in lipid catabolism while the assignment of DEGs into metabolic pathways highlighted the upregulation of those related to DNA synthesis and maintenance, confirmed by their clustering into related protein networks. Several other pathways were also identified with mixed profiles of down- and upregulated transcripts. Quantitative RT-PCR confirmed the changes in levels seen for selected mRNAs.

Conclusions: Our transcriptome-wide dataset revealed that the RIAB59 colony, found to be substantially more resistant to Bin than to the Cry48Aa/Cry49Aa toxin, developed a differential expression profile as well as metabolic features co-selected during the long-term adaptation to IAB59 and that are most likely linked to Bin resistance.