# Diversification, selective sweep, and body size in the invasive Palearctic alfalfa weevil infected with Wolbachia

Midori Tuda, Shun-Ichiro Iwase, Khadim Kébé, Julien Haran, Jiri Skuhrovec, Ehsan Sanaei, Naomichi Tsuji, Attila Podlussány, Ottó Merkl, Ahmed El-Heneidy, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Midori Tuda, Shun-Ichiro Iwase, Khadim Kébé, Julien Haran, Jiri Skuhrovec, et al. Diversification, selective sweep, and body size in the invasive Palearctic alfalfa weevil infected with Wolbachia. Scientific Reports, 2021, 11, pp.9664. 10.1038/s41598-021-88770-y. hal-03229157

HAL Id: hal-03229157

https://hal.inrae.fr/hal-03229157

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

OPEN

# Diversification, selective sweep, and body size in the invasive Palearctic alfalfa weevil infected with *Wolbachia*

Midori Tuda<fnref n="1" />,²,¹²✉, Shun-ichiro Iwase<fnref n="1" />,¹¹,¹², Khadim Kébé³,¹², Julien Haran⁴,¹², Jiri Skuhrovec⁵,¹², Ehsan Sanaei⁶, Naomichi Tsuji⁷, Attila Podlussány⁸, Ottó Merkl⁸, Ahmed H. El-Heneidy⁹ & Katsura Morimoto¹⁰

The alfalfa weevil *Hypera postica*, native to the Western Palearctic, is an invasive legume pest with two divergent mitochondrial clades in its invading regions, the Western clade and the Eastern/Egyptian clade. However, knowledge regarding the native populations is limited. The Western clade is infected with the endosymbiotic bacteria *Wolbachia* that cause cytoplasmic incompatibility in host weevils. Our aim was to elucidate the spatial genetic structure of this insect and the effect of *Wolbachia* on its population diversity. We analyzed two mitochondrial and two nuclear genes of the weevil from its native ranges. The Western clade was distributed in western/central Europe, whereas the Eastern/ Egyptian clade was distributed from the Mediterranean basin to central Asia. Intermediate mitotypes were found from the Balkans to central Asia. Most Western clade individuals in western Europe were infected with an identical *Wolbachia* strain. Mitochondrial genetic diversity of the infected individuals was minimal. The infected clades demonstrated a higher nonsynonymous/synonymous substitution rate ratio than the uninfected clades, suggesting a higher fixation of nonsynonymous mutations due to a selective sweep by *Wolbachia*. Trans-Mediterranean and within-European dispersal routes were supported. We suggest that the ancestral populations diversified by geographic isolation due to glaciations and that the diversity was reduced in the west by a recent *Wolbachia*-driven sweep(s). The intermediate clade exhibited a body size and host plant that differed from the other clades. Pros and cons of the possible use of infected-clade males to control uninfected populations are discussed.

Recent invasion events and routes of alien agricultural pests are of particular importance for the management and control of pests<fnref n="1" />⁻⁶. The knowledge of historical diversification and dispersal of agricultural pests in their native range provides insights to understand their natural and biological selective environments, including the role played by endosymbionts in pest emergence⁷⁻⁹.

The alfalfa weevil *Hypera postica* (Gyllenhal) (Coleoptera: Curculionidae: Hyperini), native to the Western Palearctic region, is a serious pest of alfalfa and other beneficial legumes in its invading territories, such as *Medicago*, *Vicia*, *Trifolium*, and *Astragalus*¹⁰⁻¹² (Nearctic, Japan, Southeast Asia, and Oceania¹³⁻¹⁵). Invading populations in the USA comprise the Western (North American) type that invaded Utah in 1904¹⁶, the Egyptian type that invaded Arizona in 1939¹⁷, and the Eastern (North American) type that invaded Maryland in 1951¹⁸.

<footnote n="1">¹Institute of Biological Control, Faculty of Agriculture, Kyushu University, Fukuoka 819-0395, Japan. ²Laboratory of Insect Natural Enemies, Department of Bioresource Sciences, Faculty of Agriculture, Kyushu University, Fukuoka, Japan. ³GRBA-BE, LE3PI Laboratory, Department of Chemical Engineering and Applied Biology, Polytechnic Higher School of Dakar, Dakar, Senegal. ⁴CBGP, Cirad, Montpellier SupAgro, INRA, IRD, Univ. Montpellier, Montpellier, France. ⁵Group Function of Invertebrate and Plant Biodiversity in Agro-Ecosystems, Crop Research Institute, Drnovska, Praha, Czech Republic. ⁶School of Biological Sciences, University of Queensland, Brisbane, Australia. ⁷Entomological Laboratory, Graduate School of Bioresource and Bioenvironmental Sciences, Kyushu University, Fukuoka, Japan. ⁸Hungarian Natural History Museum, Budapest, Hungary. ⁹Department of Biological Control, Plant Protection Research Institute, Agricultural Research Center, Giza, Egypt. ¹⁰Nata-Danchi, Fukuoka, Japan. ¹¹Present address: Research Institute of Environment, Agriculture and Fisheries, Osaka Prefecture, Japan. ¹²These authors contributed equally: Midori Tuda, Shun-ichiro Iwase, Khadim Kébé, Julien Haran and Jiri Skuhrovec. ✉email: tuda@grt.kyushu-u.ac.jp</footnote>
These types are different in their ecological, behavioral, and defensive traits (pupation site and aggregation during aestivation¹³, defensive behavior¹⁹, and encapsulation of immature endoparasitoid²⁰⁻²²) but are morphologically indistinguishable. Allozyme and mitochondrial DNA markers distinguish the Western type from the Egyptian and Eastern types²²⁻²⁴, whereas the Egyptian and Eastern types are distinguishable only by a slight difference (1–2 SNPs in *tRNA*<sup>Ser</sup>) in a mitochondrial gene sequence²²,²⁴,²⁵. The nuclear DNA polymorphism indicates that the three types share a gene pool, namely, a single species²⁵,²⁶. They also mate with each other to reproduce (but see later for incompatibility²⁷,²⁸).

*Wolbachia* (Alphaproteobacteria: Rickettsiales: Rickettsiaceae) are maternally (vertically) transmitted intracel-
lular bacteria that infect approximately 40% of insects and other arthropods alongside nematodes²⁹. These endo-
symbiotic bacteria can manipulate host reproduction via reproductive cells and the genetic mechanism of this
manipulation has recently been uncovered³⁰. Cytoplasmic incompatibility (CI), or postzygotic incompatibility
between infected and uninfected gametes, is the most commonly observed phenotype of *Wolbachia*. Theoretically,
bidirectional CI most strongly accelerates host speciation. Unidirectional CI, or postzygotic isolation of gametes
between infected males and uninfected females, can also promote host speciation. While infection by *Wolbachia*
is favored in females in populations with high *Wolbachia* prevalence, loss of *Wolbachia* can also occur through
incomplete inheritance from mothers with low *Wolbachia* density²⁹. The loss (or incomplete transmission) rate
of *Wolbachia* in insect hosts is slightly higher than the gain rate³¹. Maternal transmission and unidirectional CI
eventually reduce host mitochondrial diversity over generations. The selective sweep of mitochondria leads to a
close association between the mitochondrial clade and *Wolbachia* infection³²⁻³⁵. Several studies have discovered
that *Wolbachia* may also accelerate the fixation of nonsynonymous mutations in hosts²⁹,³⁶⁻³⁸. Various positive
fitness effects of endosymbionts on their hosts have been revealed, such as viral suppression and metabolic
provisioning³⁹,⁴⁰. The effect of endosymbionts on maternal mitochondria may also influence coevolution between
mitochondria and nuclear genomes⁴¹,⁴².

In *H. postica*, the Western clade is found to be infected by *Wolbachia* that induces unidirectional CI<sup>27,43,44</sup>.
Several invading populations of the Western clade are free of *Wolbachia*, and a cross between uninfected Western
males and Egyptian/Eastern females within these populations does produce viable offspring<sup>28</sup>. This reconfirms
that these clades, while genetically distant, remain conspecific. The CI effect between infected Western males
and uninfected Eastern females is almost perfect (only 0.1% of hybrid eggs hatch), while 29.5% of hybrid eggs
between uninfected Eastern males and infected Western females hatch<sup>27</sup>.

The presence of diverged clades and *Wolbachia* infection history in Palearctic *H. postica* in its primary range is not known to date. Here, this study aims to explore the process of selection and diversification in *H. postica* in its native range by revealing and testing mitochondrial and nuclear genetic variation geographically and phylogenetically. We also aim to test if the endosymbiont *Wolbachia* affected evolution in host weevils. The benefits and risks of the Incompatible Insect Technique⁴⁵ using infected clade males to control the uninfected clade populations are discussed.

**Results**

