Next Article in Journal
The Role of Insect Symbiotic Bacteria in Metabolizing Phytochemicals and Agrochemicals
Previous Article in Journal
Larval Crowding Did Not Enhance Adult Migration Propensity in Spodoptera frugiperda
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 13
Issue 7
10.3390/insects13070582
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Silencing Transglutaminase Genes TGase2 and TGase3 Has Infection-Dependent Effects on the Heart Rate of the Mosquito Anopheles gambiae

by
Abinaya Ramakrishnan
and
Julián F. Hillyer
*
Department of Biological Sciences, Vanderbilt University, Nashville, TN 37235, USA
*
Author to whom correspondence should be addressed.
Insects 2022, 13(7), 582; https://doi.org/10.3390/insects13070582
Submission received: 4 June 2022 / Revised: 18 June 2022 / Accepted: 22 June 2022 / Published: 26 June 2022
(This article belongs to the Section Insect Physiology, Reproduction and Development)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The immune and circulatory systems of insects are functionally integrated. An infection induces the migration of immune cells called hemocytes to the surface of the heart, where they kill pathogens around valves called ostia. In mosquitoes, a transglutaminase inhibits the infection-induced aggregation of hemocytes on the heart. Here, we studied whether transglutaminases also modify the heart contraction rate. First, we confirmed that an infection decreases the mosquito heart rate. Then, we found that disrupting transglutaminase genes has infection-dependent effects on the heart rate. Silencing TGase1 does not affect heart physiology. However, silencing TGase2 eliminates the infection-induced decrease in the heart rate, and silencing TGase3 decreases the heart rate in uninfected mosquitoes but increases the heart rate in infected mosquitoes. These experiments identify new factors that affect heart physiology in mosquitoes.

Abstract

Transglutaminases are pleiotropic enzymes that in mosquitoes participate in the formation of the mating plug and the wound-induced antimalarial response. Moreover, one transglutaminase, TGase3, negatively regulates the infection-induced aggregation of hemocytes on the heart. Given that TGase3 is an inhibitor of periostial hemocyte aggregation, we used RNAi-based gene silencing followed by intravital video imaging to scrutinize whether any of the three transglutaminases encoded in the genome of the mosquito, Anopheles gambiae, play a role in modulating the heart rate of uninfected and infected mosquitoes. Initially, we confirmed that an infection decreases the heart rate. Then, we uncovered that silencing TGase1 does not impact heart physiology, but silencing TGase2 results in a constant heart rate regardless of infection status, eliminating the infection-induced decrease in the heart rate. Finally, silencing TGase3 decreases the heart rate in uninfected mosquitoes but increases the heart rate in infected mosquitoes. We conclude that TGase2 and TGase3 modulate heart physiology and demonstrate that factors not classically associated with insect circulatory physiology are involved in the functional integration of the immune and circulatory systems of mosquitoes.

1. Introduction

Transglutaminases are calcium-dependent enzymes that crosslink proteins by catalyzing the formation of isopeptide bonds between glutamine and lysine residues [1,2]. In insects, transglutaminases are involved in cuticular morphogenesis, hemolymph coagulation, pathogen entrapment, peritrophic matrix formation, fertility, and immune system modulation [3,4,5,6,7,8]. Additionally, in the fruit fly, Drosophila melanogaster, transglutaminase is a negative regulator of the immune deficiency pathway (IMD). Specifically, transglutaminase crosslinks the N-terminal region of the IMD pathway’s transcription factor, Relish, forming a polymer that prevents its translocation into the nucleus, thereby allowing for increased bacterial proliferation following a septic infection [9]. Transglutaminase further inhibits the IMD pathway by incorporating natural polyamines into the DNA binding site of Relish [10]. However, transglutaminase does not always inhibit immunological function. Drosophila transglutaminase targets microbial surfaces, and larvae with reduced transglutaminase levels experience decreased survival [11,12]. Moreover, transglutaminase is a positive driver of immune processes in the termite, Reticulitermes chinensis, and the silk moth, Bombyx mori [13,14].
In mosquitoes, transglutaminases play roles in reproduction and immunity. The genome of the African malaria mosquito, Anopheles gambiae, encodes three transglutaminases [15]. The function of TGase1 is unknown, but TGase2 positively drives the wound-induced antimalarial response and TGase3 functions in the formation of the mating plug as well as the infection-induced aggregation of hemocytes on the heart [15,16,17,18,19]. Therefore, from an immunological perspective, TGase3 interfaces with the circulatory system.
The immune and circulatory system of mosquitoes—and insects in general—are functionally integrated [20,21,22]. An infection induces the migration of hemocytes to the surface of the heart, where they become sessile and phagocytose pathogens around heart valves called ostia. These heart-associated hemocytes, called periostial hemocytes, are advantageous because they attack circulating pathogens in locations of high hemolymph flow [23]. Although the structural mechanics of periostial hemocyte aggregation have been defined, the molecular drivers of this process are less understood. Members of the thioester containing protein family (TEP) and nimrod protein family (NIM) influence periostial hemocyte aggregation [24,25]. Moreover, the IMD pathway and the interlinked JNK pathway—as tested via the RNAi-based silencing of the positive regulators rel2, JNK1 and JNK3, and the negative regulators caspar and puckered—drive periostial hemocyte aggregation [26]. An RNAseq analysis in that same study identified transglutaminases as being more highly expressed in the heart containing periostial hemocytes than in other tissues [26], and follow-up experiments confirmed that TGase3—but not TGase1 or TGase2—is a negative regulator of periostial hemocyte aggregation during the early stages of infection, presumably by inhibiting the IMD pathway [19]. However, whether transglutaminases affect how the heart contracts remained unknown. This question is of particular interest because a systemic infection reduces the mosquito heart rate, and this is modulated by nitric oxide produced by periostial hemocytes [27].
The main circulatory organ in mosquitoes is the dorsal vessel, which is a muscular tube that traverses the length of the body and is subdivided into an aorta in the thorax and a heart in the abdomen [28,29,30,31]. Hemolymph enters the dorsal vessel through ostia that are present in every abdominal segment and is then propelled by wave-like contractions of the heart muscle [28,29]. The mosquito heart is myogenic, but its rhythmicity can be altered by both endogenous and exogenous factors [27,32,33,34,35,36]. In this study, we scrutinized whether RNAi-based silencing of transglutaminases affects the mosquito heart rate in both uninfected and infected mosquitoes. We uncovered that TGase3 silencing and, to a lesser extent TGase2 silencing, has infection-dependent effects on the mosquito heart rate.

