Hemocyte RNA-Seq analysis of Indian malarial vectors Anopheles stephensi and Anopheles culicifacies: From similarities to differences

In-depth transcriptomic analysis of Anopheles gambiae hemocytes uncovers novel genes and the oenocytoid developmental lineage

BACKGROUND

Hemocytes are immune cells that patrol the mosquito hemocoel and mediate critical cellular defense responses against pathogens. However, despite their importance, a comprehensive transcriptome of these cells was lacking because they constitute a very small fraction of the total cells in the insect, limiting the study of hemocyte differentiation and immune function.

RESULTS

In this study, an in-depth hemocyte transcriptome was built by extensive bulk RNA sequencing and assembly of hemocyte RNAs from adult A. gambiae female mosquitoes, based on approximately 2.4 billion short Illumina and about 9.4 million long PacBio high-quality reads that mapped to the A. gambiae PEST genome (P4.14 version). A total of 34,939 transcripts were annotated including 4,020 transcripts from novel genes and 20,008 novel isoforms that result from extensive differential splicing of transcripts from previously annotated genes. Most hemocyte transcripts identified (89.8%) are protein-coding while 10.2% are non-coding RNAs. The number of transcripts identified in the novel hemocyte transcriptome is twice the number in the current annotation of the A. gambiae genome (P4.14 version). Furthermore, we were able to refine the analysis of a previously published single-cell transcriptome (scRNAseq) data set by using the novel hemocyte transcriptome as a reference to re-define the hemocyte clusters and determine the path of hemocyte differentiation. Unsupervised pseudo-temporal ordering using the Tools for Single Cell Analysis software uncovered a novel putative prohemocyte precursor cell type that gives rise to prohemocytes. Pseudo-temporal ordering with the Monocle 3 software, which analyses changes in gene expression during dynamic biological processes, determined that oenocytoids derive from prohemocytes, a cell population that also gives rise to the granulocyte lineage.

CONCLUSION

A high number of mRNA splice variants are expressed in hemocytes, and they may account for the plasticity required to mount efficient responses to many different pathogens. This study highlights the importance of a comprehensive set of reference transcripts to perform robust single-cell transcriptomic data analysis of cells present in low abundance. The detailed annotation of the hemocyte transcriptome will uncover new facets of hemocyte development and function in adult dipterans and is a valuable community resource for future studies on mosquito cellular immunity.

Immune Responses and Transcriptomic Analysis of Nilaparvata lugens against Metarhizium anisopliae YTTR Mediated by Rice Ragged Stunt Virus

Plant viruses and entomopathogenic fungi (EPF) can both elicit immune responses in insects. This study was designed to clarify whether plant viruses could affect the efficacy of EPF and explore the immune responses of brown planthopper (BPH), Nilaparvata lugens, in response to different pathogen infections. In this study, a strain of Metarhizium anisopliae YTTR with high pathogenicity against BPH was selected and explored whether rice ragged stunt virus (RRSV) could affect its lethality against BPH. RNA-seq was used to detect the inner responses of BPH in response to RRSV and M. anisopliae YTTR infection. Results showed that M. anisopliae YTTR has strong lethality against BPH (RRSV-carrying and RRSV-free). RRSV invasion did not affect the susceptibility of BPH against M. anisopliae YTTR at all concentrations. At 1 × 108 spores/mL, M. anisopliae YTTR caused a cumulative mortality of 80% to BPH at 7 days post-treatment. The largest numbers of differentially expressed genes (DEGs) was obtained in BPH treated with the two pathogens than in other single pathogen treatment. In addition, KEGG enrichment analysis showed that the DEGs were mostly enriched in immune and physiological mechanisms-related pathways. Both RRSV and M. anisopliae YTTR could induce the expression changes of immune-related genes. However, most of the immune genes had varying expression patterns in different treatment. Our findings demonstrated that RRSV invasion did not have any significant effect on the pathogenicity of M. anisopliae YTTR, while the co-infection of M. anisopliae YTTR and RRSV induced more immune and physiological mechanisms -related genes' responses. In addition, the presence of RRSV could render the interplay between BPH and M. anisopliae YTTR more intricate. These findings laid a basis for further elucidating the immune response mechanisms of RRSV-mediated BPH to M. anisopliae infection.

Potential of ivermectin as an active ingredient of the attractive toxic sugar baits against the Indian malaria vectors Anopheles culicifacies and Anopheles stephensi

BACKGROUND

Attractive toxic sugar bait (ATSB) is a novel vector control tool that exploits the sugar feeding behavior of mosquitoes. The current study aims to evaluate the efficacy of ivermectin-based ATSB against insecticide susceptible and resistant strains of major Indian malaria vectors - Anopheles culicifacies and Anopheles stephensi. ATSB with different concentrations of ivermectin were tested against mosquito vectors under standard laboratory conditions.

RESULTS