**Haplotype networks and diversity.** Sequenced segments were 2001 bp; 527 bp for COI-tRNA<sup>Leu</sup>-COII, 281 bp for *Cyt b*-tRNA<sup>Ser</sup>-ND1 (*n* = 149), 801 bp for 28S (*n* = 122) and 392 bp for *EF*-1*α* (*n* = 62) (Table 1). Despite our sizable effort, PCR failed for nuclear gene segments (especially *EF*-1*α*) for a part of the specimens. The mitotype network revealed two main clades alongside intermediate variants (Fig. 1). The first group corresponded to the Eastern/Egyptian clade (Fig. 1) and consisted of diverse mitotypes with multiple connections, which contributed to a high mitochondrial genetic diversity (Table 2). This clade was widely distributed from central Asia to the Mediterranean region (Fig. 2). Within this clade, populations from the Balkan peninsula displayed high mitochondrial and nuclear genetic diversity (Table 3). The second clade corresponded to the Western clade and exhibited substantially fewer mitotypes, one dominant mitotype and rarer, closely related mitotypes in a star-shape topology, which corresponds to low genetic diversity in both gene fragments (Fig. 1, Table 2). This clade was distributed in western and central Europe, north of the Alps and Pyrenees (Fig. 2). Within this clade, 50.0% of individuals and 50.0% of populations were infected with *Wolbachia* (Fig. 2). The *Wolbachia*-infected populations demonstrated lower mitochondrial and nuclear genetic diversity than the uninfected populations (Table 2). Compared with the uninfected Eastern/Egyptian clade, the infected Western clade displayed 22 (*Cyt b*-tRNA<sup>Ser</sup>-ND1) to 82 (COI-tRNA<sup>Leu</sup>-COII) times lower mitochondrial genetic diversity (Table 2). The intermediate clade was distributed from the Balkans to central Asia (Fig. 2). The network for the nuclear fragments, 28S and *EF*-1*α*, appeared incongruent with the mitochondrial network (Figs. 1, 3), but as in mitotype variation, there was a significant difference in nuclear haplotype variation between individuals belonging to the different mitochondrial clades (*EF*-1*α*, Table 2).

**Wolbachia infection.** The Western clade individuals in northern France, the Netherlands, the Czech Repub-
lic, and Poland were infected with *Wolbachia* (Fig. 2), and all of them had identical *ftsZ* (699 bp), *coxA* (432 bp),
and *hcpA* (463 bp) sequences as reported previously⁴⁶, which corresponds to the prevailing strain *wHypera1*⁴⁷. In
contrast, Western clade individuals found in coastal southern France, Spain, Latvia, Hungary, and Croatia were
uninfected (Fig. 2). Additionally, we confirmed that *Wolbachia* was absent from all native populations studied
in the Eastern/Egyptian clade. A reconstructed phylogenetic relationship confirms that *wHypera1* belongs to
Supergroup B (Fig. 4). The *Wolbachia* strain closest to *wHypera1* to date is the one that infects the mite *Bryobia
praetiosa* (Acari: Tetranychidae)⁴⁸.

**Body size.** Elytral length was significantly different among different clades (χ²₂ = 10.96, p = 0.004) and between sexes (χ²₁ = 7.52, p = 0.006) and marginally different between infected and uninfected individuals
<table><thead><tr><th>Number on the map</th><th>Code</th><th>Collection site</th><th>n</th><th>Year</th><th>Latitude, longitude</th><th>GenBank accession<br>Mitochondrial genes</th><th>GenBank accession Nuclear<br>genes</th></tr></thead><tbody><tr><td colspan="8"><strong>Western Europe</strong></td></tr><tr><td>1</td><td>Cz</td><td>Prague, Czech Republic</td><td>5</td><td>2012</td><td>50°05'16"N,14°17'54"E</td><td>KX372573 (CO), KX372620 (CB)</td><td>KX372667 (28S), MW392102<br>(EF1a)</td></tr><tr><td>2</td><td>Ne</td><td>Amsterdam, the Netherlands</td><td>2</td><td>2014</td><td>52°21'35"N,4°57'00"E</td><td>MW393903 (CO), MW393922<br>(CB)</td><td>MW383444 (28S), MW389094<br>(EF1a)</td></tr><tr><td>3</td><td>Po</td><td>Nida Basin, Poland</td><td>2</td><td>2012</td><td>50°24'N,20°38'E</td><td>MW393912 (CO), MW393931<br>(CB)</td><td>MW383460 (28S), MW389110<br>(EF1a)</td></tr><tr><td>4</td><td>Lit</td><td>Daugavpils, Latvia</td><td>2</td><td>2007</td><td>55°52' N,26°27'E</td><td>MW393914 (CO), MW393933<br>(CB)</td><td></td></tr><tr><td colspan="8"><strong>Central Europe</strong></td></tr><tr><td>5</td><td>BuH</td><td>Budapest, Hungary</td><td>4</td><td>2012</td><td>47°28'41" N,19°01'00" E</td><td>MW393904 (CO), MW393923<br>(CB)</td><td>MW383445 (28S), MW389095<br>(EF1a)</td></tr><tr><td>6</td><td>AdH</td><td>Adyliget, Budapest, Hungary</td><td>9</td><td>2014</td><td>47°32'40" N,18°55'58" E</td><td>KX372574-75 (CO), KX372621-22 (CB)</td><td>MW383446 (28S), MW389096<br>(EF1a)</td></tr><tr><td colspan="8"><strong>France</strong></td></tr><tr><td>7</td><td>ChF</td><td>Chaussy, France</td><td>4</td><td>2013</td><td>49°07'12"N, 1°42'02"E</td><td>MW393905 (CO), MW393924<br>(CB)</td><td>MW383448 (28S), MW389098<br>(EF1a)</td></tr><tr><td>8</td><td>OrF</td><td>Orléans, France</td><td>12</td><td>2013</td><td>47°53'59"N, 1°56'24"E</td><td>MW393906,16 (CO),<br>MW393925,35 (CB)</td><td>MW383449,62 (28S),<br>MW389099,112 (EF1a)</td></tr><tr><td>9</td><td>AuF</td><td>Auradé, France</td><td>3</td><td>2014</td><td>43°33'36"N, 1°03'01"E</td><td>KX372576 (CO), KX372623 (CB)</td><td>MW383447 (28S), MW389097<br>(EF1a)</td></tr><tr><td>10</td><td>ComF</td><td>Combaillaux, France</td><td>9</td><td>2016</td><td>43°40'12"N,3°46'47"E</td><td>MW393915 (CO), MW393934<br>(CB)</td><td></td></tr><tr><td>11</td><td>AlpF</td><td>Saint-Paul-sur-Übaye, France</td><td>2</td><td>2010</td><td>44°31'12"N, 6°45'02"E</td><td>KX372579 (CO), KX372626 (CB)</td><td>MW383450 (28S), MW389100<br>(EF1a)</td></tr><tr><td>12</td><td>CorF</td><td>Casamozza, Corse Island, France</td><td>5</td><td>2016</td><td>42°30'35"N,9°26'23"E</td><td>MW393910 (CO), MW393929<br>(CB)</td><td>MW383458 (28S), MW389108<br>(EF1a)</td></tr><tr><td colspan="8"><strong>Spain</strong></td></tr><tr><td>13</td><td>Sp</td><td>La Cañada, Spain</td><td>6</td><td>2014</td><td>40°36'00"N,4°30'35"W</td><td>KX372577,82 (CO),<br>KX372624,29 (CB)</td><td>MW383452,61 (28S),<br>MW389102,11 (EF1a)</td></tr><tr><td colspan="8"><strong>Balkans and Italy</strong></td></tr><tr><td>14</td><td>Malt</td><td>Maremma, Toscana, Italy</td><td>3</td><td>2015</td><td>42°38'17"N, 11°07'29"E</td><td>MW393913 (CO), MW393932<br>(CB)</td><td></td></tr><tr><td>15</td><td>PuIt</td><td>Puglia, Italy</td><td>2</td><td>2002</td><td>41°04'N,16°26'E</td><td>KX372599-600 (CO),<br>KX372646,47 (CB)</td><td>MW383455 (28S), MW389105<br>(EF1a)</td></tr><tr><td>16</td><td>Gr</td><td>Sparti, Peloponnese, Greece</td><td>5</td><td>2005</td><td>36°51' N, 22°39' E</td><td>KX372586-87,614-16 (CO),<br>KX372633-34,61-63 (CB)</td><td></td></tr><tr><td>17</td><td>Ro</td><td>Crucea, Romania</td><td>5</td><td>2009</td><td>44°31'12" N, 28°11'59" E</td><td>KX372588,96-98,617 (CO),<br>KX372635,43-45,64 (CB)</td><td></td></tr><tr><td>18</td><td>KnBu</td><td>Knezha, Bulgaria</td><td>10</td><td>2016</td><td>43°28'48" N, 24°03'36" E</td><td>MW393917 (CO), MW393936<br>(CB)</td><td></td></tr><tr><td>19</td><td>LoBu</td><td>Lozitsa, Bulgaria</td><td>8</td><td>2016</td><td>43°34'48"N, 25°00'01"E</td><td>MW393918 (CO), MW393937<br>(CB)</td><td></td></tr><tr><td>20</td><td>Cr</td><td>Zagreb, Croatia</td><td>2</td><td>2006</td><td>45°50'35" N, 15°44'55" E</td><td>KX372578 (CO), KX372625 (CB)</td><td></td></tr><tr><td colspan="8"><strong>Africa</strong></td></tr><tr><td>21</td><td>faMo</td><td>Ouzoud falls, Morocco</td><td>10</td><td>2016</td><td>32°00'54"N,6°43'24"W</td><td>MW393909 (CO), MW393928<br>(CB)</td><td>MW383457 (28S), MW389107<br>(EF1a)</td></tr><tr><td>22</td><td>OuMo</td><td>Ouarzazate, Morocco</td><td>2</td><td>1994</td><td>30°56'N, 6°56'W</td><td>KX372589,618 (CO),<br>KX372636,65 (CB)</td><td></td></tr><tr><td>23</td><td>Li</td><td>Benghazi, Libya</td><td>5</td><td>1980</td><td>32°03'N,20°09'E</td><td>KX372603-04,13 (CO),<br>KX372650-51,60 (CB)</td><td></td></tr><tr><td>24</td><td>Eg</td><td>Sakha, Kafr El-Sheikh Governorate, Egypt</td><td>3</td><td>2013</td><td>31°05'13" N, 30°56'56" E</td><td>KX372580, MW393907 (CO),<br>KX372627, MW393926 (CB)</td><td>MW383451 (28S), MW389101<br>(EF1a)</td></tr><tr><td colspan="8"><strong>Middle East and central Asia</strong></td></tr><tr><td>25</td><td>Isr</td><td>Ga'ash, Israel</td><td>2</td><td>2014</td><td>32°13'48"N, 34°49'12"E</td><td>KX372590 (CO), KX372637 (CB)</td><td>MW383453 (28S), MW389103<br>(EF1a)</td></tr><tr><td>26</td><td>Tur</td><td>Catalan vill, Adana, Turkey</td><td>4</td><td>2002</td><td>37°15'01" N, 35°18'10" E</td><td>KX372581,601-02, MW393908<br>(CO), KX372628,48-49,<br>MW393927 (CB)</td><td></td></tr><tr><td>27</td><td>Ar</td><td>Metsamor, Armenia</td><td>2</td><td>2013</td><td>40°09'19"N, 44°07'30"E</td><td>KX372606 (CO), KX372653 (CB)</td><td>MW383456 (28S), MW389106<br>(EF1a)</td></tr><tr><td>28</td><td>AzIr</td><td>Azerbaijan, Iran</td><td>2</td><td>1999</td><td>37°56' N, 47°23' E</td><td>KX372584-85 (CO), KX372631-32 (CB)</td><td></td></tr><tr><td>29</td><td>TaIr</td><td>Taleghan, Iran</td><td>8</td><td>2015</td><td>36°12'12"N,50°51'47"E</td><td>MW393919 (CO), MW393938<br>(CB)</td><td></td></tr><tr><td colspan="8">Continued</td></tr></tbody></table>
<table><caption>**Table 1.** Sample collection information for *Hypera postica*.CO: COI-tRNA<sup>Leu</sup>-COII. CB: Cyt b-tRNA<sup>Ser</sup>-ND1.</caption><thead><tr><th>Number on the map</th><th>Code</th><th>Collection site</th><th>n</th><th>Year</th><th>Latitude, longitude</th><th>GenBank accession<br>Mitochondrial genes</th><th>GenBank accession Nuclear<br>genes</th></tr></thead><tbody><tr><td>30</td><td>HaIr</td><td>Hamedan, Iran</td><td>3</td><td>2014</td><td>34°31′16″N,48°18′43″E</td><td>MW393911 (CO), MW393930<br>(CB)</td><td>MW383459 (28S), MW389109<br>(EF1α)</td></tr><tr><td>31</td><td>FaIr</td><td>Dasht-E-Arzhan, Zagros Mts.,<br>Fars Prov., Iran</td><td>2</td><td>2000</td><td>29°34′ N, 51°56′ E</td><td>KX372608–09 (CO), KX372655–<br>56 (CB)</td><td></td></tr><tr><td>32</td><td>Tm</td><td>Anau, Turkmenistan</td><td>2</td><td>1988</td><td>37°54′N,58°30′E</td><td>MW393920 (CO), MW393939<br>(CB)</td><td></td></tr><tr><td>33</td><td>Ky</td><td>Jangy-Talap, Kyrgyzstan</td><td>4</td><td>2015</td><td>41°27′01″N, 75°01′12″E</td><td>KX372592 (CO), KX372639 (CB)</td><td>MW383454 (28S), MW389104<br>(EF1α)</td></tr></tbody></table>

