Next Article in Journal
Maternal Exposure to Ozone Modulates the Endophyte-Conferred Resistance to Aphids in Lolium multiflorum Plants
Previous Article in Journal
WingMesh: A Matlab-Based Application for Finite Element Modeling of Insect Wings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sex-Specific Wolbachia Infection Patterns in Populations of Polygraphus proximus Blandford (Coleoptera; Curculionidae: Scolytinae)

1
Laboratory of Molecular Genetics of Insects, Institute of Cytology and Genetics SB RAS, 630090 Novosibirsk, Russia
2
Laboratory of Insect Pathology, Institute of Systematics and Ecology of Animals SB RAS, 630091 Novosibirsk, Russia
3
Institute of Living Systems, Immanuel Kant Baltic Federal University, 236041 Kaliningrad, Russia
*
Author to whom correspondence should be addressed.
Insects 2020, 11(8), 547; https://doi.org/10.3390/insects11080547
Submission received: 3 August 2020 / Revised: 9 August 2020 / Accepted: 17 August 2020 / Published: 18 August 2020

Abstract

:

Simple Summary

Wolbachia bacteria are the most common symbionts of insects. These bacteria are ordinarily transmitted via oocyte cytoplasm from mother to progeny, like mitochondria, and are sporadically transmitted from one species to another. The Wolbachia symbionts have evolved to be parasitic (feminization of genetic males, male-killing, parthenogenesis, and cytoplasmic incompatibility) or/and mutualistic (increasing lifespan and fecundity, providing vitamins and nutrients, defending against viruses and parasites). Here we have studied Wolbachia infection in populations of four-eyed fir bark beetle Polygraphus proximus, which is one of the most dangerous pests of Siberian fir forests. A high rate of the only wProx Wolbachia strain in P. proximus populations was found in a vastly studied territory. Surprisingly, females were more often harboring Wolbachia than males. Besides, a comparison of the Wolbachia density in individuals has revealed that females contain much more Wolbachia symbionts than males. We suppose that the difference in infection status, as well as the difference in Wolbachia load between males and females within a population, can be found in some other Wolbachia–host associations.

Abstract

Wolbachia symbionts are maternally inherited bacteria that are widely distributed among Arthropoda hosts. Wolbachia influence their host biology in diverse ways. They may induce reproductive abnormalities, protect hosts against pathogens and parasites, or benefit hosts through metabolic provisioning. The progeny of an infected female are ordinarily infected with Wolbachia; however, Wolbachia have no future in male host progeny because they cannot transmit the symbiont to the next generation. Here, we analyze native and invasive populations of the four-eyed fir bark beetle (Polygraphus proximus) for Wolbachia prevalence and symbiont genetic diversity. This species is a dangerous pest of Siberian fir (Abies sibirica) forests. The native range of P. proximus includes the territories of the Russian Far East, Japan, Korea, and Northeast China, whereas its invasive range includes West Siberia, with further expansion westward. Surprisingly, we revealed a difference in the patterns of Wolbachia prevalence for males and females. Infection rate and Wolbachia titers were higher in females than in males. ST-533, the only haplotype of Wolbachia supergroup B, was associated with a minimum of three out of the five described mitochondrial haplotypes.

1. Introduction