2. Materials and Methods

2.1. Mosquito Rearing, Maintenance and Manipulation

Experiments were conducted on female Anopheles gambiae, Giles sensu stricto (G3 strain; Diptera: Culicidae). Mosquitoes were reared and maintained in an environmental chamber at 27 °C and 75% relative humidity, with a 12 h:12 h light–dark photoperiod [36]. Upon emergence, adults were fed ad libitum on cotton pads saturated with 10% sucrose. Prior to a manipulation, mosquitoes were anesthetized by placing them in a Petri dish resting over ice, and injections into the hemocoel were administered at the thoracic anepisternal cleft using a Nanoject III Programmable Nanoliter Injector (Drummond Scientific Company, Broomall, PA, USA).

2.2. Synthesis of Double-Stranded RNA (dsRNA) and RNA Interference (RNAi)

Double-stranded RNA was synthesized for three mosquito genes—TGase1 (AGAP009100), TGase2 (AGAP009098) and TGase3 (AGAP009099)—and one control bacterial gene—bla(ApR)—using methods previously described [19]. Specifically, gene-specific primers with T7 promoter tags (Table 1) were used to amplify a fragment of each gene using as a template cDNA from A. gambiae for transglutaminase genes or DNA from E. coli BL21(DE3) containing the pET-46 plasmid (EMD Chemicals, Gibbstown, NJ, USA) for the bacterial gene. The amplicons were subjected to agarose gel electrophoresis, and the bands of interest were excised and purified using Qiagen’s QIAquick Gel Extraction kit. The purified amplicons were then used as the template in a second PCR reaction with the same primers. The product of this reaction was purified using the QIAquick PCR Purification Kit (Qiagen, Valencia, CA, USA), and up to 1 μg of the purified product was used as the template for in vitro dsRNA synthesis using the MEGAscript T7 Kit (Applied Biosystems, Foster City, CA, USA). The dsRNA was precipitated using ethanol and re-suspended in phosphate-buffered saline (PBS). The dsRNA concentration was quantified spectrophotometrically using a BioPhotometer Plus spectrophotometer (Eppendorf AG, Hamburg, Germany), and the integrity of the dsRNA was verified by gel electrophoresis. To initiate RNAi, two-day-old female mosquitoes were injected with 300 ng of a dsRNA.

2.3. Mosquito Treatments: Uninfected and Infected with E. coli

Four days after dsRNA injection (i.e., at 6 days of adulthood), mosquitoes were subjected to one of three treatments: (i) unmanipulated (called the uninfected condition), (ii) infected with E. coli for 4 h, or (iii) infected with E. coli for 24 h. Mosquito infections were conducted using tetracycline resistant GFP-expressing E. coli (modified DH5α). These bacteria were grown overnight in Luria–Bertani’s (LB) rich nutrient medium in a 37 °C shaking incubator (New Brunswick Scientific, Edison, NJ, USA). Infection was initiated by injecting 69 nL of a dilution of the bacterial culture into the hemocoel, and the infection dose was determined by plating dilutions of the inoculum on LB agar plates and counting the colony-forming units that grew on the plates. On average, mosquitoes were infected with 15,374 GFP-E. coli. An injury group, such as a mock injection, was omitted because injury does not induce periostial hemocyte aggregation [20,22,23].

2.4. cDNA Synthesis and Quantitative RT-PCR (qPCR) for RNAi Efficiency Determination

RNA was extracted and isolated from a pool of 10 whole-body mosquitoes for each dsRNA-treatment combination. RNA was initially purified using TRIzol Reagent (Applied Biosystems, Waltham, MA, USA) and then re-purified using the RNeasy Mini Kit (Qiagen, Valencia, CA, USA). The RNA concentration was measured spectrophotometrically, and up to 5 μg of RNA was treated with RQ1 RNase-free DNase (Promega, Madison, WI, USA) and subsequently used for cDNA synthesis using Oligo(dT)20 primers and the SuperScript III First-Strand Synthesis System for RT-PCR (Applied Biosystems).
qPCR was conducted using cDNA as a template, gene-specific primers (Table 1), and Power SYBR Green PCR Master Mix (Applied Biosystems) on a CFX Connect Real-Time Detection System (Bio-Rad Laboratories, Hercules, CA, USA). The 2−ΔΔCT method was used for relative quantification, using the housekeeping gene RpS7 as the reference [36,37,38]. The mRNA fold change for each transglutaminase gene (indicative of silencing efficiency) was calculated relative to dsbla(ApR) mosquitoes of the same treatment. For each dsRNA-treatment combination, 2–3 independent trials were conducted, each with two technical replicates. These trials were conducted on some of the same mosquitoes for which heart physiology measurements were made. That is, mosquitoes were placed in TRIzol immediately after the acquisition of a video.

2.5. Measurement of Heart Contractions

Mosquitoes were placed on Sylgard 184 silicone elastomer plates (Dow Corning, Midland, MI, USA) with wings spread and ventral side down using a noninvasive procedure we have described and pictured [28,39]. After mosquitoes acclimated to room temperature, their dorsal abdomens—and the underlying heart—were viewed under bright-field trans-illumination using a Nikon SMZ1500 stereo microscope (Nikon Corporation, Tokyo, Japan) connected to a Hamamatsu ORCA-Flash 2.8 digital CMOS camera (Hamamatsu Photonics, Hamamatsu, Japan) and Nikon Advanced Research NIS-Elements software. Videos of 65 s duration were recorded, and the middle 60 s of each video were manually analyzed for heart contractions, leaving buffers of 2.5 s at the beginning and end of each video. The heart contraction rate in hertz (Hz, or contractions per second) was determined by counting the frequency of wave-like contractions of cardiac muscle, as seen through the semi-translucent abdomen [28,40]. For both the uninfected and the 24 h infection conditions, 30–40 mosquitoes were assayed per dsRNA treatment across 4 independent biological trials. For the 4 h infection condition, 15–20 mosquitoes were assayed per treatment across 2 biological trials. Some of the 24 h trials were conducted on the same mosquito batches as the 4 h trials, so some of the uninfected mosquito readings are used when reporting both the 24 h and 4 h data. Heart contraction data were analyzed by two-way ANOVA followed by Sidak’s multiple comparisons tests, using GraphPad Prism version 9.3.1. The data collected for this manuscript is contained in the Supplementary Materials (Data File S1).