Dose-response analysis of ivermectin-ATSB showed 7.8 and 19.8 ppm as 50% and 90% lethal concentration (LC₅₀ and LC₉₀ ) values for insecticide susceptible An. culicifacies. In the case of insecticide susceptible An. stephensi, the LC₉₀ value was 35 ppm which was significantly higher in comparison to the LC₉₀ for An. culicifacies. The LC₅₀ of insecticide-resistant An. culicifacies and An. stephensi were 10.6 and 15.9 ppm respectively whereas LC₉₀ values were 36.9 and 61.0. Ivermectin-ATSB resulted in 99 ± 0.8% mortality of An. culicifacies and 93 ± 3.8% mortality of An. stephensi at an ivermectin concentration of 25 ppm. In another set of experiments, the ATSB solution containing standardized dose of ivermectin was sprayed on Allysum plant and mortality of both Anopheline vectors was recorded. Here, we observed > 90% mortality in both An. stephensi and An. culicifacies.

CONCLUSION

Our study demonstrates the potential of ivermectin-based ATSB in killing Indian malaria vectors irrespective of the method of application. Further field trials with ivermectin containing ATSB may pave the way for its usage in the national vector control program. © 2022 Society of Chemical Industry.

Anopheles gambiae strain (Ag55) cultured cells originated from Anopheles coluzzii and are phagocytic with hemocyte-like gene expression

Anopheles gambiae and Anopheles coluzzii are closely related species that are predominant vectors of malaria in Africa. Recently, A. gambiae form M was renamed A. coluzzii and we now conclude on the basis of a diagnostic PCR-restriction fragment length polymorphism assay that Ag55 cells were derived from A. coluzzii. We established an Ag55 cell transcriptome, and KEGG pathway analysis showed that Ag55 cells are enriched in phagosome pathway transcripts. The Ag55 transcriptome has an abundance of specific transcripts characteristic of mosquito hemocytes. Functional E. coli bioparticle uptake experiments visualized by fluorescence microscopy and confocal microscopy and quantified by flow cytometry establish the phagocytic competence of Ag55 cells. Results from this investigation of Ag55 cell properties will guide researchers in the use and engineering of the Ag55 cell line to better enable investigations of Plasmodium, other microbes, and insecticidal toxins.

The immune deficiency and c-Jun N-terminal kinase pathways drive the functional integration of the immune and circulatory systems of mosquitoes

The immune and circulatory systems of animals are functionally integrated. In mammals, the spleen and lymph nodes filter and destroy microbes circulating in the blood and lymph, respectively. In insects, immune cells that surround the heart valves (ostia), called periostial haemocytes, destroy pathogens in the areas of the body that experience the swiftest haemolymph (blood) flow. An infection recruits additional periostial haemocytes, amplifying heart-associated immune responses. Although the structural mechanics of periostial haemocyte aggregation have been defined, the genetic factors that regulate this process remain less understood. Here, we conducted RNA sequencing in the African malaria mosquito, Anopheles gambiae, and discovered that an infection upregulates multiple components of the immune deficiency (IMD) and c-Jun N-terminal kinase (JNK) pathways in the heart with periostial haemocytes. This upregulation is greater in the heart with periostial haemocytes than in the circulating haemocytes or the entire abdomen. RNA interference-based knockdown then showed that the IMD and JNK pathways drive periostial haemocyte aggregation and alter phagocytosis and melanization on the heart, thereby demonstrating that these pathways regulate the functional integration between the immune and circulatory systems. Understanding how insects fight infection lays the foundation for novel strategies that could protect beneficial insects and harm detrimental ones.

Functional disruption of transferrin expression alters reproductive physiology in Anopheles culicifacies

BACKGROUND

Iron metabolism is crucial to maintain optimal physiological homeostasis of every organism and any alteration of the iron concentration (i.e. deficit or excess) can have adverse consequences. Transferrins are glycoproteins that play important role in iron transportation and have been widely characterized in vertebrates and insects, but poorly studied in blood-feeding mosquitoes.

RESULTS

We characterized a 2102 bp long transcript AcTrf1a with complete CDS of 1872bp, and 226bp UTR region, encoding putative transferrin homolog protein from mosquito An. culicifacies. A detailed in silico analysis predicts AcTrf1a encodes 624 amino acid (aa) long polypeptide that carries transferrin domain. AcTrf1a also showed a putative N-linked glycosylation site, a characteristic feature of most of the mammalian transferrins and certain non-blood feeding insects. Structure modelling prediction confirms the presence of an iron-binding site at the N-terminal lobe of the transferrin. Our spatial and temporal expression analysis under altered pathophysiological conditions showed that AcTrf1a is abundantly expressed in the fat-body, ovary, and its response is significantly altered (enhanced) after blood meal uptake, and exogenous bacterial challenge. Additionally, non-heme iron supplementation of FeCl3 at 1 mM concentration not only augmented the AcTrf1a transcript expression in fat-body but also enhanced the reproductive fecundity of gravid adult female mosquitoes. RNAi-mediated knockdown of AcTrf1a causes a significant reduction in fecundity, confirming the important role of transferrin in oocyte maturation.

CONCLUSION

All together our results advocate that detailed characterization of newly identified AcTrf1a transcript may help to select it as a unique target to impair the mosquito reproductive outcome.