Comparative analysis of morphospace of Neotropical Sericini (Coleoptera: Scarabaeidae): disparity in the light of species diversity and activity patterns

The striking diversity among herbivorous Scarabaeoidea is still poorly understood. Therefore, we compare in these work linear measurements of body features of Neotropical and Old World linages of Sericini chafers to uncover patterns possibly linked to foraging behavior (day or night activity) or to the success of their diversification in terms of diversity. If diversity in Sericini chafers would be linked to morphological disparity, we would expect in the much less diverse Neotropical lineage a very clear drop of morphological disparity. Comparing the two major sister clades, Old World and Neotropical Sericini, in terms of their disparity and diversity, the here elaborated distance data do not support the hypothesis that morphological disparity could be a driver for the diversity of Sericini in the Old World. The influence of the metacoxal length on morphological disparity appears to represent a possible key for a better understanding of the evolution of Sericini in the Neotropical region (including their lower diversity) and answers the question of why they are so much less diverse than their sister lineage in the Old World. The increasing length of the metacoxal plate appears to be a result of the presence of a secondary metacoxal joint in Old World Sericini, which has a crucial impact on hind leg mobility and digging behavior. General body shape and single morphological characters appear under similar general evolutionary pressure: high morphometric disparity between nocturnal and diurnal Neotropical species was also corroborated by disparity estimated from discrete morphological characters.


Introduction
Although the enormous diversity of host-specific herbivore insects is well explained by the hypothesis of an arms race between host plant and host (Ehrlich & Raven, 1964;Farrell, 1998;Mitter et al., 1991), the striking diversity among herbivorous Scarabaeoidea is still poorly understood Eberle et al., 2014) since they are prevailingly polyphagous herbivorous. While follow-up diversification based on angiosperm explosion and high survival rates are other explanations, particularly for beetle diversity (Hunt et al., 2007), these, however, do not all apply for herbivorous scarabs. In the search for explanations for scarab diversity, ecomorphological analysis is one important tool to understand the coexistence between closely related species and to uncover important evolutionary traits of the organisms (Guillerme et al., 2020;Ricklefs & Miles, 1994). Particularly, morphological disparity may unveil hidden evolutionary patterns and causalities (Deline, 2021;Deline et al., 2018;Nel et al., 2018) through the insight into the relative importance of genomic and ecological factors in macroevolution (Foote, 1997). The study of the relationship of disparity vs species diversity may reveal interesting evolutionary correlations and causalities particularly, when the objects are difficult to observe (Copilaș-Ciocianu & Sidorov, 2021;Triantis et al., 2016) or not anymore, such as in fossils (e.g., Neige, 2003;Nel et al., 2018;Ruta et al., 2013;Wan et al., 2021).
In scarabs, the use of morphological disparity, i.e., the estimate of the lineage's occupancy of a suitably defined n-dimensional morphospace (Erwin, 2007;Foote, 1997;McGhee, 1999;Wills, 2001), has been demonstrated recently in several studies. Inward et al. (2011), for example, demonstrated differences of morphospace in different functional groups of feeding habits over various assemblages of dung beetles (Scarabaeinae), even in unrelated species or groups. Eberle et al. (2014) showed that some lineages with different feeding habits may not be distinguished by morphospace in Pleurostict scarabs, while others and several phylogenetic lineages can; they also found that the morphospace coevolved with some morphological key innovations regarding locomotion, particularly in Sericini.
Sericini are herbivore pleurostict chafer beetles with more than 4000 described species and a near cosmopolitan distribution (Ahrens, 2006). It is one of the oldest extant chafer lineages Eberle et al., 2017a;McKenna et al., 2019) that had its origin in the West Gondwana during mid Cretaceous with repeated subsequent dispersals out of Africa (Eberle et al., 2017a). The oldest split of Sericini is between an Old World and a Neotropical lineage, showing an extreme contrast in species diversity: the Old World lineage is about 20 fold more diverse.
While most part of the recent taxonomic and systematic work on the tribe has focused on the Old World taxa, very little is known about the Neotropical lineage. This includes the scarce knowledge on taxonomy (Frey, 1973) as well as their diversity (i.e., number of species), distribution, and ecology. The findings of Eberle et al. (2014) were based on a few exemplar species, for which DNA data were available (Ahrens & Vogler, 2008). However, we were interested whether their findings, particularly for the Sericini, are consistent with a more complete sampling of the fauna, especially of the poorly sampled Neotropics. Therefore, we compare in this work linear measurements of body features of Neotropical and Old World linages of Sericini chafers taken routinely during taxonomic revisions to uncover patterns possibly linked to foraging behavior (day or night activity) or to the success of their diversification. If diversity in Sericini chafers would be linked to morphological disparity, we would expect in the much less diverse Neotropical lineage a very clear drop of morphological disparity. In case that disparity and diversity are not linked, patterns would be random and not significantly distinct.