($\chi^2$ = 3.77, *p* = 0.052). The elytron was longer in the intermediate variants (3.83 ± 0.110 mm, mean ± SE, *n* = 21) than in the Egyptian/Eastern (3.48 ± 0.025 mm, *n* = 64, *p* = 0.003) and Western (3.50 ± 0.057 mm, *n* = 29, *p* = 0.041) clades (no difference between the Egyptian/Eastern and Western clades, *p* = 0.934), longer in females (3.62 ± 0.042 mm, *n* = 60) than in males (3.47 ± 0.044 mm, *n* = 52), and marginally longer in uninfected individuals (3.57 ± 0.032 mm, *n* = 102) than in infected individuals (3.36 ± 0.111 mm, *n* = 11).

**Selective neutrality test and positive selection test.** For the Western clade, selective neutrality for mitochondrial segments was rejected by all indexes with minus values (*D*, *D*\*, and *F*\*), while for nuclear genes, selective neutrality was rejected by none of the indices (Table 4), suggesting recent sudden population growth after bottleneck event(s) in the mitochondrial lineage. For the Eastern/Egyptian clade, selective neutrality for both mitochondrial and nuclear segments was rejected by all indexes (Table 4), suggesting recent sudden population growth after bottleneck event(s) in the mitochondrial lineage and nuclear genes. For the intermediate clade, selective neutrality for mitochondrial segments was rejected by none of the indexes (Table 4). Selective neutrality for nuclear segments was untested because of the insufficient sample size.

A model with different $\omega$ (dN/dS) assigned for infected and uninfected clades improved the model fit most, compared to models with different $\omega$ for two or three different clades, although the improvement was nonsignificant (Table 5). The $\omega$ for infected clade was three times higher than the $\omega$ for uninfected clade even though both were < 1 (Table 5).

**Geographic history.** *Isolation by distance.* Isolation by distance (IBD) was supported for all populations (*p* = 0.003, number of pairwise comparisons *n* = 528) and for the Eastern/Egyptian clade excluding populations with intermediate mitotypes (*p* = 0.050, *n* = 153), but not for the Western clade (*p* = 0.102, *n* = 55). When the populations with intermediate mitotypes were included, the IBD within the Eastern/Egyptian clade was not supported (*p* = 0.215, *n* = 231).

**Phylogeography.** Based on available samples, the Balkan/Italian peninsulas and the Middle East are the most likely area of the origin of *H. postica*, from which the Western clade diversified via France (Fig. 5a). France is the likely area where the ancestral population was first infected with *Wolbachia* (Fig. 5a, right). We found two connections between regions that were highly supported with BF > 3: France and western Europe (BF = 7.98) and Balkan/Italy and North Africa (BF = 7.90) (BSSVS analysis, Fig. 5b).

Discussion

This study revealed the large-scale geographic distribution and genetic diversity of *H. postica* in its native range. Intermediate mitotypes with larger body sizes were found from the Balkans to central Asia. The observed reduced diversity within the Western clade is likely due to a high percentage of *Wolbachia* infection within this clade, which is known in other species. We also identified a higher substitution rate of nonsynonymous mutations, suggesting promoted (fixation of) nonsynonymous mutations in the infected Western clade. In the Western clade, recent sudden population growth after a bottleneck was suggested only for mitochondrial genes and not for nuclear genes, supporting a recent selective sweep on mitochondria by *Wolbachia* infection.

Our results demonstrated a clear pattern of geographic distribution of the two divergent mitochondrial clades across the area of study, the Eastern/Egyptian and Western clades. The low genetic variation and star-like haplotype network within the Western clade is a signature of a recent demographic expansion from a few founders⁴⁹. If the populations experienced ancient demographic bottlenecks, mitochondrial and nuclear genes are expected to have a concordant population structure. There was a weak concordance (Table 2), suggesting that these genes may have shared similar evolutionary trajectories. Two bottleneck events are likely; postglacial recolonization (see next section) and a recent mitochondrial sweep by *Wolbachia*. The former may also serve as a major driver of IBD that was supported for overall geographic populations (two clades distributed separately in the north and south with intermediate mitotypes in between). The latter likely accelerates the fixation of nonsynonymous mutations in the Western clade³³,³⁵,⁵⁰⁻⁵³.

