Next Article in Journal
Pinyon Engraver Beetle Acoustics: Stridulation Apparatus, Sound Production and Behavioral Response to Vibroacoustic Treatments in Logs
Previous Article in Journal
Species Diversity in the Leafhopper Genus Batracomorphus Lewis (Hemiptera: Cicadellidae: Iassinae) from Qinling Mountain in Shaanxi
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mitochondrial Phylogenomics of Tenthredinidae (Hymenoptera: Tenthredinoidea) Supports the Monophyly of Megabelesesinae as a Subfamily

1
College of Life Sciences, Jiangxi Normal University, Nanchang 330022, China
2
College of Forestry, Beijing Forestry University, Beijing 100083, China
3
Department of Molecular Biology and Genetics, Faculty of Science, Sivas Cumhuriyet University, Sivas 58140, Turkey
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Insects 2021, 12(6), 495; https://doi.org/10.3390/insects12060495
Submission received: 2 April 2021 / Revised: 17 May 2021 / Accepted: 18 May 2021 / Published: 26 May 2021

Abstract

:

Simple Summary

Tenthredinidae is the most speciose family of the paraphyletic ancestral grade Symphyta, including mainly phytophagous lineages. The subfamilial classification of this family has long been problematic with respect to their monophyly and/or phylogenetic placements. This article reports four complete sawfly mitogenomes of Cladiucha punctata, C. magnoliae, Megabeleses magnoliae, and M. liriodendrovorax for the first time. To investigate the mitogenome characteristics of Tenthredinidae, we also compare them with the previously reported tenthredinid mitogenomes. To explore the phylogenetic placements of these four species within this ecologically and economically important sawfly family, we perform the mitophylogenomics reconstruction and divergence time estimation analyses using a mitogenome dataset of 45 species of the superfamily Tenthredinoidea. Although these newly reported four species have currently classified in the subfamily of Allantinae, the obtained tree topology reveals the sister-group placement of Cladiucha and Megabeleses outside of Blennocampinae, Heterarthrinae, Tenthredininae, and the rest of Allantinae as a distinct lineage. In conjunction with the occurrence of morphological and molecular synapomorphic characters, apparent matches the reported divergence times of their host plants from Magnoliaceae—supporting the proposal of Megabelesesinae as a subfamily. This study represents a broad framework and valuable information for future research on this small subfamily.

Abstract

Tenthredinidae represents one of the external feeders of the most diverse superfamily, Tenthredinoidea, with diverse host plant utilization. In this study, four complete mitochondrial genomes (mitogenomes), those of Cladiucha punctata, Cladiucha magnoliae, Megabeleses magnoliae, and Megabeleses liriodendrovorax, are newly sequenced and comparatively analyzed with previously reported tenthredinid mitogenomes. The close investigation of mitogenomes and the phylogeny of Tenthredinidae leads us to the following conclusions: The subfamilial relationships and phylogenetic placements within Tenthredinidae are mostly found to be similar to the previously suggested phylogenies. However, the present phylogeny supports the monophyly of Megabelesesinae as a subfamily, with the sister-group placement of Cladiucha and Megabeleses outside of Allantinae. The occurrence of the same type of tRNA rearrangements (MQI and ANS1ERF) in the mitogenomes of Megabelesesinae species and the presence of apomorphic morphological characters also provide robust evidence for this new subfamily. The divergence and diversification times of the subfamilies appear to be directly related to colonization of the flowering plants following the Early Cretaceous. The origin time and diversification patterns of Megabelesesinae were also well matched with the divergence times of their host plants from Magnoliaceae.

1. Introduction

The exclusively herbivorous larvae of Tenthredinidae are mostly simple external feeders that feed on a wide variety of leaves, but a few are leaf miners, gall causers, or stem borers [1]. Some species of tenthredinids are economically essential pests on both deciduous and resinous trees as forest defoliators; they feed on apple, plum, pear, red currant, and gooseberry as orchard pests; and on rose and columbine as ornamental-plant pests [2,3]. The family is almost globally distributed and is more diverse in the Northern Hemisphere—in fact, its members are absent or rare in Antarctica, New Zealand, and Australia [2]. The observed distribution pattern is affected by adults who are slow fliers; this trait limits the distribution of the larvae of this diverse family on host plants. In addition to contributing to the overall distribution pattern, the host plant preferences of tenthredinid larvae have also led to an increase in speciation events by host switching, particularly in sublineages that feed on only some angiosperms [4].
This family is the largest and a considerably complex family of Symphyta, comprising 76.8% of the described species of the most diverse superfamily of ancestrally herbivorous hymenopterans, Tenthredinoidea, and approximately 65% of known Symphytan species, with 400 genera and more than 5600 extant species [5,6,7]. Currently, eight subfamilies are recognized in this family [6]: Allantinae, Athaliinae, Beldoneinae, Blennocampinae, Heterarthrinae, Nematinae, Selandriinae, and Tenthredininae. However, the subfamilial classification of Tenthredinidae has indeed been a longstanding problem, and it is known that at least some of the traditionally or recently recognized subfamilies or tribes are most likely artificial or incorrectly placed units in both morphological and molecular analyses [6,7,8,9,10,11,12,13,14,15].
Cladiucha Konow, 1902, and Megabeleses Takeuchi, 1952, are two small tenthredinid genera that feed on plants of the family Magnoliaceae (angiosperms), and most species of these genera have great economic importance, mainly in Southeast Asia [16,17,18,19,20,21,22,23]. Cladiucha is more similar in some respects to Diprionidae and some Pergidae, harboring multiple antennomeres and serrate-type antennae in the female and biramous-type antennae in the male rather than the filiform antennae usually observed in other tenthredinids [24]. This genus was first placed in the subfamily Allantinae, considering the similarity in wing venation and other structural characteristics, but generated a new tribe, Cladiuchini, for this genus [9]. Takeuchi [18] and Abe and Smith [8] also placed these genera into the subfamily Allantinae. Wei [25] also recognized the tribe Cladiuchini, but placed Cladiucha and Megabeleses in a new subfamily, Megabelesesinae Wei, 1997, along with two additional new genera, Tripidobeleses Wei, 1997, and Conobeleses Wei, 1997. This representation was followed in the new system provided by Wei and Nie [26]. However, Taeger et al. [7] placed these four genera into the subfamily Allantinae. This systematic inconsistency indicates the requirement for more detailed studies to better understand the systematic position and evolutionary history of the subfamily Megabelesesinae using comprehensive molecular data from the current approaches.
Here, four complete mitogenomes were sequenced and characterized from the family Tenthredinidae: Cladiucha punctata Wei, 2020, C. magnoliae G.R. Xiao, 1994 [not 1993], Megabeleses magnoliae Wei, 2010, M. liriodendrovorax G.R. Xiao, 1993. These mitogenomes were compared with the previously reported tenthredinid mitogenomes for a better understanding of the characteristics of the mitogenome of Tenthredinidae. We also constructed a mitogenome dataset representing eight subfamilies of the family Tenthredinidae and five families of the superfamily Tenthredinoidea by adding 29 previously and 16 newly sequenced species (Table S1). The dataset was generated for phylogenetic reconstruction of the family, applying several commonly used phylogenetic inference methods to overcome the systematic bias resulting from the non-stationarity of sequence evolution. A dated phylogeny was also rebuilt to assess the relationship between branching events and host shifts among the subfamilies.

2. Materials and Methods

2.1. Samples Analyzed

Specimens of four species were provided by Jiangxi Normal University, China, and specimen data were presented in Table S2. Total genomic DNA was isolated from the hind leg of each ethanol-preserved specimen by the DNeasy Tissue Kit (Qiagen, Hilden, Germany) according to instructions provided by the manufacturer, and the isolated genomic DNA samples were quantified using a NanoDrop (Maestrogen, Inc., Waltham, MA, USA).