Material and methods
We used six different body variables (Fig. 1), of which three comprised body ratios, and three were linear distance measurements. These measurements are taken routinely during taxonomic revisions in Sericini (e.g., Ahrens, 2004;Fabrizi et al., 2021) as they show considerable and easy to capture variation between even closely species and allow thus taxon discrimination. Furthermore they integrate quite well different body traits in a possibly functional context, such as relative size of eyes, relative length of legs, and general body proportions. We analyzed data for four available (out of five) South American genera and nineteen Old-Word genera of Sericini chafers (mainly from Asia), represented so far by 108 and 736 species, respectively (list of examined species, see Supplement Table 1). Data were taken on male holotype, syntype, or lectotype specimens of the respective species from newly measured specimens (Pacheco et al., in review a, b) and the following published taxonomic revisions Ahrens et al., 2014a, b, c;Bohacz & Ahrens 2020a, b;Chandra et al., 2021;Eberle et al., 2017b;Fabrizi et al., 2019aFabrizi et al., , b, 2021Liu et al., 2014a, b, c, d;2015aPacheco et al., 2021;Ranasinghe et al., 2020;Sreedevi et al., 2018Sreedevi et al., , 2019. In case no males were available, we used available female specimens (in a few cases in Astaena; Supplement Table 1).
All measurements were taken with the help of an ocular grid on a stereomicroscope (we had in course of the years various models), and values were converted for the distance measurements to millimeters for the different magnifications. Measurements were taken in such a way that the endpoints were in focus. Rations were calculated directly from measured values.
Obtained linear measurements were log10-transformed to render more linear relations among variables and to obtain a similar dimension of variance (Klingenberg, 1996;Ricklefs et al., 1981). Generally, the major component of variance in morphometric data sets of biological specimens is explained through size (Burnaby, 1966;Ferrario et al., 1995;Jolicoeur, 1963). However, since we consider size a crucial component of morphospace as well, we refrained to employ approaches that separate size from shape information (Eberle et al., 2014). Nevertheless, for a number of traits we used directly ratio measures, as these more directly infer some trait characteristics (e.g., the general shape of a leg) and they have been therefore used since long in taxonomic works (e.g., Ahrens, 1999Ahrens, , 2004. Patterns of morphometric covariation were analyzed with standard principal component analysis (PCA) on raw and log-normalized data. Additionally, we omitted the total body length from the analysis to reduce the effect of overweighting overall body size. Euclidean distances between species were calculated directly from PC loadings, as means of general morphological disparity being subsequentially partitioned by lineage assignment or ecological groups (i.e., day activity habits). All calculations for the analysis of morphospace were made in PAST (Hammer et al., 2001); PC axes 1 and 2 were plotted with Excel. Pearson correlation of traits was calculated in Excel.
Finally, we compared the disparity resulting from morphometric patters of the species, with the morphological disparity extracted from discrete morphological characters of the examined South American species (from Pacheco et al., in review b). For this purpose, we exported the uncorrected distances in Mesquite  and performed on these a principal coordinate analysis (Gower, 2005) in PAST.

General trait variation
While among the herbivorous scarabs, Sericini are generally smaller compared to the majority of species, the overall body size varied considerably in our data, with the length ranging from 3.6 to 14.7 mm and maximum width from 2 to 9.8 mm. Infraspecific variation was, however, not considered here, since taxonomic treatments in Sericini (see Material and methods for references) provide in general a range of observed variation rather than induvial measure assessments (mainly due to time constrains and limited taxonomic value). Among the here examined traits, the elytral length showed the strongest correlation with the total body size ( Fig. 2) (r = 0.94). The correlation between elytral length vs. body width (r = 0.91), as well as total body size vs. body width (r = 0.87) was slightly less tight. In contrast to that, none of the body ratios seemed to be very closely related to the body size (r = 0.13, − 0.05, 0.004, for BL vs. MTL/MTW, BL vs. ED/IOD, BL vs. MCL/MEL, respectively).
Diurnal species with flower-visiting habits were found to be distinctly smaller, generally smaller than 7 mm. Night active taxa resulted to be generally larger (> 7 mm). Other traits could not be revealed as particularly linked with the activity time, although diurnal species resulted to have always small eyes as well as relatively short and wide metatibiae. Night active species, which represent most of the species diversity of Sericini worldwide as well as in South America, showed a much wider variation. Neotropical lineages showed not distinctly different patterns in morphospace. However, the length of the metacoxa (in relation to the length of metepisternum) is in the mean generally lower compared to the Old World taxa, although this difference is not pronounced by a distinct gap (Fig. 3). The variation in all traits was largest in the most diverse but presumptively polyphyletic genera, Maladera and Neoserica (Ahrens & Vogler, 2008;Eberle et al., 2017a;Liu et al., 2015b).

Morphospace and disparity
The plots of the most informative first two PC axes showed two slightly different pictures for raw data and log-normalized data, with the major difference that the overlap between diurnal and nocturnal species was reduced with the log-normalized data. When the total body length was excluded from analysis ( Fig. 4B, D), results were similar to the complete data analysis. For the Neotropical genera, this led to a sharp separation between Astaena and Symmela/ Raysymmela in morphospace, while the morphospace of Symmela and Raysymmela widely overlapped. The Old World genera included in this analysis, which comprised in a large number of nocturnal lineages, had a rather pronounced overlap in morphospace (systematic specimen assignment (i.e., genus-level) not shown here in detail). After log normalization, the different traits contributed in a different way to the patterns of morphospace. Elytral length (EL) and body maximum width (EW) had the greatest impact in all analyses on the morphospace (Table 1), and for the all-raw data analysis also the total body length (BL). Interestingly, the first two had more impact in the Fig. 3 Boxplots for selected genera of the total body length (in mm), of the ratio maximal body width/ maximal elytra length (EL), of the ratio metatibial length/maximal metatibial width, of the ratio maxi-mal eye diameter/minimal interocular distance, of the ratio maximal length of metacoxa/maximal length of metepisternum log-normalized analysis than BL. Although being numerically of less impact, for the distribution of the species in morphospace also, the other variables were highly important as it was visible from the vector loadings of the traits (Fig. 4), where coxal length (i.e., R: coxa) and metatibial proportions (i.e., R: metatibia) influenced the morphospace in orthogonal direction compared to the general length measures. These vectors mainly explained the morphological divergence between Neotropical (i.e., Astaena) and Old World, while the body size (EW, EL, BL) corresponds mainly to the divergence between diurnal and nocturnal forms, for normalized data also, the size of eyes (R: eyes) became important.
The general morphological disparity was measured also employing observed distances between the species, optionally partitioned by lineage assignment or ecological groups (i.e., day activity habits). Compared to the mean and median distances between the species in morphospace, Neotropical taxa were generally above this value, while remainder Old World species were below that mean. This was the case for raw data as well as for normalized data. Maximum distances occurred, however, among the Old World species (Table 2). We found contrasting results between raw and normalized data for diurnal activity patterns: while diurnal taxa are with raw data well below mean and median values, they are above with normalized data (Table 3). When diurnal activity patterns were examined separately for Neotropical and Old World taxa (Table 4), it became evident that Old World diurnal taxa were morphologically more despair than Neotropical ones, while nocturnal taxa did not differ in their morphological disparity in relation to their origin, although maximum distances always occurred among the Old World taxa.
High morphometric disparity between nocturnal and diurnal Neotropical species was also corroborated by disparity estimated from distance data derived based on discrete morphological characters. The latter showed a clear Fig. 4 Patterns of disparity derived from morphospace: plots of PC 1 and 2 A raw data (all); B raw data (BL excluded), C log-normalized data (all); D log-normalized data (BL excluded) differentiation between the nocturnal Astaena species and the diurnal Raysymmela and Symmela species. The latter two genera, however, could not be clearly separated with the plots based on these characters (Fig. 5).

Discussion
Most Sericini are morphologically uniform, and differentiation of many species by external features is difficult or impossible (Ahrens, 2004). This extreme similarity, linked with a high degree of homoplasy of diagnostic characters, has also as a result that several of the current major genera are not monophyletic and that females are hard to classify to genus without any associated males (Ahrens, 2004;Eberle et al., 2017a, b). Therefore, morphometric methods seem to be particularly valuable for investigation of the subtle morphological disparity within this tribe.
Our results corroborated some general correlation between body shape and the period of day activity, in which the major explaining variable of divergence is total length. Body measurements and daily activity patterns were found to be in correlation with body color, which was not assessed here by appropriate analysis but by background knowledge of the senior author and observations on the here examined specimens. Many of the diurnal species, beside of having a smaller body size, show a particular bi-colored dorsal surface (with dark head and pronotum and yellow or reddish elytra with dark margin or spots), for example, some species of the genera Symmela and Raysymmela, but also the Asian Anomioserica, Microserica, and Oxyserica. This color pattern (Mallet & Joron, 1999;Pasteur, 1982) is shared with lycids (net winged beetles) (e.g., Bocek et al., 2019;Motyka et al., 2018Motyka et al., , 2021, blister beetles, or leaf beetles (Bocak & Yagi, 2010;Bocakova et al., 2016).
The metacoxal length is an important component of the observed morphological disparity (Fig. 4). It appears to represent a possible key for a better understanding of the evolution of Sericini, in general but particularly in the Neotropical region, which could answer the question, why the Neotropical lineages are so much less diverse than their sister lineage in the Old World. The increasing length of the metacoxal plate appears to be a result of the origin of a secondary metacoxal joint in Old World Sericini (Ahrens, 2006;Eberle et al., 2014) which has a crucial impact on hind leg mobility and digging behavior (Eberle et al., 2014). This hypothesis was confirmed here by all Neotropical taxa having consistently shorter metacoxa lacking the secondary metacoxal joint. The metacoxa in the latter is often even shorter than the metepisternum, while in Old World lineages, the metacoxa is always longer. However, this complex question needs further investigation, as Macrodactylini, being of similar body size as Sericini but much less diverse in the Old World, kind of replace the more rare Sericini in the new World in terms of diversity and dominance in most species habitats and assemblages (e.g., Cherman, 2011;García Atencia, 2014;García-López et al., 2010, 2012 Table 2 Euclidean distances between species partitioned by their provenience (compared to total evidence), calculated directly from PC axes from raw data ( r ) and log-normalized data ( L ) Values highlighted in italics are above the all taxa mean value  Table 3 Euclidean distances between species partitioned by their day activity pattern (compared to total evidence), calculated directly from PC axes from raw data ( r ) and log-normalized data ( L ) Values highlighted in italics are above the all taxa mean value  (Ahrens et al., 2014a, b, c, d).
Comparing Old World and Neotropical Sericini in terms of their average disparity and diversity, the here elaborated distance data (Tables 3 and 4) do not support the hypothesis that morphological disparity could be a driver for the diversity of Sericini in the Old World. The lack of clear disparity within assemblages, in contrast to those of dung beetles (Inward et al., 2011), was reason why Eberle et al. (2014) hypothesized a lack of competition among adult chafer beetles. However, as shown here, the feeding trait, such as herbivory, might not be sufficient to promote morphological disparity, especially when the herbivory is little host specific as in most scarab chafers.
In contrast to that, lineage-specific correlation between disparity and diversity was reported for many other organisms (e.g. Mindel et al., 2016;Neustupa et al., 2009;Gerovasileiou et al., 2015;Zhang et al., 2019;Tong et al., 2020Tong et al., , 2021. In some cases, correlation was evident only at higher systematic level, but not at species level (Tong et al., 2021), or among ecotypes (Tong et al., 2020). These results, however, are to be considered with caution since their measure of diversity was often the number of sampled taxa (that was partly very incomplete in respect to the true diversity of the lineage), and no phylogenetic correction was performed (Garland, 1992).
Due to the present lack of a robust molecular phylogenetic hypothesis for Neotropical Sericini (Eberle et al., 2017a), we were limited here to a rather simple comparative approach that did not consider phylogeny. A consideration of the latter would be useful in order to correct for the dependency of the signal due to the genealogical relationship of examined lineages as performed by phylogenetic contrast analysis (Garland, 1992;McPeek, 1995). Comparisons should consider therefore Table 4 Euclidean distances between species partitioned by their day activity pattern and provenience (compared to total evidence), calculated directly from PC axes from raw data ( r ) and log-normalized data ( L ) principally well-known sister clades, or functional groups (e.g., ecotypes) that are anyway not supposed to be related. Lacking phylogenetic data for the Neotropical lineages, conclusions we could make were limited. On the other hand, many Old World lineages from tree topology of Eberle et al., (2017a, b) are still in process of taxonomic discovery as well; thus, our numbers of species diversity for defined clades are in many cases rather preliminary too. Only the most recent work on Chinese Maladera species brought out more than 150 species (Fabrizi et al., 2021), which equals the total fauna of Neotropical Astaena.
Nevertheless, some initial signals of evidence regarding disparity between Old World Sericini and Neotropical ones but also between diurnal and nocturnal species became evident here. The fact, that this latter disparity was corroborated by disparity estimated from discrete morphological characters (Fig. 5), indicates that general body shape and a majority of morphological characters appear under similar general evolutionary pressure. This may point the directions for future investigations, not only for Sericini, but all evolutionary and synecological research of pleurostict chafers. Their diversification was putatively triggered by the dominance of Angiosperms (Ahrens et al., 2014a, b, c, d), and the detailed investigation of their phylogenetic relationships and morphological evolution is crucial for understanding of the processes generating extant beetle diversity.
Author contribution TP and DA designed the research. DA and MM coordinated and supervised the study. TP and DA con-ducted most of the practical analyses and drafted the manuscript. TP, DA, and MM contributed to the writing of the manuscript. All authors have read and approved the manuscript.

Availability of data and materials
The datasets used and analyzed in this study are available from the corresponding author upon request.

Declarations
Ethics approval and consent to participate Not applicable.

Competing interests The authors declare no competing interests.
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.