The asymmetric inheritance of maternal mitochondria of an infected host caused by unidirectional CI-inducing *Wolbachia* can eventually lead to a sweep, which likely explains the low mitochondrial genetic variation among infected individuals. The infected clade demonstrated accelerating nonsynonymous mutations or fixation. This result is consistent with a general trend of *Wolbachia* infected insect groups³⁸, suggesting fixation of nonsynonymous mutations in mitochondria promoted by its small effective population size under the CI-inducing
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Statistical parsimony network of mitochondrial COI-tRNA<sup>Leu</sup>-COII and Cyt b-tRNA<sup>Ser</sup>-ND1 of *Hypera postica* in its native range. Generated using TCS 1.21<sup>80</sup>.</figcaption></figure>
<table><thead><tr><th rowspan="3">Clade<br>(Wolbachia<br>infection)</th><th colspan="4">Mitochondria</th><th colspan="4">Nuclear</th></tr><tr><th colspan="2">COI-tRNA<sup>Leu</sup>-COII</th><th colspan="2">Cyt b-tRNA<sup>Ser</sup>-ND1</th><th colspan="2">28S</th><th colspan="2">EF-1α</th></tr><tr><th>n</th><th>Nucleotide<br>diversity</th><th>n</th><th>Nucleotide<br>diversity</th><th>n</th><th>Nucleotide<br>diversity</th><th>n</th><th>Nucleotide<br>diversity</th></tr></thead><tbody><tr><td>Western<br>(infected)</td><td>0.074 ± 0.150<br>(27, 6)</td><td>0.00014 ± 0.00031</td><td>a</td><td>0.148 ± 0.217<br>(27, 6)</td><td>0.00053 ± 0.00086</td><td>a</td><td>0.262 ± 0.299<br>(27, 6)</td><td>0.00033 ± 0.00042</td><td>a</td><td>3.134 ± 1.689<br>(22, 5)</td><td>0.00802 ± 0.00482</td><td>a</td></tr><tr><td>Western<br>(uninfected)</td><td>3.481 ± 1.832<br>(27, 6)</td><td>0.00661 ± 0.00387</td><td></td><td>0.729 ± 0.560<br>(27, 6)</td><td>0.00260 ± 0.00222</td><td></td><td>1.057 ± 0.724<br>(25, 6)</td><td>0.00132 ± 0.00101</td><td></td><td>6.107 ± 3.251<br>(8, 3)</td><td>0.01554 ± 0.00942</td><td></td></tr><tr><td>Egyptian/<br>Eastern<br>(uninfected)</td><td>6.057 ± 2.917<br>(74, 23)</td><td>0.01149 ± 0.00614</td><td>b</td><td>3.295 ± 1.714<br>(74, 23)</td><td>0.01173 ± 0.00676</td><td>b</td><td>1.052 ± 0.708<br>(60, 20)</td><td>0.00131 ± 0.00098</td><td>a</td><td>5.335 ± 2.646<br>(31, 12)</td><td>0.01334 ± 0.00736</td><td>b</td></tr><tr><td>Intermediate<br>(uninfected)</td><td>10.300 ± 4.897<br>(21, 4)</td><td>0.01955 ± 0.01037</td><td>c</td><td>0.733 ± 0.567<br>(21, 4)</td><td>0.00261 ± 0.00225</td><td>c</td><td>2.111 ± 1.282<br>(10, 3)</td><td>0.00264 ± 0.00181</td><td>a</td><td>0 (1, 1)</td><td>0</td><td>ab</td></tr><tr><td>Among-clade<br>variation<br>(df=2), p</td><td colspan="2">72.53%, &lt;0.00001</td><td></td><td colspan="2">87.60%, &lt;0.00001</td><td></td><td colspan="2">3.34%, 0.0168</td><td></td><td colspan="2">6.50%, 0.0544</td><td></td></tr></tbody></table>

**Table 2.** Mitochondrial and nuclear genetic diversity of the two clades (the Western and the Egyptian/ Eastern) and the intermediate clades in *Hypera postica*. n: mean number of pairwise differences; nucleotide diversity (average over loci) (mean ± SD). The genetic distance was calculated based on pairwise differences. Numbers in parentheses are sample sizes (pooled number of individuals, number of populations). The same letters within each column (gene segment) indicate no significant difference between clades (p > 0.01). For population-wise mitochondrial and nuclear diversities and distances, see Supplementary Tables S1 and S2.

<figure><img src="image_5.png" /></figure>

**Figure 2.** Geographic distribution of mitochondrial clades of *Hypera postica* in its native range. Pie chart sizes for clades are proportional to sample sizes. *Wolbachia* infection (dark gray ring: infected; white (partial) ring: uninfected) is only indicated for the Western clade individuals because none of the Egyptian/Eastern clade or the intermediate clade were infected. The background map was obtained from Fotolla, https://stock.adobe.com/jp/photos/.

*Wolbachia* infection. Furthermore, *Wolbachia* infection is advantageous for *H. postica* by enhancing resistance against its adult parasitoid *Microctonus aethiopoides*²².

In the southernmost populations of the Western clade or the geographic contact zone between the two clades, most individuals were uninfected or had lost *Wolbachia*. The imperfect maternal transmission was observed in the interclade crosses of *H. postica*²⁷; fitness costs incurred by cytoplasmic incompatibility and stochasticity during the invasion process⁴⁵,⁵⁴ may lower the *Wolbachia* infection rate. Environmental causes (e.g., extreme temperatures) may also accelerate the endosymbiont loss⁵⁵,⁵⁶. Resulting uninfected *H. postica* populations (or with lowered *Wolbachia* density) must have regained reproductive compatibility between clades and enabled crosses between the diverged clades.

The intermediate variants exhibited a large body size. Larger genitalia of the males with these mitotypes may inhibit mating with the females of other clades and promote reproductive isolation. These mitotypes also were associated with an ecological niche that differed from the niche of other clades. Bulgaria populations used *Vicia cracca* as a host plant, whereas other populations used *Medicago* and *Trifolium*. *Vicia cracca* has high contents of cyanamide⁵⁷ and canavanine⁵⁸ that are toxic to insect herbivores⁵⁹,⁶⁰.

The genetic structure of most European biota has been strongly influenced by glacial oscillations of the Holocene⁶¹,⁶², and most temperate species exhibit northward post-glacial recolonization from glacial refugia located in southern Europe through central Asia during the last glacial maxima (southern genetic richness/ northern purity⁶³; in beetles⁶⁴). In *H. postica*, we observed mtDNA differentiation for all clades and within the Eastern/Egyptian clade. Based on the estimated ancestral states in mitochondrial phylogeography and mitochondrial/nuclear genetic diversity, the Balkan and Italian peninsulas are a possible candidate for the origin of
<table><caption>**Table 3.** Mitochondrial and nuclear genetic diversity in *Hypera postica* based on geographic regions. Π (mean number of pairwise differences); nucleotide diversity (average over loci) (mean ± SD). Bolded genetic diversity indices indicate the highest diversity for each gene segment. Numbers in parentheses are sample sizes (pooled number of individuals, number of sampled populations). See Table 1 for the country codes. For population-wise mitochondrial and nuclear diversities and distances, see Supplementary Tables S1 and S2.</caption><thead><tr><th rowspan="3">Geographic region (country)</th><th colspan="4">Mitochondria</th><th colspan="4">Nuclear</th></tr><tr><th colspan="2">COI-tRNA<sup>Leu-COII</sup></th><th colspan="2">Cyt b-tRNA<sup>Ser-ND1</sup></th><th colspan="2">28S</th><th colspan="2">EF-1α</th></tr><tr><th>Π</th><th>Nucleotide diversity</th><th>Π</th><th>Nucleotide diversity</th><th>Π</th><th>Nucleotide diversity</th><th>Π</th><th>Nucleotide diversity</th></tr></thead><tbody><tr><td>Western Europe (Cz, Ne, Lit, Po)</td><td>0.51 ± 0.46 (11, 4)</td><td>0.0010 ± 0.0010</td><td>0.18 ± 0.25 (11, 4)</td><td>0.0006 ± 0.0010</td><td>0.40 ± 0.40 (10, 4)</td><td>0.0005 ± 0.0006</td><td>3.07 ± 1.85 (6, 2)</td><td>0.0078 ± 0.0055</td></tr><tr><td>France (Fr)</td><td>9.92 ± 4.65 (35, 6)</td><td>0.0188 ± 0.0098</td><td>7.06 ± 3.40 (35, 6)</td><td>0.0251 ± 0.0134</td><td>0.58 ± 0.48 (35, 6)</td><td>0.0007 ± 0.0007</td><td>4.48 ± 2.30 (21, 6)</td><td>0.0112 ± 0.0064</td></tr><tr><td>Spain (Sp)</td><td>9.13 ± 4.90 (6, 1)</td><td>0.0173 ± 0.0107</td><td><strong>8.47</strong> ± 4.57 (6, 1)</td><td><strong>0.0301</strong> ± 0.0188</td><td>0.93 ± 0.74 (6, 1)</td><td>0.0012 ± 0.0011</td><td>2.67 ± 1.65 (6, 1)</td><td>0.0068 ± 0.0049</td></tr><tr><td>Central Europe (H)</td><td>0.00 ± 0.00 (13, 2)</td><td>0.0000 ± 0.0000</td><td>0.31 ± 0.34 (13, 2)</td><td>0.0011 ± 0.0014</td><td>0.73 ± 0.58 (12, 2)</td><td>0.0009 ± 0.0008</td><td><strong>6.10</strong> ± 3.30 (7, 2)</td><td><strong>0.0155</strong> ± 0.0096</td></tr><tr><td>Balkan and Italy (Gr, Bu, Ro, Cr, It)</td><td><strong>15.53</strong> ± 7.10 (35, 7)</td><td><strong>0.0295</strong> ± 0.0150</td><td>5.08 ± 2.53 (35, 7)</td><td>0.0181 ± 0.0100</td><td><strong>1.32</strong> ± 0.85 (23, 7)</td><td><strong>0.0016</strong> ± 0.0012</td><td><strong>6.00</strong> ± 4.58 (2, 2)</td><td><strong>0.0153</strong> ± 0.0166</td></tr><tr><td>Africa (Mo, Lib, Eg)</td><td>4.23 ± 2.19 (20, 4)</td><td>0.0080 ± 0.0046</td><td>2.58 ± 1.44 (20, 4)</td><td>0.0092 ± 0.0057</td><td>1.00 ± 0.71 (16, 4)</td><td>0.0012 ± 0.0010</td><td>4.62 ± 2.42 (13, 3)</td><td>0.0117 ± 0.0069</td></tr><tr><td>Middle East and central Asia (Isr, Ir, Ar, Tur, Tm, Ky)</td><td>10.23 ± 4.81 (29, 9)</td><td>0.0194 ± 0.0102</td><td>3.08 ± 1.65 (29, 9)</td><td>0.0110 ± 0.0065</td><td>1.08 ± 0.74 (20, 8)</td><td>0.0013 ± 0.0010</td><td>5.05 ± 2.79 (7, 4)</td><td>0.0129 ± 0.0082</td></tr><tr><td>Among geographic region variation, <i>p</i></td><td>36.6%, &lt; 0.00001</td><td></td><td>63.0%, &lt; 0.00001</td><td></td><td>13.0%, &lt; 0.00001</td><td></td><td>9.2%, 0.00198</td><td></td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Statistical parsimony network of nuclear 28S and EF-1α. See Fig. 1 for colors for regions. Generated using TCS 1.21⁸⁰.</figcaption></figure>

