Phylogeny of Mesitiinae (Hymenoptera: Bethylidae): assessing their classification, character evolution and diversification

): assessing their classification, character evolution and diversification. Arthropod


Introduction
The aculeate family Bethylidae includes nearly 3,000 species of parasitoid wasps, representing one of the most diverse lineages of Chrysidoidea (Finnamore and Brothers 1993;Azevedo et al. 2018). Among bethylids, the Mesitiinae have their diversity represented by 188 species distributed into 18 genera, recorded in tropical environments of the Afrotropical, Australian, Oriental, and Palearctic zoogeographical regions. Although many bethylid fossils have been described in recent years (e.g., Azevedo and Azar 2012;Ramos et al. 2014;Colombo et al. 2021, see also Azevedo et al. 2018 for a review), no fossil species of Mesitiinae were ever described. Known host records indicate that representatives of the subfamily are parasitoids of leaf-beetle larvae of the subfamilies Clytrinae and Cryptocephalinae (Coleoptera: Chrysomelidae), which reside in close-fitting cases built of fecal material (Argaman 2003). During oviposition, these wasps may exhibit predatory habits, since the female carries the paralyzed beetle immature into preexisting soil crevices with the mandibles (Nagy 1969;Argaman 2003).

Collections
The specimens used in this study were borrowed from the following collections, with curators in parentheses:

Illustrations
The images were obtained using a Leica MZ80 Stereomicroscope attached to a Leica DFC 495 video camera and captured with LEICA LAS (Leica Application Suite V3.6.0) by Leica Microsystems (Switzerland), using a dome illumination system described by Kawada and Buffington (2016), and combined using HELICON FO-CUS (version 4.2.9). Illustrations and plates were edited for adjustments (e.g., levels, shadows/highlights).

Terminology
The terms applied to the structures follow Lanes et al. (2020) and Barbosa and Azevedo (2011), integument terminology follows Harris (1979). Abbreviation: VOL = vertex-ocular line in dorsal view.

Taxon sampling
The ingroup is composed by males of 61 species (Table 1). The species analyzed correspond to over a third of the 182 species in Mesitiinae. Species selection aimed to cover the maximum possible morphological diversity in each genus to facilitate possible taxonomic decisions.
Character definition was based on males for three main reasons: (1) the current classification by Argaman (2003) was based on males; (2) the male hypopygium and genitalia offer a range of characters not available for females; and (3) lack of conspicuous sexual dimorphism between male and female. Except for Australomesitius Barbosa & Azevedo, 2016, known only from the female, 17 out of the 18 genera currently included in the subfamily were sampled. The outgroup includes representatives of all extant subfamilies of Bethylidae (Table 1). The Bethylinae Bethylus cephalotes (Förster, 1860) was used for rooting the tree. Table 1 should be included here associated with subchapter 2.4, landscape, and maximally page-filling.

