The evolutionary history of the ancient weevil family Belidae (Coleoptera: Curculionoidea) reveals the marks of Gondwana breakup and major floristic turnovers, including the rise of angiosperms

The evolutionary history of the ancient weevil family Belidae (Coleoptera: Curculionoidea) reveals the marks of Gondwana breakup and major floristic turnovers, including the rise of angiosperms | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results The evolutionary history of the ancient weevil family Belidae (Coleoptera: Curculionoidea) reveals the marks of Gondwana breakup and major floristic turnovers, including the rise of angiosperms View ORCID Profile Xuankun Li , Adriana E. Marvaldi , Rolf G. Oberprieler , Dave Clarke , Brian D. Farrell , Andrea Sequeira , M. Silvia Ferrer , Charles O’Brien , View ORCID Profile Shayla Salzman , Seunggwan Shin , William Tang , View ORCID Profile Duane D. McKenna doi: https://doi.org/10.1101/2024.03.12.584611 Xuankun Li 1 Department of Entomology, College of Plant Protection, China Agricultural University , Beijing 100193, People’s Republic of China 2 Department of Biological Sciences, University of Memphis , Memphis, TN 38152, USA 3 Center for Biodiversity Research, University of Memphis , Memphis, TN 38152, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Xuankun Li For correspondence: xuankun.li{at}cau.edu.cn marvaldi{at}fcnym.unlp.edu.ar Adriana E. Marvaldi 4 CONICET, División Entomología, Facultad de Ciencias Naturales y Museo, Universidad Nacional de La Plata , Paseo del Bosque s/n, 1900 La Plata, Buenos Aires, Argentina 5 (AEM) CONICET, Instituto Argentino de Investigaciones de Zonas Áridas , Entomología, CC 507, 5500 Mendoza, Argentina (AEM, MSF) Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: xuankun.li{at}cau.edu.cn marvaldi{at}fcnym.unlp.edu.ar Rolf G. Oberprieler 6 CSIRO, Australian National Insect Collection , GPO Box 1700, Canberra, Australian Capital Territory 2601, Australia Find this author on Google Scholar Find this author on PubMed Search for this author on this site Dave Clarke 2 Department of Biological Sciences, University of Memphis , Memphis, TN 38152, USA 3 Center for Biodiversity Research, University of Memphis , Memphis, TN 38152, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Brian D. Farrell 7 Department of Organismic and Evolutionary Biology, Harvard University , Cambridge, MA, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Andrea Sequeira 8 Department of Biological Sciences, Wellesley College , Wellesley MA, 02481, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site M. Silvia Ferrer 5 (AEM) CONICET, Instituto Argentino de Investigaciones de Zonas Áridas , Entomología, CC 507, 5500 Mendoza, Argentina (AEM, MSF) Find this author on Google Scholar Find this author on PubMed Search for this author on this site Charles O’Brien 9 2313 West Calle Balaustre , Green Valley, Arizona 85622, USA (CO) Find this author on Google Scholar Find this author on PubMed Search for this author on this site Shayla Salzman 10 Department of Entomology, University of Georgia , Athens, GA 30602 USA (ShSa) Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Shayla Salzman Seunggwan Shin 11 School of Biological Sciences, Seoul National University , Seoul 08826, Republic of Korea Find this author on Google Scholar Find this author on PubMed Search for this author on this site William Tang 12 Florida State Collection of Arthropods, Florida Department of Agriculture - DPI , P.O. Box 147100, Gainesville, FL 32614, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Duane D. McKenna 2 Department of Biological Sciences, University of Memphis , Memphis, TN 38152, USA 3 Center for Biodiversity Research, University of Memphis , Memphis, TN 38152, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Duane D. McKenna Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract The rise of angiosperms to ecological dominance and the breakup of Gondwana during the Mesozoic marked major transitions in the evolutionary history of insect-plant interactions. To elucidate how contemporary trophic interactions were influenced by host plant shifts and palaeogeographical events, we integrated molecular data with information from the fossil record to construct a timetree for ancient phytophagous weevils of the beetle family Belidae. Our analyses indicate that crown-group Belidae originated approximately 138 Ma ago in Gondwana, associated with Pinopsida (conifer) host plants, with larvae likely developing in dead/decaying branches. Belids tracked their host plants as major plate movements occurred during Gondwana’s breakup, surviving on distant, disjunct landmasses. Some belids shifted to Angiospermae and Cycadopsida when and where conifers declined, evolving new trophic interactions, including brood-pollination mutualisms with cycads and associations with achlorophyllous parasitic angiosperms. Extant radiations of belids in the genera Rhinotia (Australian region) and Proterhinus (Hawaiian Islands) have relatively recent origins. Introduction The evolutionary interplay between plant-feeding (phytophagous) insects and vascular plants has fundamentally shaped terrestrial biodiversity for over 400 million years (e.g., Ehrlich & Raven 1964 ; Farrell et al. 1992 ; Wilf 2008 ; Misof et al. 2014 ; Swain et al. 2022 ). The fossil record has provided insights into the historical context and evolutionary processes that have shaped modern insect-plant interactions ( Labandeira & Currano 2013 ). However, (1) the degree to which contemporary trophic interactions reflect those in the past, e.g., representing relatively ancient primary associations versus more recent secondary associations, and (2) the relative roles of vicariance (including major palaeogeographical events) and dispersal in the evolution of host shifts, remain unclear. Gymnosperms (cycads, conifers and relatives), the dominant plant group during the Mesozoic Era, thrived from the Permian to the Cretaceous period, dominating most terrestrial ecosystems and forming intimate associations with many groups of phytophagous insects ( Anderson et al., 2007 ). The transition to the Cretaceous marked the rise of angiosperms (flowering plants), which gradually replaced gymnosperms in many ecological settings to become the dominant plant group ( Anderson et al., 2007 ). This transition was not only a botanical shift but also marks a critical period in the evolutionary history of insect-plant interactions. Angiosperms introduced novel ecological niches and resources, leading to the diversification of phytophagous insects and a reconfiguration of insect-plant associations (e.g., Farrell 1998 ; Wilf 2008 ; McKenna et al. 2009 , 2015 ; Swain et al. 2022 ). By the end of the Mesozoic, the continents had rifted into nearly their present forms, though their positions would continue to change. Gondwana split into South America, Africa, Australia, Antarctica and the Indian subcontinent, while Laurasia became North America and Eurasia. On account of their extraordinary taxonomic diversity and varied trophic interactions with contemporary gymnosperms and flowering plants, beetles (order Coleoptera; > 400000 described extant species) have been widely used as models to study the evolution of diversity at the insect-plant interface (e.g., Farrell 1998 ; McKenna et al. 2009 , 2019 ). The present paper addresses a significant gap in our understanding of the biogeography and evolution of host plant associations in the family Belidae (belid weevils), an ancient group of phytophagous beetles whose extant species exhibit specialized trophic associations with gymnosperms and flowering plants. By reconstructing the evolutionary history of belid weevils and examining the dynamics of their host plant associations and geographical distributions over time, this paper seeks to yield new insights into the evolution of modern insect-plant interactions. Belidae, with approximately 360 described extant species in 40 genera and two subfamilies ( Marvaldi & Ferrer, 2014 ; O’Brien & Tang, 2015 ), comprise an early branch of Curculionoidea ( Oberprieler et al., 2007 ; McKenna et al., 2009 , 2018 ; McKenna, 2011 ) ( Fig. 1 ). Belidae are indicated to have originated during the Jurassic in association with gymnosperms and subsequently colonized angiosperms, as the latter diversified and rose to ecological dominance during the Late Cretaceous and Paleogene (e.g., Farrell 1998 ; McKenna et al., 2009 ; Shin et al., 2018 ). They exhibit their highest generic diversity in the Southern Hemisphere and collectively are associated with a variety of plants, including conifers, cycads and a few families of angiosperms, including Arecaceae, Balanophoraceae, Hydnoraceae, Celastraceae, Myrtaceae and Vitaceae ( Marvaldi & Ferrer, 2014 ). The adults appear to feed mainly on stem tissues but sometimes also on pollen, some being important pollinators, and the larvae develop in the bark and woody tissues of decaying branches or twigs, gymnosperm strobili, flower buds or fruits. Neotropical Belidae in the subtribe Allocorynina develop as brood pollination mutualists in the pollen cones (or ‘male strobili’) of cycads (e.g., Tang 1987 ; Marvaldi & Ferrer, 2014 ; Salzman et al., 2020 ), a biotic interaction that has several independent origins in weevils. Download figure Open in new tab Figure 1. Phylogeny estimate of Belidae. Topology obtained from analysis of the 46t424g matrix via maximum likelihood (ML), partitioned by gene. Squares along branches are ML SH-aLRT and UFBoot values for 1–3: 1. partition by gene, 2. unpartitioned, and 3. partition by gene and codon positions. Astral support for 4, Quartet Concordance (QC) for 5, and gene concordance factor (gCF) for 6. Branches without squares indicate the highest support in all analyses. Taxa displayed on the right: a. Cyrotyphus vestitus (Agnesiotidini) (photo credit: Rolf Oberprieler), b. Pachyura australis (Pachyurini) (photo credit: Simon Grove), c. Isacanthodes ganglionicus (Belini) (photo credit: Rolf Oberprieler), d. Homalocerus sp. (Belini) (photo credit: Jeff Gruber), e. Rhinotia sp. (Belini) (photo credit: Rolf Oberprieler), f. Oxycraspedus cribricollis (Oxycorynini, Oxycraspedina) (photo credit: Adriana Marvaldi), g. Oxycorynus missionis (Oxycorynini, Oxycorynina) (photo credit: Adriana Marvaldi), h. Rhopalotria slossoni (Oxycorynini, Allocorynina) (photo credit: Shayla Salzman). The historical biogeography and evolution of host associations of Belidae have long drawn the attention of scientists ( Kuschel, 1959 ; Zimmerman, 1994 ; Anderson, 2005 ; Marvaldi et al., 2006 ). However, our understanding of belid evolution remains limited due to the lack of dated phylogeny estimates. Since the landmark study of weevil family-level phylogeny by Kuschel (1995) , several studies based on morphological characters of adults and larvae have inferred the phylogenetic relationships of Belidae and revised their tribal-level classification ( Kuschel & Leschen, 2003 ; Anderson, 2005 ; Marvaldi, 2005 ; Marvaldi et al., 2006 ), and hypotheses for generic phylogenetic relationships based on analyses of morphological data have been proposed for both the subfamilies Belinae ( Kuschel & Leschen, 2003 ) and Oxycoryninae ( Marvaldi et al., 2006 ; Anderson & Marvaldi, 2013 ). However, beyond including small numbers of exemplar taxa in higher-level studies of weevils (e.g., McKenna et al., 2009 ; Shin et al., 2018 ), the phylogeny and evolution of Belidae have not been explored using molecular data. We conducted molecular phylogenetic analyses and divergence dating of belid weevils to investigate their relationships and evolution. We also reconstructed ancestral states of host plant associations and undertook a biogeographical analysis to explore geographical patterns of diversification and the evolution of host plant organ and taxon associations. We integrated phylogenomic and Sanger data for Belidae, sampling all seven tribes and 60 % of the extant genera 1 . We performed ancestral-state reconstruction using the resulting chronogram with stochastic character mapping and event-based likelihood ancestral-area estimation. We sought to answer two main questions about the distribution and host plant associations of crown group Belidae: (1) is the development of larvae of Agnesiotidini and Pachyurini (Belinae) and Oxycraspedina (Oxycoryninae) in conifer hosts an ancient, primary association ( Farrell 1998 ) or the result of more recent, secondary colonization, and (2) how did the interplay between biogeographical process and host plant shifts influence and shape the trophic associations of belids? Results The 33t423g dataset comprised 97334 nucleotides, with 37154 parsimony-informative sites and 31.8 % missing data. The 33t424g and 46t424g datasets contained 107199 nucleotides each, with 45046 and 45157 parsimony-informative sites and 33.3 % and 51.6 % missing data, respectively. These data are summarized in Tab. S4. Flanking regions used in the present study are much shorter than those from most ultraconserved elements (9865 bp in this study, averaging 23 bp per locus versus 400 to 1000 bp per locus) and contributed less to missing data (increasing the amount by less than 1.5 %). The impact of flanking regions on the backbone (tribal-level) relationships has been tested using the reduced datasets (33t423g and 33t424g) with three different partitioning schemes. All six analyses yielded congruent tribal-level relationships and all backbone nodes were robustly supported, except the monophyly of Oxycraspedina + Aglycyderini (Fig. S1a). Analyses of datasets including flanking regions (33t434g) consistently recovered higher statistical support for the sister-group relationship between Oxycraspedina and Aglycyderini, which was further supported by the four-cluster likelihood mapping analyses (Figs S1b, c). Flanking regions were included in subsequent analyses because this increased the resolution and backbone nodal support values. Backbone topologies generated from the combined dataset (46t424g) were congruent with those from the reduced datasets ( Figs 1 , S1). The family Belidae, the subfamilies Belinae and Oxycoryninae, the tribes of Belinae (Agnesiotidini, Belini and Pachyurini) and most subtribes of Oxycoryninae (Oxycraspedina, Metrioxenina, Aglycyderina and Oxycorynina) were robustly supported as monophyletic (UFBoot ≥ 95 AND SH-aLRT ≥ 80) in the trees resulting from the three partitioning schemes ( Fig. 1 ). Tribal-level relationships were robustly recovered in Belinae, but generic relationships were only weakly supported even excluding the taxa without anchored-hybrid-enrichment data (33t423g and 33t424g) (Fig. S1). The tribe Oxycorynini in its current concept ( Marvaldi et al., 2006 ), was not recovered as a monophylum, Oxycraspedina instead forming the sister-group of Aglycyderini + Metrioxenini, and Archicorynus was found to be the sister-group of Allocorynina rather than of all other Oxycoryninae ( Anderson & Marvaldi, 2013 ). The phylogenetic positions of Metrioxenini and Archicorynini were unstable across different partitioning schemes ( Fig. 1 ). This is most likely due to insufficient data, as only Sanger data were generated for representatives of these two tribes. The phylogeny estimate generated with the partition-by-locus scheme was selected as optimal because it closely matches hypotheses based on the analysis of morphological characters ( Marvaldi et al., 2006 ). Therefore, the ML tree generated by the 46t424t dataset and partitioned by locus was used for downstream analyses. Details of node supports are summarized in Fig. S2 and Tab. S2. BEAST analyses with different tree priors/clock schemes yielded similar results, with the crown age of Belidae ranging from 160.0 to 125.6 Ma ( Tab. 1 ). The preferred BEAST analysis, applying a birth-death tree prior with 13 unlinked molecular clocks, was identified using marginal-likelihood estimation ( Tab. 1 ). Divergence-time estimation results revealed an origin of stem Belidae in the Middle Jurassic at 167.8 Ma (95 % highest posterior density = 185.9–160.0 Ma) and of crown-group Belidae in the early Lower Cretaceous at 138.5 Ma (154.9–125.6 Ma) ( Fig. 2 ). The highest likelihood among six biogeographical models tested was DIVALIKE with the unconstrained analysis (M0) (Tab. S3). Belidae were reconstructed as having a Gondwanan origin, with the two extant subfamilies Belinae and Oxycoryninae originating in the Australian and Neotropical regions, respectively ( Fig. 2 ). View this table: View inline View popup Table 1. Marginal-likelihood estimate (MLE) scores for various BEAST analyses performed for this study, and estimated ages (in Ma) for Belidae crown nodes for each tree prior/clock scheme in BEAST. Notes: SS, stepping-stone sampling marginal-likelihood estimation; PS, path-sampling marginal-likelihood estimation; median post-burn-in divergence times in millions of years (95% credibility interval). Download figure Open in new tab Figure 2. Belid timetree. Chronogram generated using BEAST and a birth-death tree prior with 13 unlinked molecular clocks; pale green diamonds on nodes indicate calibration points. Scale is in Ma. Horizontal bars depict the 95% HPD density of each estimate, cyan bars generated by birth-death tree prior, grey bars generated by Yule tree prior. Ancestral area estimation under the likelihood-based implementation of dispersal vicariance analysis and unconstrained dispersal rates (DIVALIKE, M0) performed with BioGeoBEARS. The distribution of each species is mapped to the right of the chronogram. A single most probable ancestral area is mapped at each node. Ancestral state reconstruction of belid larval host plant group and host organ usage under the ER model. The most likely ancestral host plants of Belidae were reconstructed as Pinopsida (0.53) or Angiospermae (0.43) (Fig. S2, Tab. S2). In Belinae, Pinopsida were the most likely ancestral host plants (0.79), with a shift to Angiospermae occurring in Belini (0.65) and two independent colonizations of Polypodiopsida (if indeed these are larval hosts and the two genera are not as closely related as indicated by morphology; see Discussion below). The ancestral host plants of Oxycoryninae were recovered as Angiospermae (0.76) or Pinopsida (0.21). Extant Oxycraspedina are today associated with Pinopsida and Allocorynina with Cycadopsida ( Fig. 2 ). Based on the ASR analysis, branches of the host plants were used for larval development by ancestral Belidae with high possibility (0.87), with a single shift to flowers and fruits (of Angiospermae) in Oxycorynina and at least three shifts to strobili, two to strobili of Pinopsida in Pachyurini and Oxycraspedina and one to (pollen) cones of Cycadopsida in Allocorynina ( Fig. 2 ). Discussion Phylogenetic relationships This study provides the first target-enrichment-based phylogeny estimate of the weevil family Belidae. The monophyly of both subfamilies, Belinae and Oxycoryninae, was recovered, and Aglycyderini were found deeply nested in the latter, supporting the now widely accepted two-subfamily system of belid classification (e.g., Marvaldi 2005 ). In Belinae, all three tribes (Agnesiotidini, Belini and Pachyurini) were recovered as monophyletic, supporting the classification of Kuschel & Leschen (2003) . Generic relationships in Belini were generally poorly resolved, especially when taxa with Sanger data only were included ( Figs 1 , S1). The placements of Trichophthalmus and Homalocerus are likely to be artifacts because they result from sparse and non-overlapping data (four Sanger and 13 anchored-hybrid-enrichment loci, respectively). In contrast, morphological data suggest that these two genera are closely related ( Kuschel, 1959 ; Vanin, 1976 ; Kuschel & Leschen, 2003 ) and together constitute the sister group of all Australian Belini (Fig. 195 of Kuschel & Leschen, 2003 ). In Oxycoryninae, the tribe Oxycorynini in its current concept was found to be polyphyletic, with Oxycraspedina being more closely related to Aglycyderini + Metrioxenini ( Fig. 1 ). A close relationship between Aglycyderini and Metrioxenini was also found in the phylogeny estimate derived from morphological characters ( Marvaldi et al. 2006 ). All other tribes and the subtribes were recovered as monophyletic groups, in agreement with the results of phylogenies reconstructed from morphological characters by previous authors ( Marvaldi et al., 2006 ; Anderson & Marvaldi, 2013 ). However, the placement of Archicorynus as sister-group of Allocorynina (albeit with weak support, Fig. 1 ) differs from that indicated by morphological characters, which resolved the genus as the sister-group of all other Oxycoryninae ( Anderson & Marvaldi, 2013 ). Ancestral host plants and larval feeding habits The larvae of ancestral Belidae were reconstructed as having developed endophytically in branches (not strobili) of their host plants, which is consistent with some previous hypotheses ( Marvaldi et al. 2002 ; Marvaldi 2005 ). Even though our reconstruction provides support for their ancestral hosts most likely being the Pinopsida, which is in accordance with our clade age estimates, angiosperm hosts cannot be completely ruled out, considering that crown angiosperms are now thought to have originated earlier than Belidae, in the interval between the late Permian and the latest Jurassic (256–149 Ma) ( Barba-Montoya et al., 2018 ). However, angiosperm families known to be belid hosts, such as Arecaceae, Hydnoraceae, Balanophoraceae and Fabaceae, did not diversify until the Cretaceous ( Barba-Montoya et al., 2018 ; Li et al., 2019 ). During the early evolutionary history of Belidae, conifer lineages were dominant ( Crisp & Cook, 2011 ; Anderson et al. 2007 ; Wang & Ran, 2014 ), consistent with Pinopsida being the most likely ancestral host plants of Belidae. Based on our phylogeny estimate, the host plant shift to Angiospermae may have occurred during the branching of Belini in the Middle to Late Paleogene. Possible host shifts to ferns in the Paleogene in the Neotropical genera Trichophthalmus and Homalocerus must be viewed with great reservation, as all host records of these genera are only based on the collection of adults, without any evidence of adult feeding, oviposition or larval development ( Kuschel, 1959 ; Vanin, 1976 ). As mentioned above, Homalocerus and Trichophthalmus might form the earliest-branching lineage of Belini. Therefore, the more plausible scenario for host plant shifts is a single shift to ferns from conifers during the Early to Middle Paleogene, perhaps tracking the radiation of ferns in the late Cretaceous ( Schneider et al., 2004 ; Schuettpelz & Pryer, 2009 ). In Oxycoryninae, the South American genus Oxycraspedus (the only extant genus of Oxycraspedina) retains the reconstructed ancestral host association with Pinopsida. Its current hosts, Araucariaceae, were widely distributed in Mesozoic times, and when the oxycorynines diverged by the Lower Cretaceous according to our dated phylogeny estimate, Araucaria still constituted an abundant and widely distributed host for beetles ( Sequeira et al, 2000 ). It is plausible, then, that the ancestors of the remaining oxycorynine tribes/subtribes were originally associated with these conifers before their decline by the Eocene ( Kershaw & Wagstaff, 2001 ) and that these oxycorynine lineages adapted to this decrease in conifer availability by host-shifting to various distantly related taxa of Angiospermae and Cycadopsida. As is true for weevils generally, the association of Metrioxenini with palms (Arecaceae) may date to the Upper Cretaceous, given our dating and the fossil record of the family ( Marvaldi et al, 2002 ; Matsunaga & Smith, 2021 ). However, extant Metrioxenina may only be associated with the genus Arenga ( Marvaldi et al., 2006 ) — though possible associations with other members of the tribe Caryoteae and the subfamily Coryphoideae have not been explored. The Aglycyderini evidently shifted onto angiosperms as well, but their pattern of host associations is complex and indeterminate as the larvae of most species develop in dead bark and twig tissues of a variety of plants. The larvae of Oxycorynina develop in the flowers and fruits of achlorophyllous root-parasitic plants belonging to two distantly related angiosperm families (Balanophoraceae and Hydnoraceae) ( Ferrer et al. 2011 ). The sequence of their host shifts cannot be reconstructed unequivocally; the possibilities include a parallel colonization of these taxa or a shift to Balanophoraceae first and then to Hydnoraceae or vice versa. In any case, our results suggest that these host shifts occurred in the Paleogene at c. 47.7 Ma (68.5–27.5 Ma), consistent with the estimated origins of Balanophoraceae and Hydnoraceae (crown age of Hydnoraceae at 54.7 Ma (75–37 Ma) ( Naumann et al., 2013 )). All extant species of Allocorynina are known to develop in pollen cones of cycads of the genera Dioon and Zamia ( Marvaldi et al., 2006 ; O’Brien &Tang, 2015). Our analysis recovered an initial shift to Dioon (hosts of Parallocorynus and Protocorynus ) followed by a subsequent one from Dioon to Zamia (hosts of Notorhopalotria and Rhopalotria ) ( Fig. 2 ). The shift to Dioon is indicated to have occurred between the late Cretaceous and the early Paleogene, at 58.4 Ma (67.1–49.9 Ma), which is older than the crown age of Dioon (24.6–7.5 Ma) but younger than the stem age (207.9–107.0 Ma) ( Condamine et al., 2015 ). Australian cycad weevils have also been found to have colonized cycads before their main radiation ( Hsiao et al., 2023 ). The cycad tree of life is known for long branches subtending generic-level radiations ( Nagalingum et al., 2011 ; Salas-Leiva et al., 2013 ; Condamine et al., 2015 ; Liu et al., 2022 ). The switch from Dioon to Zamia is estimated to have occurred in the middle Neogene, at 11.0 Ma (13.5– 8.6 Ma), which is in line with the crown group age of Zamia (9.5 Ma, 22.1–9.0 Ma) ( Calonje et al., 2019 ). There appears to be a pronounced conservatism in the type of tissue consumed by belid larvae, always involving parenchymatous parts of branches or reproductive structures (not pollen or seeds). The type of plant organ used for larval development shows few shifts from the ancestral state (branches): in Belinae to ovulate and/or pollen cones (or ‘female and/or male strobili’) of Pinopsida (arguably a shift; possibly just opportunistic development alongside that in twigs and branches) in some Agnesiotidini and Pachyurini, in Oxycoryninae to ovulate cones of Pinopsida in Oxycraspedus , fleshy flowers/fruits of parasitic angiosperms in Oxycorynina and pollen cones of cycads in Allocorynina (there does not appear to be any feeding on ovulate cones by Allocorynina in nature although they clearly visit them and are presumably capable of feeding on them at least in the short term; Simon et al., 2023 ). As is generally true for plant-feeding beetles ( Farrell & Sequeira, 2004 ), host taxon associations are evidently more labile than the use of tissue/organ for larval development. This is particularly noticeable even within some belid genera, whose larvae have been found developing in similar tissues (under bark and in branches) of different plant families, both between congeneric species (e.g., Sphinctobelus ) and in single species (e.g., Isacantha , Rhicnobelus ); however, there are some cases known in Belidae where the use of plant organs is apparently plastic and opportunistic (e.g. Apagobelus brevirostris (Lea) reared from both stems and cones of Araucaria ; Zimmerman (1994) ). Palaeogeographical events and host plant shifts A Gondwanan origin was recovered for Belidae at c. 138 Ma, when South America and Australia were connected via Antarctica. The two subfamilies, Belinae and Oxycoryninae, originated and diverged during the separation of East and West Gondwana ( McIntyre et al., 2017 ) and both have lineages that still develop on their ancestral host plant. In the East Gondwana clade, Belinae, the divergence time between Agnesiotis and Dicordylus was estimated at 44.7 Ma (55.9–34.5 Ma), during the Paleocene–Eocene separation of South America and Australia ( Lawver et al., 1992 ; Briggs, 1995 ; McIntyre et al., 2017 ). Thus, the evolutionary history of Dicordylus might best be explained by persisting in association with the ancestral host plant despite vicariance. The history of the other two South American genera, Homalocerus and Trichophthalmus , could potentially also be elucidated through the vicariance process (though we cannot rule out the possibility of dispersal to South America). This process may have offered an opportunity for their common ancestor to encounter the new niches of South American ferns, thereby facilitating their shift to new host plants. Such palaeogeographical events shaping host plant usage have also been documented in other phytophagous arthropods (e.g., Calatayud et al., 2016 ). Due to the lack of Afrocorynina in our sample and the weakly supported placements of Archicorynini and Metrioxenina, the evolutionary history of Oxycoryninae could not be reconstructed with high confidence. Nonetheless, a plausible scenario is that the common ancestor of Aglycyderini + Metrioxenini shifted to angiosperms in the Lower to Middle Cretaceous at 99.9 Ma (113.9–88.9 Ma) ( McIntyre et al., 2017 ), diverged in the Ethiopian region by the Middle to Upper Cretaceous and subsequently dispersed to the Palearctic and Oriental regions, the tropical Pacific islands and as far south as New Zealand (Aglycyderina). Additional taxon sampling is needed to illuminate whether this host shift promoted the dispersal of Aglycyderini + Metrioxenini. Adaptive radiations in Rhinotia and Proterhinus In Belidae, notable taxonomic diversity is observed in the genera Rhinotia and Proterhinus , with 87 and 168 described species, respectively ( Marvaldi & Ferrer, 2014 ; Brown, 2019 ). In contrast, other genera of Belidae comprise fewer than 15 known species ( Marvaldi & Ferrer, 2014 ). Host plants of Rhinotia are mostly Acacia (Fabaceae), the most species-rich plant genus in Australia, with more than 1000 endemic species ( Maslin, 2004 ; Lewis, 2005 ). Acacia is particularly dominant in arid and semi-arid areas of Australia ( Byrne et al., 2008 ). The crown age of Rhinotia is here estimated at 11.7 Ma (14.1–9.2 Ma), which postdates the origin of crown Acacia (23.9–21 Ma) ( Miller et al., 2013 ). The common ancestor of Rhinotia might have colonized Acacia in the early to middle Miocene when Australia was warm and wet ( Martin, 2006 ), and co-diversified during the aridification of Australia from 10–6 Ma ( Byrne et al., 2008 ). The origin might also have been 20 Ma earlier and involved other now extinct taxa, with Rhinotia being a surviving lineage, or it might have colonized and radiated during the aridification. A similar situation has been postulated for the thrips subfamily Phlaeothripinae, which colonized Acacia and diversified into more than 250 species in 35 genera ( McLeish et al., 2007 , 2013 ). With 159 species described from the Hawaiian Islands, Proterhinus is another example of the spectacular radiation of Hawaiian insects, along with, e.g., the nearly 1000 species of Drosophilidae (Diptera) ( O’Grady et al., 2011 ), more than 400 species of Hyposmocoma (Lepidoptera: Cosmopterigidae) ( Haines et al., 2014 ), over 110 species of Plagithmysus (Coleoptera: Cerambycidae) ( Gressitt, 1975 ) and more than 190 species of Nesophrosyne (Hemiptera: Cicadellidae) ( Bennett & O’Grady, 2013 ). Little is known about the host plants of Proterhinus , but what is known suggests that host plant ranges (e.g., in P. deceptor Perkins, P. obscurus Sharp and P. vestitus Sharp) are remarkably broader than in other species of Belidae ( Legalov, 2009 ). Such host “jumps” are typical for endemic Hawaiian phytophagous insects, e.g., Carposina (Lepidoptera: Carposinidae) ( Medeiros et al., 2016 ) and Nesosydne (Hemiptera: Delphacidae) ( Goodman et al., 2012 ), but some lineages have high host plant specificity e.g., Philodoria (Lepidoptera, Gracillariidae) ( Johns et al., 2018 ) and Nesophrosyne (Hemiptera: Cicadellidae) ( Bennett & O’Grady, 2013 ). Due to limited taxon sampling, we could not estimate the crown age of Proterhinus nor when and from where it may have arrived in Hawaii. Materials and methods Taxon sampling Forty-six taxa were included in the present study, including eight outgroups representing the chrysomeloid families Chrysomelidae and Orsodacnidae and the weevil families Anthribidae, Cimberididae and Nemonychidae. Ingroup taxon sampling spanned 38 species of Belidae in 24 genera, representing all seven tribes from both subfamilies and 60 % of extant belid genera. Genomic data from 33 taxa, 23 newly generated for this study, were used in phylogeny reconstruction. Sanger DNA sequence data (CO1, 16S, 18S and 28S) from six outgroup and 34 ingroup species were also used, including newly generated data for 28 species, and 13 species were represented by Sanger data only (Tab. S1). DNA extraction, library preparation and Illumina DNA sequencing Total genomic DNA was extracted from the legs, thoracic muscle or the whole body, depending on the size of the specimen, using the G-Biosciences OmniPrep kit (G-Biosciences, Catalog #786-136, St. Louis, MO, U.S.A.), following the manufacturer’s protocol, except that samples were incubated for 15 hours instead of 15 minutes. Final DNA extractions were eluted with 60 µL of nuclease-free water and treated with RNaseA. The remaining body parts were preserved in 95 % ethanol as vouchers. Genomic DNA QC statistics were generated for each extracted specimen using a Qubit fluorometer, and DNA quality (fragmentation/degradation and contamination with RNA) was further assessed via gel electrophoresis. The extracted DNA was fragmented by sonication with a Q800R2 Sonicator (Illumina TruSeq), using 50 µL of the DNA extractions in 0.2 mL strip tubes, targeting a modal fragment size of 350 base pairs. Genomic DNA libraries were constructed using the NEBNext Ultra II DNA Library Prep Kit (NEB #E7645L) with NEBNext Multiplex Oligos for Illumina (Dual Index Primers Sets 1 and 2) (NEB #E7600S and #E7780S), with two-sided size selection around a mode of 480 base pairs. Target enrichment through hybridization followed the myBaits Hybridization Capture for Targeted NGS (Version 5), with 65 °C chosen for the hybridization temperature. We used the published Anchored Hybrid Enrichment Coleoptera Probe set ( Haddad et al., 2018 ; Shin et al., 2018 ) and targeted 599 nuclear loci. Enriched libraries were amplified using KAPA HiFi HotStart ReadyMix. PCR cycling consisted of an initial denaturing step at 98 °C for 2 min, followed by eight cycles of denaturing at 98 °C for 20 s, annealing at 60 °C for 30 s, elongation at 72 °C for 45 s and a final elongation step at 72 °C for 5 min. The 192 enriched and multiplexed libraries were sequenced using 150-bp paired-end reads on an Illumina HiSeq Lane at Novogene Corporation Inc. (Sacramento, CA, USA). DNA isolation, PCR amplification and Sanger DNA sequencing DNA extraction and PCR amplification of Sanger data were performed at IADIZA-CONICET (Mendoza, Argentina) and Wellesley College (MA, USA). Total genomic DNA was extracted from adult voucher specimens using an adapted “salting-out” protocol (Sunnuks & Hales 1996) or the DNeasy Blood and Tissue Kit (QIAGEN, MD, USA). Tissue was processed from one to two legs or part of the thorax. The extracted DNA was stored at -20 °C. Four molecular markers (two nuclear and two mitochondrial) were used in this study: 18S rDNA (entire), 28S rDNA (regions D2, D3), 16S rDNA (regions IV, V) and COI (“barcode” or 5’ region). The primers used for amplification and sequencing of the four Sanger loci and PCR conditions are as described by Marvaldi et al. (2018) . The PCR products were purified and bi-directionally sequenced with the Sanger method, using the Sequencing Service of “Unidad de Genómica de INTA-Castelar” (Buenos Aires, Argentina) and, in Wellesley, using an ABI PRISM 3100 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Electropherograms were edited and contig-assembled using ProSeq v.2.91 ( Filatov 2002 ) and sometimes Sequencher v.5 (GeneCodes Corp.). All new sequences were deposited in GenBank under accession numbers PP832953 – PP832961 and PP840348 – PP840386 (Tab. S1). Sequence assembly and orthology prediction The dataset preparation procedure for anchored-hybrid-enrichment-targeted loci used is as outlined by Breinholt et al. (2018) . A reference set was prepared using genomic coding sequences (CDS) from nine coleopteran genomes: Anoplophora glabripennis (Motschulsky) (Cerambycidae, GCA_000390285.2), Aethina tumida Murray (Nitidulidae, GCA_001937115.1), Callosobruchus maculatus (Fabricius) (Chrysomelidae, GCA_900659725.1), Dendroctonus ponderosae (Hopkins) (Curculionidae, GCA_020466585.1), Diabrotica virgifera LeConte (Chrysomelidae, GCA_003013835.2), Gonioctena quinquepunctata (Fabricius) (Chrysomelidae, GCA_018342105.1), Leptinotarsa decemlineata Say (Chrysomelidae, GCA_000500325.2), Sitophilus oryzae (Linnaeus) (Curculionidae, GCA_002938485.2) and Tribolium castaneum (Herbst) (Tenebrionidae, GCA_000002335.3). Raw reads were assembled using an iterative baited assembly (IBA) after filtering with Trim Galore! v.0.4.0 (bioinformatics.babraham.ac.uk). Orthology was determined using the T. castaneum genome as a reference, and single-hit and genome mapping location criteria were used with NCBI Blastn ( Camacho et al., 2009 ). Cross-contamination checks were conducted with USEARCH ( Edgar, 2010 ), and sequences with > 99 % identity across different subfamilies were identified and removed. Cleaned sequences were aligned in MAFFT v.7.245 ( Katoh & Standley, 2013 ), and isoform consensuses were generated using FASconCAT-G 1.02 ( Kück & Longo, 2014 ). Following the method outlined by Teasdale et al. (2016) and Li et al. (2022a) , we used a blast-based method to extract Sanger genes from genomic sequences. In short, raw sequence data were assembled using SOAPdenovo v.2 ( Li et al., 2015 ), the four Sanger genes were identified using an all-by-all tBlastx search, reads were mapped on potential orthologous sequences using BBmap v.35.85 ( Bushnell, 2014 ), and a final consensus sequence was generated after variants were called using GATK v.4.1.1.0 ( McKenna et al., 2010 ). Extracted sequences were compared with Sanger sequencing results, and sequences with higher quality and longer reads were maintained for subsequent phylogenetic analysis. The sequences of the ribosomal markers (nuclear 18S and 28S and mitochondrial 16S) were aligned using information on the secondary structure of the arthropod rRNA genes to identify homologous positions as well as regions of ambiguous alignment to be excluded from analyses ( Gillespie et al., 2006 ; Marvaldi et al. 2009 , 2018 ). Dataset preparation Anchored-hybrid-enrichment probe sets comprise highly conserved coding probe regions (i.e., exons) and more variable, generally non-coding flanking regions (e.g., introns or intergenic regions) located on flanks of the probe region ( Lemmon et al., 2012 ; Haddad et al., 2018 ; Shin et al., 2018 ). Following the pipeline, we trimmed flanking regions with 1.5 entropy and 50 % density cutoffs at each site in the nucleotide sequence alignments ( Breinholt et al., 2018 ). AliView v1.18 ( Larsson, 2014 ) was used to manually check each nucleotide alignment to separate “flanks” from “probe regions” and ensure that the probe region was in the correct open reading frame (ORF). We used a long-branch detection protocol to investigate the possibility of external contamination, paralogous sequences and significant sequencing/assembly errors (longbranchpruner.pl available on Osiris, http://galaxy-dev.cnsi.ucsb.edu/osiris/ ). We produced maximum-likelihood (ML) gene trees from nucleotide (NT)-probe-region multiple-sequence alignments (MSAs) in IQ-TREE v.2.0.6 ( Nguyen et al., 2015 ), conducting a full model test for each gene. We pruned tip sequences that exceeded eight standard deviations from the mean tip length of the gene tree from NT MSAs. Loci with < 40 % taxon coverage were excluded. Five hundred eighty loci were assembled across 35 taxa, and 419 loci were selected for phylogenetic inference. Four Sanger genes were included, and flanking regions were concatenated and treated as a single locus. Cleaned MSAs were concatenated using Phyx v.1.1 ( Brown et al., 2017 ) to generate the dataset with 46 taxa and 424 loci for phylogenetic inference (46t424g). The 13 taxa represented only by Sanger data were excluded to evaluate the impact of missing data and flanking regions, generating two datasets, 33t424g with 33 taxa and 424 loci and 33t423g with 33 taxa and 423 loci (flanking regions excluded). Phylogenetic analyses and tests of node support We conducted ML phylogenetic analyses in IQ-TREE v2.1.3 ( Nguyen et al., 2015 ). For the 33t423g and 33t424g datasets, three partitioning schemes were used: (1) unpartitioned, (2) partitioned by locus, (3) partitioned by locus and codon position. All schemes were model-tested using ModelFinder ( Kalyaanamoorthy et al., 2017 ) as implemented in IQ-TREE. The best partitioning scheme for the 33t423g dataset was found after merging possible partitions (using the ‘-MFP+MERGE’ command) and determining the best scheme under the Bayesian information criterion (BIC). For the 33t424g and 46t424g datasets, the best schemes of the 33t423g dataset were used with the GTR+ASC model added for the flanking regions. An initial 1000 parsimony trees were generated in IQ-TREE with the command ‘-ninit 1000’, and 100 trees with the fewest steps were used to initialize the candidate set (-ntop 100), considering all possible nearest-neighbor interchanges (-allnni). These 100 trees were maintained in the candidate set during the ML tree search (-nbest 100), and unsuccessful runs were terminated after 1000 iterations (-nstop 1000). Perturbation strength was set to 0.2 (-pers 0.2), as recommended for datasets with many short sequences. We used nearest-neighbor interchange (NNI) branch swapping to improve the tree search and limit overestimating branch supports due to severe model violations (‘-bnni’ command). Node supports were computed with 1000 UFBoot replicates (‘-B’ command) ( Minh et al., 2013 ; Hoang et al., 2018 ) and SH-aLRT (‘-alrt’ command) ( Guindon et al., 2010 ). Both concatenation and gene coalescence approaches were used for tree estimation on dataset 46t424g. For concatenated analyses, partitioning schemes and parameters used were the same as those for the 33t424g and 33t423g datasets. Nodes were classified as “robustly supported” when they were recovered with support values of UFBoot ≥ 95 AND SH-aLRT ≥ 80, as “moderately supported” when UFBoot ≥ 95 OR SH-aLRT ≥ 80 and as “weakly supported” when UFBoot < 95 AND SH-aLRT < 80 ( Minh et al., 2013 ; Hoang et al., 2018 ). The best evolutionary model for coalescence analyses was found using ModelFinder Plus (‘-MFP’ command) for each gene, followed by likelihood tree searches using a partitioning scheme generated in IQ-TREE. We applied nearest-neighbor interchange (NNI) branch swapping to improve the tree search and limit overestimating branch supports due to severe model violations (‘-bnni’ command). Nodal support was computed with 1000 UFBoot replicates (‘-B’ command) ( Minh et al., 2013 ; Hoang et al., 2018 ). Multi-species-coalescent (MSC) analyses based on these single trees were calculated using ASTRAL-III v5.6.2 ( Zhang et al., 2018 ). Nodes were classified as “robustly supported” when recovered with support values ≥ 0.95. Four-cluster likelihood mapping (FcLM; Strimmer & von Haeseler, 1997 ) was performed in IQ-TREE using the 33t423g and 33t424g datasets and ModelFinder-determined partitions (partition by locus and codon position) to quantify support for the sister relationship between Allocorynina and Oxycraspedina. The four schemes used to define the four-taxon clusters are Allocorynina, Aglycyderini, Oxycraspedina and outgroups. Quartet sampling ( Pease et al., 2018 ) provided an alternate examination of branch support. Quartet sampling of internal node scores included a set of three scores: quartet concordance (QC: values near 1 being more concordant and near -1 more discordant), quartet differential (QD: the more equal the frequencies of discordant topologies, the closer to 1; 0 indicating that only one other discordant topology was found) and quartet informativeness (QI: 1 for all replicates informative, 0 for no replicates informative) ( Pease et al., 2018 ). Concordance and disagreement among genes on the selected ML tree generated by the 46t424g dataset were estimated using the gene concordance factor (gCF) implemented in IQ-TREE ( Minh et al., 2020 ), as concatenated analyses can return well-supported trees even when the level of gene incongruence is high (e.g., Jeffroy et al., 2006 ; Kumar et al., 2012 ). Nodes were classified as “robustly supported” when gCF was higher than gene discordance factor gDF1 and gDF2, “weakly supported” when gCF equaled gDF1 or gDF2 or both and “not supported” when gCF was lower than gDF1 or gDF2 (more loci supporting an alternative topology). Divergence time estimation Divergence times were estimated in a Bayesian framework using BEAST v1.10.4 ( Suchard et al., 2018 ) on the high-performance computing clusters at the University of Memphis and the China Agricultural University. SortaDate ( Smith et al., 2018 ) was used to reduce the nucleotide alignment to a computationally tractable matrix (50 loci) using combined results of clock-likeness (get_var_length.py) and least topological conflict with the species tree (get_bp_genetrees.py). Four Sanger genes were included in the BEAST analyses to ensure that all taxa were present in the reduced matrix, regardless of the SortaDate result. Three different initial partitioning strategies were used: 1) unpartitioned; 2) four partitions (three codon positions for protein-coding genes and one partition for non-protein coding genes); 3) 144 partitions by locus and nucleotide codon position. The best-fit model and partitioning schemes were selected by PartitionFinder2 ( Lanfear et al., 2017 ) using the greedy algorithm and the BIC scores across all models included in BEAST (models = beast). The reduced concatenated data matrix was imported into BEAUTi ( Drummond et al., 2012 ). Substitution and clock models were unlinked among partitions, and tree models were linked. An uncorrelated relaxed-molecular-clock model ( Drummond et al., 2006 ) and a lognormal prior were applied, two tree priors were tested for each partitioning scheme, Yule (pure birth) and birth–death, and a fixed cladogram based on the topology generated by the concatenated ML analysis was used. Four internal nodes and the root were constrained for calibration based on the fossil record of Belidae. Several supposed belid fossils have been described, but only a few could be confidently placed based on evident synapomorphies. After a careful examination following Parham et al. (2012) , five fossils were selected for calibration of four internal nodes: 1) Sinoeuglypheus daohugouensis Yu, Davis & Shih ( Yu et al., 2019 ) (together with other undescribed Daohugou specimens) for the stem of Cimberididae, with a minimum age of 160 Ma based on the age of the boundary between the Oxfordian and Callovian (ICC 2023); 2) Talbragarus averyi Oberprieler & Oberprieler ( Oberprieler & Oberprieler, 2012 ) for the stem of Rhinorhynchinae, with a minimum age of 147.3 Ma based on analysis of zircons ( Bean 2006 ); 3) Preclarusbelus vanini Santos, Mermudes & Fonseca ( Santos et al., 2007 ) and Cratonemonyx martinsnetoi Legalov ( Gratshev & Legalov, 2014 ) for the stem of Oxycoryninae, with a minimum age of 113 Ma based on the Crato Formation as Upper Aptian following Santos et al. (2011) ; 4) Pleurambus strongylus Poinar & Legalov ( Poinar & Legalov, 2014 ) for the stem of Allocorynini, with a minimum age of 15 Ma ( Iturralde-Vinent & MacPhee, 1996 , 2019 ). Fossil calibrations were introduced as minimum ages of uniform priors, and the lower margin of the estimated timing of origin of Phytophaga (195 Ma) ( McKenna et al., 2019 ) was used as a hard maximum age constraint on the calibrated nodes. Three independent analyses of each clock scheme and tree-prior combination were run to check for convergence. We evaluated the convergence and mixing of MCMC chains in Tracer version 1.6 ( Rambaut et al., 2018 ) to ensure that the effective sample sizes (ESS) exceeded 200. The resulting tree files were combined and resampled with a frequency of 100000 in LogCombiner (BEAST package) and a burn-in of 30 %. Subsampled trees were summarized as a maximum-clade-credibility tree using TreeAnnotator ( Rambaut and Drummond 2015 ), with median heights as node heights. Path sampling and stepping-stone sampling ( Xie et al., 2011 ; Baele et al., 2012 , 2013 ) were performed as part of all BEAST analyses to identify the best tree prior and clock scheme combination. Ancestral host plant reconstructions Ancestral reconstruction analyses were performed for two host plant characters on the dated phylogeny: (i) host plant higher taxa (four states): Angiospermae, Cycadopsida, Pinopsida and Polypodiopsida; (ii) host plant organ (three states): branches, strobili and flower and fruit (Tab. S1). Stochastic character mapping was conducted with the “make.simmap” command in the R package phytools ( Revell, 2012 ) with 1000 simulations. Character state transitions were assumed to have equal rates (“ER” option). Ancestral area estimation We recognize five bioregions that best account for the distribution of the sampled species in Belidae: Neotropical region (N), Australian region (A), tropical Pacific islands (I), Oriental region (O) and Palearctic + Ethiopian regions (P) (Tab. S1). We performed an event-based likelihood-ancestral-area estimation using BioGeoBEARS ( Matzke, 2014 ). Three models were used: (1) DEC (Dispersal Extinction Cladogenesis ( Ree & Smith, 2008 )), (2) DIVALIKE (a likelihood-based implementation of dispersal vicariance analysis, originally parsimony-based ( Ronquist, 1997 )), (3) BAYAREALIKE (a likelihood implementation of BayArea, originally Bayesian ( Landis et al., 2013 )). All models were also evaluated under a constrained analysis (M1), in which we considered palaeogeographical events that occurred in the past 160 Ma over four-time slices (160–125 Ma, 125–35 Ma, 35–28 Ma, 28 Ma to present; sliced by three events: separation of East and West Gondwana, separation of South America and Australia ( McIntyre et al., 2017 ), origin of Hawaii islands (Price & Clague, 2007)) and geographical distance variation for a total of six scenarios. The maximum allowed ancestral area was restricted to two. Time slices for geographical events used were those of McIntyre et al. (2017) . The Akaike Information Criterion (AIC ( Burnham & Anderson, 1998 )) and the corrected Akaike Information Criterion (AICc ( Burnham & Anderson, 2002 )) were calculated to select the best-fitting model. Figure legends Figure S1. (a) Belidae phylogenetic tree inferred from the 33t424g matrix via maximum likelihood (ML), partitioned by gene and codon position. Squares along branches are ML SH-aLRT and UFBoot values. Supporting values are based on the 33t424g matrix for 1–3 and the 33t423g matrix for 4–6. Partitioned by gene and codon position for 1 and 4, partitioned by gene for 2 and 5, and unpartitioned for 3 and 6. (b) Four-cluster likelihood mapping (FcLM) result based on the 33t424g matrix. (c) FcLM result based on the 33t423g matrix. Figure S2. Topology generated by maximum likelihood (ML) analysis of the 46t424g matrix partitioned by gene, with node numbers . Table S1. Taxa sampled in the present study, with the number of loci retained after all cleaning steps, host plant group and organ coding schemes, and distribution coding schemes. Table S2. Summary of node support generated by concatenated ML analysis, coalescence approach, gene concordance factors (gCF), and Quartet Sampling; the posterior possibility of characters generated by ancestral state reconstructions. Table S3. Results of the BioGeoBEARS analyses. Table S4. Summary of matrices analyzed. Acknowledgments This study was funded by grants from the United States National Science Foundation DEB:1355169 and 2110053 to DDM. XL was supported by the 2115 Talent Development Program of China Agricultural University. AEM acknowledges support by CONICET through grant PIP 3108. The authors acknowledge the China Agricultural University and the University of Memphis for providing high-performance computing platforms and support that contributed to the research results reported in this study. Footnotes ↵ † Deceased Results updated to show the flanking regions did not have a detrimental impact on the results of our analyses. Figures updated to fix some spelling errors and classification issues. ↵ 1 Combining Sanger sequences data with genomic data for phylogenetic inference has been demonstrated as a feasible approach to resolving deep-level relationships while adding taxa to the phylogeny for tracing the evolutionary history of characters (e.g. Zhang et al., 2016 Song et al., 2020 ; Li et al., 2022a , 2022b ). References ↵ Anderson , R.S. ( 2005 ). New Oxycoryninae from Central and South America: phylogenetic and biogeographical implications (Coleoptera: Belidae) Systematic Entomology 30 : 644 – 652 doi: 10.1111/j.1365-3113.2005.00304.x OpenUrl CrossRef ↵ Anderson , J.M. , Anderson H.M. , Cleal C.J ( 2007 ) Brief history of the gymnosperms: classification, biodiversity, phytogeography and ecology South African National Biodiversity Institute . South Africa , Pretoria , 2007 pp.x + 280 pp ↵ Anderson , R.S. , Marvaldi , A.E ( 2013 ). Finding unexpected beetles in odd places: Archicorynus kuscheli Anderson and Marvaldi, a new genus and species representing a new tribe, Archicorynini, of Oxycoryninae (Coleoptera: Belidae) from Nicaragua The Coleopterists Bulletin 61 – 71 ↵ Baele G. , Lemey P. , Bedford T. , Rambaut A. , Suchard M.A. , Alekseyenko , A.V ( 2012 ) Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty Molecular Biology and Evolution 29 : 2157 – 2167 doi: 10.1093/molbev/mss084 OpenUrl CrossRef PubMed Web of Science ↵ Baele G. , Li W.L.S. , Drummond A.J. , Suchard M.A. , Lemey P ( 2013 ) Accurate model selection of relaxed molecular clocks in Bayesian phylogenetics Molecular Biology and Evolution 30 : 239 – 243 doi: 10.1093/molbev/mss243 OpenUrl CrossRef PubMed Web of Science ↵ Barba-Montoya , J ., Dos Reis , M. , Schneider , H. , Donoghue , P. C. , Yang , Z ( 2018 ) Constraining uncertainty in the timescale of angiosperm evolution and the veracity of a Cretaceous Terrestrial Revolution New Phytologist 218 : 819 – 834 OpenUrl ↵ Bean , L. B. ( 2006 ) The leptolepid fish Cavenderichthys talbragarensis (Woodward, 1895) from the Talbragar Fish Bed (Late Jurassic) near Gulgong, New South Wales Records-Western Australian Museum 23 : 43 – 76 OpenUrl ↵ Bennett , G.M. , O’Grady , P.M ( 2013 ) Historical biogeography and ecological opportunity in the adaptive radiation of native Hawaiian leafhoppers (Cicadellidae: Nesophrosyne) Journal of Biogeography 40 : 1512 – 1523 OpenUrl ↵ Breinholt , J.W. , Earl , C. , Lemmon , A.R. , Lemmon , E.M. , Xiao , L. , Kawahara , A.Y ( 2018 ) Resolving relationships among the megadiverse butterflies and moths with a novel pipeline for anchored phylogenomics Systematic Biology 67 : 78 – 93 doi: 10.1093/sysbio/syx048 OpenUrl CrossRef ↵ Briggs , J.C. ( 1995 ) Global biogeography Elsevier Science , Amsterdam . XVII + 452 pp ↵ Brown , J.W. , Walker , J.F. , Smith , S.A ( 2017 ) Phyx: Phylogenetic tools for unix Bioinformatics 33 : 1886 – 1888 doi: 10.1093/bioinformatics/btx063 OpenUrl CrossRef ↵ Brown , S.D. ( 2019 ) A new species of Proterhinus Sharp, 1878 (Coleoptera: Curculionoidea: Belidae) from Miti’aro, Cook Islands, South Pacific Zootaxa 4664 : 481 – 496 OpenUrl ↵ Byrne , M. , Yeates , D.K. , Joseph , L. , Kearney , M. , Bowler , J. , Williams , M.A.J. , Cooper , S. , Donnellan , S.C. , Keogh , J.S. , Leys , R. , Melville , J. , Murphy , D.J. , Porch , N. , Wyrwoll , K.-H ( 2008 ) Birth of a biome: insights into the assembly and maintenance of the Australian arid zone biota Molecular Ecology 17 : 4398 – 4417 OpenUrl CrossRef PubMed Web of Science ↵ Burnham , K.P. , Anderson , D.R ( 1998 ) Practical use of the information-theoretic approach In Model selection and inference (eds KP Burnham , DR Anderson ), pp. 75 – 117 . New York, NY : Springer ↵ Burnham , K.P. , Anderson , D.R ( 2002 ) Model selection and multimodel inference: a practical information-theoretic approach , 2nd edn New York, NY : Springer . XXVI + 488 pp ↵ Bushnell , B. ( 2014 ) BBMap: A Fast, Accurate, Splice-aware Aligner (No. LBNL-7065E) Lawrence Berkeley National Lab. (LBNL) , Berkeley, CA ↵ Calatayud , J. , Hórreo , J.L. , Madrigal-Gonzalez , J. , Migeon , A. , Rodriguez , M.A. , Magalhaes , S. , Hortal , J ( 2016 ) Geography and major host evolutionary transitions shape the resource use of plant parasites Proceedings of the National Academy of Sciences 113 : 9840 – 9845 OpenUrl Abstract / FREE Full Text ↵ Calonje , M. , Meerow , A.W. , Griffith , M.P. , Salas-Leiva , D. , Vovides , A.P. , Coiro , M. , Francisco-Ortega , J ( 2019 ) A time-calibrated species tree phylogeny of the New World cycad genus Zamia L. (Zamiaceae, Cycadales) International Journal of Plant Sciences 180 : 286 – 314 OpenUrl ↵ Camacho , C. , Coulouris , G. , Avagyan , V. , Ma , N. , Papadopoulos , J. , Bealer , K. , Madden , T.L ( 2009 ) BLAST+: architecture and applications BMC Bioinformatics 10 : 421 doi: 10.1186/1471-2105-10-421 OpenUrl CrossRef PubMed ↵ Condamine , F. L. , Nagalingum , N. S. , Marshall , C. R ., Morlon , H ( 2015 ) Origin and diversification of living cycads: a cautionary tale on the impact of the branching process prior in Bayesian molecular dating BMC Evolutionary Biology 15 : 1 – 18 OpenUrl CrossRef ↵ Crisp , M. D. , Cook , L.G ( 2011 ) Cenozoic extinctions account for the low diversity of extant gymnosperms compared with angiosperms New Phytologist 192 : 997 – 1009 OpenUrl CrossRef PubMed Web of Science ↵ Drummond , A.J. , Ho , S.Y.W. , Phillips , M.J. , Rambaut , A ( 2006 ) Relaxed phylogenetics and dating with confidence PLoS Biology 4 : 699 – 710 doi: 10.1371/journal.pbio.0040088 OpenUrl CrossRef Web of Science ↵ Drummond , A.J. , Suchard , M.A. , Xie , D. , Rambaut , A ( 2012 ) Bayesian phylogenetics with BEAUti and the BEAST 1.7 Molecular Biology and Evolution 29 : 1969 – 1973 doi: 10.1093/molbev/mss075 OpenUrl CrossRef PubMed Web of Science ↵ Edgar , R.C. ( 2010 ) Search and clustering orders of magnitude faster than BLAST Bioinformatics 26 : 2460 – 2461 doi: 10.1093/bioinformatics/btq461 OpenUrl CrossRef PubMed Web of Science ↵ Ehrlich , P.R. , Raven , P.H ( 1964 ) Butterflies and plants: a study in coevolution Evolution 586 – 608 . ↵ Farrell , B.D. ( 1998 ) “Inordinate Fondness” explained: why are there So many beetles? Science 281 : 555 – 559 OpenUrl Abstract / FREE Full Text ↵ Farrell , B.D. , Mitter , C. , Futuyma , D.J ( 1992 ) Diversification at the insect-plant Interface Bioscience 42 : 34 – 42 OpenUrl CrossRef ↵ Farrell , B.D. , Sequeira , A.S ( 2004 ) Evolutionary rates in the adaptive radiation of beetles on plants Evolution 58 : 1984 – 2001 OpenUrl PubMed Web of Science ↵ Ferrer , M.S. , Marvaldi , A.E. , Sato , H.A. , Gonzalez , A.M ( 2011 ) Biological notes on two species of Oxycorynus (Coleoptera: Belidae) associated with parasitic plants of the genus Lophophytum (Balanophoraceae), and new distribution records in Argentina . Revista de la Sociedad Entomológica Argentina 70 : 351 – 355 OpenUrl ↵ Filatov , D.A. ( 2002 ) ProSeq: a software for preparation and evolutionary analysis of DNA sequence data sets Molecular Ecology Notes 2 : 621 – 624 OpenUrl CrossRef Web of Science ↵ Gillespie , J.J. , Johnston , J.S. , Cannone , J.J. , Gutell , R.R ( 2006 ) Characteristics of the nuclear (18S, 5.8 S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements Insect Molecular Biology 15 : 657 – 686 OpenUrl CrossRef PubMed Web of Science ↵ Goodman , K.R. , Welter , S.C. , Roderick , G.K ( 2012 ) Genetic divergence is decoupled from ecological diversification in the Hawaiian Nesosydne planthoppers Evolution 66 : 2798 – 2814 OpenUrl ↵ Gratshev , V.V. , Legalov , A.A. ( 2014 ) Mesozoic stage of evolution of the family Nemonychidae (Coleoptera, Curculionoidea ) Paleontological Journal 48 : 851 – 944 OpenUrl ↵ Gressitt , J.L. ( 1975 ) Evolution of the endemic Hawaiian cerambycid-beetles International Biological Program, Technical Report No. 71. U. S . 1 – 46 ↵ Guindon , S. , Dufayard , J.-F. , Lefort , V. , Anisimova , M. , Hordijk , W. , Gascuel , O ( 2010 ) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0 . Systematic Biology 59 : 307 – 321 doi: 10.1093/sysbio/syq010 OpenUrl CrossRef PubMed Web of Science ↵ Haddad , S. , Shin , S. , Lemmon , A.R. , Lemmon , E.M. , Svacha , P. , Farrell , B. , Ślipiński , A.D. , Windsor , D. , McKenna , D.D ( 2018 ) Anchored hybrid enrichment provides new insights into the phylogeny and evolution of longhorned beetles (Cerambycidae) Systematic Entomology 43 : 68 – 89 doi: 10.1111/syen.12257 OpenUrl CrossRef ↵ Haines , W.P. , Schmitz , P. , Rubinoff , D ( 2014 ) Ancient diversification of Hyposmocoma moths in Hawaii Nature Communications 5 :p. 3502 OpenUrl ↵ Hoang , D.T. , Chernomor , O. , von Haeseler , A. , Minh , B.Q. , Vinh , L.S. ( 2018 ) UFBoot2: improving the ultrafast bootstrap approximation Molecular Biology and Evolution 35 : 518 – 522 doi: 10.1093/molbev/msx281 OpenUrl CrossRef PubMed ↵ Hsiao , Y. , Oberprieler , R.G. , Zwick , A. , Zhou , Y.L. and Ślipiński , A. ( 2023 ) Museomics unveil systematics, diversity and evolution of Australian cycad-pollinating weevils Proceedings of the Royal Society B 290 :p. 20231385 OpenUrl ↵ Iturralde-Vinent , M.A. , MacPhee , R.D.E ( 1996 ) Age and paleogeographical origin of Dominican amber Science 273 : 1850 – 1852 OpenUrl Abstract / FREE Full Text ↵ Iturralde-Vinent , M.A. , MacPhee , R.D.E ( 2019 ) Remarks on the age of Dominican amber Palaeoentomology 2 : 236 – 240 OpenUrl ↵ Jeffroy , O. , Brinkmann , H. , Delsuc , F. , Philippe , H ( 2006 ) Phylogenomics: the beginning of incongruence? Trends in Genetics 22 : 225 – 231 doi: 10.1016/j.tig.2006.02.003 OpenUrl CrossRef PubMed Web of Science ↵ Johns , C.A. , Toussaint , E.F. , Breinholt , J.W. and Kawahara , A.Y ( 2018 ) Origin and macroevolution of micro-moths on sunken Hawaiian Islands Proceedings of the Royal Society B 285 :p. 20181047 OpenUrl CrossRef PubMed ↵ Kalyaanamoorthy , S. , Minh , B.Q. , Wong , T.K. , von Haeseler , A. , Jermiin , L.S ( 2017 ) ModelFinder: fast model selection for accurate phylogenetic estimates Nature Methods 14 : 587 – 589 doi: 10.1038/nmeth.4285 OpenUrl CrossRef PubMed ↵ Katoh , K. , Standley , D.M ( 2013 ) MAFFT multiple sequence alignment software version 7: improvements in performance and usability Molecular Biology and Evolution 30 : 772 – 780 doi: 10.1093/molbev/mst010 OpenUrl CrossRef PubMed Web of Science ↵ Kershaw , P. , Wagstaff , B ( 2001 ) The southern conifer family Araucariaceae: history, status, and value for paleoenvironmental reconstruction Annual Review of Ecology and Systematics 32 : 397 – 414 OpenUrl CrossRef Web of Science ↵ Kumar , S. , Filipski , A.J. , Battistuzzi , F.U. , Pond , S.L.K. , Tamura , K ( 2012 ) Statistics and truth in phylogenomics Molecular Biology and Evolution 29 : 457 – 472 doi: 10.1093/molbev/msr202 OpenUrl CrossRef PubMed Web of Science ↵ Kuschel , G. ( 1959 ) Nemonychidae, Belidae y Oxycorynidae de la fauna chilena, con algunas consideraciónes biogeográficas (Coleoptera Curculionoidea, aporte 28) Investigaciones zoológicas chilenas 5 : 229 – 271 OpenUrl ↵ Kuschel , G. ( 1995 ) A phylogenetic classification of Curculionoidea to families and subfamilies Memoirs of the Entomological Society of Washington 14 : 5 – 33 OpenUrl ↵ Kuschel , G. , Leschen , R.A.B. ( 2003 ) Appendix 1. Phylogenetic relationships of the genera of Belinae . In ‘Nemonychidae, Belidae, Brentidae (Insecta: Coleoptera: Curculionoidea)’ . (Author: G. Kuschel .) pp. 48 – 55 . Fauna of New Zealand, Volume 45 ( Manaaki Whenua Press : Lincoln, New Zealand .) OpenUrl ↵ Kück , P. , Longo , G.C ( 2014 ) FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studies Frontiers in Zoology 11 : 81 doi: 10.1186/s12983-014-0081-x OpenUrl CrossRef PubMed ↵ Labandeira , C.C. , Currano , E.D ( 2013 ) The fossil record of plant-insect dynamics Annual Review of Earth and Planetary Sciences 41 : 287 – 311 OpenUrl CrossRef GeoRef Web of Science ↵ Landis , M.J. , Matzke , N.J. , Moore , B.R. , Huelsenbeck , J.P ( 2013 ) Bayesian analysis of biogeography when the number of areas is large Systematic Biology 62 : 789 – 804 doi: 10.1093/sysbio/syt040 OpenUrl CrossRef PubMed ↵ Lanfear , R. , Frandsen , P.B. , Wright , A.M. , Senfeld , T. , Calcott , B. ( 2017 ) PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses Molecular Biology and Evolution 34 , 772 – 773 OpenUrl CrossRef PubMed ↵ Larsson A. ( 2014 ) AliView: A fast and lightweight alignment viewer and editor for large datasets Bioinformatics 30 : 3276 – 3278 doi: 10.1093/bioinformatics/btu531 OpenUrl CrossRef PubMed ↵ Lawver , L.A. , Gahagan , L.M. , Coffin , M.F ( 1992 ) The development of paleoseaways around Antarctica Anarctic Research Series vol. 56 : 7 – 30 OpenUrl ↵ Legalov , A.A. ( 2009 ) Contribution to the knowledge of the Mesozoic Curculionoidea (Coleoptera) Amurian Zoological Journal 1 : 283 – 295 OpenUrl ↵ Lemmon , A.R. , Emme , S.A. , Lemmon , E.M ( 2012 ) Anchored hybrid enrichment for massively high-throughput phylogenomics Systematic Biology 61 : 727 – 744 doi: 10.1093/sysbio/sys049 OpenUrl CrossRef PubMed ↵ Lewis , G.P. , Schrire , B. , MacKinder , B. Lock , M Lewis , G.P. ( 2005 ) Tribe Acacieae . pp. 187 – 192 In Lewis , G.P. , Schrire , B. , MacKinder , B. & Lock , M . (eds.), Legumes of the world . Kew : Royal Botanic Gardens . ↵ Li , D. , Liu , C.M. , Luo , R. , Sadakane , K. and Lam , T.W ( 2015 ) MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph Bioinformatics 31 : 1674 – 1676 OpenUrl CrossRef PubMed ↵ Li , H.T. , Yi , T.S. , Gao , L.M. , Ma , P.F. , Zhang , T. , Yang , J.B. , Gitzendanner , M.A. , Fritsch , P.W. , Cai , J. , Luo , Y. , Wang , H ( 2019 ) Origin of angiosperms and the puzzle of the Jurassic gap Nature Plants 5 : 461 – 470 OpenUrl ↵ Li , X. , Hamilton , C.A. , St Laurent , R. , Ballesteros-Mejia , L. , Markee , A. , Haxaire , J. , Rougerie , R. , Kitching , I.J. , Kawahara , A.Y ( 2022b ) A diversification relay race from Caribbean-Mesoamerica to the Andes: historical biogeography of Xylophanes hawkmoths Proceedings of the Royal Society B 289 : 20212435 doi: 10.1098/rspb.2021.2435 OpenUrl CrossRef ↵ Li , X. , St Laurent , R. , Earl , C. , Doorenweerd , C. , van Nieukerken , E.J. , Davis , D.R. , Johns , C.A. , Kawakita , A. , Kobayashi , S. , Zwick , A. , Lopez-Vaamonde , C ( 2022a ) Phylogeny of gracillariid leafumining moths: evolution of larval behaviour inferred from phylogenomic and Sanger data Cladistics 38 : 277 – 300 OpenUrl ↵ Liu , Y. , Wang , S. , Li , L. , Yang , T. , Dong , S. , Wei , T. , Wu , S. , Liu , Y. , Gong , Y. , Feng , X. , Ma , J ( 2022 ) The Cycas genome and the early evolution of seed plants Nature Plants 8 : 389 – 401 OpenUrl ↵ Martin , H.A. ( 2006 ) Cenozoic climatic changes and the development of the arid zone vegetation of Australia Journal of Arid Environments 66 : 533 – 563 OpenUrl CrossRef Web of Science ↵ Marvaldi , A.E. ( 2005 ) Larval morphology and biology of oxycorynine weevils and the higher phylogeny of Belidae (Coleoptera, Curculionoidea) Zoologica Scripta 34 : 37 – 48 OpenUrl CrossRef Web of Science ↵ Marvaldi , A.E. , Del Río , M.G. , Pereyra , V.A. , Rocamundi , N. , Lanteri , A.A. ( 2018 ) A combined molecular and morphological approach to explore the higher phylogeny of entimine weevils (Coleoptera: Curculionidae), with special reference to South American taxa Diversity 10 : 95 doi: 10.3390/d10030095 OpenUrl CrossRef ↵ Marvaldi , A.E. , Duckett , C.N. , Kjer , K.M. , Gillespie , J.J ( 2009 ) Structural alignment of 18S and 28S rDNA sequences provides insights into phylogeny of Phytophaga (Coleoptera: Curculionoidea and Chrysomeloidea) Zoologica Scripta 38 : 63 – 77 OpenUrl CrossRef Web of Science ↵ Marvaldi , A.E. , Ferrer , M.S ( 2014 ) Belidae Schoenherr , 1826 . Handbook of Zoology: Arthropoda. Insecta . ↵ Marvaldi , A.E. , Oberprieler , R.G. , Lyal , C.H.C. , Bradbury , T. , Anderson , R.S ( 2006 ) Phylogeny of the Oxycoryninae sensu lato (Coleoptera: Belidae) and evolution of host-plant associations Invertebrate Systematics 20 : 447 – 476 OpenUrl ↵ Marvaldi , A.E. , Sequeira , A.S. , O’Brien , C.W. , Farrell , B.D ( 2002 ) Molecular and morphological phylogenetics of weevils (Coleoptera, Curculionoidea): do niche shifts accompany diversification? Systematic Biology 51 : 761 – 785 OpenUrl CrossRef PubMed Web of Science ↵ Matsunaga , K.K. , Smith , S.Y ( 2021 ) Fossil palm reading: using fruits to reveal the deep roots of palm diversity American Journal of Botany 108 : 472 – 494 OpenUrl CrossRef ↵ Maslin , B.R. ( 2004 ) Classification and phylogeny of Acacia . In Evolution of ecological and behavioural diversity: Australian Acacia thrips as model organisms . Edited by Crespi BJ , Morris DC , Mound LA. Canberra : Australian Biological Resources Study and Australian National Insect Collection , CSIRO 97 – 112 ↵ Matzke , N.J. ( 2014 ) Model selection in historical biogeography reveals that founder-event speciation is a crucial process in island clades Systematic Biology 63 : 951 – 970 doi: 10.1093/sysbio/syu056 OpenUrl CrossRef PubMed ↵ McIntyre , S.R. , Lineweaver , C.H. , Groves , C.P. , Chopra , A ( 2017 ) Global biogeography since Pangaea Proceedings of the Royal Society B: Biological Sciences 284 : 20170716 OpenUrl CrossRef PubMed ↵ McKenna , D.D. ( 2011 ) Temporal lags and overlap in the diversification of weevils and flowering plants: Recent advances and prospects for additional resolution American Entomologist 57 : 54 – 55 OpenUrl CrossRef ↵ McKenna , D.D. , Clarke , D.J. , Anderson , R. , Astrin , J.J. , Brown , S. , Chamorro , L. , Davis , S.R. , De Medeiros , B. , Del Rio , M.G. , Haran , J. , Kuschel , G ( 2018 ) Morphological and molecular perspectives on the phylogeny, evolution, and classification of weevils (Coleoptera: Curculionoidea): Proceedings from the 2016 International Weevil Meeting Diversity 10 : 64 doi: 10.3390/d10030064 OpenUrl CrossRef ↵ McKenna , D.D. , Farrell , B.D. , Caterino , M.S. , Farnum , C.W. , Hawks , D.C. , Maddison , D.R. , Seago , A.E. , Short , A.E. , Newton , A.F. and Thayer , M.K ( 2015 ) Phylogeny and evolution of Staphyliniformia and Scarabaeiformia: forest litter as a stepping stone for diversification of nonphytophagous beetles Systematic Entomology 40 : 35 – 60 . OpenUrl ↵ McKenna , A. , Hanna , M. , Banks , E. , Sivachenko , A. , Cibulskis , K. , Kernytsky , A. , Garimella , K. , Altshuler , D. , Gabriel , S. , Daly , M. and DePristo , M.A. , 2010 . The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data Genome Research 20 : 1297 – 1303 OpenUrl Abstract / FREE Full Text ↵ McKenna , D.D. , Sequeira , A.S. , Marvaldi , A.E. , Farrell , B.D ( 2009 ) Temporal lags and overlap in the diversification of weevils and flowering plants Proceedings of the National Academy of Sciences 106 : 7083 – 7088 OpenUrl Abstract / FREE Full Text ↵ McKenna , D.D. , Shin , S. , Ahrens , D. , Balke , M. , Beza-Beza , C. , Clarke , D.J. , Donath , A. , Escalona , H.E. , Friedrich , F. , Letsch , H. , Liu , S ( 2019 ) The evolution and genomic basis of beetle diversity Proceedings of the National Academy of Sciences 116 : 24729 – 24737 OpenUrl Abstract / FREE Full Text ↵ McLeish , M.J. , Chapman , T.W. , Schwarz , M.P ( 2007 ) Host-driven diversification of gall-inducing Acacia thrips and the aridification of Australia BMC Biology 5 : 1 – 13 OpenUrl ↵ McLeish , M.J. , Miller , J.T. , Mound , L.A ( 2013 ) Delayed colonisation of Acacia by thrips and the timing of host-conservatism and behavioural specialisation BMC Evolutionary Biology 13 : 1 – 15 OpenUrl CrossRef ↵ Medeiros , M.J. , Bianchi , G.L. , Taschetta , L.R. , Oboyski , P.T ( 2016 ) A review of Polynesian Carposina Herrich-Schäffer (Lepidoptera: Carposinidae), with descriptions of four new species Zoological Journal of the Linnean Society 177 : 135 – 146 OpenUrl ↵ Miller , J.T. , Murphy , D.J. , Ho , S.Y. , Cantrill , D.J ., Seigler , D ( 2013 ) Comparative dating of Acacia : combining fossils and multiple phylogenies to infer ages of clades with poor fossil records Australian Journal of Botany 61 : 436 – 445 OpenUrl ↵ Minh , B.Q. , Hahn , M.W. , Lanfear , R ( 2020 ) New Methods to Calculate Concordance Factors for Phylogenomic Datasets Molecular Biology and Evolution 37 : 2727 – 2733 doi: 10.1093/molbev/msaa106 OpenUrl CrossRef PubMed ↵ Minh , B.Q. , Nguyen , M.A.T. , von Haeseler , A ( 2013 ) Ultrafast approximation for phylogenetic bootstrap Molecular Biology and Evolution 30 : 1188 – 1195 doi: 10.1093/molbev/mst024 OpenUrl CrossRef PubMed Web of Science ↵ Misof , B. , Liu , S. , Meusemann , K. , Peters , R.S. , Donath , A. , Mayer , C. , Frandsen , P.B. , Ware , J. , Flouri , T. , Beutel , R.G. , Niehuis , O ( 2014 ) Phylogenomics resolves the timing and pattern of insect evolution Science 346 : 763 – 767 OpenUrl Abstract / FREE Full Text ↵ Nagalingum , N.S. , Marshall , C.R. , Quental , T.B. , Rai , H.S. , Little , D.P. , Mathews , S ( 2011 ) Recent synchronous radiation of a living fossil Science 334 : 796 – 799 OpenUrl Abstract / FREE Full Text ↵ Naumann , J. , Salomo , K. , Der , J.P. , Wafula , E.K. , Bolin , J.F. , Maass , E. , Frenzke , L. , Samain , M.S. , Neinhuis , C. , DePamphilis , C.W. , Wanke , S ( 2013 ) Single-copy nuclear genes place haustorial Hydnoraceae within Piperales and reveal a Cretaceous origin of multiple parasitic angiosperm lineages PLoS One 8 : e79204 OpenUrl CrossRef PubMed ↵ Nguyen , L.T. , Schmidt , H.A. , von Haeseler , A. , Minh , B.Q ( 2015 ) IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies Molecular Biology and Evolution 32 : 268 – 274 doi: 10.1093/molbev/msu300 OpenUrl CrossRef PubMed ↵ Oberprieler , R. , Marvaldi , A. , Anderson , R.S ( 2007 ) Weevils, weevils, weevils everywhere . In: Linnaeus tercentenary: progress in invertebrate taxonomy Zootaxa 1668 : 491 – 520 OpenUrl ↵ Oberprieler , R.G. , Oberprieler , S.K ( 2012 ) Talbragarus averyi gen. et sp. n., the first Jurassic weevil from the southern hemisphere (Coleoptera: Curculionoidea: Nemonychidae) Zootaxa 3478 : 256 – 266 OpenUrl ↵ O’Brien , C. W. , Tang , W ( 2015 ) Revision of the New World cycad weevils of the subtribe Allocorynina, with description of two new genera and three new subgenera (Coleoptera: Belidae: Oxycoryninae) Zootaxa 3970 : 1 – 87 OpenUrl ↵ O’Grady , P.M. , Lapoint , R.T. , Bonacum , J. , Lasola , J. , Owen , E. , Wu , Y. , DeSalle , R ( 2011 ) Phylogenetic and ecological relationships of the Hawaiian Drosophila inferred by mitochondrial DNA analysis Molecular Phylogenetics and Evolution 58 : 244 – 256 OpenUrl CrossRef PubMed ↵ Parham , J.F. , Donoghue , P.C. , Bell , C.J. , Calway , T.D. , Head , J.J. , Holroyd , P.A. , Inoue , J.G. , Irmis , R.B. , Joyce , W.G. , Ksepka , D.T. , Patané , J.S ( 2012 ) Best practices for justifying fossil calibrations Systematic Biology 61 : 346 – 59 OpenUrl CrossRef PubMed ↵ Pease , J.B. , Brown , J.W. , Walker , J.F. , Hinchliff , C.E. , Smith , S.A. ( 2018 ) Quartet Sampling distinguishes lack of support from conflicting support in the green plant tree of life American Journal of Botany 105 : 385 – 403 doi: 10.1002/ajb2.1016 OpenUrl CrossRef ↵ Poinar , G.O. , Legalov , A.A ( 2014 ) Pleurambus strongylus n. gen., n. sp. (Coleoptera: Belidae) in Dominican amber Historical Biology 26 : 670 – 674 OpenUrl Price , J.P. , Clague , D.A ( 2002 ) How old is the Hawaiian biota? Geology and phylogeny suggest recent divergence Proceedings of the Royal Society of London. Series B: Biological Sciences 269 : 2429 – 2435 doi: 10.1098/rspb.2002.2175 OpenUrl CrossRef GeoRef PubMed Web of Science ↵ Rambaut , A. , Drummond , A.J ( 2015 ) TreeAnnotator v1.8.2 MCMC Output Analysis . ↵ Rambaut , A. , Drummond , A.J. , Xie , D. , Baele , G. , Suchard , A ( 2018 ) Posterior summarization in Bayesian phylogenetics using Tracer 1.7 Systematic Biology 67 : 901 – 904 OpenUrl CrossRef PubMed ↵ Ree , R.H. , Smith , S.A ( 2008 ) Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis Systematic Biology 57 : 4 – 14 doi: 10.1080/10635150701883881 OpenUrl CrossRef GeoRef PubMed Web of Science ↵ Revell , L.J. ( 2012 ) Phytools: an R package for phylogenetic comparative biology (and other things) Methods in Ecology and Evolution 3 : 217 – 223 OpenUrl ↵ Ronquist , F. ( 1997 ) Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography Systematic Biology 46 : 195 – 203 doi: 10.1093/sysbio/46.1.195 OpenUrl CrossRef Web of Science ↵ Salas-Leiva , D.E. , Meerow , A.W. , Calonje , M. , Griffith , M.P. , Francisco-Ortega , J. , Nakamura , K. , Stevenson , D.W. , Lewis , C.E. , Namoff , S ( 2013 ) Phylogeny of the cycads based on multiple single-copy nuclear genes: congruence of concatenated parsimony, likelihood and species tree inference methods Annals of Botany 112 : 1263 – 1278 OpenUrl CrossRef PubMed ↵ Salzman , S. , Crook , D. , Crall , J.D. , Hopkins , R. , Pierce , N.E ( 2020 ) An ancient push-pull pollination mechanism in cycads Science Advances 6 :p. eaay6169 OpenUrl FREE Full Text ↵ Santos , M.F.A. , Mermudes , J.R.M. , Fonseca , V.M.M. ( 2007 ) Descrição de um novo gênero e espécie de Belinae (Belidae, Curculionoidea, Coleoptera) da Formação Santana (Membro Crato, Cretáceo Inferior) da Bacia do Araripe, Nordeste do Brasil Paleontologia: Cenários de Vida 441 – 447 ↵ Santos , M.F.D.A. , Mermudes , J.R.M. , Fonseca , V.M.M.D. ( 2011 ) A specimen of Curculioninae (Curculionidae, Coleoptera) from the Lower Cretaceous, Araripe Basin, northueastern Brazil Palaeontology 54 : 807 – 814 OpenUrl GeoRef ↵ Sequeira , A.S. , Normark , B.B. and Farrell , B.D ( 2000 ) Evolutionary assembly of the conifer fauna: distinguishing ancient from recent associations in bark beetles Proceedings of the Royal Society of London. Series B: Biological Sciences 267 : 2359 – 2366 OpenUrl CrossRef GeoRef PubMed ↵ Schneider , H. , Schuettpelz , E. , Pryer , K. M. , Cranfill , R. , Magallón , S. , Lupia , R ( 2004 ) Ferns diversified in the shadow of angiosperms Nature 428 : 553 – 557 OpenUrl CrossRef PubMed Web of Science ↵ Schuettpelz , E. , Pryer , K.M ( 2009 ) Evidence for a Cenozoic radiation of ferns in an angiosperm-dominated canopy Proceedings of the National Academy of Sciences 106 : 11200 – 11205 OpenUrl Abstract / FREE Full Text ↵ Shin , S. , Clarke , D.J. , Lemmon , A.R. , Moriarty Lemmon , E. , Aitken , A.L. , Haddad , S. , Farrell , B.D. , Marvaldi , A.E. , Oberprieler , R.G. , McKenna , D.D ( 2018 ) Phylogenomic data yield new and robust insights into the phylogeny and evolution of weevils Molecular Biology and Evolution 35 : 823 – 836 OpenUrl ↵ Simon , A. , Salzman , S. , Specht , C.D. , Raguso , . A ( 2023 ) Behavior and feeding of two beetle pollinators of Zamia integrifolia (Cycadales): Rhopalotria slossoni (Coleoptera: Belidae) and Pharaxanotha floridana (Coleoptera: Erotylidae) Florida Entomologist 106 : 199 – 202 OpenUrl ↵ Smith , S.A. , Brown , J.W. , Walker , J.F ( 2018 ) So many genes, so little time: a practical approach to divergence-time estimation in the genomic era PLoS ONE , 13 : e0197433 . doi: 10.1101/114975 OpenUrl CrossRef ↵ Song , H. , Béthoux , O. , Shin , S. , Donath , A. , Letsch , H. , Liu , S. , McKenna , D.D. , Meng , G. , Misof , B. , Podsiadlowski , L. , Zhou , X ( 2020 ) Phylogenomic analysis sheds light on the evolutionary pathways towards acoustic communication in Orthoptera Nature Communications 11 : 1 – 16 OpenUrl CrossRef ↵ Strimmer , K. , von Haeseler , A ( 1997 ) Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment Proceedings of the National Academy of Sciences USA 94 : 6815 – 6819 doi: 10.1073/pnas.94.13.6815 OpenUrl Abstract / FREE Full Text ↵ Suchard , M.A. , Lemey , P. , Baele , G. , Ayres , D.L. , Drummond , A.J. , Rambaut , A. ( 2018 ) Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10 Virus Evolution 4 : vey016 doi: 10.1093/ve/vey016 OpenUrl CrossRef PubMed Sunnucks , P. , Hales , D.F ( 1996 ) Numerous transposed sequences of mitochondrial cytochrome oxidase I-II in aphids of the genus Sitobion (Hemiptera: Aphididae) Molecular Biology and Evolution 13 : 510 – 524 OpenUrl CrossRef PubMed Web of Science ↵ Swain , A. , Maccracken , S.A. , Fagan , W.F. , Labandeira , C.C ( 2022 ) Understanding the ecology of host plant–insect herbivore interactions in the fossil record through bipartite networks Paleobiology 48 : 239 – 260 . doi: 10.1017/pab.2021.20 OpenUrl CrossRef ↵ Tang , W. ( 1987 ) Insect pollination in the cycad Zamia pumila (Zamiaceae) American Journal of Botany 74 : 90 – 99 doi: 10.1002/j.1537-2197.1987.tb08582.x OpenUrl CrossRef Web of Science ↵ Teasdale , L.C. , Köhler , F. , Murray , K.D. , O’Hara , T. and Moussalli , A ( 2016 ) Identification and qualification of 500 nuclear, single-copy, orthologous genes for the Eupulmonata (Gastropoda) using transcriptome sequencing and exon capture Molecular Ecology Resources 16 : 1107 – 1123 OpenUrl ↵ Vanin , S.A. ( 1976 ) Taxonomic revision of the South American Belidae (Coleoptera) Arquivos de Zoologia 28 : 1 – 75 OpenUrl ↵ Wang , X.Q. , Ran , J.H ( 2014 ) Evolution and biogeography of gymnosperms Molecular Phylogenetics and Evolution 75 : 24 – 40 OpenUrl CrossRef PubMed ↵ Wilf , P. ( 2008 ) Insect-damaged fossil leaves record food web response to ancient climate change and extinction New Phytologist 178 : 486 – 502 OpenUrl CrossRef PubMed Web of Science ↵ Xie , W. , Lewis , P.O. , Fan , Y. , Kuo , L. , Chen , M.H ( 2011 ) Improving marginal likelihood estimation for Bayesian phylogenetic model selection Systematic Biology 60 : 150 – 160 doi: 10.1093/sysbio/syq085 OpenUrl CrossRef PubMed Web of Science ↵ Yu , Y.L. , Davis , S.R. , Shih , C.K ( 2019 ) The earliest fossil record of Belidae and its implications for the early evolution of Curculionoidea (Coleoptera) Journal of Systematic Palaeontology 17 : 2105 – 2117 OpenUrl ↵ Zhang C , Rabiee M , Sayyari E , Mirarab S. ( 2018 ) ASTRAL-III: Polynomial time species tree reconstruction from partially resolved gene trees BMC Bioinformatics 19 : 15 – 30 doi: 10.1186/s12859-018-2129-y OpenUrl CrossRef ↵ Zhang , J. , Gordon , E.R. , Forthman , M. , Hwang , W.S. , Walden , K. , Swanson , D.R. , Johnson , K.P. , Meier , R. , Weirauch , C ( 2016 ) Evolution of the assassin’s arms: insights from a phylogeny of combined transcriptomic and ribosomal DNA data (Heteroptera: Reduvioidea) Scientific Reports 6 : 22177 doi: 10.1038/srep22177 OpenUrl CrossRef ↵ Zimmerman , E.C. ( 1994 ) Australian Weevils (Coleoptera: Curculionoidea) . Volume I . Orthoceri. Anthribidae to Attelabidae. The Primitive Weevils . CSIRO Australia , Melbourne , xxxii + 741 pp OpenUrl View the discussion thread. Back to top Previous Next Posted July 11, 2024. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following The evolutionary history of the ancient weevil family Belidae (Coleoptera: Curculionoidea) reveals the marks of Gondwana breakup and major floristic turnovers, including the rise of angiosperms Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share The evolutionary history of the ancient weevil family Belidae (Coleoptera: Curculionoidea) reveals the marks of Gondwana breakup and major floristic turnovers, including the rise of angiosperms Xuankun Li , Adriana E. Marvaldi , Rolf G. Oberprieler , Dave Clarke , Brian D. Farrell , Andrea Sequeira , M. Silvia Ferrer , Charles O’Brien , Shayla Salzman , Seunggwan Shin , William Tang , Duane D. McKenna bioRxiv 2024.03.12.584611; doi: https://doi.org/10.1101/2024.03.12.584611 Share This Article: Copy Citation Tools The evolutionary history of the ancient weevil family Belidae (Coleoptera: Curculionoidea) reveals the marks of Gondwana breakup and major floristic turnovers, including the rise of angiosperms Xuankun Li , Adriana E. Marvaldi , Rolf G. Oberprieler , Dave Clarke , Brian D. Farrell , Andrea Sequeira , M. Silvia Ferrer , Charles O’Brien , Shayla Salzman , Seunggwan Shin , William Tang , Duane D. McKenna bioRxiv 2024.03.12.584611; doi: https://doi.org/10.1101/2024.03.12.584611 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Evolutionary Biology Subject Areas All Articles Animal Behavior and Cognition (7644) Biochemistry (17728) Bioengineering (13916) Bioinformatics (42037) Biophysics (21489) Cancer Biology (18637) Cell Biology (25553) Clinical Trials (138) Developmental Biology (13401) Ecology (19941) Epidemiology (2067) Evolutionary Biology (24367) Genetics (15622) Genomics (22547) Immunology (17764) Microbiology (40475) Molecular Biology (17208) Neuroscience (88747) Paleontology (667) Pathology (2842) Pharmacology and Toxicology (4834) Physiology (7659) Plant Biology (15175) Scientific Communication and Education (2047) Synthetic Biology (4304) Systems Biology (9835) Zoology (2272)