the Eastern/Egyptian clade and western Europe (France) for that of the Western clade. The primary center of genetic and species diversity of the main host *Medicago* is the Caucasus (north-western Iran and north-eastern Turkey)⁶⁵, which may also consist of the area of origin of *H. postica*.

The recent dispersal routes that include the north Mediterranean were highly supported. Anthropogenic factors may allow occasional dispersal of *H. postica* to Europe and North Africa with alfalfa traded for livestock feed (by 2,600 years ago⁶⁶). More recent international trade of alfalfa meal and pellets may continue to aid the weevil's opportunistic long-distance dispersal; France, Spain, and Italy are the major alfalfa exporters among *H. postica*'s native ranges⁶⁷.

**Conclusion.** While geographic isolation assisted continental diversification of the weevil *H. postica*, recent *Wolbachia* infection reduced diversity in a mitochondrial clade in the host weevil in western Europe. *Wolbachia*-infected males could be used as a control agent for the Incompatible Insect Technique on uninfected populations, however, the risk of heterosis in interclade crosses following accidental cure of *Wolbachia* must be assessed before application.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Bayesian consensus tree of *Wolbachia* strain based on *ftsZ*, *coxA*, and *hcpA*. Strain codes, if available, followed by host species and *Wolbachia* supergroups are shown. The *Wolbachia* strain, wHypera1, infecting *Hypera postica* in its native range, is shown in bold. Bayesian support values (posterior probabilities > 0.7) are shown near nodes. The outgroup is *Anaplasma marginale* (Alphaproteobacteria: Rickettsiales: Anaplasmataceae). Generated using MrBayes 3.2.6<sup>73</sup>. Host strains, *Wolbachia* isolates and GenBank accession numbers are listed in Supplementary Table S3.</figcaption></figure>

<table><caption>**Table 4.** Selective neutrality test results on mitochondrial and nuclear segments. *D*: Tajima's *D*; *D*\* and *F*\*: Fu and Li's *D*\* and *F*\*.</caption><thead><tr><th rowspan="2">Clade</th><th colspan="4">Mitochondria</th><th colspan="4">Nuclear</th></tr><tr><th>D (p)</th><th>D* (p)</th><th>F* (p)</th><th>n</th><th>D (p)</th><th>D* (p)</th><th>F* (p)</th><th>n</th></tr></thead><tbody><tr><td>Western</td><td>-2.389 (0)</td><td>-3.371 (0.0015)</td><td>-3.302 (0.0010)</td><td>23</td><td>-0.822 (0.152)</td><td>-0.547 (0.267)</td><td>-0.606 (0.240)</td><td>15</td></tr><tr><td>Egyptian/ Eastern</td><td>-1.510 (0.039)</td><td>-2.872 (0.0175)</td><td>-2.673 (0.018)</td><td>39</td><td>-1.491 (0.0255)</td><td>-2.007 (0.0125)</td><td>-1.959 (0.0165)</td><td>19</td></tr><tr><td>Intermediate</td><td>0.695 (0.763)</td><td>0.566 (0.707)</td><td>0.596 (0.719)</td><td>6</td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>**Table 5.** Test of positive selection on the *Hypera postica* phylogeny. Root: *Brachypera zoilus* and *H. miles*. W: Western. E: Egyptian/Eastern. Inf: infected by *Wolbachia*. Uninf: uninfected by *Wolbachia*. ΔlnL: difference in log likelihood (lnL) of each model from the model with a same single ω for W, intermediate, and E clades (i.e., the 'Root and W/intermediate/E' model as a reference model). Δdf: difference of each model in degree of freedom (df) from the model with a same single ω for W, intermediate, and E clades. The – symbol: the same ω value as the one on the left.</caption><thead><tr><th>Model (number of ω)</th><th>Root</th><th>E</th><th>Intermediate</th><th>W</th><th>InfW</th><th>ΔΔlnL</th><th>Δdf</th><th>ρ</th></tr></thead><tbody><tr><td>Root and W/intermediate/E (2)</td><td>0.0103</td><td>0.0578</td><td>–</td><td>–</td><td>–</td><td>0</td><td>0</td><td>1.000</td></tr><tr><td>Root, W and intermediate/E (3)</td><td>0.0103</td><td>0.0575</td><td>–</td><td>0.0590</td><td>–</td><td>0</td><td>1</td><td>1.000</td></tr><tr><td>Root, W/intermediate and E (3)</td><td>0.0104</td><td>0.0541</td><td>0.0640</td><td>–</td><td>–</td><td>1.50</td><td>1</td><td>0.221</td></tr><tr><td>Root, InfW and UninfW/intermediate/E (3)</td><td>0.0103</td><td>0.0571</td><td>–</td><td>–</td><td>0.1709</td><td>1.80</td><td>1</td><td>0.180</td></tr><tr><td>Root, InfW, UninfW and intermediate/E (4)</td><td>0.0103</td><td>0.0575</td><td>–</td><td>0.0554</td><td>0.1709</td><td>1.82</td><td>2</td><td>0.403</td></tr><tr><td>Root, InfW, UninfW/intermediate and E (4)</td><td>0.0104</td><td>0.0541</td><td>0.0620</td><td>–</td><td>0.1709</td><td>0.72</td><td>2</td><td>0.698</td></tr></tbody></table>

## Methods

**Sampling.** *Hypera postica* (n = 149) were obtained from 33 localities covering most of its native distribution range (Table 1). Adults were collected from cultivated and wild legume vegetation of *Medicago*<sup>47</sup>, *Trifolium* in Egypt, and *Vicia cracca* in Bulgaria, mostly during the latest decade. The samples were then stored in ethanol at 4 °C until DNA extraction.

**PCR and sequencing.** DNA was extracted from all specimens using a DNeasy Blood & Tissue kit (Qiagen, Tokyo, Japan). We amplified and sequenced two mitochondrial fragments, *COI-tRNA*<sup>Leu</sup>-COII and *Cyt b-tRNA*<sup>Ser</sup>-ND1 and two nuclear fragments, 28S and *EF-1α*. The primers used were C1-J-2797<sup>68</sup> and C2-N-3380<sup>24</sup> for the *COI-tRNA*<sup>Leu</sup>-COII fragment, CB-J-11545 and N1-N-11841<sup>68</sup> for *Cyt b-tRNA*<sup>Ser</sup>-ND1, 28S-01 and 28SR-01 for 28S<sup>69</sup>, and ef415F (5'-AACCAGAGAACATGCTCTTCTCG-3') and ef862R (5'-CTCAATTTTTTAAGT TTGTTCAATTTG-3') (this study) for *EF-1α*. PCRs were performed using GoTaq Green Master Mix (Promega, Tokyo, Japan). Cycling conditions for *COI-tRNA*<sup>Leu</sup>-COII amplification consisted of preheating at 95 °C for 2 min, followed by 38 cycles of denaturation at 94 °C for 50 s, annealing at 50 °C for 1 min, and an extension at 60 °C for 1 min. Amplification conditions for *Cyt b-tRNA*<sup>Ser</sup>-ND1 were identical, except that annealing was performed at 55 °C for 1 min. Cycling conditions for 28S amplification were as follows; preheating as above, followed by 35 cycles of denaturation at 94 °C for 30 s, annealing at 52 °C for 40 s, and an extension at 70 °C for 1 min. Those for *EF-1α* consisted of preheating as above, followed by 38 cycles of denaturation at 94 °C for 30 s, annealing at 51 °C for 40 s, and an extension at 68 °C for 1 min. Sequencing was carried out using a BigDye Terminator v3.1 Cycle Sequencing kit (Life Technologies/Applied Biosystems, Foster City, CA, USA) on a 3730 DNA Analyzer (Applied Biosystems).

**Wolbachia infection and phylogeny.** We used PCR to screen for possible *Wolbachia* infections. The *Wolbachia ftsZ* coding fragment was amplified using the primers *ftsZ-Z-f* and *ftsZ-Z-r*<sup>70</sup>. PCRs were performed using preheating as above, followed by 32 cycles at 94 °C for 40 s, 55 °C for 45 s, and 70 °C for 1 min. As a positive control, we used *Wolbachia*-infected *Callosobruchus chinensis*<sup>71</sup>. Blurred and extremely weak signals compared with the positive control were considered uninfected, which differs from a previous study<sup>47</sup>. *Wolbachia*-positive *H. postica* were further subjected to PCR and sequencing of the genes, *coxA* and *hcpA*, in addition to *ftsZ* for multilocus sequence typing of *Wolbachia*<sup>72</sup>.

For the phylogenetic reconstruction of *Wolbachia*, we used sequences from representative supergroups of *Wolbachia* (nr/nt database, Supplementary Table S1) with sequences of *Anaplasma marginale* as an outgroup. We used the GTR model, which was selected as the best fit model of nucleotide substitution by *MRBAYES*3.2.6<sup>73</sup>, based on the AICc, using *MRAIC.PL* 1.3.1<sup>74</sup>. The three gene segments were partitioned. Markov chain Monte Carlo (MCMC) simulations were performed for one million generations, with sampling conducted every 1,000 generations. The convergence of independent parallel runs was checked using *TRACER* 1.6<sup>75</sup>, and the first 25% of trees were discarded as burn-in.

**Body size.** After collecting specimens, the right elytron lengths of the samples were measured to the precision of 0.01 mm with a microscope (VH-5500, Keyence, Osaka, Japan). The sex of the samples was determined by both external and genital morphology. The effect of sex, clades (Egyptian/Eastern, Western, and intermediate), and *Wolbachia* infection on elytral lengths were tested by nonparametric Wilcoxon/Kruskal-Wallis signed rank tests. Posthoc multiple comparison was performed on the significant factor using the Steel-Dwass test. JMP 14.2.0 was used for statistical analyses.

**Selective neutrality test and positive selection test.** Selective neutrality was tested in each clade with Tajima's D, Fu and Li's D* and F*, using DNASP 6.12.03<sup>78</sup>. P values were derived from coalescent simulations with 2,000 replications. For the coalescent simulations for nuclear segments, an intermediate recombination rate was assumed. We used all mitochondrial (808 bp) or nuclear sequences (1,193 bp) of two individuals sampled per clade from each population to avoid sample size bias between populations.
<figure><img src="image_8.png" /><figcaption>a</figcaption></figure>

<figure><img src="image_9.png" /></figure>
◄**Figure 5.** (a) The maximum clade credibility (MCC) tree of *Hypera postica* resulting from a discrete Bayesian phylogeographic analysis, based on one sample per clade per country. The outgroup is *H. miles*. The nodes and branches are colored according to the most probable region (left panel; see Fig. 1 for colors indicating geographic regions and Table 1 for the country codes) and *Wolbachia* infection state (right panel; red line: infected). The numbers on the nodes indicate posterior probabilities of the most likely ancestral states [geographic region (left panel) or *Wolbachia* infection state (right panel)]. Generated using BEAST 1.10.4⁹⁰ and visualized using FigTree 1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/). (b) Map of colonization routes of *H. postica*. Thick lines: supported by BF (Bayes factor) > 7.0. Generated using SpreaD3 v0.9.7.1rc⁹² and visualized using Google Earth 7.3.0 (https://www.google.com/earth/).

The equal nonsynonymous/synonymous substitution rate ratio (dN/dS ratio, ω) between infected and unin-
fected clades and between the two major clades was tested with a phylogenetic analysis using the maximum
likelihood method (likelihood ratio test) employed by the codeml program in PAML 4.9i⁷⁹. The models with
three or more different ω for each branch were compared with a reference (basal) model with two different ω
(one for the root and the other for both the Western and Egyptian/Eastern clades). We concatenated all the open
reading frames (protein-coding fragments) and removed potential stop codons (leading to 215 codons) of the
mitochondrial sequences of two individuals randomly sampled per clade from each population to avoid sample
size bias between populations. Codons for invertebrate mitochondria were used.

**Haplotype networks and diversity.** All mitochondrial sequences were assembled using SEQUENCHER 5.0 (Gene Codes Corp, Ann Arbor, MI, USA), and we checked for the presence of pseudogenes using commonly employed methods⁸⁰,⁸¹. Statistical parsimony networks were reconstructed based on mitochondrial and nuclear fragments using TCS 1.21⁸², in which we allowed connections between haplotypes of 20 steps for mitochondrial genes and 95% for the nuclear genes, to elucidate the maximal divergence observed among haplotypes. Nucleotide diversity (average over loci) and π (mean number of pairwise differences)⁸³ were estimated for each geographic area, using ARLEQUIN 3.5.2.2⁸⁴.

**Geographic history.** *Isolation by distance.* We assessed spatial mitochondrial differentiation by testing for isolation by distance (IBD)⁸⁵. With a sweep (e.g., by *Wolbachia*) followed by rapid spread or frequent an-
thropogenic long-distance dispersal events, the IBD correlation is predicted to be weak at most. We tested if, as
predicted by IBD, pairwise geographic distances and pairwise genetic differences were positively correlated us-
ing a one-tailed Mantel test⁸⁶ based on 2,000 permutations with the ISOLDE program implemented in Genepop
4.2⁸⁷. For pairwise genetic differences, we employed corrected average pairwise differences between populations
X and Y, [π<sub>XY</sub> - (π<sub>X</sub> + π<sub>Y</sub>)/2]⁸⁸ and their *p* values were derived using ARLEQUIN.

*Phylogeography.* We estimated the historical dispersal patterns of *H. postica*, using a Bayesian discrete phylo-
geographic approach⁸⁹ with a Bayesian skyride framework implemented in the software package BEAST 1.10.4⁹⁰.
We used two mitochondrial segments (808 bp). To avoid sample size bias, we selected only one individual per
clade from a given locality but excluded intermediate mitotypes, which reduced the data set to 34 individuals,
with *H. miles* as an outgroup. We used default settings, applied the same molecular evolution model as presented
above, and used an uncorrelated relaxed clock model assuming lognormal rate distribution⁹¹. We assigned each
sequence to one of the seven geographic regions, and the symmetric exchanges between the geographic regions
throughout the entire phylogeny were modeled with the Bayesian stochastic search variable selection (BSSVS).
MCMC runs were performed for 50 million generations, sampling one tree every 25,000 generations. After
confirming the stationarity of parameter estimates using TRACER, the first 40% of trees were discarded as burn-
in, and maximum clade credibility (MCC) tree was built using TREEANNOTATORv1.10.4. As each node in each
MCMC sample is annotated with a geographic region and *Wolbachia* infection status, we assessed the certainty
of the geographic reconstruction by looking at the distribution of node states across the MCMC using FigTree
v1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/). Bayes factor (BF) values for exchange rates between each pair-
wise regions were retrieved from the log file from the BSSVS analysis using SPREAD3 v0.9.7.1rc⁹².

**Data availability**

GenBank accessions KX372573-372592, 372596-372619 and MW 393902-393920 for COI-tRNA<sup>Leu</sup>-COII, KX372620-372639, 372643-372666 and MW393921-393939 for Cyt b-tRNA<sup>Ser-ND1</sup>, KX372667 and MW383443-383462 for 28S, MW389094-389112 and 392102 for EF-1α, and MW389113-389118 for ftsZ, coxA and hcpA.

Received: 30 August 2020; Accepted: 15 April 2021

Published online: 06 May 2021

**References**

1. Bonizzoni, M. et al. On the origins of medfly invasion and expansion in Australia. *Mol. Ecol.* **13**, 3845–3855 (2004).

2. Tuda, M., Kagoshima, K., Toquenaga, Y. & Arnqvist, G. Global genetic differentiation in a cosmopolitan pest of stored beans: Effects of geography, host-plant usage and anthropogenic factors. *PLoS ONE* **9**, e106268 (2014).

3. Karsten, M., van Vuuren, B. J., Addison, P. & Terblanche, J. S. Deconstructing intercontinental invasion pathway hypotheses of the Mediterranean fruit fly (Ceratitis capitata) using a Bayesian inference approach: Are port interceptions and quarantine protocols successfully preventing new invasions?. *Divers. Distrib.* **21**, 813–825 (2015).
4. Rodriguez M. S. et al. Out of the forest: past and present range expansion of a parthenogenetic weevil pest, or how to colonize the world successfully. *Ecol. Evol.* **6**, 5431–5445 (2016).

5. Kébé, K. et al. Global phylogeography of the insect pest *Callosobruchus maculatus* (Coleoptera: Bruchinae) relates to the history of its main host *Vigna unguiculata*. *J. Biogeogr.* **44**, 2515–2526 (2017).

6. Lombaert, E. et al. Colonization history of the western corn rootworm (Diabrotica virgifera virgifera) in North America: insights from random forest ABC using microsatellite data. *Biol. Invasions* **20**, 665–677 (2018).

7. Tuda, M., Ronn, J., Buranapanichpan, S., Wasano, N. & Arnqvist, G. Evolutionary diversification of the bean beetle genus *Callosobruchus* (Coleoptera: Bruchidae): Traits associated with stored-product pest status. *Mol. Ecol.* **15**, 3541–3551 (2006).

8. Wei, S. J. et al. Population genetic structure and approximate Bayesian computation analyses reveal the southern origin and north-
ward dispersal of the oriental fruit moth *Grapholita molesta* (Lepidoptera: Tortricidae) in its native range. *Mol. Ecol.* **24**, 4094–4111
(2015).

9. Takano, S. et al. Unique clade of alphaproteobacterial endosymbionts induces complete cytoplasmic incompatibility in the coconut
beetle. *Proc. Natl. Acad. Sci. USA* **114**, 6110–6115 (2017).

10. Radcliffe, E. B. & Flanders, K. L. Biological control of alfalfa weevil in North America. *Integr. Pest Manag.* Rev. **3**, 225–242 (1998).

11. Kuwata, R., Tokuda, M., Yamaguchi, D. & Yukawa, J. Coexistence of two mitochondrial DNA haplotypes in Japanese populations of Hypera postica (Col., Curculionidae). *J. Appl. Entomol.* **129**, 191–197 (2005).

12. Skuhrovec, J. Host plants of weevils of the genus Hypera (Coleoptera: Curculionidae) occurring in the Czech Republic. *Klapalekiana* **41**, 215–255 (2005).

13. Wood, K. A., Armbrust, E. J., Bartell, D. P. & Irwin, B. J. The literature of arthropods associated with alfalfa. V. A bibliography of the alfalfa weevil, *Hypera postica* (Gyllenhal), and the Egyptian alfalfa weevil, *Hypera brunneipennis* (Boheman) (Coleoptera: Curculionidae). Illinois Agricultural Experimental Station, Special Publication, 54 (1978).

14. Kimura, H., Okumura, M. & Yoshida, T. Emergence of and recent damage by the alfalfa weevil. Shokubutsu Boeki (Plant Protec-
tion) **42**, 498–501 (in Japanese) (1988).

15. CAB International crop protection compendium. CAB International. http://www.cabicompendium.org/cpc/home.asp (2013).

16. Titus, E. G. On the life history of the alfalfa leaf-weevil. *J. Econ. Entomol.* **3**, 459–470 (1910).

17. Wehrle, L. P. The discovery of an alfalfa weevil (*Hypera brunneipennis* Boheman) in Arizona. *J. Econ. Entomol.* **33**, 119–121 (1940).

18. Poos, F. W. & Bissell, T. L. The alfalfa weevil in Maryland. *J. Econ. Entomol.* **46**, 178–179 (1953).

19. Volker, K. C. & Simpson, R. G. Behavior of alfalfa weevil larvae affecting the establishment of *Tetrastichus incertus* in Colorado. *Environ. Entomol.* **4**, 742–744 (1975).

20. Salt, G. & van den Bosch, R. The defense reactions of three species of *Hypera* (Coleoptera, Curculionidae) to an Ichneumon wasp. *J. Invertebr. Pathol.* **9**, 164–177 (1967).

21. Maund, C. M. & Hsiao, T. H. Differential encapsulation of two *Bathyplectes* parasitoids among alfalfa weevil strains, *Hypera postica* (Gyllenhal). *Can. Entomol.* **123**, 197–203 (1991).

22. Hsiao, T. H. Studies of interactions between alfalfa weevil strains, *Wolbachia* endosymbionts and parasitoids. In The ecology of agricultural pests: biochemical approaches (eds, Symondson, W. O. C. & Liddell, J. E.). 57–71 (Chapman & Hall, 1996).

23. Hsiao, T. H. & Stutz, J. M. Discrimination of alfalfa weevil strains by allozyme analysis. *Entomol. Exp. Appl.* **37**, 113–121 (1985).

24. Erney, S. J., Pruess, K. P., Danielson, S. D. & Powers, T. O. Molecular differentiation of alfalfa weevil strains (Coleoptera: Curculionidae). *Ann. Entomol. Soc. Am.* **89**, 804–811 (1996).

25. Böttger, J. A. A. Phylogenetic analysis of the alfalfa weevil complex (Coleoptera: Curculionidae) in North America. *J. Econ. Entomol.* **106**, 426–436 (2013).

26. Iwase, S., Nakahira, K., Tuda, M., Kagoshima, K. & Takagi, M. Host-plant dependent population genetics of the invading weevil *Hypera postica*. *Bull. Entomol. Res.* **105**, 92–100 (2015).

27. White, C. E., Armbrust, E. J. & Ashley, J. Cross-mating studies of eastern and western strains of alfalfa weevil. *J. Econ. Entomol.* **65**, 85–89 (1972).

28. Iwase, S. & Tani, S. New haplotype and inter-strain reproductive compatibility of *Wolbachia*-uninfected alfalfa weevil, *Hypera postica* (Coleoptera: Curculionidae), in Japan. *Entomol. Sci.* **19**, 72–76 (2016).

29. Werren, J. H. Biology of *Wolbachia*. *Annu. Rev. Entomol.* **42**, 587–609 (1997).

30. LePage, D. P. et al. Prophage WO genes recapitulate and enhance *Wolbachia*-induced cytoplasmic incompatibility. *Nature* **543**, 243–247 (2017).

31. Bailly-Bechet, M. et al. How long does *Wolbachia* remain on board? *Mol. Biol. Evol.* **34**, 1183–1193 (2017).

32. Hale, L. R. & Hoffmann, A. A. Mitochondrial DNA polymorphism and cytoplasmic incompatibility in natural populations of *Drosophila simulans*. *Evolution* **44**, 1383–1386 (1990).

33. Ballard, J. W. O. & Kreitman, M. Unravelling selection in the mitochondrial genome of *Drosophila*. *Genetics* **138**, 757–772 (1994).

34. Johnstone, R. A. & Hurst, G. D. D. Maternally inherited male-killing microorganisms may confound interpretation of mitochondrial DNA variability. *Biol. J. Linn. Soc.* **58**, 453–470 (1996).

35. Jiggins, F. M. Male-killing *Wolbachia* and mitochondrial DNA: selective sweeps, hybrid introgression and parasite population dynamics. *Genetics* **164**, 5–12 (2003).

36. Werren, J. H., Baldo, L. & Clark, M. E. *Wolbachia*: master manipulators of invertebrate biology. *Nat. Rev. Microbiol.* **6**, 741–751 (2008).

37. Shoemaker, D. D., Dyer, K. A., Ahrens, M., McAbee, K. & Jaenike, J. Decreased diversity but increased substitution rate in host mtDNA as a consequence of *Wolbachia* endosymbiont infection. *Genetics* **168**, 2049–2058 (2004).

38. Cariou, M., Duret, L. & Charlat, S. The global impact of *Wolbachia* on mitochondrial diversity and evolution. *J. Evol. Biol.* **30**, 2204–2210 (2017).

39. Teixeira, L., Ferreira, A. & Ashburner, M. The bacterial symbiont *Wolbachia* induces resistance to RNA viral infections in *Drosophila melanogaster*. *PLoS Biol.* **6**, 2753–2763 (2008).

40. Brownlie, J. C. et al. Evidence for metabolic provisioning by a common invertebrate endosymbiont, *Wolbachia pipiensis*, during periods of nutritional stress. *PLoS Pathog.* **5**, e1000368 (2009).

41. Rand, D. M., Haney, R. A. & Fry, A. J. Cytonuclear coevolution: the genomics of cooperation. *TRENDS Ecol. Evol.* **19**, 645–653 (2004).

42. Arnqvist, G. et al. The genetic architecture of metabolic rate: environment specific epistasis between mitochondrial and nuclear genes in an insect. *Evolution* **64**, 3354–3363 (2010).

43. Blickenstaff, C. C. Partial intersterility of eastern and western US strains of the alfalfa weevil. *Ann. Entomol. Soc. Am.* **58**, 523–526 (1965).

44. Hsiao, T. H. & Hsiao, C. Hybridization and cytoplasmic incompatibility among alfalfa weevil strains. *Entomol. Exp. Appl.* **37**, 155–159 (1985).

45. Laven, H. Eradication of *Culex pipiens fatigans* through cytoplasmic incompatibility. *Nature* **216**, 383–384 (1967).

46. Iwase, S. et al. Dynamics of infection with *Wolbachia* in *Hypera postica* (Coleoptera: Curculionidae) during invasion and establishment. *Biol. Invasions* **17**, 3639–3648 (2015).

47. Sanaei, E. et al. Global genetic diversity, lineage distribution and *Wolbachia* infection of the alfalfa weevil *Hypera postica* (Coleoptera: Curculionidae). *Ecol. Evol.* **9**, 9546–9563 (2019).
48. Ros, V. I. D., Fleming, V. M., Feil, E. J. & Breeuwer, J. A. J. How diverse is the genus *Wolbachia*? Multiple-gene sequencing reveals a putatively new *Wolbachia* supergroup recovered from spider mites (Acari: Tetranychidae). *Appl. Environ. Microbiol.* **75**, 1036–1043 (2009).

49. Avise, J. C. *Phylogeography: The history and formation of species* (Harvard University Press, 2000).

50. Narita, S., Nomura, M., Kato, Y. & Fukatsu, T. Genetic structure of sibling butterfly species affected by *Wolbachia* infection sweep: evolutionary and biogeographical implications. *Mol. Ecol.* **15**, 1095–1108 (2006).

51. Raychoudhury, R. et al. Phylogeography of *Nasonia vitripennis* (Hymenoptera) indicates a mitochondrial–*Wolbachia* sweep in North America. *Heredity* **104**, 318–326 (2010).

52. Jäckel, R., Mora, D. & Dobler, S. Evidence for selective sweeps by *Wolbachia* infections: phylogeny of *Altica* leaf beetles and their reproductive parasites. *Mol. Ecol.* **22**, 4241–4255 (2013).

53. Jiang, W. et al. *Wolbachia* infection status and genetic structure in natural populations of *Polytremis nascens* (Lepidoptera: Hesperiidae). *Infect. Genet. Evol.* **27**, 202–211 (2014).

54. Jansen, V. A. A., Turelli, M. & Godfray, H. C. J. Stochastic spread of *Wolbachia*. *Proc. R. Soc. Lond. B Biol. Sci.* **275**, 2769–2776 (2008).

55. Clancy, D. J. & Hoffmann, A. A. Environmental effects on cytoplasmic incompatibility and bacterial load in *Wolbachia*-infected *Drosophila simulans*. *Entomol. Exp. Appl.* **86**, 13–24 (1998).

56. Bordenstein, S. R. & Bordenstein, S. R. Temperature affects the tripartite interactions between bacteriophage WO, *Wolbachia*, and cytoplasmic incompatibility. *PLoS ONE* **6**, e29106 (2011).

57. Kamo, T. et al. Limited distribution of natural cyanamide in higher plants: Occurrence in *Vicia villosa* subsp varia, *V. cracca*, and *Robinia pseudo-acacia*. *Phytochemistry* **69**, 1166–1172 (2008).

58. Megias, C., Cortes-Giraldo, I., Giron-Calle, J., Alaiz, M. & Vioque, J. Free amino acids, including canavanine, in the seeds from 32 *Vicia* species belonging to subgenus *Vicilla*. *Biocatal. Agric. Biotechnol.* **8**, 126–129 (2016).

59. Rosenthal, G. A. & Dahlman, D. L. Incorporation of L-canavanine into proteins and the expression of its antimetabolic effects. *J. Agric. Food Chem.* **39**, 987–990 (1991).

60. Kamo, T., Tokuoka, Y. & Miyazaki, M. Quantification of canavanine, 2-aminoethanol, and cyanamide in *Aphis craccivora* and its host plants, *Robinia pseudoacacia* and *Vicia angustifolia*: Effects of these compounds on larval survivorship of *Harmonia axyridis*. *J. Chem. Ecol.* **38**, 1552–1560 (2012).

61. Hewitt, G. M. Post-glacial re-colonization of European biota. *Biol. J. Linn. Soc.* **68**, 87–112 (1999).

62. Schmitt, T. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. *Front. Zool.* **4**, 11 (2007).

63. Taberlet, P., Fumagalli, L., Wust-Saucy, A. G. & Cosson, J. F. Comparative phylogeography and postglacial colonization routes in Europe. *Mol. Ecol.* **7**, 453–464 (1998).

64. Jordal, B. H. & Kambestad, M. DNA barcoding of bark and ambrosia beetles reveals excessive NUMTs and consistent east-west divergence across Palearctic forests. *Mol. Ecol. Resour.* **14**, 7–17 (2013).

65. Quiros, C. F. & Bauchan, G. R. The genus *Medicago* and the origin of the *Medicago sativa* complex. In *Alfalfa and alfalfa improvement* (eds, Hanson, A. A., Barnes, D. K. & Hill, R. R.). 93–124 (American Society of Agronomy, Crop Science Society of America, Soil Science Society of America, 1988).

66. Small, E. *Alfalfa and Relatives: Evolution and Classification of Medicago* (NRC Research Press, 2011).

67. FAO Statistics Division. FAOSTAT: Crops and livestock products. http://www.fao.org/faostat/en/#data/TP (2017).

68. Simon, C. A. et al. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. *Ann. Entomol. Soc. Am.* **87**, 651–701 (1994).

69. Kim, C. G. et al. Pattern of morphological diversification in the *Leptocarabus* ground beetles (Coleoptera: Carabidae) as deduced from mitochondrial ND5 gene and nuclear 28S rDNA sequences. *Mol. Biol. Evol.* **17**, 137–145 (2000).

70. Holden, P. R., Brookfield, J. F. Y. & Jones, P. Cloning and characterization of an ftsZ homologue from a bacterial symbiont of *Drosophila melanogaster*. *Mol. Gen. Genet.* **240**, 213–220 (1993).

71. Kondo, N. I. et al. *Wolbachia* infections in world populations of bean beetles (Coleoptera: Chrysomelidae: Bruchinae) infesting cultivated and wild legumes. *Zool. Sci.* **28**, 501–508 (2011).

72. Baldo, L. et al. Multilocus sequence typing system for the endosymbiont *Wolbachia pipentis*. *Appl. Environ. Microbiol.* **72**, 7098–7110 (2006).

73. Huelsenbeck, J. P. & Ronquist, F. MrBayes: Bayesian inference of phylogeny. *Biometrics* **17**, 754–755 (2001).

74. Nylander, J. A. A. MrAIC.pl. Program distributed by the author. Uppsala: Evolutionary Biology Centre, Uppsala University (2004).

75. Rambaut, A. & Drummond, A. J. Tracer v1.5, http://beast.bio.ed.ac.uk/ (2009).

76. Tajima, F. Statistical methods to test for nucleotide mutation hypothesis by DNA polymorphism. *Genetics* **123**, 585–595 (1989).

77. Fu, Y.-X. & Li, W.-H. Statistical tests of neutrality of mutations. *Genetics* **133**, 693–709 (1993).

78. Librado, P. & Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. *Bioinformatics* **25**, 1451–1452 (2009).

79. Yang, Z. PAML 4: A program package for phylogenetic analysis by maximum likelihood. *Mol. Biol. Evol.* **24**, 1586–1591 (2007).

80. Song, H., Buhay, J. E., Whiting, M. F. & Crandall, K. A. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. *Proc. Natl. Acad. Sci. USA* **105**, 13486–13491 (2008).

81. Haran, J., Koutroumpa, F., Magnoux, E., Roques, A. & Roux, G. Ghost mtDNA haplotypes generated by fortuitous NUMTs can deeply disturb infra-specific genetic diversity and phylogeographic pattern. *J. Zoolog. Syst. Evol. Res.* **53**, 109–115 (2015).

82. Clement, M., Snell, Q., Walker, P., Posada, D. & Crandall, K. TCS: Estimating gene genealogies. *Parallel Distrib. Proces. Symp. Int. Proc.* **2**, 184 (2002).

83. Tajima, F. Evolutionary relationship of DNA sequences in finite populations. *Genetics* **105**, 437–460 (1983).

84. Excoffier, L. & Lischer, H. E. L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. *Mol. Ecol. Resour.* **10**, 564–567 (2010).

85. Wright, S. Isolation by distance. *Genetics* **28**, 114–138 (1943).

86. Mantel, N. The detection of disease clustering and a generalized regression approach. *Cancer Res.* **27**, 209–220 (1967).

87. Raymond, M. & Rousset, F. Genepop (version 1.2): Population genetics software for exact tests and ecumenicism. *J. Hered.* **86**, 248–249 (1995).

88. Nei, M. & Li, W. H. Mathematical model for studying genetic variation in terms of restriction endonucleases. *Proc. Natl. Acad. Sci. USA* **76**, 5269–5273 (1979).

89. Lemey, P., Rambaut, A., Drummond, A. J. & Suchard, M. A. Bayesian phylogeography finds its roots. *PLoS Comput. Biol.* **5**, e1000520 (2009).

90. Suchard, M. A. et al. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. *Virus Evol.* **4**, vey016 (2018).

91. Drummond, A. J., Suchard, M. A., Xie, D. & Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. *Mol. Biol. Evol.* **29**, 1969–1973 (2012).

92. Bielejec, F. et al. SpreaD3: Interactive visualization of spatiotemporal history and trait evolutionary processes. *Mol. Biol. Evol.* **33**, 2167–2169 (2016).
**Acknowledgements**

We thank Łukasz Kajtoch for information about the collection site in Poland, Tamás Németh and Zoltán György for guidance at the Hungarian Natural History Museum, Kumiko Kagoshima for assistance with PCR, and Anahí Espíndola for advice on phylogeographic methods.

**Author contributions**

M.T., S.I., J.H. and J.S. conceived the idea, M.T., S.I. and J.H. conducted molecular experiments, J.S., A.P., O.M., N.T., E.S., A.H.E. and K.M. performed morphological identification, M.T. and K.K. analyzed the data and M.T., J.H. and K.K. led the writing.

**Funding**

Funding was provided by Japan Society for the Promotion of Science (Grant Nos. JP23405008, JP25430194, JP26304016, JP17H04612, JP18H02207 and JP19K06840), Kyushu University (Interdisciplinary Programs in Education and Projects in Research Development (25412)) and Ministerstvo Zemědělství (Grant No. RO0418).

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-88770-y.

**Correspondence** and requests for materials should be addressed to M.T.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021