3. Results

3.1. General Experimental Design and the Efficiency of RNAi-Based Silencing of Transglutaminases

Two-day-old female mosquitoes were divided into four RNAi groups: (i) dsbla(ApR), (ii) dsTGase1, (iii) dsTGase2, and (iv) dsTGase3. Four days later (at 6 days of adulthood), each dsRNA group was divided into two treatment conditions: (i) uninfected and (ii) infected with E. coli. Each infected condition was assayed at 24 h and 4 h after treatment.
To ensure that the dsRNA injections were impacting the mRNA abundance of the target genes, we determined the efficiency of RNAi by qPCR (Figure 1). RNAi efficiency varied depending on the treatment–time combination, but on average, silencing of TGase1 was 52%, TGase2 was 54% and TGase3 was 67%. This RNAi efficiency is consistent with our earlier work silencing transglutaminases [19].

3.2. RNAi-Based Silencing of Transglutaminases Has Infection-Dependent Effects on the Heart Rate at 24 h after Treatment

We began assessing whether transglutaminase silencing affects the heart rate in uninfected and infected mosquitoes at 24 h post-treatment because we have previously observed that an infection decreases the heart rate at 1-, 3- and 5-days post-infection [27]. The heart of dsbla(ApR) uninfected mosquitoes contracted at 1.96 Hz, and infection for 24 h decreased the heart contraction rate by 9.9% (Figure 2A). This infection-induced decrease is similar to what we have previously observed following infection with bacteria [27]. Similar to dsbla(ApR) mosquitoes, infection reduced the heart contraction rate of dsTGase1 mosquitoes by 8.2%. However, the infection-induced decrease in the heart rate disappeared in dsTGase2 mosquitoes (−0.8% difference) and was reversed to an infection-induced increase of 6.4% in dsTGase3 mosquitoes.
When we instead compared the effect of injecting different types of dsRNAs within a single treatment (comparing the effect of dsRNA within the uninfected, or within 24 h after infection), additional patterns emerged (Figure 2B). The heart rate of uninfected mosquitoes that had been treated with dsbla(ApR), dsTGase1 and dsTGase2 were all similar (<0.8% difference from control), but the heart rate of uninfected dsTGase3 mosquitoes was 6% lower than that of dsbla(ApR) mosquitoes. At 24 h after infection, the heart rate of dsbla(ApR) and dsTGase1 mosquitoes was similar (2% difference), but relative to infected dsbla(ApR) mosquitoes, the heart rate of infected dsTGase2 and dsTGase3 mosquitoes was elevated by 7.3% and 11.1%, respectively.
Statistical analysis of the entire 24 h dataset revealed that infection status irrespective of the type of dsRNA affects the heart rate (p = 0.0022) but that the type of dsRNA irrespective of the infection status does not (p = 0.2199). Moreover, there was a significant interaction between infection status and dsRNA treatment (p = 0.0003), indicating that the effect of infection changes when different transglutaminases are silenced by RNAi (and vice versa). Sidak’s multiple comparisons tests then revealed that dsbla(ApR) and dsTGase1 mosquitoes display the expected infection-induced reduction in the heart rate, but that this effect disappears in dsTGase2 mosquitoes and is reversed in dsTGase3 mosquitoes. In other words, there is a dichotomy in how silencing different transglutaminases affects the heart rate in uninfected and infected mosquitoes.

3.3. The Infection-Dependent Cardiac Effect of Silencing Transglutaminases Is Recapitulated at 4 h after Treatment

After determining that transglutaminase silencing has infection-dependent effects on the mosquito heart rate at 24 h following treatment, we proceeded to test whether this effect was recapitulated at 4 h following treatment. We did so because we recently showed that silencing TGase3—but not TGase1 or TGase2—affects the infection-induced aggregation of hemocytes on the surface of the heart at 4 h after treatment [19]. The heart of dsbla(ApR) uninfected mosquitoes contracted at 1.89 Hz, and infection for 4 h decreased the heart contraction rate by 8.4% (Figure 3A). Again, this infection-induced decrease in the heart rate is similar to what we have previously observed following infection with bacteria [27], and it is also similar to what we observed here at 24 h following treatment (Figure 2). Similar to dsbla(ApR) mosquitoes, infection decreased the heart rate of dsTGase1 mosquitoes by 6.3%. However, the infection-induced decrease in the heart rate of dsTGase2 mosquitoes was attenuated to 2.6% and reversed in dsTGase3 mosquitoes; the heart rate of dsTGase3 mosquitoes increased by 8.9% at 4 h after infection relative to uninfected dsTGase3 mosquitoes.
When we compared the effect of injecting different types of dsRNAs within a single immune treatment, we found that the heart rate of dsbla(ApR) and dsTGase1 mosquitoes that were uninfected was nearly identical (−0.4% difference), as was their heart rate at 4 h after infection (1.9% difference) (Figure 3B). The heart rate of dsbla(ApR) and dsTGase2 mosquitoes that were uninfected was also nearly identical (0.4% difference), but after 4 h of infection, the heart rate of dsTGase2 mosquitoes was 6.8% higher than that of infected dsbla(ApR) mosquitoes. Finally, the heart rate of dsTGase3 uninfected mosquitoes was 6.3% lower than that of dsbla(ApR) uninfected mosquitoes, yet this phenotype reversed after infection; the heart rate at 4 h post-infection was 11.4% higher in dsTGase3 mosquitoes than in dsbla(ApR) mosquitoes.
Two-way ANOVA revealed that, across the entire 4 h dataset, infection status irrespective of the type of dsRNA does not affect the heart rate (p = 0.0825), and neither does the type of dsRNA irrespective of infection status (p = 0.2749). However, there was a significant interaction between infection status and the type of dsRNA (p < 0.0001), once again indicating that the effect of infection changes when different transglutaminases are silenced by RNAi (and vice versa). The post hoc tests revealed that dsbla(ApR) mosquitoes—and perhaps dsTGase1 mosquitoes—display the expected infection-induced reduction in the heart rate, but that this effect is lost in dsTGase2 mosquitoes and reversed in dsTGase3 mosquitoes. Overall, the data collected at 4 h and 24 h are in agreement with one another. Moreover, the expected infection-induced decrease in the heart rate is evident in control mosquitoes and mosquitoes where TGase1 was disrupted but not in mosquitoes where either TGase2 or TGase3 was disrupted.

