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.
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. 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). 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).
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.
Supplementary figures
Acknowledgements
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.
References
- New Oxycoryninae from Central and South America: phylogenetic and biogeographical implications (Coleoptera: Belidae)Systematic Entomology 30:644–652https://doi.org/10.1111/j.1365-3113.2005.00304.x
- Brief history of the gymnosperms: classification, biodiversity, phytogeography and ecologySouth African National Biodiversity Institute South Africa: Pretoria
- 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 NicaraguaThe Coleopterists Bulletin :61–71
- Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertaintyMolecular Biology and Evolution 29:2157–2167https://doi.org/10.1093/molbev/mss084
- Accurate model selection of relaxed molecular clocks in Bayesian phylogeneticsMolecular Biology and Evolution 30:239–243https://doi.org/10.1093/molbev/mss243
- Constraining uncertainty in the timescale of angiosperm evolution and the veracity of a Cretaceous Terrestrial RevolutionNew Phytologist 218:819–834
- The leptolepid fish Cavenderichthys talbragarensis (Woodward, 1895) from the Talbragar Fish Bed (Late Jurassic) near Gulgong, New South WalesRecords-Western Australian Museum 23:43–76
- Historical biogeography and ecological opportunity in the adaptive radiation of native Hawaiian leafhoppers (Cicadellidae: Nesophrosyne)Journal of Biogeography 40:1512–1523
- Resolving relationships among the megadiverse butterflies and moths with a novel pipeline for anchored phylogenomicsSystematic Biology 67:78–93https://doi.org/10.1093/sysbio/syx048
- Global biogeographyAmsterdam: Elsevier Science
- Phyx: Phylogenetic tools for unixBioinformatics 33:1886–1888https://doi.org/10.1093/bioinformatics/btx063
- A new species of Proterhinus Sharp, 1878 (Coleoptera: Curculionoidea: Belidae) from Miti’aro, Cook Islands, South PacificZootaxa 4664:481–496
- Birth of a biome: insights into the assembly and maintenance of the Australian arid zone biotaMolecular Ecology 17:4398–4417
- Practical use of the information-theoretic approachModel selection and inference New York, NY: Springer :75–117
- Model selection and multimodel inference: a practical information-theoretic approachNew York, NY: Springer
- BBMap: A Fast, Accurate, Splice-aware Aligner (No. LBNL-7065E) Lawrence Berkeley National Lab. (LBNL)Berkeley, CA
- Geography and major host evolutionary transitions shape the resource use of plant parasitesProceedings of the National Academy of Sciences 113:9840–9845
- A time-calibrated species tree phylogeny of the New World cycad genus Zamia L. (Zamiaceae, Cycadales)International Journal of Plant Sciences 180:286–314
- BLAST+: architecture and applicationsBMC Bioinformatics 10https://doi.org/10.1186/1471-2105-10-421
- Origin and diversification of living cycads: a cautionary tale on the impact of the branching process prior in Bayesian molecular datingBMC Evolutionary Biology 15:1–18
- Cenozoic extinctions account for the low diversity of extant gymnosperms compared with angiospermsNew Phytologist 192:997–1009
- Relaxed phylogenetics and dating with confidencePLoS Biology 4:699–710https://doi.org/10.1371/journal.pbio.0040088
- Bayesian phylogenetics with BEAUti and the BEAST 1.7Molecular Biology and Evolution 29:1969–1973https://doi.org/10.1093/molbev/mss075
- Search and clustering orders of magnitude faster than BLASTBioinformatics 26:2460–2461https://doi.org/10.1093/bioinformatics/btq461
- Butterflies and plants: a study in coevolutionEvolution :586–608
- “Inordinate Fondness” explained: why are there So many beetles?Science 281:555–559
- Diversification at the insect-plant InterfaceBioscience 42:34–42
- Evolutionary rates in the adaptive radiation of beetles on plantsEvolution 58:1984–2001
- Biological notes on two species of Oxycorynus (Coleoptera: Belidae) associated with parasitic plants of the genus Lophophytum (Balanophoraceae), and new distribution records in ArgentinaRevista de la Sociedad Entomológica Argentina 70:351–355
- ProSeq: a software for preparation and evolutionary analysis of DNA sequence data setsMolecular Ecology Notes 2:621–624
- 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 elementsInsect Molecular Biology 15:657–686
- Genetic divergence is decoupled from ecological diversification in the HawaiianNesosydne planthoppers Evolution 66:2798–2814
- Mesozoic stage of evolution of the family Nemonychidae (Coleoptera, Curculionoidea) Paleontological Journal 48:851–944
- Evolution of the endemic Hawaiian cerambycid-beetlesInternational Biological Program, Technical Report No. 71. U. S :1–46
- New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0Systematic Biology 59:307–321https://doi.org/10.1093/sysbio/syq010
- Anchored hybrid enrichment provides new insights into the phylogeny and evolution of longhorned beetles (Cerambycidae)Systematic Entomology 43:68–89https://doi.org/10.1111/syen.12257
- Ancient diversification of Hyposmocoma moths in HawaiiNature Communications 5
- UFBoot2: improving the ultrafast bootstrap approximationMolecular Biology and Evolution 35:518–522https://doi.org/10.1093/molbev/msx281
- Museomics unveil systematics, diversity and evolution of Australian cycad-pollinating weevilsProceedings of the Royal Society B 290
- Age and paleogeographical origin of Dominican amberScience 273:1850–1852
- Remarks on the age of Dominican amberPalaeoentomology 2:236–240
- Phylogenomics: the beginning of incongruence?Trends in Genetics 22:225–231https://doi.org/10.1016/j.tig.2006.02.003
- Origin and macroevolution of micro-moths on sunken Hawaiian IslandsProceedings of the Royal Society B 285
- ModelFinder: fast model selection for accurate phylogenetic estimatesNature Methods 14:587–589https://doi.org/10.1038/nmeth.4285
- MAFFT multiple sequence alignment software version 7: improvements in performance and usabilityMolecular Biology and Evolution 30:772–780https://doi.org/10.1093/molbev/mst010
- The southern conifer family Araucariaceae: history, status, and value for paleoenvironmental reconstructionAnnual Review of Ecology and Systematics 32:397–414
- Statistics and truth in phylogenomicsMolecular Biology and Evolution 29:457–472https://doi.org/10.1093/molbev/msr202
- 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
- A phylogenetic classification of Curculionoidea to families and subfamiliesMemoirs of the Entomological Society of Washington 14:5–33
- Appendix 1. Phylogenetic relationships of the genera of Belinae‘Nemonychidae, Belidae, Brentidae (Insecta: Coleoptera: Curculionoidea)’ Lincoln, New Zealand: Manaaki Whenua Press :48–55
- FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studiesFrontiers in Zoology 11https://doi.org/10.1186/s12983-014-0081-x
- The fossil record of plant-insect dynamicsAnnual Review of Earth and Planetary Sciences 41:287–311
- Bayesian analysis of biogeography when the number of areas is largeSystematic Biology 62:789–804https://doi.org/10.1093/sysbio/syt040
- PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analysesMolecular Biology and Evolution 34:772–773
- AliView: A fast and lightweight alignment viewer and editor for large datasetsBioinformatics 30:3276–3278https://doi.org/10.1093/bioinformatics/btu531
- The development of paleoseaways around AntarcticaAnarctic Research Series 56:7–30
- Contribution to the knowledge of the Mesozoic Curculionoidea (Coleoptera)Amurian Zoological Journal 1:283–295
- Anchored hybrid enrichment for massively high-throughput phylogenomicsSystematic Biology 61:727–744https://doi.org/10.1093/sysbio/sys049
- Legumes of the WorldLondon, United Kingdom: Royal Botanic Gardens :187–192
- MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graphBioinformatics 31:1674–1676
- Origin of angiosperms and the puzzle of the Jurassic gapNature Plants 5:461–470
- A diversification relay race from Caribbean-Mesoamerica to the Andes: historical biogeography of Xylophanes hawkmothsProceedings of the Royal Society B 289https://doi.org/10.1098/rspb.2021.2435
- Phylogeny of gracillariid leafumining moths: evolution of larval behaviour inferred from phylogenomic and Sanger dataCladistics 38:277–300
- The Cycas genome and the early evolution of seed plantsNature Plants 8:389–401
- Cenozoic climatic changes and the development of the arid zone vegetation of AustraliaJournal of Arid Environments 66:533–563
- Larval morphology and biology of oxycorynine weevils and the higher phylogeny of Belidae (Coleoptera, Curculionoidea)Zoologica Scripta 34:37–48
- A combined molecular and morphological approach to explore the higher phylogeny of entimine weevils (Coleoptera: Curculionidae), with special reference to South American taxaDiversity 10https://doi.org/10.3390/d10030095
- Structural alignment of 18S and 28S rDNA sequences provides insights into phylogeny of Phytophaga (Coleoptera: Curculionoidea and Chrysomeloidea)Zoologica Scripta 38:63–77
- Belidae SchoenherrHandbook of Zoology: Arthropoda. Insecta
- Phylogeny of the Oxycoryninae sensu lato (Coleoptera: Belidae) and evolution of host-plant associationsInvertebrate Systematics 20:447–476
- Molecular and morphological phylogenetics of weevils (Coleoptera, Curculionoidea): do niche shifts accompany diversification?Systematic Biology 51:761–785
- Fossil palm reading: using fruits to reveal the deep roots of palm diversityAmerican Journal of Botany 108:472–494
- Classification and phylogeny of AcaciaEvolution of ecological and behavioural diversity: Australian Acacia thrips as model organisms Canberra: Australian Biological Resources Study and Australian National Insect Collection :97–112
- Model selection in historical biogeography reveals that founder-event speciation is a crucial process in island cladesSystematic Biology 63:951–970https://doi.org/10.1093/sysbio/syu056
- Global biogeography since PangaeaProceedings of the Royal Society B: Biological Sciences 284
- Temporal lags and overlap in the diversification of weevils and flowering plants: Recent advances and prospects for additional resolutionAmerican Entomologist 57:54–55
- Morphological and molecular perspectives on the phylogeny, evolution, and classification of weevils (Coleoptera: Curculionoidea): Proceedings from the 2016 International Weevil MeetingDiversity 10https://doi.org/10.3390/d10030064
- Phylogeny and evolution of Staphyliniformia and Scarabaeiformia: forest litter as a stepping stone for diversification of nonphytophagous beetlesSystematic Entomology 40:35–60
- The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing dataGenome Research 20:1297–1303
- Temporal lags and overlap in the diversification of weevils and flowering plantsProceedings of the National Academy of Sciences 106:7083–7088
- The evolution and genomic basis of beetle diversityProceedings of the National Academy of Sciences 116:24729–24737
- Host-driven diversification of gall-inducing Acacia thrips and the aridification of AustraliaBMC Biology 5:1–13
- Delayed colonisation of Acacia by thrips and the timing of host-conservatism and behavioural specialisationBMC Evolutionary Biology 13:1–15
- A review of Polynesian Carposina Herrich-Schäffer (Lepidoptera: Carposinidae), with descriptions of four new speciesZoological Journal of the Linnean Society 177:135–146
- Comparative dating of Acacia: combining fossils and multiple phylogenies to infer ages of clades with poor fossil recordsAustralian Journal of Botany 61:436–445
- New Methods to Calculate Concordance Factors for Phylogenomic DatasetsMolecular Biology and Evolution 37:2727–2733https://doi.org/10.1093/molbev/msaa106
- Ultrafast approximation for phylogenetic bootstrapMolecular Biology and Evolution 30:1188–1195https://doi.org/10.1093/molbev/mst024
- Phylogenomics resolves the timing and pattern of insect evolutionScience 346:763–767
- Recent synchronous radiation of a living fossilScience 334:796–799
- Single-copy nuclear genes place haustorial Hydnoraceae within Piperales and reveal a Cretaceous origin of multiple parasitic angiosperm lineagesPLoS One 8
- IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogeniesMolecular Biology and Evolution 32:268–274https://doi.org/10.1093/molbev/msu300
- Weevils, weevils, weevils everywhereLinnaeus tercentenary: progress in invertebrate taxonomy Zootaxa 1668:491–520
- Talbragarus averyi gen. et sp. n., the first Jurassic weevil from the southern hemisphere (Coleoptera: Curculionoidea: Nemonychidae)Zootaxa 3478:256–266
- 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
- Phylogenetic and ecological relationships of the Hawaiian Drosophila inferred by mitochondrial DNA analysisMolecular Phylogenetics and Evolution 58:244–256
- ) Best practices for justifying fossil calibrationsSystematic Biology 61:346–59
- Quartet Sampling distinguishes lack of support from conflicting support in the green plant tree of lifeAmerican Journal of Botany 105:385–403https://doi.org/10.1002/ajb2.1016
- Pleurambus strongylus n. gen., n. sp. (Coleoptera: Belidae) in Dominican amberHistorical Biology 26:670–674
- How old is the Hawaiian biota? Geology and phylogeny suggest recent divergenceProceedings of the Royal Society of London. Series B: Biological Sciences 269:2429–2435https://doi.org/10.1098/rspb.2002.2175
- TreeAnnotator v1.8.2MCMC Output Analysis
- Posterior summarization in Bayesian phylogenetics using Tracer 1.7Systematic Biology 67:901–904
- Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesisSystematic Biology 57:4–14https://doi.org/10.1080/10635150701883881
- Phytools: an R package for phylogenetic comparative biology (and other things)Methods in Ecology and Evolution 3:217–223
- Dispersal-vicariance analysis: a new approach to the quantification of historical biogeographySystematic Biology 46:195–203https://doi.org/10.1093/sysbio/46.1.195
- Phylogeny of the cycads based on multiple single-copy nuclear genes: congruence of concatenated parsimony, likelihood and species tree inference methodsAnnals of Botany 112:1263–1278
- An ancient push-pull pollination mechanism in cycadsScience Advances 6
- 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 BrasilPaleontologia: Cenários de Vida :441–447
- A specimen of Curculioninae (Curculionidae, Coleoptera) from the Lower Cretaceous, Araripe Basin, northueastern BrazilPalaeontology 54:807–814
- Evolutionary assembly of the conifer fauna: distinguishing ancient from recent associations in bark beetlesProceedings of the Royal Society of London. Series B: Biological Sciences 267:2359–2366
- Ferns diversified in the shadow of angiospermsNature 428:553–557
- Evidence for a Cenozoic radiation of ferns in an angiosperm-dominated canopyProceedings of the National Academy of Sciences 106:11200–11205
- Phylogenomic data yield new and robust insights into the phylogeny and evolution of weevilsMolecular Biology and Evolution 35:823–836
- 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
- So many genes, so little time: a practical approach to divergence-time estimation in the genomic eraPLoS ONE 13https://doi.org/10.1101/114975
- Phylogenomic analysis sheds light on the evolutionary pathways towards acoustic communication in OrthopteraNature Communications 11:1–16
- Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignmentProceedings of the National Academy of Sciences USA 94:6815–6819https://doi.org/10.1073/pnas.94.13.6815
- Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10Virus Evolution 4https://doi.org/10.1093/ve/vey016
- Numerous transposed sequences of mitochondrial cytochrome oxidase I-II in aphids of the genus Sitobion (Hemiptera: Aphididae)Molecular Biology and Evolution 13:510–524
- Understanding the ecology of host plant–insect herbivore interactions in the fossil record through bipartite networksPaleobiology 48:239–260https://doi.org/10.1017/pab.2021.20
- Insect pollination in the cycad Zamia pumila (Zamiaceae)American Journal of Botany 74:90–99https://doi.org/10.1002/j.1537-2197.1987.tb08582.x
- Identification and qualification of 500 nuclear, single-copy, orthologous genes for the Eupulmonata (Gastropoda) using transcriptome sequencing and exon captureMolecular Ecology Resources 16:1107–1123
- Taxonomic revision of the South American Belidae (Coleoptera)Arquivos de Zoologia 28:1–75
- Evolution and biogeography of gymnospermsMolecular Phylogenetics and Evolution 75:24–40
- Insect-damaged fossil leaves record food web response to ancient climate change and extinctionNew Phytologist 178:486–502
- Improving marginal likelihood estimation for Bayesian phylogenetic model selectionSystematic Biology 60:150–160https://doi.org/10.1093/sysbio/syq085
- The earliest fossil record of Belidae and its implications for the early evolution of Curculionoidea (Coleoptera)Journal of Systematic Palaeontology 17:2105–2117
- ASTRAL-III: Polynomial time species tree reconstruction from partially resolved gene treesBMC Bioinformatics 19:15–30https://doi.org/10.1186/s12859-018-2129-y
- Evolution of the assassin’s arms: insights from a phylogeny of combined transcriptomic and ribosomal DNA data (Heteroptera: Reduvioidea)Scientific Reports 6https://doi.org/10.1038/srep22177
- Australian Weevils (Coleoptera: Curculionoidea)Orthoceri. Anthribidae to Attelabidae. The Primitive Weevils Melbourne: CSIRO Australia
Article and author information
Author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
Copyright
© 2024, Li et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 324
- downloads
- 24
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.