Bacteria of the Wolbachia genus are the most successfully inherited symbionts that inhabit numerous arthropod and nematode hosts [1,2,3]. Wolbachia are vertically transmitted through host generations via oocytes [4], and are sporadically transmitted, horizontally, to non-related hosts [5,6,7,8,9]. These symbionts may affect their host biology in different ways, from manipulating host reproduction to increasing host fitness. Thus, Wolbachia increase lifespan and fecundity [10,11,12], protect against viruses and parasites [13,14,15,16], provide vitamins and other nutrients [17,18], and suppress some mutations in their hosts [19,20,21]. They also influence host reproduction by inducing feminization of genetic males, male killing, parthenogenesis, or cytoplasmic incompatibility (CI) [22]. These phenomena promote the spread of infection in a host population. The first three phenomena lead to a shift in the sex ratio towards females in a host population [22], whereas CI is manifested by embryonic mortality in the progeny of symbiont-free females and Wolbachia-harboring males. Therefore, Wolbachia in males commonly have no future and, in the case of CI, only function in controlling the spread of the symbiont in a host population.
We studied Wolbachia infection in native and invasive populations of the four-eyed fir bark beetle Polygraphus proximus Blandford, 1894 (Coleoptera; Curculionidae: Scolytinae). This species is one of the most dangerous pests of Siberian fir (Abies sibirica) forests. Since 2014, P. proximus has been listed in the European and Mediterranean Plant Protection Organization (EPPO) Alert List (https://gd.eppo.int/taxon/POLGPR). The native range of P. proximus includes the territories of the Russian Far East, Japan, Korea, and Northeast China [23,24,25]. The secondary (invasive) range covers the territories of West Siberia [24,25], with further expansion occurring westward (Figure 1).
Wolbachia infections have been reported in Scolytinae species, particularly in the following genera: Coccotrypes, Euwallacea, Hypothenemus, Ips, Pityogenes, Taphrorychus, Xyleborus, and Xylosandrus [26,27,28,29,30]. The infection rates vary, but are above 70% in most cases. The majority of Wolbachia isolates belong to supergroup A, with a few belonging to supergroup B.
Here, we first report on Wolbachia in the four-eyed fir bark beetle and determine infection rates of the population, including the infection prevalence according to sex. To do so, we carried out Wolbachia titer by qPCR and genotyped Wolbachia isolates using a multilocus sequence typing (MLST) protocol [31]. In addition, the barcode region of mitochondrial DNA was sequenced for infected and uninfected individuals collected from different localities, because Wolbachia symbionts and mitochondria are inherited together, according to the host’s maternal lineage, and can therefore form cytotypes. A cytotype results from the mitotype and the genotype of an infection, meaning cytotypic analysis can reveal the evolutionary history of Wolbachia infection [6,32,33,34,35].

2. Materials and Methods

2.1. Insect Collection

Adult specimens of P. proximus were collected in four regions of Russia: Sakhalin province (native range) (N = 60), Tomsk province (N = 271), Udmurtia (N = 38), and Krasnoyarsk Krai (N = 48) (invasive territories) (Figure 1 and Table 1). In Tomsk province and Udmurtia, the beetles were collected in their parental nests, whereas in Krasnoyarsk they were collected from forest products. In Sakhalin, overwintered beetles and those intruding into the bark of a host tree were collected. The samples were stored in air-dry conditions or fixed in 96% ethanol. The primary study was designed to analyze the infection status of each individual regardless of beetle sex. However, in the course of the study, we had to change the design to specifically include sexes as defined according to forehead morphology [36]. Finally, 7 out of 13 samples were subdivided by sex, namely one from Sakhalin province (beetles forming pairs), one from Udmurtia, two from Tomsk province, and three from Krasnoyarsk krai.

2.2. DNA Extraction

The total DNA was individually extracted from whole specimens. The insects were homogenized in 200 µL of DB extraction buffer (10 mM Tris-HCl (pH 8.0), 25 mM EDTA, 0.5% SDS, 0.1 M NaCl) with 1 µL of proteinase K (AppliChem) and incubated at 56 °C for 1 h. Total DNA was purified with NaOAc, precipitated, and dissolved in 100 µL of double-distilled water.

2.3. PCR and Sequencing

All DNA samples were tested by PCR with HCO/LCO primers for the cytochrome c oxidase subunit I (COI) gene [37]. Wolbachia infection status was determined by conventional PCR using the W-SpecF/R primer set [38] and coxAf1/r1 [31], and additionally confirmed for random samples with primer sets hcpAf1/r1 and fbpAf1/r1 [31]. Wolbachia-negative samples and controversial samples were checked by nested PCR with ftsZunif [39] for the first round and ftsZf1/r1 for the second, and/or fbpAf2/r2 and fbpAf1/r1, respectively. The PCR reactions were performed using BioMaster HS-Taq PCR (2×) (BiolabMix, Novosibirsk, Russia) with a 20 µL volume. The thermocycler protocol was 95 °C for 5 min, followed by 35 cycles of conventional PCR, and 15 + 30 cycles of nested PCR at 95 °C for 15 s, 55 °C for 40 s (1 min for HCO/LCO), and 72 °C for 40 s, plus a final elongation step of 72 °C for 3 min. To genotype a Wolbachia isolate, we amplified five MLST genes [31]. The amplicons (COI and MLST genes) were purified with exonuclease (ExoI) enzyme (New England Biolabs), and sequenced using a BrightDye Terminator Cycle Sequencing Kit (Nimagen). Full MLST profiles of four Wolbachia isolates were deposited in the PubMLST (https://pubmlst.org/wolbachia/) database, and all MLST alleles and fragments of the COI gene were deposited in the GenBank database (Table S1).

2.4. Quantitative Real-Time PCR Assay

We designed a qPCR assay to measure the relative number of Wolbachia genomes per beetle genome. A single copy gene per Wolbachia genome fbpA, and single copy gene per beetle genome Inx2 were used to estimate the Wolbachia load in DNA samples. For a nuclear reference, the 750 bp region of the Inx2 gene was amplified by nested PCR according to Che et al. [40], and then sequenced for confirmation of the target product (deposited in GenBank MT786531). Primers for qPCR were designed using the NCBI primer-BLAST tool (https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?LINK_LOC=BlastHome). The primer sets Inx2prox (F 5′-AGTGGGTGTGCTTCGTATTGT-3′, R 5′-GTCCGCCTTACAGTCCTCG-3′), and fbpAprox (F 5′-CTGCGTGGCTGTTGGGTTTA-3′, R 5′-GCGCAGCATAGGCAATAACA-3′) amplified 143 bp and 188 bp fragments of Inx2 and fbpA loci, respectively. Each qPCR mixture of 20 µL volume included 10 µL of ×2 qPCR HS-Taq mix (Biolabmix, Novosibirsk, Russia), 0.3 µM of each primer, and a 1 µL DNA sample that was in the range of 10–40 ng. All qPCRs were performed in duplicate on a CFX 96 Real-Time PCR Detection System (Bio-Rad Laboratories Inc., Hercules, CA, USA). Conditions of the qPCR were optimized by gradient PCR; we used 95 °C for 5 min, followed by 50 cycles of 95 °C for 10 s, and 63 °C for 30 s. Melt curves were generated to ensure specificity of the PCR product (65–95 °C with a step of 0.5 °C). Efficiencies of amplification and R2 were determined using a 5-fold dilution of a DNA sample for both loci (Inx2prox; E 98.21%, R2 0.9987; fbpAprox: E 101.50%, R2 0.9989).

2.5. Data Analysis

All statistical calculations were performed using Minitab 17.1.0 software (Minitab Inc., State College, PA, USA). The sequences were analyzed using FinchTV v. 1.4.0 (Geospiza Inc., Seattle, WA, USA). The alignments were generated using the MUSCLE algorithm [41]. For phylogenetic analysis of P. proximus mtDNA, the 620 bp region of the COI gene was used (38 sequences). To identify mtDNA haplogroups, we used the data from Kononov et al. [42]. Phylogenetic trees were reconstructed by the maximum likelihood algorithm in MEGA 6.0 [43], using the Tamura 3-parameter model for the COI gene and the GTR model for MLST data. For the MLST data analysis, we also used four Wolbachia isolates of Scolytinae species (Seqence Types (STs): 138, 207, 208, 210), seven STs (1, 8, 9, 35, 62, 19, 90) as Wolbachia supergroup references, and three STs (184, 303, 482) that had similar haplotype characteristics to ST-533 (Table S2; SM1). Each allele of ST-533 was independently analyzed for supergroup clustering to check the inter-supergroup recombination of Wolbachia haplotypes [8,9]. Wolbachia load was measured as a relative copy number (fbpA/Inx2) = 2 Cq(Inx2)−Cq(fbpA). Because the efficiencies of both primer sets were close to 100%, the amplification differences were not taken into account in calculations.

3. Results

3.1. Wolbachia Infection Rates and Wolbachia Titer

Wolbachia prevalence in populations of P. proximus was initially estimated by conventional PCR, followed by symbiont negative DNA samples being rechecked by nested PCR (Table S3). Finally, we estimated Wolbachia titer by real time quantitative PCR across local collections.

3.1.1. Conventional PCR

The Wolbachia symbionts were found in 12 out of 13 P. proximus samples (Table 2). The only Wolbachia-free sample was Tomarinsky-2, Sakhalin, with a sample size of N = 8. The infection rates varied among sampling sites from 39% to 77% and were not statistically different between the studied regions (Pearson’s chi-square, p = 0.237).
Surprisingly, we found a great difference in Wolbachia prevalence for males and females, with 6 out of 13 samples being partially or totally subdivided based on sex. The proportion of infected females was significantly higher than males: across all collections p < 0.001; in localities, for three populations p < 0.001; for two populations p < 0.05; and for one p = 0.138.

3.1.2. Nested PCR

Because a significant difference in Wolbachia rates within populations based on the sex of individuals was observed, we tried to reexamine all Wolbachia negative DNA samples by nested PCR. In total, we analyzed 189 samples, some of which were checked for two loci. As a result, 69 new infected samples were found. Previously seemingly Wolbachia-free samples from Tomarinsky-2 appeared to be nearly totally infected (7 out of 8). Taking into account the results based on both conventional and nested PCR, the infection rates varied from 48% to 92%. The infection rates in Udmurtia and Sakhalin were higher than in Tomsk and Krasnoyarsk (Fisher exact test for each of the four comparisons, p < 0.05 and p < 0.001).
Overall differences between infected females and males remained significant (Fisher exact test p < 0.01); however, in terms of localities, significant differences were only found in Tomarinsky-1 (p < 0.05) and Bakcharsky (p < 0.01) (Table 2).

3.1.3. Estimation of Wolbachia Titer by qPCR

The values of Wolbachia genomes per beetle genome were measured by real-time quantitative PCR. Here, we used 27 conventional-positive and 28 nested-positive samples across all regions; among them were 13 females, 20 males, and 22 samples with undetermined sex. Variation of the Wolbachia copy number (fbpA/Inx2) was in the range of 2.9 × 10−4 to 1.73. Our results for the conventional and nested PCR were in a good agreement with the estimation of Wolbachia load by qPCR. The threshold of the two PCR approaches was about 0.018–0.020 (Figure 2). The differences between the most heavily and lowest infected females was 22-fold, and the difference was 864-fold in males. The median of the Wolbachia titer was 0.1792 in females and 0.0035 in males, which was a significant difference (U Mann–Whitney p < 0.0001). In particular, there were no females below the threshold, whereas there were seven males (detected by conventional PCR) above it. In two cases, Wolbachia infection was successfully detected by conventional PCR, although the infection titers were ~0.002 and 0.005.

3.2. Characterization of Wolbachia Isolates and Mitochondrial DNA

Complete MLST profiles were obtained for four Wolbachia isolates (SM-1) from the Tomsk, Krasnoyarsk and Sakhalin regions; all were identical and corresponded to the ST-533 haplotype. Alleles ftsZ-263 and fbpA-461 were first characterized in the PubMLST database; then, we checked them in the GenBank database using the BLASTn tool and confirmed their uniqueness. For this reason, we characterized additional isolates from populations of Udmurtia, Tomsk, and Krasnoyarsk (IDs 1940–1943; SM-1) with regard to the ftsZ and fbpA genes and did not detect any other alleles. The concatenated sequence of ST-533 and each MLST allele was clustered into supergroup B (Figure 3 and Table S2). The ST-482 haplotype was the closest to ST-533, and there were three substitutions in the ftsZ and fbpA loci. Unfortunately, the host of ST-482 was not indicated in the PubMLST database (an isolate was not registered). The Wolbachia haplotype of Pityogenes chalcographus (Scolytinae) was also clustered into supergroup B, with a total of 72 substitutions from ST-533. Wolbachia infection of Leptopilina clavipes (parasitoid wasp) had the same hcpA-148 as ST-533 and ST-482; however, its haplotype (ST-303) was genetically closer to the isolate of P. chalcographus than that of P. proximus. The ftsZ-91 allele of Wolbachia from Hylyphantes graminicola (spider) differed from ftsZ-263 of wProx by only two substitutions, although other alleles were found to have much greater differences.
As Wolbachia infection was uniform across a vast territory, we analyzed a 620 bp barcode region of the COI gene to characterize maternal inheritance in more detail. In particular, we aimed to determine the mtDNA haplotype of infected and uninfected P. proximus specimens at all sampling sites. As references, we used data from Kononov et al. [42], where five mtDNA haplogroups have been described. The analyzed mtDNA of infected and uninfected samples of the Sakhalin and Tomsk collections belonged to haplogroup I (Figure 4). In Udmurtia they belonged to haplogroups I and II, and in Krasnoyarsk haplogroup IV was detected. Thus, we did not find dependence between the Wolbachia infection status and mtDNA haplogroup.

4. Discussion

Wolbachia symbionts in the subfamily Scolytinae have previously been found in 11 species of 8 genera [26,27,28,29,30], and many of these species are characterized by high infection rates (70–100%). Here, we first provided data on Wolbachia infection in populations of P. proximus in four distant regions of its native and invasive ranges. The overall infection rate in the studied populations was 68% (95% confidence interval, 63–73%). Surprisingly, the infection rate of females was higher than in males (79.7% vs. 61.8%; Fisher exact test p = 0.003). Moreover, most of the infected males were characterized by significantly lower Wolbachia titers than the females. Genetic diversity of Wolbachia in P. proximus populations indicated only one strain in the vast studied territory, which is associated with at least three out of five mitochondrial haplogroups.
Maternal transmission of Wolbachia implies that the symbiont has no future in males; therefore, it seemed logical that Wolbachia could be absent in males. We have not found other reports of sex-specific Wolbachia distribution in other populations, as is the case with single infected P. proximus. However, Arai et al. [44] reported, in an experiment on subdividing triple infection, a 79% transmission rate for the wHm-b strain in males of Homona magnanima (Lepidoptera: Tortricidae) versus a rate of 100% in daughters. Dutton and Sinkins [45], in double infected Aedes albopictus, found wAlbA loss in males. Further, Tortosa et al. [46] confirmed a loss of the wAlbA strain and demonstrated that wild females were double infected, whereas significant proportions of wild males were only single infected. There is a prediction of reduced expression of cytoplasmic incompatibility (CI) in males through selection of the host genome [47], and Wolbachia loss in males could be considered [46]. Therefore, if Wolbachia maintenance in a host population is based on CI, CI repression due to Wolbachia loss in males should lead to a loss of infection in a whole population. We do not have sufficient information to infer the effect(s) of Wolbachia on P. proximus biology, although the Wolbachia-inducing phenomena that lead to a sex ratio bias should be excluded, because we observed the primary sex ratio to be 1:1. There is a question as to whether Wolbachia in P. proximus is driven by CI and therefore will be lost in future, or high Wolbachia rates maintained by mutualistic effects for females, as, for instance, has been shown for Wolbachia-infected females of the coffee berry borer Hypothenemus hampei (Scolytinae), which produce more eggs than tetracycline-treated females [12].
The uninfected P. proximus males are likely the result of Wolbachia loss in ontogenesis. Most infected males are characterized by low Wolbachia density over all of the studied regions. Similar results were found for the wAlbA strain of Ae. albopictus [45,46]. In Scolytine, cases of low Wolbachia titers were found in populations of Pityogenes chalcographus [29]; however, there was no comparison between males and females. Therefore, there are few examples of sex-specific differences in Wolbachia titers, whereas there are a number of low symbiont titers in other species [29,48,49,50,51,52,53,54,55].

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/11/8/547/s1. Table S1. GenBank accession numbers of Wolbachia MLST genes and mtDNA COI gene used for phylogenetic analysis. Table S2. MLST profiles of Wolbachia isolates. Table S3. Data of conventional, nested PCR, and Wolbachia titer. Supplementary Material 1. Fasta file of concatenated sequences of five Wolbachia MLST genes used in this study.

Author Contributions

Y.I., R.B. and I.K. conceived and designed the research; I.K. collected P. proximus specimens and conducted morphology analysis; R.B., A.R., Y.I. and M.D. conducted the molecular genetic analysis; R.B. and Y.I. wrote the draft and created the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the State Budgeted Project No. 0324-2019-0041, the Russian Foundation for Basic Research (RFBR) No. 18-316-00099, No 18-34-20060, and No. 19-04-00983.

Acknowledgments

The authors express sincere gratitude to E. N. Akulov for the collection, and for kindly providing insects from Krasnoyarsk Krai; A. Fedorov and N. M. Debkov for samples from the Republic of Udmurtia; and to three anonymous reviewers for valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Taylor, M.J.; Bandi, C.; Hoerauf, A. Wolbachia. Bacterial endosymbionts of filarial nematodes. Adv. Parasitol. 2005, 60, 245–284. [Google Scholar] [CrossRef] [PubMed]
  2. Haegeman, A.; Vanholme, B.; Jacob, J.; Vandekerckhove, T.T.; Claeys, M.; Borgonie, G.; Gheysen, G. An endosymbiotic bacterium in a plant-parasitic nematode: Member of a new Wolbachia supergroup. Int. J. Parasitol. 2009, 39, 1045–1054. [Google Scholar] [CrossRef] [PubMed]
  3. Zug, R.; Hammerstein, P. Still a host of hosts for Wolbachia: Analysis of recent data suggests that 40% of terrestrial arthropod species are infected. PLoS ONE 2012, 7, e38544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Werren, J.H. Biology of Wolbachia. Annu. Rev. Entomol. 1997, 42, 587–609. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Vavre, F.; Fleury, F.; Lepetit, D.; Fouillet, P.; Boulétreau, M. Phylogenetic evidence for horizontal transmission of Wolbachia in host-parasitoid associations. Mol. Biol. Evol. 1999, 16, 1711–1723. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Baldo, L.; Ayoub, N.A.; Hayashi, C.Y.; Russell, J.A.; Stahlhut, J.K.; Werren, J.H. Insight into the routes of Wolbachia invasion: High levels of horizontal transfer in the spider genus Agelenopsis revealed by Wolbachia strain and mitochondrial DNA diversity. Mol. Ecol. 2008, 17, 557–569. [Google Scholar] [CrossRef]
  7. Ahmed, M.Z.; Breinholt, J.W.; Kawahara, A.Y. Evidence for common horizontal transmission of Wolbachia among butterflies and moths. BMC Evol. Biol. 2016, 16, 118. [Google Scholar] [CrossRef] [Green Version]
  8. Ilinsky, Y.; Kosterin, O.E. Molecular diversity of Wolbachia in Lepidoptera: Prevalent allelic content and high recombination of MLST genes. Mol. Phylogenet. Evol. 2017, 109, 164–179. [Google Scholar] [CrossRef]
  9. Shaikevich, E.; Bogacheva, A.; Rakova, V.; Ganushkina, L.; Ilinsky, Y. Wolbachia symbionts in mosquitoes: Intra-and intersupergroup recombinations, horizontal transmission and evolution. Mol. Phylogenet. Evol. 2019, 134, 24–34. [Google Scholar] [CrossRef]
  10. De Barro, P.J.; Hart, P.J. Antibiotic curing of parthenogenesis in Eretmocerus mundus (Australian parthenogenic form). Entomol. Exp. Appl. 2001, 99, 225–230. [Google Scholar] [CrossRef]
  11. Dong, P.; Wang, J.J.; Hu, F.; Jia, F.X. Influence of Wolbachia infection on the fitness of the stored-product pest Liposcelis tricolor (Psocoptera: Liposcelididae). J. Econ. Entomol. 2007, 100, 1476–1481. [Google Scholar] [CrossRef]
  12. Mariño, Y.; Verle Rodrigues, J.; Bayman, P. Wolbachia affects reproduction and population dynamics of the coffee berry borer (Hypothenemus hampei): Implications for biological control. Insects 2017, 8, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Hedges, L.M.; Brownlie, J.C.; O'Neill, S.L.; Johnson, K.N. Wolbachia and virus protection in insects. Science 2008, 322, 702. [Google Scholar] [CrossRef] [PubMed]
  14. Teixeira, L.; Ferreira, Á.; Ashburner, M. The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster. PLoS Biol. 2008, 6, e1000002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Cattel, J.; Martinez, J.; Jiggins, F.; Mouton, L.; Gibert, P. Wolbachia-mediated protection against viruses in the invasive pest Drosophila suzukii. Insect Mol. Biol. 2016, 25, 595–603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Van Nouhuys, S.; Kohonen, M.; Duplouy, A. Wolbachia increases the susceptibility of a parasitoid wasp to hyperparasitism. J. Exp. Biol. 2016, 219, 2984–2990. [Google Scholar] [CrossRef] [Green Version]
  17. Hosokawa, T.; Koga, R.; Kikuchi, Y.; Meng, X.Y.; Fukatsu, T. Wolbachia as a bacteriocyte-associated nutritional mutualist. Proc. Natl. Acad. Sci. USA 2010, 107, 769–774. [Google Scholar] [CrossRef] [Green Version]
  18. Nikoh, N.; Hosokawa, T.; Moriyama, M.; Oshima, K.; Hattori, M.; Fukatsu, T. Evolutionary origin of insect–Wolbachia nutritional mutualism. Proc. Natl. Acad. Sci. USA 2014, 111, 10257–10262. [Google Scholar] [CrossRef] [Green Version]
  19. Starr, D.J.; Cline, T.W. A host–parasite interaction rescues Drosophila oogenesis defects. Nature 2002, 418, 76. [Google Scholar] [CrossRef]
  20. Clark, M.E.; Anderson, C.L.; Cande, J.; Karr, T.L. Widespread prevalence of Wolbachia in laboratory stocks and the implications for Drosophila research. Genetics 2005, 170, 1667–1675. [Google Scholar] [CrossRef] [Green Version]
  21. Ote, M.; Ueyama, M.; Yamamoto, D. Wolbachia protein TomO targets nanos mRNA and restores germ stem cells in Drosophila sex-lethal mutants. Curr. Biol. 2016, 26, 2223–2232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Werren, J.H.; Baldo, L.; Clark, M.E. Wolbachia: Master manipulators of invertebrate biology. Nat. Rev. Microbiol. 2008, 6, 741. [Google Scholar] [CrossRef] [PubMed]
  23. Tokuda, M.; Shoubu, M.; Yamaguchi, D.; Yukawa, J. Defoliation and dieback of Abies firma (Pinaceae) trees caused by Parendaeus abietinus (Coleoptera: Curculionidae) and Polygraphus proximus (Coleoptera: Scolytidae) on Mount Unzen, Japan. Appl. Entomol. Zool. 2008, 43, 1–10. [Google Scholar] [CrossRef] [Green Version]
  24. Baranchikov, Y.; Akulov, E.; Astapenko, S. Bark beetle Polygraphus proximus: A new aggressive Far Eastern invader on Abies species in Siberia and European Russia. In Proceedings of the 21st US Department of Agriculture Interagency Research Forum on Invasive Species 2010, Annapolis, MD, USA, 12–15 January 2010; pp. 64–65. [Google Scholar]
  25. Krivets, S.A.; Bisirova, E.M.; Kerchev, I.A.; Pats, E.N.; Chernova, N.A. Transformation of taiga ecosystems in the Western Siberian invasion focus of four-eyed fir bark beetle Polygraphus proximus Blandford (Coleoptera: Curculionidae, Scolytinae). Russ. J. Biol. Inv. 2015, 6, 94–108. [Google Scholar] [CrossRef]
  26. Stauffer, C.; Van Meer, M.M.M.; Riegler, M. The presence of the protobacteria Wolbachia in European Ips typographus (Col., Scolytidae) populations and the consequences for genetic data. Proc. Ger. Soc. Gen. Appl. Entomol. 1997, 11, 709–711. [Google Scholar]
  27. Vega, F.E.; Benavides, P.; Stuart, J.A.; O’Neill, S.L. Wolbachia infection in the coffee berry borer (Coleoptera: Scolytidae). Ann. Entomol. Soc. Am. 2002, 95, 374–378. [Google Scholar] [CrossRef] [Green Version]
  28. Zchori-Fein, E.; Borad, C.; Harari, A.R. Oogenesis in the date stone beetle, Coccotrypes dactyliperda, depends on symbiotic bacteria. Physiol. Entomol. 2006, 31, 164–169. [Google Scholar] [CrossRef]
  29. Arthofer, W.; Markus, R.; Dimitrios, N.A.; Christian, S. Evidence for low-titre infections in insect symbiosis: Wolbachia in the bark beetle Pityogenes chalcographus (Coleoptera, Scolytinae). Environ. Microbiol. 2009, 11, 1923–1933. [Google Scholar] [CrossRef]
  30. Kawasaki, Y.; Schuler, H.; Stauffer, C.; Lakatos, F.; Kajimura, H. Wolbachia endosymbionts in haplodiploid and diploid scolytine beetles (Coleoptera: Curculionidae: Scolytinae). Environ. Microbiol. Rep. 2016, 8, 680–688. [Google Scholar] [CrossRef]
  31. Baldo, L.; Hotopp, J.C.D.; Jolley, K.A.; Bordenstein, S.R.; Biber, S.A.; Choudhury, R.R.; Hayashi, C.; Maiden, M.C.J.; Tettelin, H.; Werren, J.H. Multilocus sequence typing system for the endosymbiont Wolbachia pipientis. Appl. Environ. Microbiol. 2006, 72, 7098–7110. [Google Scholar] [CrossRef] [Green Version]
  32. Hurst, G.D.; Jiggins, F.M.; von der Schulenburg, J.H.G.; Bertrand, D.; West, S.A.; Goriacheva, I.I.; Zakharov, I.A.; Werren, J.H.; Stouthamer, R.; Majerus, M.E.N. Male–killing Wolbachia in two species of insect. Proc. R. Soc. Lond. B Biol. Sci. 1999, 266, 735–740. [Google Scholar] [CrossRef] [Green Version]
  33. Richardson, M.F.; Weinert, L.A.; Welch, J.J.; Linheiro, R.S.; Magwire, M.M.; Jiggins, F.M.; Bergman, C.M. Population genomics of the Wolbachia endosymbiont in Drosophila melanogaster. PLoS Genet. 2012, 8, e1003129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Ilinsky, Y. Coevolution of Drosophila melanogaster mtDNA and Wolbachia genotypes. PLoS ONE 2013, 8, e54373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Bykov, R.A.; Yudina, M.A.; Gruntenko, N.E.; Zakharov, I.K.; Voloshina, M.A.; Melashchenko, E.S.; Danilova, M.V.; Mazunin, I.O.; Ilinsky, Y.Y. Prevalence and genetic diversity of Wolbachia endosymbiont and mtDNA in Palearctic populations of Drosophila melanogaster. BMC Evol. Biol. 2019, 19, 48. [Google Scholar] [CrossRef] [Green Version]
  36. Kerchev, I.A. On monogyny of the four-eyed fir bark beetle Polygraphus proximus Blandf. (Coleoptera, Curculionidae: Scolytinae) and its reproductive behavior. Entomol. Rev. 2014, 94, 1059–1066. [Google Scholar] [CrossRef]
  37. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 1994, 3, 294–299. [Google Scholar]
  38. Werren, J.H.; Windsor, D.M. Wolbachia infection frequencies in insects: Evidence of a global equilibrium? Proc. R. Soc. Lond. B Biol. Sci. 2000, 267, 1277–1285. [Google Scholar] [CrossRef] [Green Version]
  39. Lo, N.; Casiraghi, M.; Salati, E.; Bazzocchi, C.; Bandi, C. How many Wolbachia supergroups exist? Mol. Biol. Evol. 2002, 19, 341–346. [Google Scholar] [CrossRef] [Green Version]
  40. Che, L.H.; Zhang, S.Q.; Li, Y.; Liang, D.; Pang, H.; Ślipiński, A.; Zhang, P. Genome-wide survey of nuclear protein-coding markers for beetle phylogenetics and their application in resolving both deep and shallow-level divergences. Mol. Ecol. Resour. 2017, 17, 1342–1358. [Google Scholar] [CrossRef]
  41. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  42. Kononov, A.; Ustyantsev, K.; Blinov, A.; Fet, V.; Baranchikov, Y.N. Genetic diversity of aboriginal and invasive populations of four-eyed fir bark beetle Polygraphus proximus Blandford (Coleoptera, Curculionidae, Scolytinae). Agric. For. Entomol. 2016, 18, 294–301. [Google Scholar] [CrossRef]
  43. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Arai, H.; Hirano, T.; Akizuki, N.; Abe, A.; Nakai, M.; Kunimi, Y.; Inoue, M.N. Multiple infection and reproductive manipulations of Wolbachia in Homona magnanima (Lepidoptera: Tortricidae). Microb. Ecol. 2019, 77, 257–266. [Google Scholar] [CrossRef] [PubMed]
  45. Dutton, T.J.; Sinkins, S.P. Strain-specific quantification of Wolbachia density in Aedes albopictus and effects of larval rearing conditions. Insect Mol. Biol. 2004, 13, 317–322. [Google Scholar] [CrossRef] [PubMed]
  46. Tortosa, P.; Charlat, S.; Labbe, P.; Dehecq, J.S.; Barré, H.; Weill, M. Wolbachia age-sex-specific density in Aedes albopictus: A host evolutionary response to cytoplasmic incompatibility? PLoS ONE 2010, 5, e9700. [Google Scholar] [CrossRef] [Green Version]
  47. Koehncke, A.; Telschow, A.; Werren, J.H.; Hammerstein, P. Life and death of an influential passenger: Wolbachia and the evolution of CI-modifiers by their hosts. PLoS ONE 2009, 4, e4425. [Google Scholar] [CrossRef] [Green Version]
  48. Jeyaprakash, A.; Hoy, M.A. Long PCR improves Wolbachia DNA amplification: Wsp sequences found in 76% of sixty-three arthropod species. Insect Mol. Biol. 2000, 9, 393–405. [Google Scholar] [CrossRef]
  49. Arthofer, W.; Riegler, M.; Schneider, D.; Krammer, M.; Miller, W.J.; Stauffer, C. Hidden Wolbachia diversity in field populations of the European cherry fruit fly, Rhagoletis cerasi (Diptera, Tephritidae). Mol. Ecol. 2009, 18, 3816–3830. [Google Scholar] [CrossRef]
  50. Hughes, G.L.; Allsopp, P.G.; Brumbley, S.M.; Woolfit, M.; McGraw, E.A.; O’Neill, S.L. Variable infection frequency and high diversity of multiple strains of Wolbachia pipientis in Perkinsiella planthoppers. Appl. Environ. Microbiol. 2011, 77, 2165–2168. [Google Scholar] [CrossRef] [Green Version]
  51. Schuler, H.; Arthofer, W.; Krumböck, S.; Bertheau, C.; Stauffer, C. Wolbachia infection in the Walnut-husk fly Rhagoletis completa Cresson 1929 (Diptera: Tephritidae). Mitt. Dtsch. Ges. Allg. Angew. Entomol. 2012, 18, 243–245. [Google Scholar]
  52. Mee, P.T.; Weeks, A.R.; Walker, P.J.; Hoffmann, A.A.; Duchemin, J.B. Detection of low-level Cardinium and Wolbachia infections in Culicoides. Appl. Environ. Microbiol. 2015, 81, 6177–6188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Schneider, D.I.; Parker, A.G.; Abd-alla, A.M.; Miller, W.J. High-sensitivity detection of cryptic Wolbachia in the African tsetse fly (Glossina spp.). BMC Microbiol. 2018, 18, 140. [Google Scholar] [CrossRef] [PubMed]
  54. Uldig, C.K.S.; Wenzel, M.; Fernandez-Begne, P.; Bubnell, J.; Aquadro, C.F. The persistence of low-titer Wolbachia Pipientis Infection in Antibiotic-Treated Drosophila. TAGC 2020. Available online: https://tagc2020.figshare.com/articles/poster/The_persistence_of_low-titer_Wolbachia_pipientis_infection_in_antibiotic-treated_Drosophila/12150591/1 (accessed on 12 August 2020). [CrossRef]
  55. Yudina, M.A.; Bykov, R.A.; Kotti, B.K.; Vysochina, N.P.; Stakheev, V.V.; Broshkov, A.D.; Zakharov, I.K.; Ilinsky, Y.Y. Wolbachia Infection in Flea Populations (Insecta: Siphonaptera). Biol. Bull. Rev. 2019, 9, 403–411. [Google Scholar] [CrossRef]
Figure 1. Map of the native and invasive ranges of Polygraphus proximus and the range of Abies sibirica. Sampling sites for P. proximus are indicated (the same numbers in Table 1 and Table 2): 1—Tomarinsky-1; 2—Tomarinsky-2; 3—Yuzhno-Sakhalinsk (Sakhalin province); 4—Bolshemurtinsky district; 5—Emelianovsk forestry; 6—Kozulsk forestry (Krasnoyarsk krai); 7—Tomsky; 8—Krivosheinsky; 9—Molchanovsky; 10—Verkhneketsky; 11—Chainsky; 12—Bakcharsky (Tomsk province); 13—Malopurginsky district (Udmurtia).
Figure 1. Map of the native and invasive ranges of Polygraphus proximus and the range of Abies sibirica. Sampling sites for P. proximus are indicated (the same numbers in Table 1 and Table 2): 1—Tomarinsky-1; 2—Tomarinsky-2; 3—Yuzhno-Sakhalinsk (Sakhalin province); 4—Bolshemurtinsky district; 5—Emelianovsk forestry; 6—Kozulsk forestry (Krasnoyarsk krai); 7—Tomsky; 8—Krivosheinsky; 9—Molchanovsky; 10—Verkhneketsky; 11—Chainsky; 12—Bakcharsky (Tomsk province); 13—Malopurginsky district (Udmurtia).
Insects 11 00547 g001
Figure 2. Wolbachia titer of P. proximus estimated by quantitative real time PCR.
Figure 2. Wolbachia titer of P. proximus estimated by quantitative real time PCR.
Insects 11 00547 g002
Figure 3. The maximum likelihood (ML) phylogenetic tree of Wolbachia isolates was reconstructed based on concatenated sequences of five multilocus sequence typing (MLST) genes using the GTR model of nucleotide replacement. The sequence type (ST), ID (id) numbers from the PubMLST database, and GenBank accession numbers, with host species, supergroups, and bootstrap values higher than 75 (1000 replicates), are indicated. Wolbachia isolates of the Scolytinae species are underlined. Studied isolates of ST-533 are in bold.
Figure 3. The maximum likelihood (ML) phylogenetic tree of Wolbachia isolates was reconstructed based on concatenated sequences of five multilocus sequence typing (MLST) genes using the GTR model of nucleotide replacement. The sequence type (ST), ID (id) numbers from the PubMLST database, and GenBank accession numbers, with host species, supergroups, and bootstrap values higher than 75 (1000 replicates), are indicated. Wolbachia isolates of the Scolytinae species are underlined. Studied isolates of ST-533 are in bold.
Insects 11 00547 g003
Figure 4. The ML phylogenetic tree of P. proximus was reconstructed with a Tamura 3-parameter model of nucleotide replacement based on sequences of a 620 bp region of the COI gene. Regions of collection, GenBank accession numbers, and bootstrap values higher than 75 (1000 replicates) are indicated. Samples of the study are indicated in bold. Wolbachia-infected samples are indicated with an asterisk (*). Haplogroups of mtDNA follow Kononov et al. [42].
Figure 4. The ML phylogenetic tree of P. proximus was reconstructed with a Tamura 3-parameter model of nucleotide replacement based on sequences of a 620 bp region of the COI gene. Regions of collection, GenBank accession numbers, and bootstrap values higher than 75 (1000 replicates) are indicated. Samples of the study are indicated in bold. Wolbachia-infected samples are indicated with an asterisk (*). Haplogroups of mtDNA follow Kononov et al. [42].
Insects 11 00547 g004
Table 1. Collection of Polygraphus proximus.
Table 1. Collection of Polygraphus proximus.
Locality No.Sampling Sites, Coordinates (Where Available), and Date of CollectionNSex
Unknown
Sakhalin province (in total)60232413
1Tomarinsky-1, 48°29′22.2″ N 142°01′49.7″ E, 8 June 201837-2413
2Tomarinsky-2, 48°30′26.6″ N 142°00′43.9″ E, 8 June 201888--
3Yuzhno-Sakhalinsk, 46°57′50.8″ N 142°45′16.7″ E, 11 June 20181515--
Krasnoyarsk krai (in total)48-3414
4Bolshemurtinsky district, 25 December 20192-11
5Emelianovsk forestry, December 201916-88
6Kozulsk forestry, 12–17 February 202030-255
Tomsk province (in total)2711486360
7Tomsky, 56°27′55.0″ N 85°06′46.0″ E, 28 August 20183333--
8Krivosheinsky, 57°24′27.0″ N 83°55′16.0″ E, 7 August 20182626--
9Molchanovsky, 57°29′34.0″ N 84°16′27.0″ E, 18 August 20183131--
10Verkhneketsky, 58°23′30.2″ N 84°06′35.8″ E, 29 August 201989-4643
11Chainsky, 57°47′17.0″ N 82°12′34.0″ E, 10 August 20183434--
12Bakcharsky, 57°16′35.0″ N 81°30′18.0″ E, 15 August 201858241717
Udmurtia
13Malopurginsky district, 56°38′40.0″ N 53°05′57.6″ E, 11 October 20193851716
Total:417176138103
Table 2. Wolbachia infection in P. proximus.
Table 2. Wolbachia infection in P. proximus.
Locality No.Region, LocalityConventional PCRConventional and Nested PCR
Nuninf/N♀inf/N♂infTotal %inf (95%CI)Nuninf/N♀inf/N♂infTotal %inf (95%CI)
Sakhalin province (in total)8/23/458.33 (44.88–70.93)20/23/986.67 (75.41–94.06)
1Tomarinsky-10/23/472.97 (55.88–86.21)0/23/986.49 (71.23–95.46)
2Tomarinsky-20/0/0-7/0/087.5 (47.35–99.68)
3Yuzhno-Sakhalinsk8/0/053.33 (26.59–78.73)13/0/086.67 (59.54–98.34)
Krasnoyarsk Krai (in total)0/25/460.42 (45.27–74.23)0/25/664.58 (49.46–77.84)
4Bolshemurtinsky district0/1/0-0/1/0-
5Emelianovsk forestry0/8/368.75 (41.34–88.98)0/8/581.25 (54.35–95.95)
6Kozulsk forestry0/16/156.67 (37.43–74.54)0/16/156.67 (37.43–74.54)
Tomsk province (in total)74/41/1547.97 (41.89–54.1)87/45/3461.26 (55.17–67.09)
7Tomsky20/0/060.61 (42.14–77.09)22/0/066.67 (48.17–82.04)
8Krivosheinsky20/0/076.92 (56.35–91.03)21/0/080.77 (60.65–93.45)
9Molchanovsky12/0/038.71 (21.85–57.81)15/0/048.39 (30.15–66.94)
10Verkhneketsky1/24/1241.57 (31.21–52.51)1/28/2459.55 (48.62–69.83)
11Chainsky15/0/044.12 (27.19–62.11)19/0/055.88 (37.89–72.82)
12Bakcharsky6/17/344.83 (31.74–58.46)9/17/1062.07 (48.37–74.49)
Udmurtia
13Malopurginsky district3/17/155.26 (38.30–71.38)4/17/1492.11 (78.62–98.34)
Total:85/106/2451.56 (46.65–56.45)111/110/6368.11 (63.40–72.56)
NuninfWolbachia-infected samples with undetermined sex; NinfWolbachia-infected females; NinfWolbachia-infected males.

Share and Cite

MDPI and ACS Style

Bykov, R.; Kerchev, I.; Demenkova, M.; Ryabinin, A.; Ilinsky, Y. Sex-Specific Wolbachia Infection Patterns in Populations of Polygraphus proximus Blandford (Coleoptera; Curculionidae: Scolytinae). Insects 2020, 11, 547. https://doi.org/10.3390/insects11080547

AMA Style

Bykov R, Kerchev I, Demenkova M, Ryabinin A, Ilinsky Y. Sex-Specific Wolbachia Infection Patterns in Populations of Polygraphus proximus Blandford (Coleoptera; Curculionidae: Scolytinae). Insects. 2020; 11(8):547. https://doi.org/10.3390/insects11080547

Chicago/Turabian Style

Bykov, Roman, Ivan Kerchev, Marya Demenkova, Artem Ryabinin, and Yury Ilinsky. 2020. "Sex-Specific Wolbachia Infection Patterns in Populations of Polygraphus proximus Blandford (Coleoptera; Curculionidae: Scolytinae)" Insects 11, no. 8: 547. https://doi.org/10.3390/insects11080547

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