4. Discussion

Infection decreases the heart rate, and this is modulated by nitric oxide produced by heart-associated hemocytes [27,41,42]. Although multiple factors affect heart contraction dynamics [43,44,45], whether regulators of periostial hemocyte aggregation affect the heart rate remained unknown. In the present study, we confirmed that an infection decreases the heart rate and uncovered that TGase3—and to a lesser extent TGase2—affect the heart rate in an infection-dependent manner. Specifically, silencing TGase1 has no impact on heart physiology and recapitulates that an infection decreases the heart rate. Silencing TGase2 leads to a constant heart rate regardless of whether a mosquito is infected, with the heart rate always being similar to that of an uninfected mosquito; in other words, the infection-induced decrease in the heart rate does not occur. Finally, silencing TGase3 decreases the heart rate in uninfected mosquitoes but increases the heart rate in infected mosquitoes (Figure 4).
Prior to the initiation of this study, one of our working hypotheses was that neither TGase1 nor TGase2 plays a role in heart physiology, and we based this hypothesis on our earlier finding that these genes do not affect periostial hemocyte aggregation [19]. Our other working hypothesis was that TGase3 does not play a role in heart physiology in uninfected mosquitoes but dampens the infection-induced reduction in the mosquito heart rate; we based this hypothesis on our earlier finding that TGase3 does not affect periostial hemocyte aggregation in uninfected mosquitoes but negatively affects this process following infection [19]. The hypothesis on TGase1 is supported by the data but the hypotheses on TGase2 and TGase3 are not. The phenotypes observed are difficult to explain, but below, we present several plausible explanations, which are not mutually exclusive.
A possible explanation as to why the heart rate of TGase3- or TGase2-silenced mosquitoes increases following infection (relative to infected mosquitoes injected with control dsRNA) pertains to the infection-induced deposition of melanin on the heart. Melanization is a multifunctional biochemical cascade driven by phenoloxidase and other enzymes that is essential in the immune response against bacteria and parasites [46,47]. Transglutaminase and prophenoloxidase are both produced by hemocytes, and hemolymph clot formation involves both transglutaminase and phenoloxidase [7,48,49]. Moreover, activation of the melanization cascade in mosquitoes leads to the accumulation of melanin in the periostial regions of the heart [20,23,25,50,51]. We previously observed that silencing TGase3, and to a lesser extent TGase2, leads to a large increase in melanin deposition within the periostial regions of infected mosquitoes [19]. This melanin could partially block hemolymph entry into the heart via the ostia, and to compensate for the decrease in stroke volume, the heart of these mosquitoes may contract faster. A problem with this explanation is that it does not explain the heart rate phenotype for TGase3-silenced uninfected mosquitoes, or TGase3- or TGase2-silenced mosquitoes that had been infected for 4 h; there is negligible melanin on the heart of those mosquitoes. However, the observation of a more tepid phenotype when silencing TGase2 when compared to TGase3 could be due to TGase3 being more highly expressed in the heart with periostial hemocytes than TGase2, as uncovered in our RNAseq screen.
Another potential explanation pertains to downstream effects of the IMD pathway. The genome of D. melanogaster encodes one transglutaminase, and it is a negative regulator of the IMD pathway [9,10]. In mosquitoes, TGase3 is a negative regulator of periostial hemocyte aggregation [19], and because the IMD pathway is a positive regulator of periostial hemocyte aggregation [26], we hypothesized that TGase3 impacts periostial hemocyte aggregation by inhibiting the IMD pathway. The IMD pathway produces an array of antimicrobial peptides [52,53,54], and therefore, removing an inhibitor of the IMD pathway may increase antimicrobial peptide production. More antimicrobial peptides may quell an infection more quickly, which would decrease the need for nitric oxide. Nitric oxide is an antimicrobial produced by many tissues—including periostial hemocytes—and also reduces the heart rate [27,55,56,57,58], and so, a reduction in nitric oxide may eliminate or dampen the infection-induced decrease in the heart rate. Again, this explanation would be consistent with our data from infected mosquitoes, but it would not explain the phenotype observed in uninfected mosquitoes with silenced TGase3.
Other possibilities rely on comparisons with vertebrate animals. For example, transglutaminase causes age-associated vascular stiffening and dysfunction in mice, resulting in a loss of vascular compliance as this mammal ages [59]. It is possible that silencing mosquito TGase3 reduces heart stiffening, thereby increasing cardiac efficiency and reducing the heart rate in uninfected mosquitoes. Although this does not explain the phenotype observed in TGase3-silenced mosquitoes that were infected, it is in agreement with the progressive heart rate increase that takes place as a mosquito ages [32].

5. Conclusions

Here, we analyzed the entire transglutaminase gene family of A. gambiae in the context of the circulatory system. We conclude that TGase2 and TGase3 modulate heart physiology in ways that depend on the infection status of the mosquito and demonstrate that factors not classically associated with circulatory physiology are involved in the functional integration of the immune and circulatory systems of mosquitoes.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/insects13070582/s1, Data File S1.

Author Contributions

Conceptualization, J.F.H.; Methodology, J.F.H.; Validation, A.R. and J.F.H.; Formal analysis, A.R. and J.F.H.; Investigation, A.R.; Resources, J.F.H.; Data curation, A.R. and J.F.H.; Writing—original draft preparation, A.R. and J.F.H.; Writing—review and editing, A.R. and J.F.H.; Visualization, A.R. and J.F.H.; Supervision, J.F.H.; Project administration, J.F.H.; Funding acquisition, J.F.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by U.S. National Science Foundation grant IOS-1949145.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data are contained within the article and the Supplementary Material.

Acknowledgments

Experimental support by Tania Estévez-Lao and useful discussions with Cole Meier, Jordyn Sanner and Lindsay Martin are greatly appreciated.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Eckert, R.L.; Kaartinen, M.T.; Nurminskaya, M.; Belkin, A.M.; Colak, G.; Johnson, G.V.; Mehta, K. Transglutaminase regulation of cell function. Physiol. Rev. 2014, 94, 383–417. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Lorand, L.; Graham, R.M. Transglutaminases: Crosslinking enzymes with pleiotropic functions. Nat. Rev. Mol. Cell Biol. 2003, 4, 140–156. [Google Scholar] [CrossRef] [PubMed]
  3. Clark, K.D. Insect hemolymph immune complexes. Subcell Biochem. 2020, 94, 123–161. [Google Scholar] [CrossRef] [PubMed]
  4. Dziedziech, A.; Shivankar, S.; Theopold, U. Drosophila melanogaster responses against entomopathogenic nematodes: Focus on hemolymph clots. Insects 2020, 11, 62. [Google Scholar] [CrossRef] [Green Version]
  5. Baxter, R.H.; Contet, A.; Krueger, K. Arthropod innate immune systems and vector-borne diseases. Biochemistry 2017, 56, 907–918. [Google Scholar] [CrossRef] [Green Version]
  6. Shibata, T.; Kawabata, S.I. Pluripotency and a secretion mechanism of Drosophila transglutaminase. J. Biochem. 2018, 163, 165–176. [Google Scholar] [CrossRef]
  7. Schmid, M.R.; Dziedziech, A.; Arefin, B.; Kienzle, T.; Wang, Z.; Akhter, M.; Berka, J.; Theopold, U. Insect hemolymph coagulation: Kinetics of classically and non-classically secreted clotting factors. Insect Biochem. Mol. Biol. 2019, 109, 63–71. [Google Scholar] [CrossRef]
  8. Lindgren, M.; Riazi, R.; Lesch, C.; Wilhelmsson, C.; Theopold, U.; Dushay, M.S. Fondue and transglutaminase in the Drosophila larval clot. J. Insect Physiol. 2008, 54, 586–592. [Google Scholar] [CrossRef]
  9. Shibata, T.; Sekihara, S.; Fujikawa, T.; Miyaji, R.; Maki, K.; Ishihara, T.; Koshiba, T.; Kawabata, S. Transglutaminase-catalyzed protein-protein cross-linking suppresses the activity of the NF-kappaB-like transcription factor relish. Sci. Signal. 2013, 6, ra61. [Google Scholar] [CrossRef]
  10. Maki, K.; Shibata, T.; Kawabata, S.I. Transglutaminase-catalyzed incorporation of polyamines masks the DNA-binding region of the transcription factor Relish. J. Biol. Chem. 2017, 292, 6369–6380. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, Z.; Wilhelmsson, C.; Hyrsl, P.; Loof, T.G.; Dobes, P.; Klupp, M.; Loseva, O.; Morgelin, M.; Ikle, J.; Cripps, R.M.; et al. Pathogen entrapment by transglutaminase—A conserved early innate immune mechanism. PLoS Pathog. 2010, 6, e1000763. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Sekihara, S.; Shibata, T.; Hyakkendani, M.; Kawabata, S.I. RNA interference directed against the transglutaminase gene triggers dysbiosis of gut microbiota in Drosophila. J. Biol. Chem. 2016, 291, 25077–25087. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zhao, X.; Liu, L.; Zhou, W.; Cai, Q.; Huang, Q. Roles of selenoprotein T and transglutaminase in active immunization against entomopathogenic fungi in the termite Reticulitermes chinensis. J. Insect Physiol. 2020, 125, 104085. [Google Scholar] [CrossRef] [PubMed]
  14. Phillips, D.R.; Clark, K.D. Bombyx mori and Aedes aegypti form multi-functional immune complexes that integrate pattern recognition, melanization, coagulants, and hemocyte recruitment. PLoS ONE 2017, 12, e0171447. [Google Scholar] [CrossRef] [Green Version]
  15. Rogers, D.W.; Baldini, F.; Battaglia, F.; Panico, M.; Dell, A.; Morris, H.R.; Catteruccia, F. Transglutaminase-mediated semen coagulation controls sperm storage in the malaria mosquito. PLoS Biol. 2009, 7, e1000272. [Google Scholar] [CrossRef] [Green Version]
  16. Nsango, S.E.; Pompon, J.; Xie, T.; Rademacher, A.; Fraiture, M.; Thoma, M.; Awono-Ambene, P.H.; Moyou, R.S.; Morlais, I.; Levashina, E.A. AP-1/Fos-TGase2 axis mediates wounding-induced Plasmodium falciparum killing in Anopheles gambiae. J. Biol. Chem. 2013, 288, 16145–16154. [Google Scholar] [CrossRef] [Green Version]
  17. Silveira, H.; Gabriel, A.; Ramos, S.; Palma, J.; Felix, R.; Custodio, A.; Collins, L.V. CpG-containing oligodeoxynucleotides increases resistance of Anopheles mosquitoes to Plasmodium infection. Insect Biochem. Mol. Biol. 2012, 42, 758–765. [Google Scholar] [CrossRef]
  18. Le, B.V.; Nguyen, J.B.; Logarajah, S.; Wang, B.; Marcus, J.; Williams, H.P.; Catteruccia, F.; Baxter, R.H. Characterization of Anopheles gambiae transglutaminase 3 (AgTG3) and its native substrate plugin. J. Biol. Chem. 2013, 288, 4844–4853. [Google Scholar] [CrossRef] [Green Version]
  19. Yan, Y.; Ramakrishnan, A.; Estevez-Lao, T.Y.; Hillyer, J.F. Transglutaminase 3 negatively regulates immune responses on the heart of the mosquito, Anopheles gambiae. Sci. Rep. 2022, 12, 6715. [Google Scholar] [CrossRef]
  20. King, J.G.; Hillyer, J.F. Infection-induced interaction between the mosquito circulatory and immune systems. PLoS Pathog. 2012, 8, e1003058. [Google Scholar] [CrossRef] [Green Version]
  21. Yan, Y.; Hillyer, J.F. The immune and circulatory systems are functionally integrated across insect evolution. Sci. Adv. 2020, 6, eabb3164. [Google Scholar] [CrossRef] [PubMed]
  22. King, J.G.; Hillyer, J.F. Spatial and temporal in vivo analysis of circulating and sessile immune cells in mosquitoes: Hemocyte mitosis following infection. BMC Biol. 2013, 11, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sigle, L.T.; Hillyer, J.F. Mosquito hemocytes preferentially aggregate and phagocytose pathogens in the periostial regions of the heart that experience the most hemolymph flow. Dev. Comp. Immunol. 2016, 55, 90–101. [Google Scholar] [CrossRef] [Green Version]
  24. Sigle, L.T.; Hillyer, J.F. Eater and draper are involved in the periostial haemocyte immune response in the mosquito Anopheles gambiae. Insect Mol. Biol. 2018, 27, 429–438. [Google Scholar] [CrossRef] [PubMed]
  25. Yan, Y.; Hillyer, J.F. Complement-like proteins TEP1, TEP3 and TEP4 are positive regulators of periostial hemocyte aggregation in the mosquito Anopheles gambiae. Insect Biochem. Mol. Biol. 2019, 107, 1–9. [Google Scholar] [CrossRef] [PubMed]
  26. Yan, Y.; Sigle, L.T.; Rinker, D.C.; Estevez-Lao, T.Y.; Capra, J.A.; Hillyer, J.F. The IMD and JNK pathways drive the functional integration of the immune and circulatory systems of mosquitoes. bioRxiv 2022. [Google Scholar] [CrossRef]
  27. Estevez-Lao, T.Y.; Sigle, L.T.; Gomez, S.N.; Hillyer, J.F. Nitric oxide produced by periostial hemocytes modulates the bacterial infection-induced reduction of the mosquito heart rate. J. Exp. Biol. 2020, 223, jeb225821. [Google Scholar] [CrossRef]
  28. Glenn, J.D.; King, J.G.; Hillyer, J.F. Structural mechanics of the mosquito heart and its function in bidirectional hemolymph transport. J. Exp. Biol. 2010, 213, 541–550. [Google Scholar] [CrossRef] [Green Version]
  29. Sigle, L.T.; Hillyer, J.F. Structural and functional characterization of the contractile aorta and associated hemocytes of the mosquito Anopheles gambiae. J. Exp. Biol. 2018, 221, jeb181107. [Google Scholar] [CrossRef] [Green Version]
  30. Barbosa da Silva, H.; Godoy, R.S.M.; Martins, G.F. The basic plan of the adult heart Is conserved across different species of adult mosquitoes, but the morphology of heart-associated tissues varies. J. Med. Entomol. 2019, 56, 984–996. [Google Scholar] [CrossRef]
  31. Leodido, A.C.M.; Ramalho-Ortigao, M.; Martins, G.F. The ultrastructure of the Aedes aegypti heart. Arthropod. Struct. Dev. 2013, 42, 539–550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Doran, C.R.; Estevez-Lao, T.Y.; Hillyer, J.F. Mosquito aging modulates the heart rate and the proportional directionality of heart contractions. J. Insect. Physiol. 2017, 101, 47–56. [Google Scholar] [CrossRef]
  33. Hillyer, J.F.; Estevez-Lao, T.Y.; Mirzai, H.E. The neurotransmitters serotonin and glutamate accelerate the heart rate of the mosquito Anopheles gambiae. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 2015, 188, 49–57. [Google Scholar] [CrossRef]
  34. Ellison, H.E.; Estevez-Lao, T.Y.; Murphree, C.S.; Hillyer, J.F. Deprivation of both sucrose and water reduces the mosquito heart contraction rate while increasing the expression of nitric oxide synthase. J. Insect Physiol. 2015, 74, 1–9. [Google Scholar] [CrossRef] [PubMed]
  35. Hillyer, J.F.; Estevez-Lao, T.Y.; de la Parte, L.E. Myotropic effects of FMRFamide containing peptides on the heart of the mosquito Anopheles gambiae. Gen. Comp. Endocrinol. 2014, 202, 15–25. [Google Scholar] [CrossRef]
  36. Estevez-Lao, T.Y.; Boyce, D.S.; Honegger, H.W.; Hillyer, J.F. Cardioacceleratory function of the neurohormone CCAP in the mosquito Anopheles gambiae. J. Exp. Biol. 2013, 216, 601–613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Coggins, S.A.; Estevez-Lao, T.Y.; Hillyer, J.F. Increased survivorship following bacterial infection by the mosquito Aedes aegypti as compared to Anopheles gambiae correlates with increased transcriptional induction of antimicrobial peptides. Dev. Comp. Immunol. 2012, 37, 390–401. [Google Scholar] [CrossRef] [PubMed]
  38. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  39. Andereck, J.W.; King, J.G.; Hillyer, J.F. Contraction of the ventral abdomen potentiates extracardiac retrograde hemolymph propulsion in the mosquito hemocoel. PLoS ONE 2010, 5, e12943. [Google Scholar] [CrossRef] [Green Version]
  40. Hillyer, J.F.; Estevez-Lao, T.Y.; Funkhouser, L.J.; Aluoch, V.A. Anopheles gambiae corazonin: Gene structure, expression and effect on mosquito heart physiology. Insect Mol. Biol. 2012, 21, 343–355. [Google Scholar] [CrossRef]
  41. Da Silva, R.; da Silva, S.R.; Lange, A.B. The regulation of cardiac activity by nitric oxide (NO) in the Vietnamese stick insect, Baculum extradentatum. Cell Signal 2012, 24, 1344–1350. [Google Scholar] [CrossRef] [PubMed]
  42. Broderick, K.E.; Feala, J.; McCulloch, A.; Paternostro, G.; Sharma, V.S.; Pilz, R.B.; Boss, G.R. The nitric oxide scavenger cobinamide profoundly improves survival in a Drosophila melanogaster model of bacterial sepsis. FASEB J. 2006, 20, 1865–1873. [Google Scholar] [CrossRef] [PubMed]
  43. Chowanski, S.; Lubawy, J.; Urbanski, A.; Rosinski, G. Cardioregulatory functions of neuropeptides and peptide hormones in insects. Protein Pept. Lett. 2016, 23, 913–931. [Google Scholar] [CrossRef]
  44. Hillyer, J.F. Insect heart rhythmicity is modulated by evolutionarily conserved neuropeptides and neurotransmitters. Curr. Opin. Insect Sci. 2018, 29, 41–48. [Google Scholar] [CrossRef] [PubMed]
  45. Hillyer, J.F.; Pass, G. The insect circulatory system: Structure, function, and evolution. Annu. Rev. Entomol. 2020, 65, 121–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Christensen, B.M.; Li, J.; Chen, C.C.; Nappi, A.J. Melanization immune responses in mosquito vectors. Trends Parasitol. 2005, 21, 192–199. [Google Scholar] [CrossRef]
  47. Whitten, M.M.A.; Coates, C.J. Re-evaluation of insect melanogenesis research: Views from the dark side. Pigment. Cell Melanoma Res. 2017, 30, 386–401. [Google Scholar] [CrossRef]
  48. Dziedziech, A.; Schmid, M.; Arefin, B.; Kienzle, T.; Krautz, R.; Theopold, U. Data on Drosophila clots and hemocyte morphologies using GFP-tagged secretory proteins: Prophenoloxidase and transglutaminase. Data Brief 2019, 25, 104229. [Google Scholar] [CrossRef]
  49. Khalili, D.; Kalcher, C.; Baumgartner, S.; Theopold, U. Anti-fibrotic activity of an antimicrobial peptide in a Drosophila model. J. Innate Immun. 2021, 13, 376–390. [Google Scholar] [CrossRef]
  50. Sousa, G.L.; Bishnoi, R.; Baxter, R.H.G.; Povelones, M. The CLIP-domain serine protease CLIPC9 regulates melanization downstream of SPCLIP1, CLIPA8, and CLIPA28 in the malaria vector Anopheles gambiae. PLoS Pathog. 2020, 16, e1008985. [Google Scholar] [CrossRef]
  51. Yassine, H.; Kamareddine, L.; Chamat, S.; Christophides, G.K.; Osta, M.A. A serine protease homolog negatively regulates TEP1 consumption in systemic infections of the malaria vector Anopheles gambiae. J. Innate Immun. 2014, 6, 806–818. [Google Scholar] [CrossRef] [PubMed]
  52. Hillyer, J.F. Insect immunology and hematopoiesis. Dev. Comp. Immunol. 2016, 58, 102–118. [Google Scholar] [CrossRef] [Green Version]
  53. Myllymaki, H.; Valanne, S.; Ramet, M. The Drosophila IMD signaling pathway. J. Immunol. 2014, 192, 3455–3462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kleino, A.; Silverman, N. The Drosophila IMD pathway in the activation of the humoral immune response. Dev. Comp. Immunol. 2014, 42, 25–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Hillyer, J.F.; Estevez-Lao, T.Y. Nitric oxide is an essential component of the hemocyte-mediated mosquito immune response against bacteria. Dev. Comp. Immunol. 2010, 34, 141–149. [Google Scholar] [CrossRef] [PubMed]
  56. Herrera-Ortiz, A.; Martinez-Barnetche, J.; Smit, N.; Rodriguez, M.H.; Lanz-Mendoza, H. The effect of nitric oxide and hydrogen peroxide in the activation of the systemic immune response of Anopheles albimanus infected with Plasmodium berghei. Dev. Comp. Immunol. 2011, 35, 44–50. [Google Scholar] [CrossRef]
  57. Luckhart, S.; Vodovotz, Y.; Cui, L.; Rosenberg, R. The mosquito Anopheles stephensi limits malaria parasite development with inducible synthesis of nitric oxide. Proc. Natl. Acad. Sci. USA 1998, 95, 5700–5705. [Google Scholar] [CrossRef] [Green Version]
  58. Powers, J.C.; Turangan, R.; Joosse, B.A.; Hillyer, J.F. Adult mosquitoes infected with bacteria early in life have stronger antimicrobial responses and more hemocytes after reinfection later in life. Insects 2020, 11, 331. [Google Scholar] [CrossRef]
  59. Wang, H.; Chen, J.; Jandu, S.; Melucci, S.; Savage, W.; Nandakumar, K.; Kang, S.K.; Barreto-Ortiz, S.; Poe, A.; Rastogi, S.; et al. Probing tissue transglutaminase mediated vascular smooth muscle cell aging using a novel transamidation-deficient Tgm2-C277S mouse model. Cell Death Discov. 2021, 7, 197. [Google Scholar] [CrossRef]
Insects 13 00582 g001 550
Figure 1. RNAi-based knockdown efficiency of TGase1, TGase2 and TGase3 in uninfected and infected mosquitoes. Graphs display the relative mRNA abundance in mosquitoes treated with bla(ApR), TGase1, TGase2 or TGase3 dsRNA. Mosquitoes were uninfected or had been infected with E. coli for 4 or 24 h, and mRNA abundance was calculated relative to the bla(ApR) for each treatment. Column heights mark the mean, and whiskers indicate the standard error of the mean (S.E.M.).
Figure 1. RNAi-based knockdown efficiency of TGase1, TGase2 and TGase3 in uninfected and infected mosquitoes. Graphs display the relative mRNA abundance in mosquitoes treated with bla(ApR), TGase1, TGase2 or TGase3 dsRNA. Mosquitoes were uninfected or had been infected with E. coli for 4 or 24 h, and mRNA abundance was calculated relative to the bla(ApR) for each treatment. Column heights mark the mean, and whiskers indicate the standard error of the mean (S.E.M.).
Insects 13 00582 g001
Insects 13 00582 g002 550
Figure 2. Effect of RNAi-based silencing of TGase1, TGase2 or TGase3 on the heart rate at 24 h after treatment. Graphs show the average heart contraction rates in mosquitoes treated with bla(ApR), TGase1, TGase2 or TGase3 dsRNA, and were either uninfected or infected with E. coli for 24 h. The same data are presented in two configurations: panel (A) groups the data by dsRNA treatment whereas panel (B) groups the data by infection status. Data were analyzed by two-way ANOVA followed by Sidak’s post hoc test. Column heights mark the means, and whiskers show the S.E.M. Asterisks (*) indicate post hoc p < 0.05.
Figure 2. Effect of RNAi-based silencing of TGase1, TGase2 or TGase3 on the heart rate at 24 h after treatment. Graphs show the average heart contraction rates in mosquitoes treated with bla(ApR), TGase1, TGase2 or TGase3 dsRNA, and were either uninfected or infected with E. coli for 24 h. The same data are presented in two configurations: panel (A) groups the data by dsRNA treatment whereas panel (B) groups the data by infection status. Data were analyzed by two-way ANOVA followed by Sidak’s post hoc test. Column heights mark the means, and whiskers show the S.E.M. Asterisks (*) indicate post hoc p < 0.05.
Insects 13 00582 g002
Insects 13 00582 g003 550
Figure 3. Effect of RNAi-based silencing of TGase1, TGase2 or TGase3 on the heart rate at 4 h after treatment. Graphs show the average heart contraction rates in mosquitoes treated with bla(ApR), TGase1, TGase2 or TGase3 dsRNA, and were either uninfected or infected with E. coli for 4 h. The same data are presented in two configurations: panel (A) groups the data by dsRNA treatment whereas panel (B) groups the data by infection status. Data were analyzed by two-way ANOVA followed by Sidak’s post hoc test. Column heights mark the means, and whiskers show the S.E.M. Asterisks (*) indicate post hoc p < 0.05.
Figure 3. Effect of RNAi-based silencing of TGase1, TGase2 or TGase3 on the heart rate at 4 h after treatment. Graphs show the average heart contraction rates in mosquitoes treated with bla(ApR), TGase1, TGase2 or TGase3 dsRNA, and were either uninfected or infected with E. coli for 4 h. The same data are presented in two configurations: panel (A) groups the data by dsRNA treatment whereas panel (B) groups the data by infection status. Data were analyzed by two-way ANOVA followed by Sidak’s post hoc test. Column heights mark the means, and whiskers show the S.E.M. Asterisks (*) indicate post hoc p < 0.05.
Insects 13 00582 g003
Insects 13 00582 g004 550
Figure 4. Diagrammatic representation of the mosquito circulatory system and the process of periostial hemocyte aggregation, together with the effect of silencing transglutaminase genes on the heart contraction rate of uninfected and infected mosquitoes.
Figure 4. Diagrammatic representation of the mosquito circulatory system and the process of periostial hemocyte aggregation, together with the effect of silencing transglutaminase genes on the heart contraction rate of uninfected and infected mosquitoes.
Insects 13 00582 g004
Table
Table 1. Gene names, gene IDs, and primers used in this study.
Table 1. Gene names, gene IDs, and primers used in this study.
GeneVectorBase
ID a		ApplicationSequence (Forward and Reverse) bAmplicon (bp) c
TranscriptGenomic
RpS7AGAP010592qPCRGACGGATCCCAGCTGATAAA132281
GTTCTCTGGGAATTCGAACG
TGase1AGAP009100qPCRCTGCACAAGGGACTGTTCCA191259
AACGCCAAAAAGCCATCCAC
TGase2AGAP009098qPCRCGGTGGACGCTGACTATCAA225297
GTACTGGCCGAGCTTCCATT
TGase3AGAP009099qPCRTACAGCAGCCAGCGGTTTAG236236
ATATCGCGCCCAGTGTAGTC
bla(ApR)(Bacterial gene)RNAiTAATACGACTCACTATAGGGCCGAGCGCAGAAGTGGT214214
TAATACGACTCACTATAGGGAACCGGAGCTGAATGAA
TGase1AGAP009100RNAiTAATACGACTCACTATAGGGCATTCCGGTTAATCAGT361433
TAATACGACTCACTATAGGGCGTAGTCGATTGTAAGA
TGase2AGAP009098RNAiTAATACGACTCACTATAGGGTCAGAGCTGTCTAACAAA490490
TAATACGACTCACTATAGGCGTACCGCTCAACTCC
TGase3AGAP009099RNAiTAATACGACTCACTATAGGGAAAACCTTCCACACGTC501501
TAATACGACTCACTATAGGGTTGAACAGCACAAACAA
a Vectorbase IDs were obtained from the AgamP4 assembly in www.vectorbase.org (exception: bla(ApR)) (website last accessed 1 June 2022). b Underlined sequences are specific to the T7 RNA polymerase promoter sites needed for dsRNA synthesis. c Amplicon sizes are based on the sequences in Vectorbase. For dsRNA primers, amplicon sizes include the T7 promoter sequence tags.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ramakrishnan, A.; Hillyer, J.F. Silencing Transglutaminase Genes TGase2 and TGase3 Has Infection-Dependent Effects on the Heart Rate of the Mosquito Anopheles gambiae. Insects 2022, 13, 582. https://doi.org/10.3390/insects13070582

AMA Style

Ramakrishnan A, Hillyer JF. Silencing Transglutaminase Genes TGase2 and TGase3 Has Infection-Dependent Effects on the Heart Rate of the Mosquito Anopheles gambiae. Insects. 2022; 13(7):582. https://doi.org/10.3390/insects13070582

Chicago/Turabian Style

Ramakrishnan, Abinaya, and Julián F. Hillyer. 2022. "Silencing Transglutaminase Genes TGase2 and TGase3 Has Infection-Dependent Effects on the Heart Rate of the Mosquito Anopheles gambiae" Insects 13, no. 7: 582. https://doi.org/10.3390/insects13070582

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop