- Open Access
Precopulatory acoustic interactions of the New World malaria vector Anopheles albimanus (Diptera: Culicidae)
Parasites & Vectors volume 12, Article number: 386 (2019)
Anopheles albimanus is a malaria vector in Central America, northern South America and the Caribbean. Although a public health threat, An. albimanus precopulatory mating behaviors are unknown. Acoustics play important roles in mosquito communication, where flight tones allow males to detect and attract potential mates. The importance of sound in precopulatory interactions has been demonstrated in Toxorhynchites brevipalpis, Aedes aegypti, Culex quinquefasciatus and Anopheles gambiae; convergence in a shared harmonic of the wing beat frequency (WBF) during courtship is thought to increase the chance of copulation. To our knowledge, An. albimanus precopulatory acoustic behaviors have not been described to date. Here, we characterized An. albimanus (i) male and female flight tones; (ii) male–female precopulatory acoustic interactions under tethered and free flight conditions; and (iii) male-male acoustic interactions during free flight.
We found significant increases in the WBFs of both sexes in free flight compared to when tethered. We observed harmonic convergence between 79% of tethered couples. In free flight, we identified a female-specific behavior that predicts mate rejection during male mating attempts: females increase their WBFs significantly faster during mate rejection compared to a successful copulation. This behavior consistently occurred during mate rejection regardless of prior mating attempts (from the same or differing male). During group flight, males of An. albimanus displayed two distinct flying behaviors: random flight and a swarm-like, patterned flight, each associated with distinct acoustic characteristics. In the transition from random to patterned flight, males converged their WBFs and significantly decreased flight area, male-male proximity and the periodicity of their trajectories.
We show that tethering of An. albimanus results in major acoustic differences compared to free flight. We identify a female-specific behavior that predicts mate rejection during male mating attempts in this species and show that male groups in free flight display distinct flying patterns with unique audio and visual characteristics. This study shows that An. albimanus display acoustic features identified in other mosquito species, further suggesting that acoustic interactions provide worthwhile targets for mosquito intervention strategies. Our results provide compelling evidence for swarming in this species and suggests that acoustic signaling is important for this behavior.
Malaria is caused by Plasmodium parasites that are transmitted through the bites of infected Anopheles mosquito vectors and represents a major disease that afflicts the inhabitants of tropical countries. Although most lethal malaria cases are concentrated in sub-Saharan Africa , numerous cases of the disease occur in South America . In Colombia, an estimated 12 million people live in malaria endemic areas and are at risk of infection . In 2017, ~53,000 cases of malaria were reported in Colombia , with the majority originating in areas near the Atlantic and Pacific coasts. Despite the implementation of various control methods, which include indoor and outdoor insecticide spraying and the use of insecticide impregnated bednets, malaria continues to be a major public health issue in Colombia [5, 6].
Anopheles albimanus (subgenus Nyssorhynchus) is a vector of malaria throughout Central America, the northern portion of South America and the Caribbean . In Colombia, An. albimanus is one of three primary vectors of malaria . Despite its status as an important vector in the Americas, there exists a large gap in our knowledge of the behavior, ecology and biology of this species, making the development and/or improvement of control methods difficult. Much of our understanding of anopheline reproduction is inferred from studies in An. gambiae (subgenus Cellia), whose post-mating reproductive biology differs substantially from An. albimanus .
Mosquito vectorial capacity is dependent on biological factors that include population density, female blood-feeding rate and female longevity , each related to successful reproduction [8, 10, 11]. Thus, processes of reproduction offer promising targets to control mosquito vectors and, by extension, the diseases they spread [12,13,14]. Contemporary mosquito control methods include the mass release of transgenic or sterile males to suppress native mosquito populations [15, 16]. Therefore, mass-released males need to compete with wild males for mating opportunities, which require males to detect and successfully interact with females. In Aedes [17, 18], Culex  and Anopheles species [20, 21], acoustic signals are used by males to locate and recognize females. Males detect nearby flying females of the same species by the sound of their wing beat frequency (WBF) [20, 22,23,24], pursue them, and appear to acoustically interact with them at a local scale [25, 26] prior to copulation. It has also been suggested that females use acoustic signals as an indicator of male fitness [27, 28].
In several Anopheles species, the initiation of reproductive behavior begins when males form a swarm . In African vectors such as An. funestus, An. gambiae and An. arabiensis, swarm formation begins when a few males begin to fly simultaneously [30,31,32,33]; additional males rapidly join the flight, increasing the swarm size and density [29, 34]. At the initiation of swarm formation, males fly in a non-specific circular motion before forming a more tightly patterned, cohesive group . Once formed, males either mate with females that penetrate the swarm or depart the swarm to mate with a female flying nearby . To date, acoustic interactions between swarming males have not been described for Anopheles species. However, cohesive acoustic behavior between multiple tethered Ae. aegypti males has been observed : male groups exhibit similar flight tones, a phenomenon more evident in larger groups , although both frequency convergence and divergence behaviors are observed between pairs of tethered Ae. aegypti males in close proximity , showing how variable such interactions can be. Frequency divergence has also been reported between tethered Cx. quinquefasciatus males . Characterization of male-male acoustic interactions may be key in understanding swarm formation and swarm cohesion.
Studies describing acoustic interactions between opposite-sex mosquitoes are more numerous. In several species, female flight tones act as a mating call that attracts males of the same species [23, 38, 39]. In experiments with tethered Ae. aegypti , Cx. quinquefasciatus  and An. gambiae [20, 21], males and females modulate their WBFs to match each other during courtship, a phenomenon known as harmonic convergence. Such convergence does not happen in the WBF of the flight tones but instead in a shared harmonic [17, 21]. An additional acoustic behavior, characterized by the rapid frequency modulation (RFM) of the WBF, was first detected in Cx. quinquefasciatus  and later described in An. coluzzi , An. gambiae  and Ae. aegypti . RFM is observed after a steep increase of the male WBF when in close proximity to a tethered female or an artificial tone source [23, 38, 39]. In all species studied, RFM occurs prior to copulation. In Ae. aegypti, RFM coincides with initial contact of a female by a courting male and ends prior to copulation or separation of the pair . Moreover, RFM and harmonic convergence appear to be independent events [25, 26, 40].
To date, flight acoustics and male–female pre-mating acoustic interactions have not been investigated in An. albimanus. Furthermore, it is unknown if males of this species swarm in the wild. To gain insight into An. albimanus precopulatory acoustic behaviors, we first characterized flight acoustics in this species by determining male and female WBFs when tethered and in free flight. We observed a significant increase in the WBF of each sex when in free flight compared to when tethered. We next characterized male–female acoustic interactions when tethered and in free flight during courtship. We found that 79% of tethered couples display harmonic convergence. Using criteria established for tethered couples, we did not detect harmonic convergence in free fight, possibly due to the brevity of mating attempts (1.38 ± 0.32 s) and the absence of paired flight after couples made contact (when harmonic convergence likely occurs; ). However, we observed that during mate rejection, females increase their WBFs significantly faster than those that mate. Finally, using audio and visual analysis, we show that male groups in free flight change flight trajectories and match flight tones during a stereotypic, patterned flight, which could be suggestive of swarm-like behavior in An. albimanus. Our results represent important first steps toward dissecting the subtle, pre-mating behavioral interactions that occur in this important malaria vector.
Anopheles albimanus were obtained from the colony maintained at Programa de Estudios y Control de Enfermedades Tropicales (PECET), Universidad de Antioquia, Colombia. These mosquitoes were originally collected in Santa Rosa, Bolivar, Colombia, and have been kept in colony since 1995. Mosquitoes were reared and maintained in a walk-in climate chamber at 27 °C and 80% relative humidity (RH) under a 12:12 h light:dark photoperiod. Eggs were collected in 473-ml cups filled with dH2O. Upon hatching, larvae were individually transferred to trays (~60 per tray) containing 0.5 l of dH2O. Larvae were fed daily with 20 mg of fish food (Tetracolor; Tetra, Melle, Germany). Pupae were individually transferred to 15-ml vials to ensure virginity. Upon eclosion, adults were separated by sex, transferred to 4-l plastic cages, and held in same-sex groups until experiments commenced. Mosquitoes had access to 20% sucrose ad libitum. We used 4–6-day-old adults in all experiments.
Audio recording setup
Recordings were performed in a soundproof chamber that contained two particle velocity microphones (NR-23158-000, Knowles; Itasca, IL, USA). Microphone signals were amplified and digitalized by an USB audio interface (M-Track Quad Four Channel Audio; M-Audio, Cumberland, USA). Mosquito flight tones were processed using Matlab (R2016a, The Mathworks Inc, Natick, USA) at a sample rate of 11025 Hz/24 bits.
Males and females were recorded to determine their individual WBF distribution. Virgin mosquitoes were anesthetized on ice and tethered to a human hair using a cyanoacrylate-based adhesive (Super Pega Infinita, Medellin, Colombia). The hair was attached on one end to an insect pin while the opposite end was attached to the pronotum (dorsal section of the thorax) as in . Before recording, the position of the hair was checked to ensure natural wing movement. To record flight tones, mosquitoes were placed ~1 cm above a single particle velocity microphone. When necessary, flight was triggered by blowing or by gently moving the legs of the individual.
Using tethered individuals as described above, male–female interactions (couples) were recorded to determine if An. albimanus display the precopulatory behavior of harmonic convergence . One mosquito was placed 1 cm above one microphone and kept stationary. A tethered individual of the opposite sex, placed ~1 cm above a second microphone, was brought in and out of the near field of the first, stationary mosquito for intervals of 10 s with 5 s of rest between intervals. After repeating 3 times, the positions of the individuals were switched. Wing length of males and females were measured as previously described  to estimate individual size. All experiments were performed at 26 ± 2 °C. Tone analysis was performed using recordings that had > 10 s of continuous flight.
Since mosquito flight tone has been described as a harmonic sound , a sinusoidal-harmonic model was used to extract the WBF information of the recorded signals. From the spectrogram [fast Fourier transform (FFT) length of 4096 points, hamming window of 80 ms and 50% overlapping], we extracted the WBF by tracking the lower frequency peak of the spectrogram that persisted for at least 1 s and fulfilled two conditions: between consecutive windows, the frequency and power did not change more than ± 20 Hz and ± 10 dB FS, respectively. For individual flights, we analyzed WBF mean and WBF range of both sexes. For male–female interactions, we identified harmonic convergence as in Aldersley et al. . Comparing the extent of frequency modulation and frequency distribution prior to and during interactions, we identified events where there was an active modulation of the WBF in response to the conspecific flight tone. To label convergence events as positive or negative, we used a time-varying fundamental frequency ratio between females and males (♀WBF/♂WBF). For instance, a ratio of 0.667 means that the female’s third harmonic is converging with the males second. We set a positive convergence as an event of duration greater than 1 s, in which the convergence ratio was within an interval of tolerance ± 1.5%. The time between individual introduction and the onset of convergence was noted as latency . Analysis was performed using Matlab (R2016a, The Mathworks Inc, Natick, USA)
Audio and visual recording setup
Audio and video of free-flying mosquitoes were simultaneously recorded in an experimental arena (20 × 20 × 25 cm transparent plastic box) placed within a sound-proof chamber. The experimental arena contained four electret microphones (Knowles FG-23329-C05), three placed on the side and back walls and one on the bottom of the box. Microphone signals were amplified by an operational amplifier (INA128; Texas Instruments, Dallas, USA) and digitalized by a 16-channel digital-analog converter (779676-01; National Instruments, Austin, USA) at a sample rate of 11025/16 bits. For video recording, a webcam (HD Pro Webcam C920; Logitech, Lausanne, Switzerland) was placed 30 cm in front of the arena to follow mosquito movement.
Three types of recordings were performed. (i) We recorded individual males and females inside the arena (n = 30 for each sex). (ii) Male–female acoustic interactions were recorded by introducing two males into the arena followed by one female. Flight tones were continuously recorded until copulation was observed or the individuals stopped flying. Mosquitoes were removed after a successful copulation. Successful copulations were verified by dissecting females to determine insemination status (i.e. presence/absence of sperm in the spermathecae). A mating attempt was counted as successful if a couple locked genitalia for more than 5 s and sperm was transferred to the female. An attempt was considered a rejection when the male and female genitalia did not come into contact or the female did not allow the male to lock genitalia. As females often rejected courting males, we recorded 167 total mating attempts from 42 different females. Recordings from 9 females were discarded, as it was not possible to clearly identify the frequency components of each sex. (iii) Acoustic interactions of male groups were recorded by introducing 8 mosquitoes into the arena 5 h before recording started. Although 8 mosquitoes were used in every experiment, only groups of 2, 3 and 4 mosquitoes flying simultaneously were analyzed. For each experiment, recordings (2.5 h in length) were segmented into 10 min sections. When a specific flying behavior was observed visually, unique fragments of 5 s were used in our analysis (for the respective number of males flying). All recordings were performed at 26 ± 2 °C between 16:00 and 18:30, corresponding to the 2 h prior to sunset.
Audio signal analysis
Signals from each of the four microphones were used to extract frequency components of the flight tones. By obtaining the cross spectrum among the microphones (FFT length of 4096, window length 0.08 s, 50% overlapping) , we generated a high order spectrogram which was used to analyze the signals continuously, regardless of the distance between mosquitoes and microphones. Depending on the interaction analyzed, we extracted different frequency features from the spectrograms. (i) For male and female individual flight, the WBF distribution and range was determined. (ii) For male–female interactions, we visually selected relevant steps of mating and analyzed the WBF ratio (♀WBF/♂WBF) 1–1.5 s prior to a mating attempt, the rate of WBF increase during a mating interaction and the response latency of the WBF increase of the female following that of the male. The change in frequency (ΔF = Fmax − F0) was measured as the difference between the initial frequency at the onset of the increase (F0) and the maximum frequency reached (Fmax). For females, however, the change in frequency was divided by two, as measurements where extracted from the second harmonic (ΔF = (Fmax − F0)/2). The duration of the increase was measured by using the time at the same two points (Δt = tmax − t0). Finally, we calculated the rate of increase (ΔF/Δt) of each sex. (iii) For male groups, we measured the average difference between the highest and the lowest WBFs, referred hereafter as frequency difference.
Audio recording from the microphone integrated to the camera was used to synchronize audio and video. We used a 1 kHz pure tone to synchronize the camera and electret microphones. Once audio and video were synchronized, we extracted 5 s video fragments from recordings in which the specific behavior was observed. Each fragment was segmented into individual frames and each frame converted to grayscale and adjusted through the equalization of its histogram to improve the contrast of the image. Finally, for each frame, X and Y coordinates of each mosquito were obtained manually to reconstruct flight trajectories using Matlab (R2016a, The Mathworks Inc, Natick, USA). Trajectories were represented by ellipse due to the natural flight of the individuals, following previous methodologies . In brief, ellipse area was calculated to represent the area covered by the movement of each mosquito. The center of the ellipse corresponded to the mean position of each mosquito trajectory (the X and Y coordinates) during the 5 s analyzed. The horizontal and vertical diameters of each ellipse were calculated as 4 times the standard deviation of each trajectory point of each individual. The average distance of the centroids of each ellipse in relation to the others was calculated to assess the extent of the aggregation of the group. Finally, the loop period of each trajectory was calculated to assess the periodicity of the movement, as in Gibson .
For tethered mosquitoes, we performed a linear regression analysis to determine the relationship between wing length (as a size index) and WBF. To compare frequency distributions and ranges of tethered and free-flying individuals, we performed a t-test and a factor analysis. As factors, we evaluated sex and recording method. Residuals of t-tests and factorial models were tested for normality, homogeneity of variance and independence using Shapiro-Wilks, Bartlett and Durbin-Watson tests, respectively. Aspects influencing male–female interactions during free flight trials were assessed using the non-parametric Mann–Whitney U-test. The frequency difference broadcast by male groups and the differences between the characteristics of flight trajectories were assessed using a Mann–Whitney-U test. Results are presented as mean ± SEM. Sample sizes are indicated for each experiment.
Flight tones of tethered and free-flying An. albimanus males and females
To characterize An. albimanus flight tones, virgin males and females were individually recorded while tethered or in free flight. When tethered, we observed a WBF of 524.11 ± 63.73 Hz (mean ± SEM, n = 78) for males and 368.91 ± 34.28 Hz (n = 116) for females (Fig. 1a). There was a significant difference in the WBF between males and females (t-test: t(192) = − 21.94, P < 0.01). By measuring the wing length of males (2.78 ± 0.02 mm) and females (2.89 ± 0.02 mm), we did not find evidence of a relationship between individual size and WBF (Linear regression; males: R2 < 0.01, Pearson’s r(78) = 0.06, P = 0.58; females: R2 = 0.03, Pearson’s r(116) = − 0.20, P = 0.05; Additional file 1: Figure S1a). Thus, size did not appear to significantly influence the WBF in tethered individuals of either sex.
Because various studies of insect flight have reported that WBFs of tethered individuals differ compared to individuals in free flight [4, 25, 44, 45], we sought to determine how An. albimanus acoustics might differ in this regard by examining mosquitoes in flight. Males and females were individually introduced into a soundproof box and their flight tones recorded. During free flight, we observed a WBF of 713.76 ± 43.16 Hz (n = 30) for males, and 430.34 ± 28.48 Hz (n = 30) for females (Fig. 1a). A significant increase in flight tone was observed for individuals in free flight compared to tethered individuals (t-test; males: t(106) = − 9.03, P < 0.01; females: t(144) = − 15.00, P < 0.01). Similar to tethered individuals, we did not observe a strong relationship between mosquito size (males: 2.83 ± 2.02 mm; females: 3.05 ± 220.53 mm) and WBF (Linear regression; males: R2 = 0.12, Pearson’s r(45) = 0.35, P = 0.02; females: R2 < 0.01, Pearson’s r(54) = 0.03, P = 0.81; Additional file 1: Figure S1b).
When analyzing each spectrogram, we noted that the tones of individuals were characterized by larger modulations of the WBF when in free flight compared to when tethered (Fig. 1b). To evaluate the effect of tethering on frequency modulation, we randomly selected a sample of 30 tethered males and 30 tethered females to compare them with 30 males and 30 females in free flight by performing a factorial ANOVA, using the recording method and sex as factors. Sex (factorial ANOVA: F(1, 116)=13.1, P < 0.01) and tethering (factorial ANOVA: F(1, 116) = 37.57, P < 0.01) each affected the range of modulation, and an interaction between these two factors was observed (factorial ANOVA: F(1, 116) = 4.66, P = 0.03). Thus, when males and females are tethered, the range of their WBF modulation is similar (males: 53.33 ± 7.89 Hz; females: 45.25 ± 3.14 Hz) but is significantly increased (males: 99.10 ± 5.75 Hz; females: 67.17 ± 4.67 Hz) during free flight (Fig. 1c).
Harmonic convergence between tethered An. albimanus males and females
Since harmonic convergence has been previously described using tethered mosquitoes of other species [17, 19, 20], we characterized harmonic convergence using tethered male and female An. albimanus. We observed harmonic convergence in 41 of 52 couples (78.8%). By analyzing the ratio of the female WBF with respect to the males’, we observed two distinct groups where convergence was observed: in 21 couples (40.4%), the female’s third harmonic converged with the male’s second (Fig. 2, left panel), while in 20 cases (38.5%), the female’s fourth harmonic converged with the male’s third (Fig. 2, right panel). The average time of convergence was 2.00 ± 0.81 s and the response latency ranged from 4 to 13 s. The remaining couples did not match harmonics. We did not detect a relationship between the ratio of couple size (♂size/♀size) and the ratio of convergence (t-test: t(40) = 0.68, P = 0.49): couples where the female’s third harmonic converged with the male’s second had a size ratio of 0.954 ± 0.015; couples where the female’s fourth harmonic converged with the male’s third had a size ratio of 0.942 ± 0.007.
Male–female precopulatory acoustic behaviors during free flight
To characterize mosquito flight tones prior to copulation, we examined male–female acoustic behaviors during mating attempts. Male mating attempts were determined by visual and audio analysis: a male “chasing” a female [25, 40] was classified as the onset of a mating attempt. From 33 different females, 134 mating attempts were analyzed. Mating attempts resulted in either a copulation (i.e. the female was inseminated) or a mating rejection; 12 of 33 females mated. The duration of a mating interaction was brief (1.38 ± 0.32 s; n = 33). Ten of the mated females rejected a male at least once, while two mated after the first mating attempt. All mated females copulated within 1–5 min after introduction to the arena. Females that did not mate in our assays rejected males 1–9 times, with females that rejected once (n = 3) not flying the remainder of the assay. To compare mating attempts that ended in a successful copulation or a rejection, we examined audio recordings of the 12 interactions that resulted in a successful copulation and randomly selected one mating attempt from each of the 21 females that did not mate in order to identify male- or female-specific acoustic behaviors that occur immediately prior to An. albimanus copulation or mate rejection.
We consistently observed that all mating attempts, regardless of outcome, were characterized by an increase in male WBF followed by an increase in female WBF (Fig. 3a; Additional file 2: Video S1). The increase in the WBF of the male coincided with the onset of the female chase, while that of the female coincided with her flying away from the approaching male. After increasing their WBFs, males and females rapidly modulated their frequency until couples departed or stopped flying. Although these oscillations appear to be similar to those reported in other species [25, 26, 40], it was not always possible to observe this behavior in every mating attempt as 25% of the couples that copulated landed upon making contact. The time interval of rapid frequency oscillation after reaching maximum WBF ranged between 0.34–1.38 s and 0.18–1.56 s for males and females, respectively.
Using the criteria for tethered pairs in Aldersley et al. , we did not detect harmonic convergence in free-flying pairs. The brevity of the mating interaction, the WBF range of individuals in free flight, and that several couples ceased flying upon contact suggest that our methods were not sufficient to detect this phenomenon. However, we analyzed the latency of the response of the female, the female/male WBF ratio prior the interaction and the WBF increase during the mating attempt to determine if these characteristics influenced the outcome of a mating attempt. The latency of the females’ WBF increase in response to that of the male did not significantly differ between interactions that resulted in a copulation or a rejection (copulation: 0.375 ± 0.074 s; rejection: 0.542 ± 0.047 s; Mann Whitney U-test: U(12, 21) = 78, Z = − 1.79, P = 0.07). The female/male frequency ratio average (♀WBF/♂WBF) 1–1.5 s prior to a copulation, before any observable male–female interaction occurred, was similar for mating attempts that resulted in copulation (♀WBF/♂WBF: 0.595 ± 0.024) or rejection (♀WBF/♂WBF: 0.596 ± 0.024; Mann–Whitney U-test: U(12, 21) = 125, Z = − 0.03, P = 0.97) and did not influence mating outcome. Finally, to analyze the WBF increase during the mating attempt, we measured the change in WBF (ΔF) and the time interval between the initial WBF increase and the maximum frequency reached (Δt), allowing us to calculate the WBF rate of increase (ΔF/Δt) of both sexes during this interaction (Table 1, Fig. 3a). In males, we found that the rate of WBF increase was similar for mating attempts that resulted in a copulation or rejection (Mann–Whitney U-test: U(12, 21) = 109, Z = 0.63, P = 0.52) (Table 1, Fig. 3b). While females also increased their WBFs regardless of mating outcome, the rate of the frequency increase was significantly greater prior to a mating rejection compared to a copulation (Table 1, Fig. 3b) (Mann–Whitney U-test: U(12, 21) = 24, Z = − 3.81, P < 0.01), owing to differences in the extent of the frequency increase (Mann–Whitney U-test: U(12, 21) = 58, Z = − 2.54, P = 0.01) and the time interval of the increase (Mann–Whitney U-test: U(12, 21) = 40.5, Z = 3.19, P < 0.01) (Table 1).
As we detected differences in the female rate of WBF increase between mate rejection or copulation subsequent to a mating attempt, we further analyzed this behavior. In females that rejected a male more than once (n = 18), we performed a paired comparison between two different mating attempts of the same female and found no differences in the rate of WBF increase (Wilcoxon matched pairs test: Z = 0.97, P = 0.32) (Additional file 1: Figure S2a). Thus, females consistently increased their WBFs during rejections. However, in females that rejected and then copulated with a male (n = 10), we performed a paired comparison between the copulation and a rejection of the same female and found significant differences in the rate of WBF increase between outcomes (Wilcoxon matched pairs test: Z = 2.84, P < 0.01) (Additional file 1: Figure S2b): the rate of increase was lower for mating attempts that resulted in copulation compared with a rejection. Thus, in our assays, the increasing rate of the female WBF predicted copulation or rejection. However, there are likely additional behavioral and physiological factors that influence mating.
Male-male acoustic behaviors during free flight
We next analyzed acoustic and flying behaviors among groups of males. We introduced eight males into a soundproof chamber and made an audio record of their flight tones and a visual record of their flight characteristics. As only a portion of the introduced males flew at any given time, we have differing numbers of recordings of two, three and four males flying simultaneously. For our analysis, we selected 5 s intervals in which flying behavior was clear to analyze male flight trajectories using an ellipse , determining the distance of each individual in relation to the rest of the group and following the X and Y coordinates of each individual over this time frame (see “Methods”). To assess flight acoustics of male groups, we analyzed the average frequency difference of each group; the difference between the highest and lowest WBFs observed in the spectrogram represented the frequency bandwidth of the interaction among males during the 5 s segments analyzed.
Two distinct flying behaviors were apparent among groups of males flying simultaneously: random flight and a swarm-like, patterned flight (Fig. 4a, b; Additional file 3: Video S2), referred to hereafter as “patterned flight”. During random flight, male WBFs diverged or were distinct (Fig. 4c, d, bottom panels), but converged once patterned flight was initiated, reducing the frequency differences of the male flight tones (two males: Fig. 4c; four males: Fig. 4d; Additional file 3: Video S2). By examining flight behavior in 5 s intervals, we observed a significant change in the frequency difference between males in random flight and those in patterned flight. This change was consistent among groups of two (Mann–Whitney U-test: U(8, 8) = 0, Z = − 3.36, P < 0.01), three (Mann–Whitney U-test: U(5, 5) = 0, Z = − 2.61, P < 0.01) and four (Mann–Whitney U-test: U(5, 5) = 0, Z = − 2.61, P < 0.01). The difference in WBF of groups in patterned flight was always smaller than groups in random flight (Fig. 4e). The mean frequency difference of groups of two, three and four males in patterned flight was 38.20 ± 3.33, 52.60 ± 4.01 and 68.21 ± 7.11 Hz, respectively. On the other hand, the mean frequency difference of groups of two, three and four males in random flight was 110.88 ± 19.91, 177.05 ± 31.16 and 167.65 ± 19.88 Hz, respectively.
During patterned flight, we observed that An. albimanus fly in an ellipsoid pattern around a central position, similar to what has been reported in Cx. quinquefasciatus . Therefore, we quantified visual differences between random and patterned flight (where we observed acoustic interactions) by examining the trajectory of individual mosquitoes and analyzing the ellipsoid shape formed while flying . We found significant differences in the characteristics of male trajectories between the two types of flight observed: the ellipse area representing random flight is significantly larger than the area representing patterned flight [random flight: (1.55 ± 0.29) × 105 px2; patterned flight: (0.49 ± 0.06) × 105 px2; Mann–Whitney U-test: U(18, 18) = 58, Z = 3.24, P < 0.01]. We further measured the distance (in pixels) of each mosquito in relation to the rest of the group by calculating the average distance among the centroids of each ellipse in comparison with the others. We found that mosquitoes flying randomly are significantly less aggregated than males in patterned flight (random flight: 244.66 ± 21.03 px; patterned flight: 95.79 ± 11.15 px; Mann–Whitney U-test: U(18, 18) = 69, Z = 2.89, P < 0.01). Finally, by separately analyzing the X and Y trajectories of males in flight, we observed characteristic loops (i.e. periodic movement) similar to those defined in previous studies . Mosquitoes in patterned flight display loops with significantly shorter periods (random flight: 2.63 ± 0.39 s; patterned flight: 0.90 ± 0.08 s; Mann–Whitney U-test: U(18, 18) = 33.5, Z = 4.03, P < 0.01) (Fig. 4f). Taken together, these data show that patterned flight is associated with distinctive acoustic behaviors.
Mosquito flight tones differ by species and are used to develop novel trapping tools  and recognition algorithms , highlighting the importance of flight tone characterization of disease vectors. While there is ample information about Ae. aegypti, Cx. quinquefasciatus and An. gambiae acoustic behavior prior to and during mating, there is no information about acoustic behavior of one the most important malarial vectors in the Americas, An albimanus. Furthermore, most acoustic characterization of mosquito vectors has been performed using immobilized specimens. Here, we characterized acoustic behavior of An. albimanus. We observed significant differences in the flight acoustics between free-flying and tethered An. albimanus, describe acoustic interactions between free-flying and tethered male–female couples, and identify a relationship between male-male acoustic interactions and the coordination of movement that occurs when males transition from random flight to patterned flight.
Tethered An. albimanus exhibit lower WBFs than other reported mosquito species (Ae. aegypti [4, 48], An coluzzi , An. gambiae [21, 22]). While WBFs vary among mosquito species , they can be affected by factors that include ambient temperature [4, 50], humidity and age . We observed major acoustic differences between tethered and free-flying mosquitoes: (i) the mean WBF of free-flying An. albimanus is higher for both sexes; (ii) the mean WBF ratio between females and males (♀WBF/♂WBF) changes (tethered = 0.71, free flight = 0.61); and (iii) free-flying mosquitoes modulate their WBF to a greater extent than tethered individuals. Tethered insects have been used to study locomotion [52, 53], migration [52, 53] and wing movement ; these studies have demonstrated that tethering leads to unnatural flight behaviors and distorted wing strokes. Tethered individuals are restricted to a horizontal flight path, do not produce the equivalent lift required to achieve normal flight, and do not support their own body mass, factors that may give an inaccurate view of natural flight . As flight tones are linked to motor function , our results show that restricting motion alters two flight tone characteristics: the WBF and the ability to modulate it. Although a decrease in the WBF linked to tethering has been reported for Ae. aegypti , Cx. quinquefasciatus , midges  and locusts , Villareal et al.  found no difference in WBF of tethered and free-flying Ae. aegypti females when assays were performed at the same temperature. Although our assays were performed at the same temperature, we found that tethering impacted An. albimanus WBFs.
Harmonic convergence, or the synchronization of frequencies between two specific harmonics during male–female interactions, has been described in Ae. aegypti [17, 18], Cx. quinquefasciatus  and An. gambiae [20, 21] using tethered individuals. This behavior appears to be related to the male auditory strategy used to track females  and, more recently, has been suggested to be a manifestation of male–female coordination prior to formation of the copula position . We observed harmonic convergence in tethered An. albimanus. The time of convergence and the latency of the response were similar to An. gambiae  and the convergence events we detected have been observed in other species [18, 21, 37]. When examining male–female acoustic interactions during a mating attempt in free flight, we were unable to observe harmonic convergence using the criteria of Aldersley et al. . Harmonic convergence presumably occurs rapidly during a mating attempt, resembling the early stages of Ae. aegypti precopulatory behavior (i.e. from the male approach to the end of RFM) . It is possible that convergence similar to what we observed under tethered conditions occurs after male–female contact, as has been demonstrated in Ae. aegypti . Thus, we might not have detected convergence in our assays as An. albimanus couples often did not exhibit long paired flight after making contact.
However, we observed a common precopulatory behavior in free flight: a rapid increase in male WBF immediately followed by a rapid increase in female WBF. A similar behavior is seen in Cx. quinquefasciatus , An. coluzzii and An. gambiae , where males increase their WBF in response to female flight tone. The rate of frequency increase in males was similar regardless of mating outcome, while a more rapid increase in female frequency resulted in mate rejection. It is possible that the frequency increase results from controlled flight to reach or escape a potential mate. Results from the Culex  and Anopheles  auditory system suggest that to locate females, males use the difference between their own frequency and that of the female  rather than the female WBF itself. The rapid change in female WBF during mate rejection immediately modifies the male–female WBF difference. Remaining within the optimal auditory sensitivity range may be a critical characteristic to reach and copulate with a female.
Male swarming behavior appears to be an obligatory feature of copulation for some Anopheles species [35, 58]. However, swarming behaviors have not been described for Latin American species. Furthermore, male-male interactions during swarming, and female mate selection within a swarm, are not well understood. Acoustic signaling might be an important factor during these interactions, although studies of male-male and male–female interactions in free flight are few. In our assays, An. albimanus males performed two identifiable types of flying behavior regardless of the number of males assessed: random and patterned flight. During random flight, males flew in large trajectories within the entire experimental arena. Once patterned flight was initiated, the large trajectories immediately gave way to a form of flight composed of small loops in a specific region of the arena. Interestingly, males in patterned flight acoustically interacted within a narrower band of frequency, linking flight pattern and acoustic behavior. This result is likely related to the clustering of male flight tones of closely located, tethered Ae. aegypti males . During swarming, males must coordinate their flight patterns and recognize females that enter the swarm. It has been proposed that flight tones are used to coordinate movement during group flight , dividing into frequency clusters to reduce acoustic interference .
Linking acoustic and flying behaviors under natural, unrestricted conditions will provide important information regarding behaviors during courtship and ultimately, female mate selection. Many of the male–female acoustic interactions we observed during mating attempts in free flight have been reported in other mosquito species, suggesting that such interactions are common during courtship in mosquito species and supporting the idea that acoustic interactions are crucial for successful copulation. This study also reveals a connection between male acoustics and flight characteristics, giving insight into the importance of male-male acoustic interaction during patterned flight and potentially giving clues about swarm formation and cohesion. Furthermore, An. albimanus males modulate their flight tones in the presence of other males. While this behavior likely entails additional sensory cues, targeting acoustic male interactions may disrupt essential reproductive behaviors of this, and potentially other non-swarming species. Control programmes that rely on the mass release of laboratory-reared mosquitoes [60,61,62,63] need to consider acoustic behavior, since the results presented here, as well as those of previous reports, demonstrate that different mosquito species display similar precopulatory acoustic behaviors, highlighting the importance of acoustics in mating success [27, 28]. Characterization of An. albimanus acoustic flight behaviors, an important Latin American malaria vector, is a step toward understanding male-male and male–female interactions prior to copulation and will aid in the improvement of vector surveillance and control programmes in areas affected by this species.
Availability of data and materials
The datasets generated during and/or analyzed during the present study are available from the corresponding authors upon reasonable request.
wing beat frequency
rapid frequency modulation
Hay SI, Guerra CA, Tatem AJ, Noor AM, Snow RW. The global distribution and population at risk of malaria: past, present, and future. Lancet Infect Dis. 2004;4:327–36.
Recht J, Siqueira AM, Monteiro WM, Herrera SM, Herrera S, Lacerda MVG. Malaria in Brazil, Colombia, Peru and Venezuela: current challenges in malaria control and elimination. Malar J. 2017;16:273.
Padilla JC, Lizarazo FE, Murillo OL, Mendigana FA, Pachon E, Vera MJ. Transmission scenarios of major vector-borne diseases in Colombia, 1990–2016. Biomedica. 2017;37:27–40.
Arthur BJ, Emr KS, Wyttenbach RA, Hoy RR. Mosquito (Aedes aegypti) flight tones: frequency, harmonicity, spherical spreading, and phase relationships. J Acoust Soc Am. 2014;135:933–41.
Instituto Nacional de Salud. Boletín epidemiológico semanal. 2019. https://www.ins.gov.co/buscador-eventos/BoletinEpidemiologico/2019%20Bolet%C3%ADn%20epidemiol%C3%B3gico%20semana%204.pdf. Accessed 30 May 2019.
Blair S. Retos para la eliminación de la malaria en Colombia: un problema de saber o de poder. Biomedica. 2012;32:31–148.
Montoya-Lerma J, Solarte YA, Giraldo-Calderon GI, Quinones ML, Ruiz-Lopez F, Wilkerson RC, et al. Malaria vector species in Colombia: a review. Mem Inst Oswaldo Cruz. 2011;106(Suppl. 1):223–38.
Mitchell SN, Kakani EG, South A, Howell PI, Waterhouse RM, Catteruccia F. Mosquito biology. Evolution of sexual traits influencing vectorial capacity in anopheline mosquitoes. Science. 2015;347:985–8.
Kramer LD, Ciota AT. Dissecting vectorial capacity for mosquito-borne viruses. Curr Opin Virol. 2015;15:112–8.
Villarreal SM, Pitcher S, Helinski MEH, Johnson L, Wolfner MF, Harrington LC. Male contributions during mating increase female survival in the disease vector mosquito Aedes aegypti. J Insect Physiol. 2018;108:1–9.
Helinski MEH, Harrington LC. Male mating history and body size influence female fecundity and longevity of the dengue vector Aedes aegypti. J Med Entomol. 2011;48:202–11.
Childs LM, Cai FY, Kakani EG, Mitchell SN, Paton D, Gabrieli P, et al. Disrupting mosquito reproduction and parasite development for malaria control. PLoS Pathog. 2016;12:e1006060.
Mitchell SN, Catteruccia F. Anopheline reproductive biology: Impacts on vectorial capacity and potential avenues for malaria control. Cold Spring Harb Perspect Med. 2017;7:a025593.
Diabate A, Tripet F. Targeting male mosquito mating behaviour for malaria control. Parasit Vectors. 2015;8:347.
Alphey L, Benedict M, Bellini R, Clark GG, Dame DA, Service MW, et al. Sterile-insect methods for control of mosquito-borne diseases: an analysis. Vector Borne Zoonotic Dis. 2010;10:295–311.
Benelli G, Jeffries CL, Walker T. Biological control of mosquito vectors: past, present, and future. Insects. 2016;7:E52.
Cator LJ, Arthur BJ, Harrington LC, Hoy RR. Harmonic convergence in the love songs of the dengue vector mosquito. Science. 2009;323:1077–9.
Cator LJ, Arthur BJ, Ponlawat A, Harrington LC. Behavioral observations and sound recordings of free-flight mating swarms of Ae. aegypti (Diptera: Culicidae) in Thailand. J Med Entomol. 2011;48:941–6.
Warren B, Gibson G, Russell IJ. Sex Recognition through midflight mating duets in Culex mosquitoes is mediated by acoustic distortion. Curr Biol. 2009;19:485–91.
Pennetier C, Warren B, Dabiré KR, Russell IJ, Gibson G. “Singing on the wing” as a mechanism for species recognition in the malarial mosquito Anopheles gambiae. Curr Biol. 2010;20:131–6.
Cator LJ, Ng’Habi KR, Hoy RR, Harrington LC. Sizing up a mate: variation in production and response to acoustic signals in Anopheles gambiae. Behav Ecol. 2010;21:1033–9.
Gibson G, Warren B, Russell IJ. Humming in tune: sex and species recognition by mosquitoes on the wing. J Assoc Res Otolaryngol. 2010;11:527–40.
Belton P. Attraction of male mosquitoes to sound. J Am Mosq Control Assoc. 1994;10:297–301.
Wishart G, van Sickle GR, Riordan DF. Orientation of the males of Aedes aegypti (L.) (Diptera: Culicidae) to sound. Can Entomol. 1962;94:613–26.
Simoes PMV, Ingham RA, Gibson G, Russell IJ. A role for acoustic distortion in novel rapid frequency modulation behaviour in free-flying male mosquitoes. J Exp Biol. 2016;219:2039–47.
Aldersley A, Cator LJ. Female resistance and harmonic convergence influence male mating success in Aedes aegypti. Sci Rep. 2019;9:2145.
Cator LJ, Harrington LC. The harmonic convergence of fathers predicts the mating success of sons in Aedes aegypti. Anim Behav. 2011;82:627–33.
Cator LJ, Zanti Z. Size, sounds and sex: interactions between body size and harmonic convergence signals determine mating success in Aedes aegypti. Parasit Vectors. 2016;9:622.
Zawada JW, Dahan-Moss YL, Muleba M, Dabire RK, Maiga H, Venter N, et al. Molecular and physiological analysis of Anopheles funestus swarms in Nchelenge. Zambia. Malar J. 2018;17:49.
Diabate A, Yaro AS, Dao A, Diallo M, Huestis DL, Lehmann T. Spatial distribution and male mating success of Anopheles gambiae swarms. BMC Evol Biol. 2011;11:184.
Charlwood JD, Thompson R, Madsen H. Observations on the swarming and mating behaviour of Anopheles funestus from southern Mozambique. Malar J. 2003;2:2.
Hassan MM, Zain HM, Basheer MA, Elhaj HEF, Sayed El BB. Swarming and mating behavior of male Anopheles arabiensis Patton (Diptera: Culicidae) in an area of the Sterile Insect Technique Project in Dongola, northern Sudan. Acta Trop. 2014;132(Suppl):S64–9.
Manoukis NC, Diabate A, Abdoulaye A, Diallo M, Dao A, Yaro AS, et al. Structure and dynamics of male swarms of Anopheles gambiae. J Med Entomol. 2009;46:227–35.
Marchand RP. Field observations on swarming and mating in Anopheles gambiae mosquitoes in Tanzania. Neth J Zool. 1983;34:367–87.
Howell PI, Knols BGJ. Male mating biology. Malar J. 2009;8(Suppl. 2):S8.
Aldersley A, Champneys A, Homer M, Bode NWF, Robert D. Emergent acoustic order in arrays of mosquitoes. Curr Biol. 2017;27:R1208–10.
Aldersley A, Champneys A, Homer M, Robert D. Quantitative analysis of harmonic convergence in mosquito auditory interactions. J R Soc Interface. 2016;13:20151007.
Charlwood JD, Jones MDR. Mating behaviour in the mosquito, Anopheles gambiae s.l. save. I. Close range and contact behaviour. Physiol Entomol. 1979;4:111–20.
Gibson G, Russell I. Flying in tune: sexual recognition in mosquitoes. Curr Biol. 2006;16:1311–6.
Simoes PMV, Gibson G, Russell IJ. Pre-copula acoustic behaviour of males in the malarial mosquitoes Anopheles coluzzii and Anopheles gambiae s.s. does not contribute to reproductive isolation. J Exp Biol. 2017;220:379–85.
Heuvel MJ. The effect of rearing temperature on the wing length, thorax length and ovariole number of the adult mosquito, Aedes aegypti (L.). Trans R Entomol Soc Lond. 2009;115:197–216.
Swami A, Mendel JM, Nikias C. Higher-order spectral analysis toolbox: MATLAB user’s guide. Natick: The Math Works Inc.; 2001.
Gibson G. Swarming behaviour of the mosquito Culex pipiens quinquefasciatus: a quantitative analysis. Physiol Entomol. 1985;10:283–96.
de Silva P, Nutter B, Bernal XE. Use of acoustic signals in mating in an eavesdropping frog-biting midge. Anim Behav. 2015;103:45–51.
Kutsch W, Stevenson P. Time-correlated flights of juvenile and mature locusts: a comparison between free and tethered animals. J Insect Physiol. 1981;27:455–9.
Johnson BJ, Rohde BB, Zeak N, Staunton KM, Prachar T, Ritchie SA. A low-cost, battery-powered acoustic trap for surveilling male Aedes aegypti during rear-and-release operations. PLoS ONE. 2018;13:e0201709.
Mukundarajan H, Hol FJH, Castillo EA, Newby C, Prakash M. Using mobile phones as acoustic sensors for high-throughput mosquito surveillance. eLife. 2017;6:e27854.
Gopfert MC, Briegel H, Robert D. Mosquito hearing: sound-induced antennal vibrations in male and female Aedes aegypti. J Exp Biol. 1999;202:2727–38.
Kahn MC, Celestin W, Offenhauser W. Recording of sounds produced by certain disease-carrying mosquitoes. Science. 1945;101:335–6.
Villarreal SM, Winokur O, Harrington L. The impact of temperature and body size on fundamental flight tone variation in the mosquito vector Aedes aegypti (Diptera: Culicidae): implications for acoustic lures. J Med Entomol. 2017;54:1116–21.
Costello RA. Effects of environmental and physiological factors on the acoustic behavior of Aedes aegypti (L.) (Diptera: Culicidae). 1974. https://core.ac.uk/download/pdf/56367578.pdf. Accessed 11 July 2019.
Minter M, Pearson A, Lim KS, Wilson K, Chapman JW, Jones CM. The tethered flight technique as a tool for studying life-history strategies associated with migration in insects. Ecol Entomol. 2018;43:397–411.
Ribak G, Barkan S, Soroker V. The aerodynamics of flight in an insect flight-mill. PLoS ONE. 2017;12:e0186441.
Thomas ALR, Taylor GK, Srygley RB, Nudds RL, Bomphrey RJ. Dragonfly flight: free-flight and tethered flow visualizations reveal a diverse array of unsteady lift-generating mechanisms, controlled primarily via angle of attack. J Exp Biol. 2004;207:4299–323.
Dudley R. The biomechanics of insect flight: form, function, evolution. Ann Entomol Soc Am. 2014;93:1195–6.
Bomphrey RJ, Nakata T, Phillips N, Walker SM. Smart wing rotation and trailing-edge vortices enable high frequency mosquito flight. Nature. 2017;544:92–5.
Simoes PMV, Ingham R, Gibson G, Russell IJ. Masking of an auditory behaviour reveals how male mosquitoes use distortion to detect females. Proc Biol Sci. 2018;285:20171862.
Yuval B. Mating systems of blood-feeding flies. Annu Rev Entomol. 2006;51:413–40.
Lapshin DN, Vorontsov DD. Frequency tuning of individual auditory receptors in female mosquitoes (Diptera, Culicidae). J Insect Physiol. 2013;59:828–39.
Carvalho DO, McKemey AR, Garziera L, Lacroix R, Donnelly CA, Alphey L, et al. Suppression of a field population of Aedes aegypti in Brazil by sustained release of transgenic male mosquitoes. PLoS Negl Trop Dis. 2015;9:e0003864.
Flores HA, O’Neill SL. Controlling vector-borne diseases by releasing modified mosquitoes. Nat Rev Microbiol. 2018;16:508–18.
O’Neill SL. The use of Wolbachia by the World Mosquito Program to interrupt transmission of Aedes aegypti transmitted viruses. Adv Exp Med Biol. 2018;1062:355–60.
Servick K. Winged warriors. Science. 2016;354:164–7.
We would like to thank the Programa de Estudio y Control de Enfermedades Tropicales (PECET, Universidad de Antioquia) for generously providing the mosquitoes used in this study and technical support; SISTEMIC for computing support; Laura Harrington, Ethan Degner, Garrett League and Miguel Toro for helpful comments on the manuscript; and J. Francisco Vargas B. and Freddy Ruiz-Lopez for technical advice.
This study was supported by COLCIENCIAS Grants CD 325671250804 (CT 409-2016) (to CAP), 647 (to support HPS), COLCIENCIAS, Universidad de Antioquia and the Max Planck Society cooperation Grant 566-1 (2014) (to FWA) and Universidad de Antioquia Grant CPT1410 (to support VV).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The relationship between size and wingbeat frequency of a tethered and b free-flying mosquitoes recorded at 26 ± 2 °C. Linear regression between wing length and frequency of tethered males (blue, tethered: R2 < 0.01, Pearson’s r = 0.06, P = 0.58; free flight: R2 = 0.12, Pearson’s r = 0.35, P = 0.02) and females (red, tethered: R2 = 0.03, Pearson’s r = − 0.20, P = 0.05; free flight: R2 < 0.01, Pearson’s r = 0.03, P = 0.81). Although our rearing method did not intend to generate individuals of different sizes, the determination coefficient (R2) and Pearson’s correlation coefficient (r) shows that in males and females there is a weak or a non-significant relationship between size and wingbeat frequency. Figure S2. Paired comparison between a two rejections and b a rejection and a successful copulation of the same female. While there is no significant difference between rejections (Wilcoxon matched pairs test: Z = 0.97, P = 0.32), there is a significant decrease of the WBF rate of increase (Wilcoxon matched pairs test: Z = 2.84, P < 0.01) between rejections and successful copulations.
Video of male–female mating attempts: rejections and successful copulation. The left panel shows the experimental arena and the reconstruction of the flight trajectories of the male (in blue) and the female (in red). Right panels indicate the time variation of X (top) and Y (middle) coordinates of the flight trajectories and the spectrogram (bottom) of the fundamental frequency.
Comparison between random and patterned flight for groups of two and four mosquitoes. The left panel shows the experimental arena and the reconstruction of the flight trajectories of males. Right panels indicate the time variation of X (top) and Y (middle) coordinates of the flight trajectories and the spectrogram (bottom) of the fundamental frequency.
About this article
Cite this article
Pantoja-Sánchez, H., Gomez, S., Velez, V. et al. Precopulatory acoustic interactions of the New World malaria vector Anopheles albimanus (Diptera: Culicidae). Parasites Vectors 12, 386 (2019). https://doi.org/10.1186/s13071-019-3648-8
- Harmonic convergence