2.2. Mitogenome Sequencing, Annotation, and Analyses

The total genomic DNA extracts of four species were pooled and sequenced by the Illumina HiSeq 4000 next-generation sequencing (NGS) platform using 150-bp paired-end reads, conducted at Shanghai Majorbio Bio-pharm Technology Co., Ltd., China. Raw NGS reads of each species were firstly considered for quality control by FastQC v0.11.9 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc) (accessed on 10 January 2021). The FASTQ files were imported into Geneious R11 (Biomatters Ltd, Auckland, New Zealand) [27], and the reads were then trimmed with BBDuk as implemented in Geneious R11. Duplicate reads and reads < 100 bp were filtered out to eliminate low-quality scores. The high-quality reads were then used to construct the mitogenome sequences of the four species. De novo assemblies of sequences were performed using MIRA assembler implemented in Geneious R11, and these assemblies were then mapped using the reference mitogenomes with ‘medium-low sensitivity’ parameters for each species. The assemblies constructed by these approaches were aligned, manually compared, and finally compiled into a single contig for each species. tRNA genes were found by their presumed secondary structure and anticodon sequence using the tRNAscan-SE server [28] and DOGMA [29] with the mito/chloroplast genetic code under the default search options. The boundaries and locations of the protein-coding genes (PCGs) and rRNA genes were identified by comparing the reported tenthredinid homologous gene sequences by ORF Finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html) (accessed on 10 January 2021) and by homology-based BLAST searches in GenBank. The precise ends of rRNA genes were predicted from the boundaries of the neighboring tRNA genes. The rRNA secondary structures were inferred by comparison with those of other reported sawfly species [30,31]. These predicted structures were visualized by VARNA v3-93 [32] and RNAviz 2.0.3 [33] with reference to the results of the CRW site [34]. Intergenic spacers and overlapping regions were inferred manually. The sequences of these four mitogenomes were deposited in GenBank under accession numbers MT295305-MT295306 and MW255939-MW255940. Summary statistics on the base composition, nucleotide substitution, and codon usage were analyzed by MEGA v6.0 [35]. Strand asymmetries were calculated using the following formulae: AT skew = (A − T)/(A + T) and GC skew = (G − C)/(G + C) [36].

2.3. Phylogenetic Analyses

Alignment and model selection. Phylogenetic analyses were carried out using a mitogenome dataset for the superfamily Tenthredinoidea: Thirty species (16 previously and 14 newly sequenced) from Tenthredinidae, representing eight known subfamilies, and 15 species (13 previously and two newly sequenced) from five other families, representing Argidae (three species), Cimbicidae (nine species), Diprionidae (one species), Heptamelidae (one species), and Pergidae (one species). One species from the superfamily Xyeloidea was included as an outgroup (Table S1). Mitogenomes of 29 previously reported tenthredinoid species and one outgroup were retrieved from the NCBI annotated database using Geneious R11. Each PCG was aligned under codon-based multiple alignments with ClustalW [37] as implemented in MEGA v6.0. The RNA genes were aligned as DNA using the MAFFT algorithm [38] as implemented in Geneious R9. Each obtained alignment was then concatenated using SequenceMatrix v.1.7.8 [39]. The optimal partition strategy and best fit evolutionary model of each partition were designated by PartitionFinder v1.1.1 [40], applying the Bayesian information criterion (BIC) and the ‘greedy’ algorithm based on branch lengths estimated as ‘unlinked’. The data blocks were defined by genes and codons to generate an input configuration file with 63 (with all codon positions) and 50 (without 3rd codon positions) partitions. The best partitioning schemes and related models were used in all subsequent phylogenetic analyses (Table S2).
Test of substitution saturation. The genetic saturation levels of different codon positions and genes were calculated by correlation coefficient analysis implemented in R core packages [41], comparing the distances estimated by applying the best fit evolutionary model GTR + G + I proposed by jModelTest2 with the uncorrected p-distances [42]. The distance values were measured with PAUP v4.0b10 [43].
Phylogenetic inference. The phylogenetic reconstructions were performed with four different datasets: (i) The 13 PCGs with all codon positions plus 22 tRNAs and two rRNAs (P123RNA) (16,719 bp); (ii) the 13 PCGs with the 1st and 2nd codon positions plus 22 tRNAs and two rRNAs (P12RNA) (12,811 bp); (iii) the PCGs with all codon positions excluding three genes plus 22 tRNAs and two rRNAs (P123RNAexc3genes) (15,871 bp), and (iv) the PCGs with the 1st and 2nd codon positions excluding three genes plus 22 tRNAs and two rRNAs (P12RNAexc3genes) (12,150 bp). The fourth dataset was generated by considering the result of the substitution saturation test, which indicated lower degrees of correlation between all codon positions of PCGs and the first and second codon positions of atp8, nad4l, and nad6 (Table S3). Maximum likelihood (ML) and Bayesian inference (BI) approaches were used for each dataset to infer whether the datasets were sensitive to phylogenetic inference methods. ML analyses were conducted in RAXML v8.0.9 [44] as implemented in Geneious R11, using the suggested GTR + G + I substitution model by PartitionFinder v1.1.1 for each nucleotide partition under 1000 bootstrap replicates. BI analyses were performed for unlinked branch lengths of each partition scheme under MrBayes v3.2.2 [45]. Distribution of posterior parameters were estimated in two independent runs with four Markov chains (three cold, one heated) based on 10 million generations and sampling every 1000 generations. The log-likelihood files produced by each run were assessed considering ESS values equal to or greater than 200 for all priors using Tracer v1.7 [46]. After the assessment, the first 25% of trees sampled in each run were scrapped as burn-in, and a majority-rule consensus tree (BI tree) was constructed from the remaining trees.
Divergence time estimation. The divergence times among the tenthredinids were estimated in BEAST v1.8.3 utilizing the dataset of P123RNA [47]. The uncorrelated relaxed lognormal clock was preferred because the dataset was found to be “not clock-like” (likelihood ratio test: −ln + c 210,126.327, −ln−c 209,813.739, d.f. = 15, p = 5.14E–92). Yule speciation tree priors were applied to model rate variation among species [48]. To estimate divergence times, a second calibration point was used: The age of (i) the split between the clade of Argidae-Pergidae and all remaining Tenthredinoidea (a mean of 175 Ma and ranging from 150–200 Ma) and (ii) the Argidae and Pergidae split (a mean of 135 Ma and ranging from 110–160 Ma) based on previous studies [4,49,50]. The run was performed using 100 million generations with samples taken every 10,000 generations, and the result was evaluated to assess convergence and confirm effective sample sizes (ESS > 200) in Tracer v1.6. The maximum clade credibility of the trees was quantified using TreeAnnotator v1.8.3 [47] after discarding 25% of the samples as burn-in. The obtained tree was then visualized in FigTree v1.4.2 [51].

3. Results

3.1. Mitogenome Architectures and Nucleotide Compositions