Characters
A total of 112 characters (Appendix 1) were analyzed. Many of them were taken from descriptions in Móczár  (1970,1971), Nagy (1969Nagy ( , 1972, and Argaman (2003); additionally, new characters are proposed here for the first time.

Character matrix
The character matrix (Table S1) was produced using DELTA software (Dallwitz et al. 1993). All characters were treated as unordered. Inapplicable characters were coded as "?".

Parsimony analyses
The searches for the most parsimonious trees were carried out in TNT version 1.5 (Goloboff et al. 2016, using the Ratchet, Sectorial Searches and Tree-Fusing searching strategies (Goloboff 1999, Nixon, 1999. Parameters were as follows: collapsing rules selected for TBR; random seed set to 0; Sectorial Search in default mode; 200 iterations of Ratchet; 20 cycles for Drift; 10 rounds for Tree Fusing. It has been argued that results based on characters properly weighted are to be preferred over those with all characters equally weighted (Farris 1969, Goloboff 1993, Goloboff et al. 2008a. Implied weighting is the most widely used method for attributing different weights during tree search, as it is independent of previous analyses and weighting schemes unlike, for example, successive weighting (e.g., Farris 1969). The weighting against homoplasy under implied weighting is related to a constant k -the lower the value of k, the higher the strength against homoplasy Goloboff et al. (2008b). Here, we used the TNT script setk.run, written by Salvador Arias (Instituto Miguel Lillo, San Miguel de Tucuman, Argentina), to calculate the value of k. The script returned a value of k = 11.674805 for our data set.

Bayesian analyses
Bayesian analyses were conducted in MRBAYES 3.2.7 (Ronquist et al. 2012). We used the Mk model to morphological data, with correction for ascertainment bias (lset coding = variable), since autapomorphic characters were included. We first conducted an analysis without partitioning the original matrix, accounting for among-character rate heterogeneity using a discrete Gamma distribution with four rate categories (lset rates = gamma) and the prior on branch lengths described by an exponential distribution with scale parameter = 10 (prset brlenspr = Unconstrained:Exp(10)). We also conducted a similar analysis partitioning characters according to their degree of homoplasy. For this purpose, we retrieved homoplasy scores from implied-weighting analyses in TNT (see above). These values, derived from Goloboff's measure of homoplasy, are normalized between 0 and 1, with the lowest value representing no homoplasy (Goloboff et al. 2008b). Branch lengths were maintained linked among partitions, and site-specific rates within partitions were not considered, as suggested by Rosa et al. (2019). MCMC analyses ran for 5,000,000 generations, sampling every 1,000, with four chains, and two independent runs. Convergence was assessed with Tracer 1.6 (Rambaut et al. 2018). Trees shown are majority-rule consensus trees (Contype = Allcompat).

Bayesian analyses
Bayesian analyses largely corroborated the backbone of the relationships retrieved in parsimony (Figs S2, S3). Results from unpartitioned and partitioned analyses dif-fered. In both analyses the genus Anaylax was not recovered as monophyletic, with Anaylax simplicitus being recovered as a distinct lineage relative to other species of the genus included in the present account. The unpartitioned analysis recovered Incertosulcus krombeini as sister group to Metrionotus, while in the analysis using partitioning by homoplasy score it was recovered nested within Metrio notus. However, the posterior probability of the clade Metrionotus + I. krombeini was very low in both cases (i.e., < 0.42). Both analyses also recovered different taxa as the sister group to all other mesitiine lineages: Moczariella centenaria in the unpartitioned analysis and the genus Bradepyris in the partitioned analysis.

Taxonomic Accounts
The interpretation of topologies obtained allowed us to propose 17 nomenclatural changes: two new genera, one genus synonymy, three revalidations in species status, and 11 new specific combinations (Figs 5-6; Table 2, S2). Because a recent review for diagnostic characteristics for Bethylidae genera was published by in Azevedo et al. (2018), we describe here only the diagnostic characteristics for the new genera proposed and the changes for the genera reinterpreted in this work. Incertosulcus krombeini Móczár, 1970 Incertosulcus krombeini Móczár, 1970 Brachymesitius krombeini (Móczár, 1970) comb. nov.

New genera
Etymology. The name Hadesmesitius, masculine, is a combination of the "Hades", the Greek mythology god that has a forked weapon with the same shape of the hypopygium in this genus, which is diagnostic for the group, and the name "Mesitius", the type genus of Mesitiinae. Description. Wings: hyaline. Head: As long as wide; malar space as long as VOL, convergent; clypeus with median lobe rounded, median clypeal carina inclined; antenna with pubescence sparse and short; pedicel cylindrical, first flagellomere as long as pedicel, flagellomeres short; eye very small, without pubescence; frons foveolate, with frontal carina; ocelli very small; anterior ocellus crossing supra-ocular line; dorsal and ventral half of occipital carina low. Pronotum: Dorsal pronotal area shorter than wide, foveolate, with humeral angle rounded, side straight, anterior margin outcurved, posterior margin incurved, median pronotal line absent; mesoscutum coriaceous, median mesonotal sulcus absent, notaulus present and narrow; mesoscutellum not touching the metapectal-propodeal disc; metapectal-propodeal disc as long as its half width, metapostnotal median carina complete, with longitudinal ridge between metapostnotal median carina and metapostnotal-propodeal carina, posterior propodeal projection very short and thick; spiracle shape elliptical; propodeal declivity areolate, with median and lateral carinae; lateral surface of metapectal-propodeal complex areolate, without carinae. Wings: Hind wing with four distal hamuli. Metasoma dorsal and ventral region of terga III-VI polished; hypopygium bilobate, spiculum short, with lobate and short branch, wider than long, lateral margins convergent, with lateral anterior projection.
Etymology. The name Brachymesitius, masculine, is a combination of the names "brachy", from the Greek "short", and refers to the reduced size of structures, such as eye size, flagellomeres, ocelli, length of dorsal pronotal area, posterior propodeal projection, and hypopygium, which are diagnostic for the group, and the name "Mesitius", the type genus of Mesitiinae. Móczár, 1970, now Brachymesitius krombeini (Móczár, 1970 stat. rev. et comb. nov., removed from the synonym of Parvoculus indicus Kieffer, 1905.

Pycnomesitius Móczár, 1971
Remarks. This genus is characterized by the head as long as wide, the anteromesoscutum without median mesonotal line, the posterior propodeal projection short, and the metasomal tergum II densely punctured (Azevedo et al. 2018). However, it was found to be polyphyletic . In order to solve this problem, Sulcomesitius wahisi Móczár, 1984 is herein transferred from Sulcomesitius to Pycnomesitius, P. wahisi (Móczár, 1984) comb. nov.

Sulcomesitius Móczár, 1970
Remarks. This genus is characterized by the malar space as long as vertex-ocular line, convergent anteriorly, in front view, the anteromesoscutum with median mesonotal line well impressed, the forewing with nebulous Cu and A veins, the hypopygium with branches lobate and long, and the male genitalia with dorsal paramere S-shaped, ventral paramere narrower than dorsal (Azevedo et al. 2018). However, it was found to be polyphyletic (Figs 4-6). The same applies to Sulcomesitius nepalensis, which was recovered as sister group to a clade formed by four species of Metrionotus and Brachymesitius krombeini in the partitioned Bayesian analysis and recovered nested within species of Metrionotus in the unpartitioned analysis; and Gerbekas laoensis, which was recovered as sister group to a clade formed by four species of Sulcomesitius in the partitioned Bayesian analysis and as single clade in the unpartitioned analysis. In both cases, the support for such groupings is low, indicated by posterior probability values below 0.4. To solve this problem, Sulcomesitius nepalensis Móczár, 1986 is herein transferred to Metrio notus, M. nepalensis (Móczár, 1986) comb. nov. and Gerbekas laoensis Móczár, 1975 is herein transferred from Gerbekas to Sulcomesitius, S. laoensis (Móczár, 1975) comb. nov.
The topologies obtained allowed the identification of morphological characters which potentially played important roles during the diversification of Mesitiinae, including sculpture of frons (#24); sculpture of dorsal pronotal area (#33); presence of median pronotal line (#36); presence of posterior propodeal projection (#63); length of hypopygium (#94); shape of posterior hypopygeal margin (#96); and length of hypopygium branches (#98). These characters are unique to the subfamily and allow us to hypothesize about their evolution. These hypotheses are largely based on convergent characteristics shared between Mesitiinae and Chrysidinae (Chrysididae) (Argaman, 2003).

Integumental adaptations
From the 19 genera proposed for Mesitiinae, 15 exhibit roughly sculptured frons (character #24, state 1) and pronotal area (character #33, state 1), with foveolate patterns (Figs 1A, D), while only two genera completely lack these features. Mesitiinae attack beetle larvae of Chrysomelidae (Coleoptera) (Argaman 2003), the author found them living into ant nets, hence the thick and robust integument of Mesitiinae is presumably associated with the lifestyle of hosts. Michener (2000) postulated that the rough sculpturation (lamellae, carinae and foveolation) and projections could be related to the strengthening of the integument, providing defensive mechanisms for vulnerable areas such as the neck, base of metasoma, and other membranous regions in kleptoparasitic bees. Similar integumental sculpturation is observed in Nyssonini (Crabronidae), a group of apoid wasps that also exhibit kleptoparasitic behavior (Bohart and Menke 1976) and Chrysidinae (Chrysididae), which attack bees and aculeate wasps (Kimsey 1992). The same and convergent features can be observed in Mutillidae (Ronchetti and Polidori 2020). Thus, the integumental thickening in Mesitiinae seems to be associated with defense against their aculeate hosts, as mentioned by Lucena and Almeida (2022?) for Chrysidinae.
Cryptocephalini and Clytrini (Cryptocephalinae) leaf beetles have close association with ant nests. The larval stages remain in the ant nest in a positive interaction. The cocoon brought by the beetle mother is carried by the ants into the nest to complete its development (Agrain et al. 2015). Thus, to reach their beetle hosts, the Mesitiinae need to enter the ant nests. The convergent behavior among all taxa above is that all of them have dangerous hosts (bees and ants), as cited Thus, dense foveolation and thick integument could be associated with defenses.

Relation between pronotal structure and head movements
The median pronotal line (Character #36, state 1) characteristic of many mesitiine lineages (Fig. 1D) is associated with the pronotum-postoccipital muscle (Vilhelmsen et al. 2010), which has its origin at the internal ridge associated with this impression. This ridge could increase the anchorage insertion point allowing stronger contractions. The muscle is the pronotal elevator of the head, so the increase of power for this muscle allows more possibilities for head movements. Nagy (1968) described the parasitoid behavior of females of Mesitiinae and recorded that they steal pre-pupal beetles from the ant nests using their mandibles. Therefore, the increased range of head movements could be adaptative in the context of the female parasitoid behavior.

Propodeal adaptions
Among bethylids, the posterior projections on the propodeum ( Fig. 2A) are exclusive for Mesitiinae (Character #63, state 1). However, it is absent in Anaylax, Astromesitius, Bradepyris, Clytrovorus, Hadesmesitius gen. nov., and Moczariella. Argaman (2003) argued that this posterior projection could facilitate the opening of the cocoon wall during adult emergence, but this was never confirmed for Mesitiinae species. Perhaps a more plausible hypothesis is that the posterior propodeal projection is associated with defense against ants, with the projection protecting the base of metasoma, preventing damage to the petiole. On the other hand, the musculature could indicate another adaptation associated with this structure. The muscle T1-S/T2 has its origin at the posterior corner of the metapectal-propodeal complex and inserts at anterior margin of second metasomal segment. We dissected some Mesitiinae specimens with this projection and observed that the T1-S/T2 muscle has its origin inside the projection, thus increasing the anchorage insertion point of this muscle and giving it more strength. Additionally, this muscle is related to the sternal torsion of the metasoma (Mikó et al. 2007).
All mesitiine wasps have the second metasomal segment longer than the others, an exclusive feature for the subfamily (Barbosa and Azevedo 2011), which was also recovered herein as a synapomorphy for Mesitiinae. The great degree of metasomal segment modification, mainly the length of the third segment longer than others, is associated with oviposition, copulation and defense. Moreover, an additional muscle row was recorded in association with this segment expansion (Kimsey 1992), which is an additional anchorage point for the T1-S/T2 muscle into the posterior propodeal projection.

Hypopygium modifications
There are several shapes of hypopygium exclusive to Mesitiinae, which are described by seven characters in the present analyses (Characters #94 to #100).
Clade A (Fig. 4 and 6) includes seven genera, comprising 113 described species, which represent 63.8% of the total diversity of the subfamily. In this clade are also included the largest species of Mesitiinae.
Muscles located at the base of male genitalia are responsible for movements such as protraction as well as copulation, being inserted at the anterior region of the hypopygium, including the spiculum and anterolateral apodeme. Therefore, the contraction and relaxation between the genitalia base and the spiculum provides the movement of genitalia structures. Thus, the shape and size of the spiculum have direct association with insertion of muscles in the genitalia, affecting the kind and potential of its movements that is, with more and diversified muscle insertions, structures will be capable of performing more complex movements.
The modification of the length of the hypopygial branches is associated with the deformation of the hypopygium, which results from the contraction of the muscles. More muscles inserted at a longer spiculum promote a higher degree of hypopygium deformation, hence the long branches are associated with long median indentation, providing the hypopygium with an area of deformation, giving the structure more flexibility. This is also associated with a wide spiculum. On the other hand, short branches are associated with a simple acute spiculum, since the muscle contraction provides less deformation to the hypopygium, without the need of a deformation area.
The hypopygium shape is associated with muscle insertion and hence it could provide specific functions and adaptations for each genus. Schulmeister (2003) recognized these muscles and named them as "a", "b" and "c". These muscles originate in the gonocondyle in the cupula and inserts at the spiculum and laterally at the ninth sternite (= hypopygium). The cupula is attached to the male genitalia base, and the muscles among these sclerites promote some movement of the genitalia, thus the muscles gonocondyle-spiculum (a), laterally of gonocondyle-spiculum (b), and gonocondyle-laterally ninth sternite (c), have indirect action in male genitalia action (Schulmeister 2003). These muscles expose the male genitalia by the elevation of the basal margin of the cupula (Boudinot 2013), thus probably a larger insertion point could increase the torque movement.

Distribution and biogeography
Presently, Mesitiinae are known from warm regions of the Old World, encompassing all of its four zoogeographical regions: Afrotropical (including Madagascar), Australian, Oriental, and Palearctic (Fig. 7). Unfortunately, Australomesitius mirus Barbosa & Azevedo, 2006, the single species of the family ever recorded in Australia, could not be included in this analysis. However, its position can be inferred based on the shape of the apex of median clypeal carina, arched [in profile], the presence of fusion between sublateral and inner discal carina of propodeal disc, and the orientation of inner discal carina of propodeal disc not parallel with median carina. Observing the distribution patterns among Mesitiinae lineages, there is an apparent association of the early-diverging lineages (e.g., the genera Bradepyris and Moczariella) with the Palearctic region. Additionally, thirteen of nineteen Mesitiinae genera have species recorded from the Palearctic region. From the remaining six genera, Australomesitius is endemic of Australia and Pilomesitius is endemic of Madagascar, while the other four are recorded in the Oriental and Afrotropical regions.
Based on this pattern, it may be that the early diversification of Mesitiinae occurred in the Palearctic region, with lineages progressively occupying the adjacent Oriental and Afrotropical regions and, later, Australia and Madagascar. Occupation of Madagascar likely occurred several times independently within the subfamily, while only a single lineage was able to reach Australia. In the future, a dated phylogeny of the family allied to the exploration of its relationships among other bethylid lineages and discoveries regarding fossil history may provide valuable evidence for a detailed approach on biogeography and diversification of Mesitiinae.

Conclusions
The present study is the most comprehensive cladistic treatment focusing on Mesitiinae tribal classification and character evolution, and the first to treat a large representative group and more accurate classification for each tribe.
Triglenusini were recovered as paraphyletic, and Domonkosini, Heterocoeliini and Mesitiini as polyphyletic, showing the classification in Argaman (2003) as unsup-ported. Morphological characters previously used in the former studies (Nagy 1969(Nagy , 1972Móczár 1970Móczár , 1971 were shown to be inconsistent regarding the monophyly of tribes. Thus, we corroborate the elimination of tribal treatment, following Azevedo et al. (2018).

Authors' contributions
DNB. planned, prepared, and designed the study. MH. and AL. supervised the study. DNB. performed the photography. DNB and MH. Performed the cladistic analyzes, AL. performed the Bayesian analyzes.

Competing interests
The authors have declared that no competing interests exist.

Acknowledgments
We thank the curators of the museums listed for providing the material required for this study, Géllert Púskás and Sándor Csósz for hosting DNB at HNHM, Celso O. Azevedo for being the main advisor to DNB during this study and facilitating access to the laboratory structure; thanks also to Geane Lanes, Ricardo Kawada, Fernando Noll, and Yuri Leite for careful revisions of the manuscript; to the Ernest Mayr Grant