Mating Rhythm of T. castaneum
During ten days after emergence, daily mating rates of the same geographical population varied significantly within different ages of T. castaneum (one-way ANOVA or Kruskal-Wallis test: P < 0.05), and daily mating rates on the same day were not significantly different among six populations of beetles (one-way ANOVA or Kruskal-Wallis test: P > 0.05). Furtherly, daily mating rates varied significantly within different ages across six populations of beetles (two-way ANOVA: F = 89.99; df = 9, 240; P < 0.01), and no significant variations in the daily mating rhythm were found among six populations of beetles (two-way ANOVA: F = 0.32; df = 5, 240; P = 0.90). Six populations of T. castaneum exhibited an identical daily mating rhythm. Beetles started to mate once they emerged. The daily mating frequency kept increasing with the increment of age, and reached a peak on the fourth day (52.0–58.0%), and then declined to some extent, and maintained at a relatively high level till the end of the test (31.0–45.0%). It was also found that T. castaneum mated more than once in its lifetime and the average mating frequency of an adult was about 3.5 during ten days after emergence (Fig. 1).
As the mating rate of T. castaneum was very low during the first two days after emergence, bihourly mating frequencies of 3 to 10-day-old beetles were collected to analyze the diurnal mating rhythm to avoid the disturbance of random mating. As for 3 to 10-day-old beetles, they could mate at any time. Bihourly mating rates of the same population of beetles varied significantly within different mating intervals (one-way ANOVA or Kruskal-Wallis test: P < 0.05). Bihourly mating rates of the same mating interval were not significantly different among six populations of beetles (one-way ANOVA or Kruskal-Wallis test: P > 0.05). Furtherly, bihourly mating rates of beetles varied significantly within 12 mating intervals across six populations (two-way ANOVA: F = 10.18; df = 11, 288; P < 0.01), and no significant variations in the diurnal mating rhythm were found among six populations (two-way ANOVA: F = 0.48; df = 5, 288; P = 0.79). The mating frequency reached the maximum during 16:00–20:00 and the bihourly mating rates ranged 3.9–4.9% (Fig. 2).
Mating Preference of T. castaneum
In the same geographical population, the attraction rates of females to males ranged 48.7–52.7% during 7–10 days after emergence, which were similar to random choice (50%). The attraction rate of males to females ranged 51.3–85.3% (Fig. 3). Male attraction was significantly higher than the female attraction in the populations of HBWH, SHJD, and HNZZ (one-way ANOVA: P < 0.05), while it was not found in the populations of ZJJX, GA1, and GDST (one-way ANOVA: P > 0.05). Furtherly, male attraction was significantly higher than female attraction across six populations (two-way ANOVA: F = 94.76; df = 1, 24; P < 0.01). This suggests that in the process of orientating a mating partner, male T. castaneum was more like the pheromone sender and females searched for their sexual partners with the aid of the sexual pheromone. Furthermore, the attraction to the opposite sex varied significantly across six populations of beetles (two-way ANOVA: F = 11.79; df = 5, 24; P < 0.01). The attraction between two genders from the HBWH population was highest sequentially followed by beetles from the populations of SHJD, GA1, HNZZ, ZJJX, and GDST. This indicates that there was detectable divergence in the chemical communication among six populations of beetles.
When the female preference was measured, females from six populations all preferred males to the cracked wheat and the average attraction rates of males to females ranged 51.3–92.7%. Males from six populations showed significantly different attractiveness to the same population of females (one-way ANOVA: P < 0.05), and females from six populations showed similar preferences to the same population of males (one-way ANOVA: P > 0.05). Furtherly, male attraction was significantly different among six populations (two-way ANOVA: F = 176.62; df = 5, 72; P < 0.01) and males from the HBWH population were most attractive to females sequentially followed by males from the populations of SHJD, HNZZ, GA1, GDST, and ZJJX (Fig. 4). Conversely, the female preference was not significantly different among six populations (two-way ANOVA: F = 1.21; df = 5, 72; P = 0.31). This suggests that the population divergence in the mate choice could be attributed to changes in the male pheromone production rather than changes in the female pheromone perception.
As six populations of females exhibited a similar mate choice, females from the population of HBWH were used as the odor receiver in the following odor choice tests. As the attraction of males varied significantly within the populations of HBWH, HNZZ, and ZJJX, these three populations were used in the following male volatile analysis test.
Constitution Of Male Volatiles
Volatiles from male T. castaneum mainly consisted of methyl-1,4-benzoquinone, DMD, 1-tetradecene, and 1-pentadecene, all of which accounted for > 98% of the whole volatiles measured by the area normalization method (Table S3). So standard chemicals of these four volatiles were further applied in the odor choice test.
Among all experimental designs, only DMD exhibited olfactory attractiveness to females at a moderate release rate and methyl-1,4-benzoquinone even exhibited repellency to females at the highest release rate (χ2 test, P < 0.05) (Fig. 5). The attraction of four male volatiles to females was significantly different and DMD showed the strongest attractiveness to females (two-way ANOVA: F = 419.01; df = 3, 40; P < 0.01). The release rate of volatiles significantly influenced the choice of females (two-way ANOVA: F = 156.93; df = 4, 40; P < 0.01) and the moderate release rate (0.1 and 1 mM) of volatiles could attract more females.
Furthermore, the DMD production by males was significantly different among six populations during 7–10 days after emergence (one-way ANOVA: F = 46.18; df = 5, 24; P < 0.01). The DMD production by individual males from the population of HBWH was highest sequentially followed by males from the populations of SHJD, HNZZ, GA1, GDST, and ZJJX (Fig. 6), which was positively correlated with male attraction in the mate choice test. This corroborated that variations in male DMD production played an important role in the population divergence of the mate choice.
Effect Of Gut Bacteria On Host Pheromone Communication
Six culturable gut bacteria were isolated from male T. castaneum from the geographical population of HBWH. They were identified as Bacillus cereus, C. freundii, Enterococcus faecium, E. coli, Klebsiella pneumoniae, and Staphylococcus aureus through the sequence analysis of the 16S rDNA gene (Table S4). The quantitative analysis indicated that the load of gut bacteria community in T. castaneum males were effectively removed by the antibiotic treatment and the abundance of certain gut bacteria were restored to some extent after reinoculation of a certain bacterial strain to axenic beetles (Table S5).
The antibiotic treatment significantly decreased the olfactory attraction of males to females and gnotobiotic reinoculation of gut bacteria had significant effects on male attraction (one-way ANOVA: F = 54.36; df = 7, 16; P < 0.01) (Fig. 7A). Two bacterial strains, C. freundii and E. coli significantly enhanced the attraction of males to females, while other four bacterial strains could not. Furthermore, the antibiotic treatment significantly decreased the male pheromone production, and gnotobiotic reinoculation with some bacterial strains significantly rescued the male DMD production (one-way ANOVA: F = 731.46; df = 7, 32; P < 0.01) (Fig. 7B). In particular, the male DMD production was almost restored to the normal level after gnotobiotic reinoculation with C. freundii and E. coli, which was consistent with their performance in male attraction. This suggests that two bacterial strains could mediate the chemical communication in T. castaneum by interfering with the male pheromone production. Although E. faecium also could enhance the male DMD production, no difference in male attraction was found between the gnotobiotic group and the axenic group, which suggests that females did not respond to this magnitude of difference in DMD production.