The complete mitogenomes of four species representing two genera were sequenced and characterized: C. magnoliae (15,761 bp), C. punctata (16,187 bp), M. magnoliae (16,219 bp), and M. liriodendrovorax (15,466 bp) (Figure 1, Table S4). Each mitogenome consisted of a typical set of 37 genes: Thirteen PCGs, twenty-two tRNAs, two rRNAs, and an A + T rich region. The observed length variation among these mitogenomes was primarily due to variation in the A + T-rich region (Table S4). The genes were mostly located on the majority J strand, except for four PCGs (nad1, nad4, nad4l, and nad5), two rRNAs, and eight tRNAs. The mitochondrial gene orders of the four species were entirely conserved (Figure 1). However, their arrangements were slightly different from those of the inferred ancestral pancrustacean mitogenome (Daphnia pulex), with the occurrence of several tRNA gene rearrangements. The first rearrangement event was found in the IQM gene cluster; here, the trnQ gene was inverted, and the trnM and trnI genes were transposed, moving upstream and downstream of the trnQ gene, respectively (arranged as MQI, Figure 1). The second event was a reverse transposition of the trnR gene upstream of the trnF gene in the ARNS1EF gene cluster (arranged as ANS1ERF, Figure 1).
The nucleotide compositions of these mitogenomes were biased toward A and T, ranging from 80.7% (M. liriodendrovorax) to 83.1% (C. punctata and C. magnolia) (Table S5). Excluding the A + T-rich regions, the highest AT content was found in rRNA genes of C. punctata (85.5%), and the lowest was observed in the PCGs of M. liriodendrovorax (79.5%). The skew metrics in the whole mitogenomes were different with respect to the AT skew, but similar in terms of the GC skew. A more pronounced A bias was observed in Megabeleses, with an average AT skew of 0.033, than in Cladiucha (0.013). However, a deviation was found in the PCGs: T-skewed in all species (with an average AT skew of −0.113) and G-skewed in C. magnoliae (0.006), M. liriodendrovorax (0.005), and M. magnoliae (0.026). This deviation probably results from T bias in the 2nd codon position and G bias in the 1st codon position. All PCGs and rRNA genes were similar in length (Table S4). However, a variation was observed in cox2, with six codon differences between M. magnoliae (225 aa) and the other species (227 aa). The estimated initiation codons commonly encoded either isoleucine or methionine (ATN), and termination codons were conserved, except for cob in two Megabeleses species and nad4 in two Cladiucha species and M. magnoliae, which all use partial termination codons (T-). A significant correlation was observed between nucleotide composition and codon preference (Figure 2, Table S6). UUA-Leu had the highest relative synonymous codon usage (RSCU), with an average value of 3.53, and all the remaining codons with RSCU greater than 2.00 had T or, particularly, A in the third codon position (Figure 2).
All the predicted tRNAs had standard anticodons, ranging in size from 60 bp (trnV in M. liriodendrovorax) to 71 bp (trnK in Megabeleses, trnT in C. magnoliae) and folding into canonical clover-leaf structures, except for all trnS1 in all four species and trnV in two Megabeleses species (Figure 3). The identical location of the rRNA genes was found in all species, with rrnS between the trnV and A + T-rich regions and rrnL positioned between trnL1 and trnV (Figure 1, Table S4). The putative secondary structures of rrnL and rrnS genes in four species were consistent with the models suggested for other insects, consisting of 49 helices from five structural domains (domain III is absent as in other arthropods) (Figure 4 and Figure 5). Variation was observed in the total length of intergenic regions among the species. The total length of the intergenic regions was 123 bp (at 19 different locations) in M. liriodendrovorax, 253 bp (at 17 locations) in C. magnoliae, 324 bp (at 17 locations) in C. punctata, and 461 bp (at 19 locations) in M. magnoliae (Table S4). Homology searches for these intergenic regions showed no significant similarity with any identified nucleotide sequence. The length and location of overlaps between adjacent genes were similar in both Cladiucha (between 1 and 7 bp) and Megabeleses (between 1 and 8 bp) (Table S4).

3.2. Phylogeny of Tenthredinidae

Eight phylogenetic reconstructions (four datasets by two inference approaches) recovered four different tree topologies (Figure S1 and Table 1). The obtained topologies were sensitive to both inference approach and dataset differences, and the support values were mostly lower for ML trees than for the BI of the corresponding dataset (Figure 6 and Figure S1). The exclusion of third codon positions did not affect the topology of the recovered phylogenetic tree, but the tree was sensitive to the saturated PCGs (atp8, nad4l, and nad6; Table S2). The most frequently obtained tree was recovered from the datasets of P123RNA and P12RNA, with both ML and BI approaches finding the same topology. This tree topology was also obtained from the dataset of P123RNAexc3genes applying just the ML inference method. The recovered tree supported the monophyly of Pergidae, Argidae, Diprionidae, and Cimbicidae, while a paraphyletic relationship was observed in Tenthredinidae, with the placement of Athaliinae (Figure 6). Argidae and Pergidae formed a monophylum that was sister to the remaining taxa. Next, Athalia (currently classified in the tenthredinid subfamily Allantinae) was separated from a large clade consisting of the representatives of Diprionidae, Cimbicidae, Heptamelidae, and Tenthredinidae with strong nodal support [posterior probability (PP) = 1.00 and ML = 100%] (Figure 6). A sister-group relationship was supported between Diprionidae and Cimbicidae (PP = 1.00, ML = 73%), and this clade was also identified as sister to Heptamelidae + Tenthredinidae (PP = 1.00, ML = 78%). At the family level, the topology recovered P12RNAexc3genes under both inferences and P123RNAexc3genes under BI supported a relationship of (Heptamelidae + (Diprionidae + Cimbicidae)) + Tenthredinidae (PP = 1.00, ML = 61%) (Figure S1, Table 1). Within Tenthredinidae s. str. (except for a rogue taxon Athalia), all analyses supported the monophyly of the subfamilies Nematinae, Selandriinae, and Tenthredininae with high support, while Allantinae, Heterarthrinae, and Blennocampinae were found to be paraphyletic (Figure 6 and Figure S1). Nematinae was robustly placed as the most basal tenthredinid subfamily with strong nodal support (PP = 1.00, ML = 100%), except for the topology recovered under P12RNAexc3genes under BI, which was obtained a sister-group relationship between Nematinae and Selandriinae. After the divergence of Nematinae, Selandriinae was recovered as a sister to all remaining tenthredinids (Figure 6 and Figure S1). Next, the four newly sequenced species representing Cladiucha and Megabeleses (currently in Allantinae) formed a distinct lineage with high support values (PP = 1.00, ML = 98%). They were sisters to the representatives of Blennocampinae, Heterarthrinae, Tenthredininae, and the rest of Allantinae (Figure 6 and Figure S1). For the remaining four subfamilies, a sister-group relationship between Heterarthrinae, Blennocampinae, and Hemibeleses tianmunicus (Allantinae) forming a paraphyletic grade, and Tenthredininae + Allantinae s. str with strong nodal supports (PP = 1.00, ML = 78%), except for the topology recovered under P12RNAexc3genes under ML, which Tenthredininae had a basal placement (ML = 88%) (Figure 6 and Figure S1, Table 1).

3.3. Divergence Time Estimations at the Subfamily Level

The inferred divergence and diversification time estimates of tenthredinids were based on the dataset of P123RNA and a secondary-calibrated relaxed molecular clock (Figure 7). According to the dated tree, the stem-group age of Tenthredinidae was estimated to be 114.46 Ma (94.10–134.71) in the Early Cretaceous, which also corresponds to the split of Heptamelidae and Tenthredinidae. Diversification of Tenthredinidae occurred between the end of the Early Cretaceous and the late Miocene (Figure 7). The oldest division in the tenthredinids is the subfamily Nematinae, with an estimated age of 104.39 (85.63–122.66) Ma. The stem age of the Cladiucha and Megabeleses lineages was dated at 95.13 (83.64–118.87) Ma in the Early–Late Cretaceous boundary, while the split between these genera was estimated to have occurred at the beginning of the Eocene at 53.34 (39.92–67.67) Ma (Figure 7). Within this lineage, the divergence time of M. magnoliae and M. liriodendrovorax corresponded to the Oligocene [31.76 Ma, (21.57–42.62 Ma)], while the split of C. magnoliae and C. punctata was the most recent and was dated to 10.98 (6.97–15.49) Ma, corresponding to the late Miocene.

4. Discussion

4.1. Mitogenome Organization of Cladiucha and Megabeleses

Four complete and twelve nearly complete mitogenomes of Tenthredinidae have been reported so far [14,53,54,55,56,57,58,59,60,61]. Here, the complete mitogenomes of four species representing Cladiucha and Megabeleses were sequenced for the first time, and their mitogenome characteristics were summarized at the subfamily level by comparing them with the published and/or unpublished mitogenomes representing the family. All sequenced genomes are highly conserved in size, nucleotide content, codon usage, and secondary structures of tRNAs and rRNAs (Figure 2, Figure 3, Figure 4 and Figure 5, Tables S4 and S5). Their mitogenome sizes fit well within the range reported for tenthredinid mitogenomes, as well as for Symphytan mitogenomes, ranging from 15,108 bp in Hemathlophorus brevigenatus (MW632125 in GenBank) to 20,370 bp in Trachelus iudaicus [62]. Similar to the reported rearrangement hotspot points in hymenopteran mitogenomes [63,64], the most rearrangements were frequently observed in the IQM, ARNS1EF, and WCY gene clusters (Figure 1), possibly related to illicit priming of mitochondrial replication and/or illegitimate intramitochondrial recombination [65,66]. tRNA rearrangements have been most frequently reported; the mitogenomes of the derived hymenopteran suborder Apocrita appear to be more prone to rearrangements than Symphytan mitogenomes, but the observed pattern in tenthredinid mitogenomes, as well as other reported Symphytan mitogenomes, supports the frequent rearrangement of tRNA genes in Symphyta [31,60,62,67,68,69]. Here, the same type of architecture of tRNA rearrangement in the IQM (arranged as MQI) and ARNS1EF (arranged as ANS1ERF) gene clusters in the sequenced mitogenomes of both Cladiucha and Megabeleses species may provide further evidence of molecular synapomorphy for this lineage (Figure 1). These derived gene rearrangements were also not found in other tenthredinid taxa, indicating that they occurred during the divergence of this lineage, and therefore, supported their recognition as a subfamily [25].

4.2. Phylogenetic Placement of Cladiucha and Megabeleses

Among the known Symphytan families, the internal phylogeny of Tenthredinidae has attracted much attention because of the unstable phylogenetic placement of the subfamilies, as well as their monophyly within the family depending on the combination or absence of enough morphological characters, relatively narrow range of taxonomic coverage and a limited number of easily accessible sequence datasets [10,12]. To elucidate relationships within the family, previous phylogenetic studies based on morphological characters and/or recent comprehensive molecular analyses have provided several conclusions. (i) The genus Athalia was moved to the subfamily Athaliinae rather than Allantinae [10], but its taxonomic status remains unstable. This conclusion was mainly based on the placement of Athalia, which was placed as either sister to all other members of the family [10,70,71] or outside of the tenthredinids, with changing placement among the families Cimbicidae, Diprionidae, Heptamelidae, and Tenthredinidae [49,72,73,74,75]. (ii) The genus Heptamelus or the tribe Heptamelini should be excluded from Selandriinae [12]. Moreover, the genera Heptamelus, Carinoscutum, Parahemitaxonus, and Pseudoheptamelus have recently been recognized as members of a small new family, Heptamelidae [6,10]. (iii) The monophyly of Nematinae and Tenthredininae was supported in most studies, but the obtained topologies can also introduce uncertainty in the placement of the remaining traditional subfamilies [10,12,49].
Here, the internal phylogeny of the family was produced with a mitogenome dataset, including 13 PCGs (three saturated genes were excluded in some analyses, Figure S1), two rRNA genes, and 22 tRNA genes (15,365 bp in length), but with a relatively low number of species (30 species representing eight subfamilies). Despite the fact that the present phylogeny suggests a complex evolutionary history, the recovered branching pattern with high nodal support was mostly compatible with previously reported phylogenies (Figure 6) [10,12,49,72,73,75]. The phylogenetic placement of the genus Athalia, outside of Diprionidae, Cimbicidae, Heptamelidae, and Tenthredinidae, obviously did not support its current classification under the subfamily Athaliinae and required reassessment to provide a solid foundation of its taxonomic position. The sister-group relationship between Heptamelidae and Tenthredinidae was supported in most of the analyses (Figure 6), as described by Boevé et al. [70]. At the same time, the placement of Heptamelidae as a sister taxon to a clade comprising Diprionidae, Cimbicidae, and Tenthredinidae was also recovered by Malm and Nyman [10] (Table 1, Figure S1). Within Tenthredinidae, the present phylogeny is broadly congruent with previously reported relationships, supporting the monophyly of the subfamilies Nematinae, Selandriininae, and Tenthredininae (Figure 6) [10,11,49]. However, the formation of the paraphyletic grade comprising the representatives of Heterearthrinae, Blennocampinae, and Hemibeleses (Allatinae) indicates that on the one hand, these subfamilies need to revision in their current classifications, but on the other hand, this suggestion would be premature to recommend a new subfamilial classification for Tenthredinidae until a more comprehensive phylogeny becomes available, including more genera and/or species representatives of these subfamilies. But even so, Cladiucha and Megabeleses were consistently identified as monophyletic and outside of Allantinae (Figure 6 and Table 1). The synapomorphy of Megabelesesinae is as follows: The female lance has a long and broad membranous lobe and approaches the upper 0.3–0.5 of the lancet (Figure 7). Another possible synapomorphic character is the vein 1m-cu being very close to the base of vein Rs. Some shared morphological characters among the species of these genera, but not in the species of Allantinae, also support Megabelesesinae as a plesiomorphic: (i) The mandibles are simple and bidentate, (ii) the head is not distinctly enlarged behind the eyes, and (iii) the anterior lobe of the pronotum is very narrow. Their phylogenetic placement as a distinct lineage, in conjunction with (i) the morphological synapomorphy discussed above; (ii) the shared molecular features (occurrence of the same type of tRNA rearrangement, arranged as MQI and ANS1ERF in their sequenced mitogenomes; Figure 1, Table S4); and (iii) host specialists feeding on Magnoliaceae, may provide evidence supporting the proposal of the subfamily Megabelesesinae.

4.3. Divergences of Tenthredinid Subfamilies through Time

The stem group age and crown age of Tenthredinidae are estimated to be older than those reported by Peters et al. [50], but broadly congruent with previously reported dated trees [4,49,73]. A close association between the host preferences and divergence times of the subfamilies was observed in the dated tree (Figure 7), supporting the suggestion that the diversification rate increments in tenthredinids can be directly related to dominant colonization of flowering plants (angiosperms) following their initial appearance by the Early Cretaceous between approximately 130 and 100 Ma [4,76,77]. The estimated crown ages of the subfamilies also suggest that the major divergence within the family corresponds to the Cenomanian/Turonian (estimated at approximately 95 Ma) and Cretaceous–Paleogene boundaries (estimated at approximately 70 Ma). These important geological intervals, recognized as severe biotic crises in the history of life [78], played an unignorable role in influencing the floral compositions, with minor changes in the diversification of nonflowering seed plants (gymnosperms), significant extinction events in spore-bearing plants (pteridophytes) and rapid increases in the origin and diversification rates of angiosperms [77]. The sudden and rapid changes in the global biodiversity patterns during these periods most likely triggered diversification and colonization of new host plants for tenthredinid sawflies [4].
The present dated tree also suggests that the divergence and diversification times of Cladiucha and Megabeleses nearly coincide with the divergence times of their host plants (Figure 7). The members of these genera constituting a distinct lineage, Megabelesesinae, appear to be more special with a relatively narrow host range, feeding on the members of the family Magnoliaceae [20,24,79]. Based on fossil records and evidence from molecular phylogenetic studies, this family has generally been considered one of the relatively primitive families of flowering plants with woody members [80,81]. The origin of Magnoliaceae and split time between two subfamilies—namely, Magnolioideae and Liriodendroideae, are estimated to correspond to the Early–Late Cretaceous boundary, occurring approximately 93.5 Ma [80], 100 Ma [82], or 113 Ma [83]. However, the crown ages of Magnolioideae and Liriodendroideae are dated to approximately 42 Ma or 54 Ma and 28 Ma or 15 Ma, respectively [80,82]. The subfamily and species divergence of the family Magnoliaceae mostly matches the divergence of Cladiucha and Megabeleses, as well as their species diversification (Figure 7). However, concomitant radiation does not necessarily imply a reciprocally causative evolutionary scenario—at best, our data do not exclude the hypothesis that speciation of herbivorous insects was triggered by host shifts depending on the radiation of their host plants.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects12060495/s1, Figure S1: Phylogenetic tree recovered from the dataset of P12RNAexc3genes using (a) ML, (b) BI; of P123RNAexc3genes using (c) ML and (d) BI; of P12RNA using (e) ML and (f) BI. An optimal partitioning scheme was selected with PartitionFinder. The posterior probabilities and/or bootstrap values were shown in the trees, Table S1: Summary information of mitogenomes used in phylogenetic analyses, Table S2: Best partitioning scheme and model selected by PartitionFinder for phylogenetic analyses, Table S3: Regression of the pairwise distances of the different genes and codon positions in the Megabelesesinae species, Table S4: Mitogenome organizations of Megabelesesinae species, Table S5: Nucleotide composition of the mitogenomes of Megabelesesinae species, Table S6: Relative synonymous codon usage (RSCU) of the Megabelesesinae mitogenomes.

Author Contributions

Conceptualization, M.W.; methodology, M.W.; validation, D.W.; formal analysis, S.J., Ö.D., and M.B.; investigation, G.N.; resources, M.W.; data curation, S.J.; writing-original draft, S.J., and Ö.D.; writing-review and editing, G.N. and E.M.K.; visualization, G.N.; supervision, M.W.; project administration, M.W.; funding acquisition, G.N. and M.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (31501885, 31970447).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in Science Data Bank at DOI: 10.11922/sciencedb.00702.

Acknowledgments

The members of Lab of Insect Systematics and Evolutionary Biology (LISEB) from Jiangxi Normal University and Evolutionary Bioinformatics Research Groups (EBRG) from Sivas Cumhuriyet University are thanked for their contributions to laboratory work. We thank the anonymous reviewers for their careful reading and many constructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jervis, M.; Vilhelmsen, L. Mouthpart evolution in adults of the basal, “symphytan”, hymenopteran lineages. Biol. J. Linn. Soc. 2000, 70, 121–146. [Google Scholar] [CrossRef]
  2. Boevé, J.L. Sawflies (Hymenoptera: Tenthredinidae). In Encyclopedia of Entomology; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2004; pp. 1952–1954. [Google Scholar]
  3. Robinson, W.H. Hymenoptera. In Handbook of Urban Insects and Arachnids; Cambridge University Press: Cambridge, UK, 2005; pp. 224–290. [Google Scholar]
  4. Nyman, T.; Onstein, R.E.; Silvestro, D.; Wutke, S.; Taeger, A.; Wahlberg, N.; Blank, S.M.; Malm, T. The early wasp plucks the flower: Disparate extant diversity of sawfly superfamilies (Hymenoptera: ‘Symphyta’) may reflect asynchronous switching to angiosperm hosts. Biol. J. Linn. Soc. 2019, 128, 1–19. [Google Scholar] [CrossRef]
  5. Goulet, H.; Hubert, J.T. Hymenoptera of the World: An Identification Guide to Families; Agriculture Canada Publication: Ottawa, ON, Canada, 1993; ISBN 0660149338. [Google Scholar]
  6. Taeger, A.; Liston, A.D.; Prous, M.; Groll, E.K.; Gehroldt, T.; Blank, S.M. EcatSym—Electronic World Catalog of Symphyta (Insecta, Hymenoptera); Program Version 5.0 (19 December 2018), Data Version 40 (23 September 2018); Senckenberg Deutsches Entomologisches Institut (SDEI): Müncheberg, Germany; Available online: https://sdei.de/ecatsym/ (accessed on 15 February 2020).
  7. Taeger, A.; Blank, S.M.; Liston, A.D. World catalog of Symphyta (Hymenoptera). Zootaxa 2010, 2580, 1–1064. [Google Scholar] [CrossRef] [Green Version]
  8. Abe, M.; Smith, D.R. The genus-group names of Symphyta (Hymenoptera) and their type species. Esakia 1991, 31, 1–115. [Google Scholar] [CrossRef]
  9. Benson, R.B. On the classification of sawflies (Hymenoptera: Symphyta). Trans. R. Entomol. Soc. Lond. 1938, 87, 353–384. [Google Scholar] [CrossRef]
  10. Malm, T.; Nyman, T. Phylogeny of the symphytan grade of Hymenoptera: New pieces into the old jigsaw(fly) puzzle. Cladistics 2015, 31, 1–17. [Google Scholar] [CrossRef]
  11. Nyman, T.; Zinojev, A.G.; Vikberg, V.; Farrell, B.D. Molecular phylogeny of the sawfly subfamily Nematinae (Hymenoptera: Tenthredinidae). Syst. Entomol. 2006, 31, 569–583. [Google Scholar] [CrossRef]
  12. Vilhelmsen, L. Morphological phylogenetics of the Tenthredinidae (Insecta: Hymenoptera). Invertebr. Syst. 2015, 29, 164. [Google Scholar] [CrossRef]
  13. Wei, M.; Niu, G. Review of Anhoplocampa Wei (Hymenoptera, Tenthredinidae), with description of a new species and a new combination. Zookeys 2011, 159, 81–90. [Google Scholar] [CrossRef] [Green Version]
  14. Niu, G.; Zhang, Y.; Li, Z.; Wei, M. Characterization of the mitochondrial genome of Analcellicampa xanthosoma gen. et sp. nov. (Hymenoptera: Tenthredinidae). PeerJ 2019, 7, e6866. [Google Scholar] [CrossRef] [Green Version]
  15. Lacourt, J. Contribution à une révision mondiale de la sous-famille des Tenthredininae (Hymenoptera: Tenthredinidae). Annales de la Société Entomologique de France (N. S.) 1996, 32, 363–402. [Google Scholar]
  16. Boquan, W.; Xunru, A.; Qiong, P.; Zongrong, L.; Wei, W. Preliminary investigation of the pests of Magnolia officinalis. Plant Prot. 2011, 37, 132–134. [Google Scholar]
  17. Ding, Y.Z.; Zheng, H.S. Study on the sawfly Megabeleses crassitarsis (Hymenoptera: Tenthredinidae, Allantinae). Sci. Silvae Sin. 1999, 35, 68–71. [Google Scholar]
  18. Takeuchi, K. A Generic Classification of the Japanese Tenthredinidae (Hymenoptera: Symphyta); Issued in Celebration of the Sixtieth Birthday of Kichizo Takeuchi by his Friends: Kyoto, Japan, 1952. [Google Scholar]
  19. Wei, M. A new species of Cladiucha Konow (Hymenoptera, Tenthredinidae) from China with a key to sawfly species damaging plants of Magnoliaceae. Acta Zootaxonomica Sin. 2010, 35, 635–640. [Google Scholar]
  20. Wei, M. Revision of Megabeleses Takeuchi (Hymenoptera, Tenthredinidae) with description of two new species from China. Zootaxa 2010, 2729, 36–50. [Google Scholar] [CrossRef] [Green Version]
  21. Niu, G.; Zhao, H.; Wei, M. Review of Cladiucha Konow (Hymenoptera: Tenthredinidae) with description of a new species from China. Entomotaxonomia 2020, 42, 311–318. [Google Scholar] [CrossRef]
  22. Sundukov, Y.N. A new species of the genus Cladiucha Konow, 1902 (Hymenoptera: Tenthredinidae) from Vietnam. Far East. Entomol. 2020, 407, 25–28. [Google Scholar] [CrossRef]
  23. Konow, F.W. Neue Blattwespen. (Hym.). J. Syst. Hymenopterol. Dipterol. 1902, 2, 384–390. [Google Scholar]
  24. Smith, D.R. Rediscovery of Cladiucha insolita Konow (Hymenoptera, Tenthredinidae), description of the male and intraspecific variation. J. Hymenopt. Res. 2017, 60, 173–179. [Google Scholar] [CrossRef] [Green Version]
  25. Wei, M. A new subfamily and two new genera of Tenthredinidae (Hymenoptera: Tenthredinomorpha). Entomotaxonomia 1997, 19, 69–76. [Google Scholar]
  26. Wei, M.; Nie, H. Generic list of Tenthredinoidea s. str. (Hymenoptera) in new systematic arrangement with synonyms and distribution data. J. Cent. For. Coll. 1998, 18, 23–31. [Google Scholar]
  27. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  28. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef]
  29. Wyman, S.K.; Jansen, R.K.; Boore, J.L. Automatic annotation of organellar genomes with DOGMA. Bioinformatics 2004, 20, 3252–3255. [Google Scholar] [CrossRef] [Green Version]
  30. Korkmaz, E.M.; Doğan, Ö.; Budak, M.; Başıbüyük, H.H. Two nearly complete mitogenomes of wheat stem borers, Cephus pygmeus (L.) and Cephus sareptanus Dovnar-Zapolskij (Hymenoptera: Cephidae): An unusual elongation of rrnS gene. Gene 2015, 558, 254–264. [Google Scholar] [CrossRef]
  31. Niu, G.; Korkmaz, E.M.; Doğan, Ö.; Zhang, Y.; Aydemir, M.N.; Budak, M.; Du, S.; Başıbüyük, H.H.; Wei, M. The first mitogenomes of the superfamily Pamphilioidea (Hymenoptera: Symphyta): Mitogenome architecture and phylogenetic inference. Int. J. Biol. Macromol. 2019, 124, 185–199. [Google Scholar] [CrossRef]
  32. Darty, K.; Denise, A.; Ponty, Y. VARNA: Interactive drawing and editing of the RNA secondary structure. Bioinformatics 2009, 25, 1974–1975. [Google Scholar] [CrossRef] [Green Version]
  33. De Rijk, P.; Wuyts, J.; de Wachter, R. RnaViz 2: An improved representation of RNA secondary structure. Bioinformatics 2003, 19, 299–300. [Google Scholar] [CrossRef] [Green Version]
  34. Cannone, J.J.; Subramanian, S.; Schnare, M.N.; Collett, J.R.; D’Souza, L.M.; Du, Y.; Feng, B.; Lin, N.; Madabusi, L.V.; Müller, K.M.; et al. The comparative RNA web (CRW) site: An online database of comparative sequence and structure information for ribosomal, intron, and other RNAs. BMC Bioinform. 2002, 3, 2. [Google Scholar] [CrossRef] [Green Version]
  35. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  36. Perna, N.T.; Kocher, T.D. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J. Mol. Evol. 1995, 41, 353–358. [Google Scholar] [CrossRef] [PubMed]
  37. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The CLUSTAL X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25, 4876–4882. [Google Scholar] [CrossRef] [Green Version]
  38. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  39. Vaidya, G.; Lohman, D.J.; Meier, R. SequenceMatrix: Concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics 2011, 27, 171–180. [Google Scholar] [CrossRef]
  40. Lanfear, R.; Calcott, B.; Ho, S.Y.W.; Guindon, S. PartitionFinder: Combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol. Biol. Evol. 2012, 29, 1695–1701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; ISBN 3900051070. [Google Scholar]
  42. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef] [Green Version]
  43. Swofford, D.L.; PAUP. Phylogenetic analysis using parsimony (and other methods). Version 4. Sinauer Assoc. 2002, 294–307. [Google Scholar] [CrossRef] [Green Version]
  44. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [PubMed]
  45. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. Mrbayes 3.2: Efficient bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [Green Version]
  46. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior summarisation in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 2018, syy032. [Google Scholar] [CrossRef] [Green Version]
  47. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [Green Version]
  48. Drummond, A.J.; Ho, S.Y.W.; Phillips, M.J.; Rambaut, A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4, 699–710. [Google Scholar] [CrossRef]
  49. O’Reilly, J.E.; dos Reis, M.; Donoghue, P.C.J. Dating tips for divergence-time estimation. Trends Genet. 2015, 31, 637–650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Peters, R.S.; Krogmann, L.; Mayer, C.; Donath, A.; Gunkel, S.; Meusemann, K.; Kozlov, A.; Podsiadlowski, L.; Petersen, M.; Lanfear, R.; et al. Evolutionary history of the Hymenoptera. Curr. Biol. 2017, 27, 1–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Rambaut, A. FigTree v1.4.2: Tree Drawing Tool; Institute of Evolutionary Biology, University of Edinburgh: Edinburgh, UK, 2014. [Google Scholar]
  52. Gillespie, J.J.; Johnston, J.S.; Cannone, J.J.; Gutell, R.R. Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): Structure, organization, and retrotransposable elements. Insect Mol. Biol. 2006, 15, 657–686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wei, S.-J.; Niu, F.-F.; Du, B.-Z. Rearrangement of trnQ-trnM in the mitochondrial genome of Allantus luctifer (Smith) (Hymenoptera: Tenthredinidae). Mitochondrial DNA 2016, 27, 856–858. [Google Scholar] [CrossRef]
  54. He, H.; Niu, G.; Zhang, B.; Wei, M. The complete mitochondrial genome of Athalia proxima (Hymenoptera: Tenthredinidae) and phylogenetic analysis. Mitochondrial DNA Part B 2019, 4, 3868–3869. [Google Scholar] [CrossRef] [Green Version]
  55. Li, Y.; Wei, M.; Liu, J.; Niu, G. Characterization of the mitochondrial genome of Eutomostethus vegetus Konow, 1898 (Hymenoptera: Tenthredinidae) and phylogenetic analysis. Mitochondrial DNA Part B 2020, 5, 3051–3052. [Google Scholar] [CrossRef]
  56. Wu, R.; Wei, M.; Liu, M.; Niu, G. Advancement in sequencing the mitochondrial genome of Birmella discoidalisa Wei, 1994 (Hymenoptera: Tenthredinidae) and the phylogenetic classification of Fenusini. Mitochondrial DNA Part B 2019, 4, 4100–4101. [Google Scholar] [CrossRef]
  57. Wu, D.; Wang, H.; Wei, M.; Niu, G. The nearly complete mitochondrial genome of Colochela zhongi Wei, 2016 (Hymenoptera: Tenthredinidae) and phylogenetic analysis. Mitochondrial DNA Part B 2020, 5, 3359–3360. [Google Scholar] [CrossRef]
  58. Luo, X.; Wei, M.; Niu, G. Nearly complete mitochondrial genome of Siobla xizangensis Xiao, Huang & Zhou, 1988 (Hymenoptera: Tenthredinidae) and phylogenetic analysis. Mitochondrial DNA Part B 2019, 4, 4102–4103. [Google Scholar] [CrossRef] [PubMed]
  59. Wei, S.J.; Wu, Q.L.; Liu, W. Sequencing and characterization of the Monocellicampa pruni (Hymenoptera: Tenthredinidae) mitochondrial genome. Mitochondrial DNA 2013. [Google Scholar] [CrossRef]
  60. Song, S.-N.; Tang, P.; Wei, S.-J.; Chen, X.-X. Comparative and phylogenetic analysis of the mitochondrial genomes in basal hymenopterans. Sci. Rep. 2016, 6. [Google Scholar] [CrossRef] [Green Version]
  61. Song, S.-N.; Wang, Z.-H.; Li, Y.; Wei, S.-J.; Chen, X.-X. The mitochondrial genome of Tenthredo tienmushana (Takeuchi) and a related phylogenetic analysis of the sawflies (Insecta: Hymenoptera). Mitochondrial DNA Part A 2016, 27, 2860–2861. [Google Scholar] [CrossRef] [PubMed]
  62. Korkmaz, E.M.; Aydemir, H.B.; Temel, B.; Budak, M.; Başıbüyük, H.H. Mitogenome evolution in Cephini (Hymenoptera: Cephidae): Evidence for parallel adaptive evolution. Biochem. Syst. Ecol. 2017, 71, 137–146. [Google Scholar] [CrossRef]
  63. Dowton, M.; Cameron, S.L.; Dowavic, J.I.; Austin, A.D.; Whiting, M.F. Characterization of 67 mitochondrial tRNA gene rearrangements in the Hymenoptera suggests that mitochondrial tRNA gene position is selectively neutral. Mol. Biol. Evol. 2009, 26, 1607–1617. [Google Scholar] [CrossRef] [Green Version]
  64. Dowton, M.; Austin, A.D. Evolutionary dynamics of a mitochondrial rearrangement “hot spot” in the Hymenoptera. Mol. Biol. Evol. 1999, 16, 298–309. [Google Scholar] [CrossRef] [Green Version]
  65. Dowton, M.; Castro, L.R.; Campbell, S.L.; Bargon, S.D.; Austin, A.D. Frequent mitochondrial gene rearrangements at the hymenopteran nad3-nad5 junction. J. Mol. Evol. 2003, 56, 517–526. [Google Scholar] [CrossRef]
  66. Kong, X.; Dong, X.; Zhang, Y.; Shi, W.; Wang, Z.; Yu, Z. A novel rearrangement in the mitochondrial genome of tongue sole, Cynoglossus semilaevis: Control region translocation and a tRNA gene inversion. Genome 2009, 52, 975–984. [Google Scholar] [CrossRef]
  67. Doğan, Ö.; Korkmaz, E.M. Nearly complete mitogenome of hairy sawfly, Corynis lateralis (Brullé, 1832) (Hymenoptera: Cimbicidae): Rearrangements in the IQM and ARNS1EF gene clusters. Genetica 2017, 145, 341–350. [Google Scholar] [CrossRef] [PubMed]
  68. Korkmaz, E.M.; Doğan, Ö.; Durel, B.S.; Temel Altun, B.; Budak, M.; Başıbüyük, H.H. Mitogenome organization and evolutionary history of the subfamily Cephinae (Hymenoptera: Cephidae). Syst. Entomol. 2018, 43, 606–618. [Google Scholar] [CrossRef]
  69. Aydemir, M.N.; Korkmaz, E.M. Comparative mitogenomics of Hymenoptera reveals evolutionary differences in structure and composition. Int. J. Biol. Macromol. 2020. [Google Scholar] [CrossRef]
  70. Boevé, J.-L.; Blank, S.M.; Meijer, G.; Nyman, T. Invertebrate and avian predators as drivers of chemical defensive strategies in tenthredinid sawflies. BMC Evol. Biol. 2013, 13, 198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Schulmeister, S.; Wheeler, W.C.; Carpenter, J.M. Simultaneous analysis of the basal lineages of Hymenoptera (Insecta) using sensitivity analysis. Cladistics 2002, 18, 455–484. [Google Scholar] [CrossRef]
  72. Heraty, J.; Ronquist, F.; Carpenter, J.M.; Hawks, D.; Schulmeister, S.; Dowling, A.P.; Murray, D.; Munro, J.; Wheeler, W.C.; Schiff, N.; et al. Evolution of the hymenopteran megaradiation. Mol. Phylogenet. Evol. 2011, 60, 73–88. [Google Scholar] [CrossRef] [Green Version]
  73. Ronquist, F.; Klopfstein, S.; Vilhelmsen, L.; Schulmeister, S.; Murray, D.L.; Rasnitsyn, A.P. A total-evidence approach to dating with fossils, applied to the early radiation of the Hymenoptera. Syst. Biol. 2012, 61, 973–999. [Google Scholar] [CrossRef] [Green Version]
  74. Sharkey, M.J.; Carpenter, J.M.; Vilhelmsen, L.; Heraty, J.; Liljeblad, J.; Dowling, A.P.G.; Schulmeister, S.; Murray, D.; Deans, A.R.; Ronquist, F.; et al. Phylogenetic relationships among superfamilies of Hymenoptera. Cladistics 2012, 28, 80–112. [Google Scholar] [CrossRef] [Green Version]
  75. Schulmeister, S. Morphology and evolution of the tarsal plantulae in Hymenoptera (Insecta), focussing on the basal lineages. Zool. Scr. 2003, 32, 153–172. [Google Scholar] [CrossRef]
  76. Crane, P.R.; Friis, E.M.; Pedersen, K.R. The origin and early diversification of angiosperms. Nature 1995, 374, 27–33. [Google Scholar] [CrossRef]
  77. Silvestro, D.; Cascales-Miñana, B.; Bacon, C.D.; Antonelli, A. Revisiting the origin and diversification of vascular plants through a comprehensive Bayesian analysis of the fossil record. New Phytol. 2015, 207, 425–436. [Google Scholar] [CrossRef]
  78. Monnet, C. The Cenomanian–Turonian boundary mass extinction (Late Cretaceous): New insights from ammonoid biodiversity patterns of Europe, Tunisia and the Western Interior (North America). Palaeogeogr. Palaeoclimatol. Palaeoecol. 2009, 282, 88–104. [Google Scholar] [CrossRef]
  79. Xiao, G. Redescription of the genus Cladiucha (Hymenoptera: Tenthredinidae) and descriptions of two new species from China. J. Beijing For. Univ. 1994, 3, 15–22, English edition. [Google Scholar]
  80. Nie, Z.L.; Wen, J.; Azuma, H.; Qiu, Y.L.; Sun, H.; Meng, Y.; Sun, W.B.; Zimmer, E.A. Phylogenetic and biogeographic complexity of Magnoliaceae in the Northern Hemisphere inferred from three nuclear data sets. Mol. Phylogenet. Evol. 2008. [Google Scholar] [CrossRef]
  81. Massoni, J.; Couvreur, T.L.; Sauquet, H. Five major shifts of diversification through the long evolutionary history of Magnoliidae (angiosperms). BMC Evol. Biol. 2015, 15, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Azuma, H.; García-Franco, J.G.; Rico-Gray, V.; Thien, L.B. Molecular phylogeny of the Magnoliaceae: The biogeography of tropical and temperate disjunctions. Am. J. Bot. 2001, 88, 2275–2285. [Google Scholar] [CrossRef] [PubMed]
  83. 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.; et al. Origin of angiosperms and the puzzle of the Jurassic gap. Nat. Plants 2019, 5, 461–470. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Mitochondrial genome organization of Megabelesesinae with reference to the ancestral type of insect mitochondrial genomes. Genes transcribed from the J- and N-strands are shown in green and orange, respectively. The A + T-rich region is indicated by blue, and tRNA genes are labeled by their single-letter amino acid code.
Figure 1. Mitochondrial genome organization of Megabelesesinae with reference to the ancestral type of insect mitochondrial genomes. Genes transcribed from the J- and N-strands are shown in green and orange, respectively. The A + T-rich region is indicated by blue, and tRNA genes are labeled by their single-letter amino acid code.
Insects 12 00495 g001
Figure 2. Relative synonymous codon usage (RSCU) of the Megabelesesinae mitogenomes. The RSCU values are shown as the accumulative bar diagrams; from left to right, column shows Megabeleses magnoliae, M. liriodendrovorax, Cladiucha punctata, and C. magnoliae, respecitively. The accumulative bar diagrams of codon families are provided on the x-axis.
Figure 2. Relative synonymous codon usage (RSCU) of the Megabelesesinae mitogenomes. The RSCU values are shown as the accumulative bar diagrams; from left to right, column shows Megabeleses magnoliae, M. liriodendrovorax, Cladiucha punctata, and C. magnoliae, respecitively. The accumulative bar diagrams of codon families are provided on the x-axis.
Insects 12 00495 g002
Figure 3. Predicted secondary structures of 22 tRNA genes of Cladiucha. Dashes indicate Watson-Crick base pairs, and dots indicate G-U base pairing.
Figure 3. Predicted secondary structures of 22 tRNA genes of Cladiucha. Dashes indicate Watson-Crick base pairs, and dots indicate G-U base pairing.
Insects 12 00495 g003
Figure 4. Predicted rrnS secondary structure in the Cladiucha mitochondrial genome. The numbering of helices follows Gillespie et al. [52]. Roman numerals refer to domain names. Tertiary inter-actions and base triples are connected by continuous lines. C. magnoliae as a basemap and base change among Cladiucha species are presented in circles with red (C. magnoliae) and pink (C. punctata) colors.
Figure 4. Predicted rrnS secondary structure in the Cladiucha mitochondrial genome. The numbering of helices follows Gillespie et al. [52]. Roman numerals refer to domain names. Tertiary inter-actions and base triples are connected by continuous lines. C. magnoliae as a basemap and base change among Cladiucha species are presented in circles with red (C. magnoliae) and pink (C. punctata) colors.
Insects 12 00495 g004
Figure 5. Predicted rrnL secondary structure in the Cladiucha mitochondrial genome. The numbering of helices follows Gillespie et al. [52]. Roman numerals refer to domain names. Tertiary inter-actions and base triples are connected by continuous lines. C. magnoliae as a basemap and base change among Cladiucha species are presented in circles with red (C. magnoliae) and pink (C. punctata) colors.
Figure 5. Predicted rrnL secondary structure in the Cladiucha mitochondrial genome. The numbering of helices follows Gillespie et al. [52]. Roman numerals refer to domain names. Tertiary inter-actions and base triples are connected by continuous lines. C. magnoliae as a basemap and base change among Cladiucha species are presented in circles with red (C. magnoliae) and pink (C. punctata) colors.
Insects 12 00495 g005
Figure 6. Phylogenetic tree recovered using the dataset P123RNA using Bayesian inference and maximum likelihood methods. Only support values <1.0 (posterior probabilities) and <100% (bootstraps) are shown.
Figure 6. Phylogenetic tree recovered using the dataset P123RNA using Bayesian inference and maximum likelihood methods. Only support values <1.0 (posterior probabilities) and <100% (bootstraps) are shown.
Insects 12 00495 g006
Figure 7. Dated phylogeny constructed with P123RNA dataset in BEAST v1.8.3. The axis on the bottom refers to million years and shows the geological time. The blue bars on the nodes represent 95% of the high posterior density of divergence times obtained from the MCMC tree. The divergence times of each node obtained from the MCMC tree analysis are written in black on the nodes. The yellow shaded area delineates the arising of Tenthredinoidea, while the green and purple shaded areas delineate the arising of Megabelesesinae and Cladiucha, respectively, depending on the minimum or maximum ages. The traits in green and purple boxes are their synapomorphies, respectively. The phylogenetic trees given at the bottom are the angiosperm tree and the Magnolia tree, respectively.
Figure 7. Dated phylogeny constructed with P123RNA dataset in BEAST v1.8.3. The axis on the bottom refers to million years and shows the geological time. The blue bars on the nodes represent 95% of the high posterior density of divergence times obtained from the MCMC tree. The divergence times of each node obtained from the MCMC tree analysis are written in black on the nodes. The yellow shaded area delineates the arising of Tenthredinoidea, while the green and purple shaded areas delineate the arising of Megabelesesinae and Cladiucha, respectively, depending on the minimum or maximum ages. The traits in green and purple boxes are their synapomorphies, respectively. The phylogenetic trees given at the bottom are the angiosperm tree and the Magnolia tree, respectively.
Insects 12 00495 g007
Table 1. Summary of the phylogenetic relationships recovered by different datasets and inference approaches.
Table 1. Summary of the phylogenetic relationships recovered by different datasets and inference approaches.
Datasets Within TenthredinoideaWithin Tenthredinidae
Inference MethodsBIMLBIML
P123RNA(P + A) + (Ath + ((D + C) + (H + T)))(P + A) + (Ath + (D + C) + (H + T)))N + (S + (M + ((Ht + (Al + B)) + (Ht + B)) + (Th + Al))))N + (S + (M + ((Ht + (Al + B)) + (Ht + B)) + (Th + Al))))
P12RNA(P + A) + (Ath + ((D + C) + (H + T)))(P + A) + (Ath + ((D + C) + (H + T)))N + (S + (M + ((Ht + (Al + B)) + (Ht + B)) + (Th + Al))))N + (S + (M + ((Ht + (Al + B)) + (Ht + B)) + (Th + Al))))
P123RNAexc3genes(P + A) + (Ath + ((H + (D + C)) + T))(P + A) + (Ath + ((D + C) + (H + T)))N + (S + (M + ((Ht + (Al + B)) + (Ht + B)) + (Th + Al))))N + (S + (M + ((Ht + (Al + B)) + (Ht + B)) + (Th + Al))))
P12RNAexc3genes(P + A) + (Ath + ((H + (D + C)) + T))(P + A) + (Ath + ((H + (D + C)) + T))N + (S + (M + (((Ht + (Al + B)) + (Ht + B)) + (Th + Al)))N + (S + (M + (Th + (((Ht + (Al + B)) + (Ht + B)) + Al))))
ML, maximum likelihood; BI, Bayesian inference; P123, all codon positions of PCGs; P12, 1st and 2nd codon positions of PCGs; RNA, nucleotide sequences of rrnS, rrnL and tRNA genes; ‘exc3genes’ indicates that nad4l, nad6 and atp8 excluded from analyses; P, Pergidae; A, Argidae; Ath, Athaliinae; D, Diprionidae; C, Cimbicidae; H, Heptamelidae; T, Tenthredinidae; N, Nematinae; S, Selandriinae; M, Megabelesinae; Ht, Heterarthrinae; Al, Allantinae; B, Blennocampinae; Th, Tenthredininae. Please see Figure S1 for details.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Niu, G.; Jiang, S.; Doğan, Ö.; Korkmaz, E.M.; Budak, M.; Wu, D.; Wei, M. Mitochondrial Phylogenomics of Tenthredinidae (Hymenoptera: Tenthredinoidea) Supports the Monophyly of Megabelesesinae as a Subfamily. Insects 2021, 12, 495. https://doi.org/10.3390/insects12060495

AMA Style

Niu G, Jiang S, Doğan Ö, Korkmaz EM, Budak M, Wu D, Wei M. Mitochondrial Phylogenomics of Tenthredinidae (Hymenoptera: Tenthredinoidea) Supports the Monophyly of Megabelesesinae as a Subfamily. Insects. 2021; 12(6):495. https://doi.org/10.3390/insects12060495

Chicago/Turabian Style

Niu, Gengyun, Sijia Jiang, Özgül Doğan, Ertan Mahir Korkmaz, Mahir Budak, Duo Wu, and Meicai Wei. 2021. "Mitochondrial Phylogenomics of Tenthredinidae (Hymenoptera: Tenthredinoidea) Supports the Monophyly of Megabelesesinae as a Subfamily" Insects 12, no. 6: 495. https://doi.org/10.3390/insects12060495

APA Style

Niu, G., Jiang, S., Doğan, Ö., Korkmaz, E. M., Budak, M., Wu, D., & Wei, M. (2021). Mitochondrial Phylogenomics of Tenthredinidae (Hymenoptera: Tenthredinoidea) Supports the Monophyly of Megabelesesinae as a Subfamily. Insects, 12(6), 495. https://doi.org/10.3390/insects12060495

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop