Next Article in Journal
Fungi Present in the Organic and Mineral Layers of Six Broad-Leaved Tree Plantations as Assessed by the Plate Dilution Method
Next Article in Special Issue
Systematics, Ecology and Taxonomy of Collembola: Introduction to the Special Issue
Previous Article in Journal
Use of Human Dominated Landscape as Connectivity Corridors among Fragmented Habitats for Wild Asian Elephants (Elephas maximus) in the Eastern Part of Thailand
Previous Article in Special Issue
A New Genus of Sminthurididae (Collembola, Symphypleona) from Brazil, with Notes on the Systematics of the Family
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Evolution of Collembola Higher Taxa (Arthropoda, Hexapoda) Based on Mitogenome Data †

by
Bruno Cavalcante Bellini
1,*,
Feng Zhang
2,
Paolla Gabryelle Cavalcante de Souza
1,
Renata Clicia dos Santos-Costa
1,
Gleyce da Silva Medeiros
1 and
Nerivânia Nunes Godeiro
3,*
1
Department of Botany and Zoology, Biosciences Center, Federal University of Rio Grande do Norte (UFRN), Highway BR-101, Lagoa Nova, Campus Universitario, Natal 59072-970, RN, Brazil
2
Department of Entomology, College of Plant Protection, Nanjing Agricultural University, Nanjing 210095, China
3
Shanghai Natural History Museum, Shanghai Science and Technology Museum, Shanghai 200041, China
*
Authors to whom correspondence should be addressed.
This published work and the nomenclatural acts it contains have been registered in ZooBank, the online registration system for the ICZN (International Code of Zoological Nomenclature). The LSID (Life Science Identifier) for this publication is: urn:lsid:zoobank.org:pub:E44A878B-F6AB-4BC0-929B-E23695F2DC43.
Diversity 2023, 15(1), 7; https://doi.org/10.3390/d15010007
Submission received: 18 November 2022 / Revised: 15 December 2022 / Accepted: 16 December 2022 / Published: 21 December 2022
(This article belongs to the Special Issue Systematics, Ecology and Taxonomy of Collembola)

Abstract

:
Mitogenomes represent useful tools for investigating the phylogeny of many metazoan clades. Regarding Collembola, the use of mitogenomics has already shown promising results, but few published works include sufficient taxon sampling to study its evolution and systematics on a broader scale. Here, we present a phylogenetic study based on the mitogenomes of 124 species from 24 subfamilies, 16 families, and four orders—one of the most comprehensive datasets used in a molecular study of Collembola evolution to date—and compare our results with the trees from recently published papers and traditional systematic hypotheses. Our main analysis supported the validity of the four orders and the clustering of Poduromorpha with Entomobryomorpha (the traditional Arthropleona). Our data also supported the split of Symphypleona s. str. into the Appendiciphora and Sminthuridida suborders, and the division of the Neelipleona into two subfamilies: Neelinae and Neelidinae subfam. nov. On the other hand, the traditional Symphypleona s. lat., Isotomoidea, and all the Isotomidae subfamilies were refuted by our analyses, indicating a need for a systematic revision of the latter family. Though our results are endorsed by many traditional and recent systematic findings, we highlight a need for additional mitogenomic data for some key taxa and the inclusion of nuclear markers to resolve some residual problematic relationships.

1. Introduction

Springtails (Collembola) are tiny arthropods, with the majority of the species smaller than 5 mm [1]. Because of this limited size, they are unknown to most people. Even so, springtails are among the most widely distributed and abundant extant terrestrial animals [2]. Also, Collembola represents the richest lineage among the early diverging (“apterygotan”) hexapod lineages, gathering more species than Protura, Diplura, Archaeognatha, and Zygentoma combined. Compared to other non-Holometabola hexapods, this clade is only smaller in number of species than the Hemiptera and Orthoptera insect orders [3]. Even with more than 9400 described species, many others are yet to be described [2,4].
The evolution and systematics of the Collembola main lineages have been debated in several revisions for more than a century. Many traditional propositions to organize the supraspecific taxa and infer their evolution have been discussed by different authors. Such hypotheses were based on the external morphology of adult specimens, first instar and post-embryonic development patterns, and even internal structures. However, they mostly resulted in discrepant propositions, showing little consensus compared to each other [5,6,7,8,9,10,11,12,13,14,15,16,17,18]. The largest internal groups of springtails, nowadays considered as orders, were also continuously reorganized. They were first suggested as only two main lineages: the Arthropleona and Symphypleona sensu Börner, 1901 [5,6], going up to five orders, including the ambiguous Metaxypleona [9,16]. More recently, different studies agreed that there are four main well-delimited lineages: Neelipleona, Symphypleona, Entomobryomorpha, and Poduromorpha [19,20], but even so, some species and genera show so many autapomorphies that their placements question the validity and boundaries of such orders [8,16,21].
With the advent of molecular tools and modern bioinformatics to investigate biological evolution, many traditional and more recent views on the systematics of the Collembola were put to the test. Taxonomic evolution, affinities, boundaries, and validity were investigated at different scales, confirming or refuting previous hypotheses. Starting from a smaller taxonomical scope, intra- and interspecific levels of genetic variation were investigated to provide data on species delimitation, population structuring, and the existence of cryptic species [4,22,23,24,25,26,27,28,29,30,31,32,33]. Based on similar methodologies, species–groups affinities and validity were also analyzed [34,35,36,37]. The placement and relationships of genera within subfamilies, as well as their monophyletic status, were also evaluated using molecular-based approaches [38,39,40,41]. Likewise, the limits, affinities, and main diagnostic features of subfamilies, families, and superfamilies were the main subject of many recent studies, especially among the Entomobryomorpha [42,43,44,45,46,47]. Finally, on a broader scale, the relationships of the current four Collembola orders were also recently investigated, but so far they show very little consensus regarding their affinities [27,28,35,46,48,49,50,51,52,53,54,55,56,57].
Mitogenomes are useful tools for studying the phylogeny of many metazoan lineages. They can provide a large set of useful comparative data for arthropod systematics, especially when confronted with single locus markers [55,58,59,60]. Regarding Collembola, the use of mitogenomics already showed promising results, but most of the published papers were limited to a small sample size and/or focused on a specific higher taxon [41,52,53,55,57,60,61,62,63,64]. The main exception is the recent study of Cucini et al. [56], which provided a wider view of Collembola phylogenetics based on a large dataset of mitogenome sequences, using first and second codon positions in the analyses. In this work, many important observations on the use of mitogenomes for phylogenetic investigation and Collembola gene orders were provided, but the authors opted to discuss their systematic findings in light of recent molecular studies, avoiding in-depth comparisons of their results with traditional hypotheses or looking for morphological evidence to support their results. Additionally, some incongruence between traditional systematics and the obtained results was observed, like the finding of polyphyletic Tomoceridae, Neelidae, and Entomobryomorpha.
In the current study, we aimed to provide a large-scale investigation of the main clades within Collembola based on mitogenomes, using a broad set of species, genera, subfamilies, and families of its four orders. In contrast to Cucini et al. [56], we used amino acid sequences in our analyses, similarly to Sun et al. [53]. Such efforts resulted in one of the most representative phylogenies for the entire class so far, facilitating our ability to investigate the affinities of higher taxa and, at the same time, confirming previously published hypotheses concerning smaller clades. Additionally, we provided a comprehensive review of the main previous systematic hypotheses based on morphology and molecular data and how they compare with our results.

2. Materials and Methods

2.1. Species and Mitogenomic Data Acquisition

We surveyed the NCBI database in May 2022 and retrieved mitogenomic data for 123 collembolan, two dipluran, and two proturan species for analyses. Our dataset gathers data published in 38 studies, with 12 sequences never used in mitogenomic phylogenetic studies before (Ascocyrtus cinctus Schäffer, Desoria tigrina (Tullberg), Lepidosira calolepsis (Börner), Plutomurus gul (Yosii), Pogonognathelus flavensis (Tullberg), Pogonognathelus longicornis (Müller), Pseudobourletiella spinata (MacGillivray), Seira boneti Denis, Sminthurides aquaticus (Bourlet), Tomocerus maximus (Liu, Hou and Li), Tomocerus nigrus Sun, Liang and Huang, and Tomocerus vulgaris (Tullberg)). Only one species newly sequenced was included in our dataset: Seira boneti. The paper with its mitogenome announcement and description is under preparation.
A few mitogenomes were previously discarded from our dataset due to the absence of some coding genes or incorrect phylogenomic placement (possible identification errors). Only mitochondrial protein coding genes were used for the analyses. In total, the ingroup included 124 taxa of Collembola: three species representing the Neelipleona, nine Symphypleona, 15 Poduromorpha, and 97 Entomobryomorpha, totaling 24 sampled subfamilies and 16 families. Four outgroups, two species of Protura and two of Diplura, were chosen based on molecular studies focused on basal hexapods, which supported a closer relationship between these three groups [49,50,65,66,67].
The detailed classification information and voucher numbers of the 128 species analyzed in this study are listed in Table S1.

2.2. Phylogenetic Analyses

Firstly, we organized matrices of nucleotide data (with and without the third codon) and ran a priori analyses. Nucleotide data were investigated with different alignment and trimming tools/parameters. At the beginning of each run, IQTree performs a composition chi-square test to verify the homogeneity of the character composition in the alignment. This test showed that 87 sequences failed to achieve the optimal value of chi2, suggesting possible inconsistencies with our nucleotide matrices. The resulting ML trees of nucleotide data, with or without the third codon, failed to support Collembola (with the outgroup clustered together with Actaletidae using the partitioned dataset), or to support more than two orders at the same time (with Neelipleona and Entomobryomorpha as polyphyletic taxa in the same topology using the unpartitioned dataset). Because of this, we dismissed nucleotide data from our final analyses.
TransDecoder v5.5.0 [68] was used to translate the nucleotide sequences of the 13 protein-coding genes (PCGs) of each species into amino acids. MAGUS strategy [69], employing MAFFT [70], was used to align each PCG independently. BMGE v1.12 [71] performed the trimming, with the defaults. PhyKIT v1.9.0 [72] was used to concatenate the genes and generate the matrix and partition scheme. Maximum Likelihood (ML) analyses were performed with IQTREE v2.0.7 [73], 1000 ultrafast bootstrap [74], and SH-aLRT replicates. The “edge-linked-proportional partition model with separate substitution models and separate rates across sites” method was used for two ML partitioned analyses: one tree was constructed using ModelFinder [75] to select the best-fitting substitution model for each gene partition (Table S2), and another tree was made with the mixture model (EX-EHO). The resultant tree of the mixture model was used as a guide for the unpartitioned method with the posterior mean site frequency (PMSF) model option “-m mtART + C60 + FO + R” [76]. Additionally, the traditional unpartitioned method was tested with the best model (mtZOA + F + R10) selected by ModelFinder. The Bayesian Inference (BI) tree was reconstructed using Phylobayes-MPI v1.8 [77], unpartitioned dataset, default model CAT+GTR with discrete gamma (four categories). Two chains were run until the likelihood had satisfactorily converged (maxdiff < 0.3). The consensus tree was generated using a burn-in of 1000 trees and subsampling every 10 trees. Maximum Parsimony (MP) analysis was performed using MPBoot v1.1.0 [78], with 1000 bootstrap replicates. All resulting phylogenies were visualized and first edited in FigTree v1.3.1 [79]. Concerning ML and MP analyses, we considered clades with 95–100 of bootstrap as resolved (with high node support), 90–94.9 as acceptable, and lower than 90 as unresolved. For BI analysis, any clade with less than a 0.95 posterior probability was considered as unresolved.
A priori, we considered the different types of analyses, partition schemes, and models as equally capable of solving the phylogeny. To choose the main topology, we defined the following criteria: (1) all four orders should be recovered as monophyletic groups; (2) superfamilies and families with strong morphological evidence supporting their validities should also be recovered, as the Neanuroidea, Entomobryoidea, Sminthuroidea, Sminthurididae, and Neelidae; (3) since most sampled taxa belong to Entomobryoidea, the higher Entomobryidae (the group with elongate fourth abdominal segment) should also be recovered, following the findings of [43,44,45,60]; and (4) an overall higher node support for the main relationships, above the subfamily level, should be recovered (only bootstrap for ML analyses [74]). Due to the high number of terminals, we manually reconstructed all trees highlighting the sampled springtails higher taxa (main subfamilies, families, superfamilies, suborders, and orders, and a few genera). To better illustrate these branches, we drew potential representatives of each lineage loosely based on the photograph repository of [1] using CorelDraw 2021, except for Spinactaletes Soto-Adames (Actaletidae), which was adapted from [80], and the depictions of the orders, adapted from [81].

2.3. Tree Topology Tests

Using the RELL approximation method with 10,000 replicates [74], we performed tree topology tests on the constraining monophyly of some of the suprageneric taxa of Collembola. The following six hypotheses were proposed: (1) best tree without any constraints (Figure 1); (2) Monophyly of Orchesellidae (including Orchesellinae and Heteromurinae as independent taxa) and Paronellidae; (3) Monophyly of Hypogastruridae; (4) (Sminthuridae + Sminthurididae) + Bourletiellidae; (5) Sminthurididae + Sminthuroidea; (6) (Bourletiellidae + Dicyrtomidae) + (Sminthuridae + Sminthurididae). Kishino-Hasegawa (KH), Shimodaira-Hasegawa (SH), weighted KH (WKH) and weighted SH (WSH), expected likelihood weight (ELW), and approximately unbiased (AU) tests were performed in IQTREE.

3. Results

3.1. Matrices and Trees

The data matrix used for the phylogenetic analyses included 128 taxa and 3003 amino acid sites representing the 13 mitochondrial protein-coding genes. Six trees were generated based on different models or phylogenetic inference methods. Among our analyses, the ML tree, unpartitioned, under the model mtZOA + F + R10 (best model suggested by Model Finder), better represented the currently accepted classifications of the higher taxa of Collembola in topology and overall node support, as discussed ahead. This phylogeny was chosen as our main tree and is represented in Figure 1 and Figure S6 (the latter with detailed node support values). The other five trees are depicted in the Supplementary Materials as Figures S1–S5. A synthesis of the different analyses, parameters, tree codes, and related figures is presented in Table 1.

3.2. Phylogeny of the Orders

The ordinal relationships of all obtained trees are summarized in Figure 2. Our main phylogeny recovered the four Collembola orders as monophyletic taxa, with the following topology: Neelipleona + (Symphypleona + (Entomobryomorpha + Poduromorpha)) (Figure 1 and Figure 2A). Two other ML trees also reached this same topology concerning the orders (ML_2 and ML_4), while ML_3 resolved (Neelipleona + Symphypleona) + (Entomobryomorpha + Poduromorpha) (Figure 2A,B and Figures S1–S3). Our BI and MP trees did not recover the monophyly of all orders, with Neelipleona as an ingroup of Poduromorpha in the former, and Actaletidae outside of Entomobryomorpha in the latter (Figure 2C,D, Figures S4 and S5).
Entomobryomorpha + Poduromorpha, the traditional Arthropleona grouping, was strongly supported by our main phylogeny. The bootstrap support for Symphypleona as its sister group was lower (91), slightly below of a 95 optimal threshold, while the SH-aLRT support for this relationship was only 58.1. On the other hand, the Neelipleona appeared as the most basal order of Collembola, with absolute node support (Figure 1, Figure 2A and Figure S6). Similar results were obtained for ML_4 (Figure S3), while ML_2 reached the same order topology but with higher support for the Arthropleona + Symphypleona clade (Figure S1). The ML_3 tree also recovered the Arthropleona with bootstrap support = 95 and SH-aLRT = 89.7, while the clade Symphypleona + Neelipleona was recovered with very low support levels (Figure S2). The Bayesian Inference tree found a very low posterior probability value for Neelipleona as an ingroup of Poduromorpha (Figure S4), while the Entomobryomorpha lacking Actaletidae was recovered with a bootstrap value of only 59 in the MP analysis (Figure S5).
Our main tree supported all orders with bootstrap values above 95, except for Neelipleona, with 94 (Figure 1 and Figure S6). The Maximum Likelihood analyses using other models achieved a higher node support for Neelipleona (Figures S1–S3), while BI and MP reached a posterior probability value of 0.63, and parsimony bootstrap of 80, respectively (Figures S4 and S5).

3.3. Entomobryomorpha

The phylogenies of the Entomobryomorpha higher taxa are summarized in Figure 3. Our main tree supported the topology Tomoceridae + (Actaletidae + (Isotomidae + Entomobryoidea)), as well as our BI analysis (Figure 1, Figure 3A and Figure S4). Similar trees were obtained by other ML models, but with different internal organizations for the Entomobryoidea (Figure 3B,D and Figures S1–S3), while the MP tree recovered the Entomobryomorpha without Actaletidae (Figure 3E and Figure S5). All trees supported Tomoceridae as the most basal branch of Entomobryomorpha, and did not support the Isotomoidea sensu Soto-Adames et al. [18], clustering Isotomidae and Actaletidae, and all subfamilies of Isotomidae, with Pachyotominae (Paranurophorus Denis), Proisotominae (Proisotoma Börner and Subsitoma Stach), Anurophorinae (Anurophorus Nicolet, Cryptopygus Willem, Folsomia Willem, and Tetracanthella Schött) mixed and as ingroups of Isotominae (the other isotomid genera) (Figure 1 and Figures S1–S5).
The ML_1 and BI trees found Entomobryidae s. lat. topology as Entomobryinae s. lat. + (Paronellinae + (Lepidocyrtinae + Seirinae)), with Entomobryinae s. lat. gathering the Salininae (Akabosia Kinoshita, Callyntrura Börner, and Salina MacGillivray) and the Entomobryinae sensu Zhang et al. [44]. The same trees did not support the Orchesellidae and Orchesellinae, while Heteromurinae appeared as an ingroup of part of Orchesellinae, and Nothobryinae was recovered as the most basal branch of Entomobryoidea (Figure 1, Figure 3A, and Figure S4). Similar results were found by ML_4, but with Callyntrura as an ingroup of Paronellinae (Figure 3D and Figure S3). The ML_2, ML_3 and MP trees recovered the Orchesellidae, but failed to support the Orchesellinae as well (Figure 3B,C,E, Figures S1, S2 and S5). Tree topology tests aimed at investigating the Entomobryoidea inner relationships also rejected, with full confidence, both Orchesellidae and the traditional Paronellidae in the same topology. Such tests also refuted independent Orchesellinae and Heteromurinae (Table S3). The main phylogeny support levels were high or acceptable for the majority of the main Entomobryoidea branches, with the most important exceptions of the Entomobryinae s. lat. cluster, with a bootstrap support of 85 and an SH-aLRT of 67.3, and Entomobryidae s. lat. + Orchesella Templeton, with a bootstrap support of 79 and an SH-aLRT of 80.1 (Figure 1 and Figure S6). The Bayesian Inference analysis, which found the most similar topology to the main tree, recovered Entomobryidae s. lat. with very high node support, but the clade Paronellinae + (Lepidocyrtinae + Seirinae) was supported by a posterior probability value of only 0.65 (Figure S4).

3.4. Poduromorpha

The internal relationships of the Poduromorpha higher clades of all trees are summarized in Figure 4. All ML analyses supported the following internal topology for the order: Onychiuridae + (Ceratophysella Börner + (Gomphiocephalus Carpenter + (Poduridae + Neanuroidea))) (Figure 1, Figure 4A and Figures S1–S3). A close topology was obtained by MP analysis, but the internal organization of the Onychiuridae and Neanuroidea subfamilies was different from ML results (Figure 4A,C and Figure S5). Bayesian Inference recovered the Neelidae as an ingroup of Poduromorpha, and as the sister group of the Onychiuridae (Figure 4B and Figure S4). All analyses dismissed the Hypogastruridae, represented by Ceratophysella and Gomphiocephalus, and included the Brachystomellidae as an ingroup of the Neanuridae, represented by the Frieseinae, Pseudachorutinae and Neanurinae subfamilies (Figure 1, Figure 4, and Figures S1–S5). Hypogastruridae was also rejected by most of the tree topology tests (Table S3).
All internal nodes of Poduromorpha had high support in our main analysis, with the exception of the SH-aLRT value supporting the Onychiurinae (represented by Allonychiurus Yoshii, Orthonychiurus Stach, and Thalassaphorura Bagnall), a result similar to the ML_2 tree (Figure 1, Figures S1 and S6). Low node support values were found by BI clustering the Neelidae, Neelidae + Onychiuridae, and Poduromorpha + Neelipleona clades (Figure S4). The Maximum Parsimony analysis had mixed bootstrap values, with some internal nodes of Neanuroidea and Onychiuridae with low support (Figure S5).

3.5. Symphypleona and Neelipleona

The relationships of the Symphypleona families and Neelipleona genera of all the obtained trees are summarized in Figure 5 and Figure 6, respectively. Concerning the Symphypleona, in all analyses the sampled families were recovered as monophyletic independent groups (Figure 1, Figure 5 and Figures S1–S5). All ML models reached the same results for the order, with the Sminthurididae (suborder Sminthuridida sensu Sánchez-García and Engel [82]) as the sister group of the suborder Appendiciphora sensu Bretfeld [17]. Within the latter group, Dicyrtomidae was found to be the sister group of the Sminthuroidea (Sminthuridae + Bourletiellidae) (Figure 1, Figure 5A and Figures S1–S3). Alternative hypotheses which disregarded the suborders of Symphypleona were recovered by BI: Dicyrtomidae + (Sminthurididae + (Sminthuridae + Bourletiellidae)); and MP: (Dicyrtomidae + Bourletiellidae) + (Sminthuridae + Sminthurididae) (Figure 5B,C, Figures S4 and S5). The main tree showed high node support values for all the internal branches of Symphypleona, with the exception of an SH-aLRT value of 66.2 for the Sminthuroidea clade (Figure 1 and Figure S6). Other ML analyses also achieved lower SH-aLRT support for this node, as well as for other internal groups of Appendiciphora (Figures S1–S3). Tree topology tests better supported two competing hypotheses: (1) Appendiciphora + Sminthuridida (the unconstrained ML_1 tree); and (6) (Bourletiellidae + Dicyrtomidae) + (Sminthuridae + Sminthurididae) (found by MP tree) (Table S3).
Regarding the Neelipleona (Neelidae), all trees recovered the same topology for the sampled genera: Neelides Caroli + (Neelus Folsom + Megalothorax Willem) (Figure 1, Figure 6, and Figures S1–S5). The support values were mostly high for both nodes within the order in all ML analyses, except for an SH-aLRT support of 58.6 in the basal branch of ML_1 tree (Figure 1 and Figure S6). Bayesian Inference found a low value of posterior probability for the node containing the Neelipleona within Poduromorpha (Figure S4), while in the MP analysis, the bootstrap support for the order was 80 (Figure S5).
Finally, Figure 7 translates the full view of the Collembola higher taxa relationships obtained in our main ML phylogeny, ML_1 (Figure 1).

4. Discussion

4.1. Phylogeny of the Orders

The validity and affinities of the current four Collembola orders have been subjects of many systematic revisions and phylogenetic studies. Although some previous published trees, as well as our BI and MP analyses, did not recover the monophyly of all four orders, the internal and at least part of the external morphology advocate that they are distinct monophyletic lineages [10,12,13,15,17,18,83,84]. Even the puzzling Podura Linnæus, Mackenziella Hammer, and Actaletidae, which were previously allocated in the distinct order Metaxypleona by Salmon [9], have morphological traits which assign them to the Poduromorpha, Symphypleona, and Entomobryomorpha, respectively [18,21,84,85]. Molecular phylogenetic studies based on mitogenomes or other markers also support the current systematic placing of Podura and Actaletidae, while the placement of Mackenziella has yet to be evaluated with genetic data [27,28,35,46,48,49,50,51,52,53,55,56,63,86].
Assuming the four orders are monophyletic taxa, a more important question arises: what are their relationships? A traditional view provided by Börner proposed the higher clade Arthropleona to gather the Poduromorpha and Entomobryomorpha, and Symphypleona s. lat. to gather the current Symphypleona s. str. and Neelipleona [5,6]. Such division was mainly based on body shape and the presence or absence of clear abdominal segments. Concerning the Arthropleona, more recently published trees supported this clade to some extent [27,35,48,49,50,51,53], while others refuted it [28,46,54,55,56,57,87]. Even our results did not reach a consensus for the relationships of these orders, with all ML trees supporting the Arthropleona, while the BI and MP trees rejected this clade (Figure 2). At first, it may be tempting to consider the Arthropleona as a monophyletic taxon compared to the very dissonant trunk morphologies of Symphypleona s. str. and Neelipleona. However, the body tagmosis of Poduromorpha and Entomobryomorpha is likely a plesiomorphy and potentially has no phylogenetic signal to group them, as a segmented ancestral state for the Collembola is expected when observing the relationships of basal hexapods and the segmented Protura and Diplura bauplans [65,66,67]. In this scenario, for now, there is no clear, striking morphological trait that could be assigned as a synapomorphy of Arthropleona, whereas the loss of the pronotum in Entomobryomorpha is even shared with Neelipleona and Symphypleona s. str. (except for Mackenziellidae), and the first instar tibiotarsal chaetotaxy may point to a closer relationship between the Entomobryoidea and Symphypleona [1,54,88]. So, the morphological support for this higher clade is, for now, unfulfilling. If the Arthropleona clade is valid, the synapomorphies of the group are likely to be more discrete and can be related to the chaetotaxy or even to the internal morphology, as studied by Cassagnau [12].
Contrastingly, the clade Symphypleona s. lat. has less support based on more recent molecular-based phylogenetic studies, and it was recovered only by Xiong et al. [51] and Schnider et al. [28], while it was refuted by most of the published trees, including those phylogenies relying on mitogenomes [27,35,46,50,54,55,56,57]. The findings of Sun et al. [53] were conflicting regarding this clade, with the phylogeny based on mitogenomes partitioned dataset of nucleotide sequences dismissing the Symphypleona s. lat., whereas the tree inferred from amino acid sequences supporting it, but with a bootstrap value of 73.9. In our analyses, only the ML_3 tree recovered Symphypleona s. str. + Neelipleona, but with very low node support (Figure 2 and Figure S2). Bretfeld [17,83] suggested three apomorphies for the Symphypleona s. lat.: thoracic and first abdominal segments fused to some extent, resulting in a globular body; the presence of neosminthuroid chaetae on the parafurcal area (furca basis), which were posteriorly lost in different internal lineages of Symphypleona s. str. and Neelipleona, and the presence of a gutterlike mucro. Other morphological traits, however, may refute a closer relationship between the two orders. The spherical trunk seen in Neelipleona is likely not homologous to the one of Symphypleona s. str., as, in the former, the large abdomen is mostly constituted by thorax II and III, with a reduction of abdominal segments in size and chaetotaxy (at least in most genera), while the Symphypleona s. str. have a large abdomen mostly formed by abdominal segments I–IV, holding a more complex chaetotaxy compared to the thorax [28,35,83,89,90,91]. Other features of the trunk chaetotaxy are also not clearly comparable between both orders, like the presence of papillate long bothriotricha on the large and small abdomens of Symphypleona s. str., which already emerge in the first instar. In Neelipleona, if these are present, they are short discrete non-papillate bothriotricha, clearly not homologous compared to the latter [28,35,90,91]. The complex chaetotaxy of the small abdomen seen in the Symphypleona s. str., even during the first instar, does not match the strongly reduced chaetotaxy of the Neelipleona [28,83,91,92,93]. The tibiotarsal chaetotaxy of Neelidae is also strongly reduced compared to the Symphypleona first instar, and it may be more in line with the Poduromorpha [28,35,91,93,94,95,96]. Finally, the reduced antennae, long coxae, tenaculum shape, and midgut morphology of the Neelipleona advocate they may not be related to the Symphypleona s str. as well [1,83,91].
Our main tree topology, Neelipleona + (Symphypleona + (Entomobryomorpha + Poduromorpha)) (Figure 1 and Figure 2), was also recovered by Gao et al. [50], and was similar to Sun et al. [53] tree based on mitogenomic partitioned nucleotide sequences. The disparity of tree topologies regarding Collembola orders in recent studies is high and does not point to any clear consensus [56]. We believe the inclusion of large sets of nuclear markers would improve the results obtained by the use of mitogenomes or few isolated mitochondrial markers in this case. Nuclear DNA can be more resilient to mutation compared to mitochondrial genes, making the latter markers more fitting to investigate a more recent evolution [97]. In this case, deeper relationships within the Collembola phylogeny could be better investigated with the inclusion of high amounts of nuclear DNA, which hold the potential to unveil more solid data to evaluate the validity and affinities of the current four orders.

4.2. Entomobryomorpha

The Entomobryoidea higher clades have been investigated in many recent phylogenetic studies, which provided further grounds for their evolutionary affinities and insights into the phylogenetical signal of morphological traits on which traditional systematics relied on [38,39,40,41,42,43,44,45,56,60,61,62,63,98,99]. For instance, body scales have emerged more than once within the clade and have little or no phylogenetic signal within the Entomobryinae s. lat [43]. Similarly, the straight smooth dens has also emerged more than once within the superfamily, and the traditional Paronellidae sensu Soto-Adames et al. [18] is likely a polyphyletic group [39,43,44], as supported by our ML, BI, and MP analyses and tree topology tests. On the other hand, trunk chaetotaxy patterns show more reliable attributes to group the higher Entomobryoidea, especially concerning the sensilla [14,43,44,45].
The results of our ML_1 and BI trees (Figure 1, Figure 3A, and Figure S4) are similar to recently published studies. Our data support the Seirinae as the sister group of Lepidocyrtinae, similarly to Godeiro et al. [41,60,63] trees based on mitogenomes, but differently from the BI tree of Cucini et al. [56] based on the first and second codon positions. Other ML and MP trees also retrieved the same topology, reinforcing the validity of this clade (Figure 3B–E, Figures S1–S3 and S5). Morphology also better supports this relationship rather than the alternative hypothesis of Lepidocyrtinae + (Entomobryinae s. lat. + Seirinae) [14,44,100]. Within Seirinae, our data indicate that at least some oriental species of Seira Lubbock (Seira sanloemensis Godeiro and Cipola and S. boneti Denis) may represent an independent clade of Lepidocyrtinus Börner + (Seira + Tyrannoseira Bellini and Zeppelini), similar to the results obtained by Godeiro et al. [61]. Within the Lepidocyrtinae, we reached a very similar internal topology compared to Godeiro et al. [60], which was also based in mitogenomes. We found the Paronellinae closer to the Lepidocyrtinae, which was expected (see [44]), but as the sister group of Seirinae + Lepidocyrtinae, differently from other studies, such as those from Zhang et al. [39,43,44]. Nevertheless, our results better solved the position of the Paronellinae within Entomobryoidea when compared to other mitogenomic phylogenies [53,56,60], matching at some level morphology and the obtained topology based on a larger dataset of Paronellinae, Seirinae, and Lepidocyrtinae. If our results are confirmed in future studies, this would mean the resemblance of the Paronellinae to the Lepidocyrtinae, like the reduction of dorsal macroquetotaxy and the same sensillar pattern [14,44], are possibly due to the plesiomorphies of this clade, and the more complex dorsal macrochaetotaxy of Seirinae was achieved posteriorly within the lineage.
The Entomobryinae s. lat. was recovered by our ML_1 (Figure 1) and BI (Figure S4) trees with high node support and gather genera with a wide range of dorsal chaetotaxy and furca morphologies. This clade clusters unscaled and scaled genera, with highly variable dorsal main chaetotaxy, ranging from a polymacrochaetotic coverage (more common) to a reduced dorsal macrochaetotaxy [39,43,44]. Within this branch, the dens varies from crenulate to smooth, while the mucro morphology is somewhat similar to the Paronellinae in lineages like the Salininae sensu Zhang et al. [44], Zhuqinia Zhang, Ma and Greenslade, and Paronellides Schött, or more on par with the Entomobryidae s. str [43,44]. The main morphological feature grouping these very distinct taxa is the trunk sensillar pattern of 2, 2|1, 2, 2 from the mesothorax to the third abdominal segment. However, even this potential synapomorphy is secondarily modified in some internal branches, like the Salininae [44]. As in many other recent papers based on mitogenomes or other markers, our data could not clearly resolve some internal relationships within the subfamily, as the lower node support of our ML_1 and BI trees pointed out.
Our analyses could not solve the internal relationships of the Orchesellidae, with the ML_1, ML_4, and BI trees supporting a polyphyletic condition of the family, with some branches having low node support (Figure 1 and Figure 3A,D, Figures S3, S4 and S6). Our tree topology tests rejected the family as well (Table S3). The trees which recovered the Orchesellidae, ML_2, ML_3, and MP, could not separate the Heteromurinae from the Orchesellinae (Figure 3B,C,E, Figures S1, S2 and S5), a result also endorsed by our tree topology tests (Table S3). A similar problem was found in the BI tree of Godeiro et al. [60], which was also based on mitogenomes. Apparently, the lack of nuclear genes in the analyses obscures the relationships of the Orchesellidae, which were better solved in studies based on a few mitochondrial and nuclear markers combined, like in Zhang et al. [38,39,42,43] and Nunes et al. [40]. This may represent a similar issue to the one discussed in the previous topic, in which the absence of more mutation-resilient markers prevents a clearer understanding of deeper nodes within the phylogeny. Even so, it is worth noting that there is no clear synapomorphy to circumscribe the Orchesellidae, and the reduced abdomen and the presence of a postantennal organ in some genera are likely plesiomorphic features of the Entomobryoidea. On the other hand, morphology, to some extent, advocates that Orchesellinae and Heteromurinae are distinct lineages [14,60,98], a result not supported by our findings. We believe our limited dataset for the Orchesellidae and the mitogenomes alone are insufficient to clearly unveil the affinities and validities of the basal Entomobryoidea suprageneric taxa.
All our ML, BI, and MP trees did not support the current systematics of the Isotomoidea in two main lines: firstly, the internal division of the Isotomidae subfamilies as presented in Bellinger et al. [1] or as proposed by Potapov [101]; secondly, the positioning of Actaletidae next to Isotomidae, following Soto-Adames et al. [18]. The essential revision provided by Potapov [101] suggested three subfamilies for the Isotomidae: Pachyotominae, Anurophorinae, and Isotominae, disregarding Proisotominae due to its similarities with Anurophorinae. The author also highlighted the absence of “absolutely strict differences” between Anurophorinae and Isotominae and the absence of phylogenetic grounds to support his classification. Contrarily, the Bellinger et al. [1] database kept Proisotominae as a valid subfamily, separating it from Anurophorinae based on the number of body sensilla. The morphology within Isotomidae is remarkably variable, and the boundaries among the different subfamilies, genera, and subgenera may be difficult to determine in some cases [18,101]. Our data point to Pachyotominae, Anurophorinae, and Proisotominae as being internal groups of Isotominae. Moreover, Proisotoma and Subsitoma, the representative Proisotominae taxa in our analyses, were recovered apart from each other, within a clade gathering all sampled Anurophorinae (Anurophorus, Cryptopygus, Folsomia, and Tetracanthella) together with Paranurophorus, the sole Pachyotominae we sampled (Figure 1 and Figures S1–S5). Although such results support Potapov’s [101] view of merging Anurophorinae with Proisotominae, they did not sustain Pachyotominae as a full subfamily. Additionally, the clade mixing Pachyotominae, Anurophorinae, and Proisotominae genera is not independent as it emerges inside a higher group with Pseudisotoma Handschin, Semicerura Maynard, Parisotoma Bagnall, and Folsomotoma Bagnall (all Isotominae genera) in its base (Figure 1 and Figures S1–S5). Similar results were also obtained by some other studies dealing with Isotomidae taxa, with at least the Anurophorinae as an ingroup of a paraphyletic Isotominae [46,55,56,63]. Our data are not definitive by any means, as some of the internal nodes of Isotomidae had low support in our analyses (Figure 1 and Figures S1–S6). However, the high number of isotomid taxa used in this study is unmatched compared to the previously cited papers, providing some level of robustness to our findings. Such results point to a need for an urgent systematic/phylogenetic revision of the Isotomidae, including a review of the features used to delimit its subfamilies.
In contrast to the uncertainties of the internal systematics of the Isotomidae, all of our phylogenies supported it as a valid independent family, refuting the hypothesis of Yosii [8] and the findings of D’Haese [48], Schneider et al. [27,28], and Schneider and D’Haese [35]; and as the sister group of the Entomobryoidea (Figure 3), as hypothesized by Szeptycki [14] and was recovered in many recent published trees [46,47,51,53,55,56,57,63]. Our analyses also place Actaletidae apart from Isotomoidea, differently from Cassagnau [12], Massoud [13], and Soto-Adames et al. [18]. Our ML and BI trees supported Actaletidae as the sister group of Isotomidae + Entomobryoidea, suggesting the similarities between Actaletidae and Isotomidae, like the presence and shape of the postantennal organ, the absence of a trochanteral organ, body scales and head bothriotricha, plus other chaetotaxic features [18,80], may actually be due to the resemblance of both families with the ancestral of the branch Actaletidae + (Isotomidae + Entomobryoidea). In order to endorse this hypothesis, basal Entomobryoidea, like part of Nothobryinae, also share a similar postantennal organ and absence of scales with the Actaletidae and Isotomidae, while having a reduced number of metatrochanteral spines [102,103]. Previously, only the recent study of Godeiro et al. [63] (based on mitogenomes) used the Actaletidae in molecular phylogeny, reaching a topology different from our results: Actaletidae + (Tomoceridae + (Isotomidae + Entomobryoidea)). The Entomobryomorpha dataset of the previous study was limited compared to ours, and the overall high node support retrieved in our ML and BI analyses combined to the morphological evidence of a potentially closer relationship of Isotomidae and Actaletidae, suggest the clade Tomoceridae + (Actaletidae + (Isotomidae + Entomobryoidea)) is more plausible.
Our ML, MP, and BI trees found the Tomoceroidea to be the most basal branch of Entomobryomorpha (Figure 3). Such results are on par with the findings of Sun et al. [53] tree based on partitioned nucleotide sequences and Yu et al. [47], but contrast with Sun et al. [53] tree inferred from amino acid sequences, Xiong et al. [51], and Yu et al. [46]. Our dataset is limited regarding the Tomoceroidea, with only three genera of Tomocerinae (Tomoceridae) (Table S1). With the absence of Oncopoduridae and Lepidophorellinae representatives, we cannot clearly compare our results with the previously cited papers. Nevertheless, none of our trees clustered the Tomoceridae with Poduromorpha, and found the sampled Tomocerinae as an ingroup of Entomobryomorpha with overall high node support (with the exception of MP analysis, Figure S5), as endorsed by the morphology [18].

4.3. Poduromorpha

The main results obtained by the ML analyses (Figure 4) are remarkably similar to other published phylogenies based on mitogenomes [52,53,55,56,57]. These studies, as well as our data, support the Onychiuroidea apart from the clade gathering the Hypogastruridae, Poduridae, and Neanuroidea lineages. Sun et al. [53] ML tree inferred from amino acid sequences and Cucini et al. [56] BI tree excluding the third codon position recovered the Onychiuroidea as a valid superfamily, gathering Onychiuridae and Tullbergiidae. On the other hand, Sun et al. [53] ML tree based on partitioned nucleotide sequences and Leo et al. [55] BI tree based only on the first and second codon positions did not support the clade Onychiuridae + Tullbergiidae, while Carapelli et al. [52], Ma et al. [57], and our study only sampled the Onychiuridae. Most of these studies also refuted the Hypogastruridae in the same way we found in our ML, BI, and MP analyses and the majority of our tree topology tests, putting apart Ceratophysella and Gomphiocephalus, with the exception of Sun et al. [53] ML tree based on partitioned nucleotide sequences and Carapelli et al. [52], which only sampled Gomphiocephalus. Additionally, in most phylogenies, including ours, Poduridae was recovered as the sister group of the Neanuroidea, with the exception of Ma et al. [57], which did not sample Podura. Finally, Cucini et al. [56] and our ML and BI trees achieved the same internal topology for the Neanuroidea, with Neanurinae as the sister group of Pseudachorutinae, and Brachystomellidae as an ingroup of Neanuridae and the sister group of Frieseinae.
On the other hand, other phylogenetic studies based on fewer genes recovered mixed results. We found D’Haese’s [48] tree achieved a quite similar topology compared to our ML phylogenies, but a straightforward comparison would be unwise since D’Haese’s dataset was remarkably better represented than ours. Nevertheless, this study found the Onychiuroidea + Triacanthella Schäffer apart from the clade gathering most of the polyphyletic Hypogastruridae plus Poduridae and Neanuroidea. An unexpected basal position for Triacanthella was obtained by Luan et al. [49] and Xiong et al. [51]. These data combined support that the actual genus placement should be revised in future systematics studies of the Poduromorpha. Luan et al. [49] also found the topology Hypogastruridae in part + (Poduridae + Neanuroidea); however, their dataset for the Poduromorpha was limited to five species. Greenslade et al. [86] found the Onychiuroidea in part as the sister group of the branch clustering the paraphyletic Hypogastruidae, Poduridae, and Neanuroidea; however, in this study, the sister group of the latter was Gomphiocephalus, and the Poduridae appeared mixed with Hypogastruridae lineages. Contrarily, Schneider et al. [27,28] trees found Poduridae as an ingroup of Neanuroidea, and Hypogastruridae as a monophyletic family; Xiong et al. [51] recovered the topology Neanuroidea + (Tullbergiidae + (Poduridae + Hypogastruridae in part)), and Yu et al. [46] found Hypogastruridae in part + ((Triacanthella + Odontellidae) + (Poduridae + Neanuroidea)).
At least part of our findings is also supported by morphology. D’Haese’s [87] detailed morphology-based phylogeny also pointed out for a non-monophyletic Hypogastruridae, and the Poduridae closer to the Neanuroidea. The developed chewing mouthparts seen among the Hypogastruridae are undoubtedly a plesiomorphy shared by most lineages of the other orders of springtails and supposedly do not have any phylogenetical signal within the Poduromorpha. Conversely, the modified mouthparts of the Neanuroidea, with the reduction or complete loss of mandibles, notable modifications on the maxillae capitulum, and variable elongations of the oral cone, are clearly synapomorphies of derived taxa. Still, it is not clear if the Odontellidae belongs to the Neanuroidea, as some morphological traits have supported [87,104], or if its peculiar maxillae lacking the cardo point to an independent path of mouthparts modification within the Onychiuroidea, or even of an independent basal Poduromorpha branch, as supported by some molecular studies [27,28,35,48,51,86]. Our findings also endorse that Brachystomellidae should be considered a subfamily of Neanuridae, following Massoud’s [104] view. A similar placement for Brachystomellidae was also obtained by D’Haese [48] and Cucini et al. [56]. The limited dataset of Neanuroidea in our analyses, with the absence of representatives of Caputanurininae, Morulininae, and Uchidanurinae, prevents us from providing further comments on the systematics of the superfamily. We are aware that Dr. D’Haese’s team is currently working on a large phylogenetic study of the Neanuroidea, which will likely be more conclusive about the validity and internal relationships of its subfamilies.

4.4. Symphypleona and Neelipleona

Recent advances in the systematics of Collembola using molecular markers provided different views for the internal relationships of the Symphypleona s. str. families and the order’s validity. For example, D’Haese [48] and Luan et al. [49] found the Symphypleona s. str. to be the paraphyletic basal stock of all springtails, a result not endorsed by morphology or traditional systematics of Collembola [6,10,15,17,20,83,105]. On the other hand, more recent studies recovered the Symphypleona s. str. as a valid order, but with different internal topologies. Using five sampled species, Xiong et al. [51] found the ML/BI tree: Sminthurididae + (Bourletiellidae + Sminthuridae); Schneider et al. [27,28] and Schneider and D’Haese [35], with nine sample taxa, found the MP/ML tree: (Katiannidae + Arrhopalitidae) + (Dicyrtomidae in part + (Dicyrtoma Bourlet + (Bourletiellidae + Sminthuridae))); using six taxa, Yu et al. [46] found the BI/ML tree (Bourletiellidae + (Katiannidae + Sminthuridae)) + (Sminthurididae in part + (Sminthurides Börner + Dicyrtomidae)); while Sun et al. [54], based on four species, found the ML tree: Arrhopalitidae + (Sminthurididae + (Sminthuridae + Bourletiellidae)).
Concerning the use of mitogenomes to investigate Symphypleona s. str. phylogeny, the results were also divergent. Leo et al. [55] and Ma et al. [57], with four and three species, respectively, and based on first and second codon positions, reached BI and ML trees with the same topology: Dicyrtomidae + (Sminthuridae + Bourletiellidae). Using seven sampled species and based on protein-coding genes, Nardi et al. [64] reached the BI tree Sminthurididae + (Sminthuridae + (Bourletiellidae + Dicyrtomidae)). Sun et al. [53], with six species and based on partitioned nucleotide sequences, found the ML tree Sminthuridae in part + (Dicyrtomidae + ((Katiannidae + Sminthurididae) + (Sminthurus Latreille + Bourletiellidae))), while using amino acid sequences recovered (Dicyrtomidae + (Sminthurididae + Katiannidae)) + (Bourletiellidae + Sminthuridae)). The more recent study of Cucini et al. [56], based on nine species and on first and second codon positions, reached the BI tree Bourletiellidae/Dicyrtomidae + (Sminthuridae + (Katiannidae + Sminthurididae)).
All the above discrepant results are due to different types and analysis parameters (similar to our conflicting ML, BI, and MP trees), mismatched sampled taxa, and/or molecular markers. It is also worth noting that, so far, no representative molecular phylogeny of the Symphypleona s. str. has been published, with datasets reaching at most nine species from few genera, as in Schneider et al. [27,28], Cucini et al. [56], and in our study. In this scenario, all conclusions about the internal evolution of the Symphypleona s. str. based on molecular data should be taken as preliminary. Nevertheless, our ML trees (Figure 5) support previous hypotheses of the internal relationships of the order. The suborders Sminthuridida and Appendiciphora Bretfeld [17] sensu Sánchez-García and Engel [82] were recovered with high node support in all ML trees, while our tree topology tests recovered mixed results, both supporting Appendiciphora and Sminthuridida, or refuting them as independent taxa (Table S3). The morphology strongly supports the two suborders, as the Sminthurididae and Mackenziellidae (Sminthuridida) share the apomorphic antennal clasper of the males and the plesiomorphic short and spherical ventral tube sacs and absence of the subanal appendages; at the same time, the Appendiciphora, gathering all the other families of Symphypleona s. str., share the derived long ventral tube sacs and the subanal appendage of the females [10,17,21,83]. The absence of any Katiannoidea taxa (Katiannidae, Spinothecidae, Arrhopalitidae, and Collophoridae) in the analyses prevented us from testing Börner’s [6,105] original view of a closer relationship between Sminthurididae and Katiannoidea, which was endorsed by the results of Sun et al. [53] and Cucini et al. [56]. Additionally, Mackenziella psocodes Hammer, the sole species of Mackenziellidae, has never been included in any molecular phylogeny so far. Such an omission prevents a clearer view of its relationships with other Collembola lineages, as its morphology is unique, and its systematic position within different orders was disputed by several authors [8,16,21,106,107].
Our ML trees also supported the Sminthuroidea sensu Bretfeld [108], a superfamily that gathers the Sminthuridae and Bourletiellidae. Such lineages are remarkably similar, with some taxa presenting overlapping morphologies between them [109]. Contrarily, our ML and BI trees did not support Sminthurinae and Sphyrothecinae, as Lipothrix Börner, a representative of the latter, was found inside the former (Figure 1 and Figures S1–S4). At this time, with only three species and genera of Sminthuridae sampled, such results are insufficient to propose the fusion or regrouping of these subfamilies. Additionally, our analyses did not sample any Songhaicinae taxon, the third subfamily of Sminthuridae [110].
We could not compare our results concerning the Neelipleona with other studies investigating molecular phylogenies of the group. In most of them, the dataset was restricted to a single species [46,50,51,53], or the presence of a species without an assigned genus obscured the comparison between the topologies [55,56]. The detailed trees of Schneider et al. [27,28] and Schneider and D’Haese [35] only sampled two genera, Neelus and Megalothorax, and focused mainly on the latter internal lineages. The same sampled genera limitation was observed in the phylogeny of Ma et al. [57]. Although our data are also limited by sample size, the use of the three main genera of Neelipleona in the analyses provided us with an important insight on the internal evolution of the order. The sensory fields are unique structures without morphological parallel in other Collembola taxa, likely representing an autapomorphy of the Neelipleona [89]. Some authors consider Neelides as the sole genus of Neelidae lacking sensory fields [83,93], while others provided evidence that they are present but are rudimentary and discrete [89,91]. In all other genera of the family, including Neelus and Megalothorax, such structures are very well delimited and hold specialized inner chaetae [83,91,93]. In our analysis, Neelides was recovered as the most basal branch of the order, with Megalothorax + Neelus as its sister group (Figure 6). Based on this finding, we believe the sensory fields present in Neelides (if homologous to the ones seen in other genera [91]) represent a more basal state of this character in Neelipleona internal evolution, while the well-delimited sensory fields are a derived state, a synapomorphic trait of the higher taxa. This hypothesis, supported by our data, allows us to suggest the subdivision of Neelidae into two subfamilies: Neelinae Handlirsch, 1929 [83], including the sampled Neelus and Megalothorax, the largest genera of the order [1], plus Zelandothorax Delamare Deboutteville & Massoud, Spinaethorax Papáč and Palacios-Vargas and Acanthoneelidus Bretfeld & Griegel; and Neelidinae subfam. nov. Bellini, Godeiro and Zhang, represented by Neelides species. This proposition is included in Figure 7.

5. Conclusions

Our results advocate that mitogenomes are highly suitable for the study of Collembola phylogeny, supporting many previously traditional and more recent hypotheses regarding the relationships between its main internal lineages. For instance, our main tree supported the validity of the four current orders, the superorder Arthropleona, the Entomobryoidea as the sister group of Isotomidae, the Neanuroidea (putting Brachystomellidae inside Neanuridae), the split of Symphypleona s. str. into the Appendiciphora and Sminthuridida suborders, the Sminthuroidea clustering Sminthuridae and Bourletiellidae, and Neelides as a basal branch of Neelipleona, allowing us to divide the Neelidae into two subfamilies: Neelinae and Neelidinae subfam. nov. Contrarily, our main tree refuted the superorder Symphypleona s. lat., Paronellidae sensu Soto-Adames et al. [18], all subfamilies of Isotomidae, Hypogastruridae, and Poduridae as a basal lineage of Poduromorpha. Though this study presents one of the most taxonomically and genetically comprehensive springtail phylogenies to date, it also highlights a great need for the inclusion of additional representative taxa and nuclear markers to better resolve the evolutionary relationships among some persistently problematic springtail lineages.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d15010007/s1. Figure S1: Phylogeny of Collembola inferred from Maximum Likelihood, partitioned dataset, mixture model (ML_2). Node labels show SH-aLRT support and bootstrap values, respectively; Figure S2: Phylogeny of Collembola inferred from Maximum Likelihood, unpartitioned dataset, model mtART + C60 + FO + R (ML_3). Node labels show SH-aLRT support and bootstrap values, respectively; Figure S3: Phylogeny of Collembola inferred from Maximum Likelihood, partitioned dataset detailed in Table S2 (ML_4). Node labels show SH-aLRT support and bootstrap values, respectively; Figure S4: Phylogeny of Collembola inferred from Bayesian Inference (BI), default model CAT+GTR with four categories. Node labels show posterior probability support; Figure S5: Phylogeny of Collembola inferred from Maximum Parsimony (MP), ‘*’ marks Actaletidae outside the Entomobryomorpha. Node labels show bootstrap support; Figure S6: Phylogeny of Collembola inferred from Maximum Likelihood, unpartitioned dataset, substitution model mtZOA + F + R10 (ML_1). Node labels show raw SH-aLRT support and bootstrap values, respectively. Branches of each Collembola order are represented in different colors. Table S1: Taxonomic information, GenBank accession numbers, country and information source of all analyzed species; Table S2: Partitioning scheme and substitution models selected by Model Finder used for partitioned maximum likelihood analyses. Table S3: Tree topology tests of non monophyletic groups of Collembola higher taxa.

Author Contributions

Conceptualization, B.C.B. and N.N.G.; methodology, N.N.G.; software, B.C.B. and N.N.G.; validation, B.C.B., N.N.G. and F.Z.; formal analysis, B.C.B. and N.N.G.; investigation, B.C.B. and N.N.G.; resources, B.C.B., N.N.G. and F.Z.; data curation, N.N.G., B.C.B., P.G.C.d.S., R.C.d.S.-C. and G.d.S.M.; writing—original draft preparation, B.C.B., N.N.G. and F.Z.; writing—review and editing, B.C.B., N.N.G., P.G.C.d.S., R.C.d.S.-C., G.d.S.M. and F.Z.; visualization, B.C.B. and N.N.G.; supervision, B.C.B. and N.N.G.; project administration, B.C.B.; funding acquisition, B.C.B., N.N.G. and F.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Council for Scientific and Technological Development (CNPq), grant number 309114/2021-7 (B.C.B. research project); the Coordination for the Improvement of Higher Education Personnel (CAPES), grant numbers 156828/2016-01 (R.C.d.S.-C. Ph.D. scholarship) and 001 (G.d.S.M. Ph.D. scholarship); Fundação de Apoio à Pesquisa do Estado da Paraíba (FAPESQ/PB), grant number 07/2021-SEECT (P.G.C.d.S. Ph.D. scholarship); Research Foundation of Shanghai Science and Technology Museum and Postdoctoral fund of Haibo Program of Pudong New Area in 2021 (N.N.G. postdoctoral grants); and the National Natural Science Foundation of China, grant numbers 31970434 and 32270470 (F.Z. research projects).

Institutional Review Board Statement

Ethical review and approval were waived for this study since we used the public NCBI database to gather the studied molecular samples.

Data Availability Statement

All data is contained within the article and the Supplementary Materials.

Acknowledgments

We would like to thank the anonymous reviewers for carefully revising the manuscript and providing ideas to improve it.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bellinger, P.F.; Christiansen, K.A.; Janssens, F. Checklist of the Collembola of the World. Available online: http://www.collembola.org (accessed on 16 November 2022).
  2. Potapov, A.; Bellini, B.C.; Chown, S.L.; Deharveng, L.; Janssens, F.; Kovác, L.; Kuznetsova, N.; Ponge, J.-F.; Potapov, M.; Querner, P.; et al. Towards a global synthesis of Collembola knowledge—Challenges and potential solutions. Soil Org. 2020, 92, 161–188. [Google Scholar] [CrossRef]
  3. Zhang, Z.-Q. Animal biodiversity: An outline of higher-level classification and survey of taxonomic richness. Zootaxa 2011, 3148, 3–6. [Google Scholar] [CrossRef]
  4. Cicconardi, F.; Fanciulli, P.P.; Emerson, B.C. Collembola, the biological species concept and the underestimation of global species richness. Mol. Ecol. 2013, 22, 5382–5396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Börner, C. Neue Collembolenformen und zur Nomenclatur der Collembola Lubbock. Zool. Anz. 1901, 24, 696–712. [Google Scholar]
  6. Börner, C. Das System der Collembolen nebst Beschreibung neuer Collembolen des Hamburger Naturhistorischen Museums. Mitt. Nat. Mus. Hambg. 1906, 23, 147–188. [Google Scholar]
  7. Paclt, J. Biologie der primär flügellosen Insekten. Jena 1956, 6, 1–258. [Google Scholar]
  8. Yosii, R. Phylogenetische bedeutung der chaetotaxie bei den Collembolen. Contrib. Biol. Lab. Kyoto Univ. 1961, 12, 1–37. [Google Scholar]
  9. Salmon, J.T. An Index to the Collembola. R. Soc. N. Z. Bull. 1964, 7, 1–144. [Google Scholar]
  10. Richards, W.R. Generic classification, evolution, and biogeography of the Sminthuridae of the world (Collembola). Mem. Entomol. Soc. Can. 1968, 53, 3–54. [Google Scholar] [CrossRef]
  11. Uchida, H. Tentative key to the Japanese genera of Collembola, in relation to the world genera of this order I. Sci. Rep. Hirosaki Univ. 1971, 18, 63–76. [Google Scholar]
  12. Cassagnau, P. La phylogénie des Collemboles à la lumière des structures endocrines rétrocérébrales. I Symp. Inter. Zoolfilogenia 1971, 1, 333–349. [Google Scholar]
  13. Massoud, Z. Essai de synthèse sur la phylogénie des Collemboles. Rev. Ecol. Biol. Sol. 1976, 13, 241–252. [Google Scholar]
  14. Szeptycki, A. Chaetotaxy of the Entomobryidae and Its Phylogenetical Significance. Morpho-Systematic Studies on Collembola; Polska Akademia Nauk: Kraków, Poland, 1979; Volume IV, pp. 1–219. [Google Scholar]
  15. Betsch, J.M. Éléments pour une monographie des Collemboles Symplyplêones (Hexapodes, Aptérygotes). Mém. Mus. Natl. Hist. Nat. Sér. A Zool. 1980, 116, 1–227. [Google Scholar]
  16. Moen, P.; Ellis, W.N. Morphology and Taxonomic Position of Podura aquatica (Collembola). Entomol. Gen. 1984, 9, 193–204. [Google Scholar] [CrossRef]
  17. Bretfeld, G. Phylogenetic systematics of the higher taxa of Symphypleona Börner, 1901 (Insecta, Entognatha, Collembola). Proc. 2nd Intern. Sem. Apterygota 1986, 1, 302–311. [Google Scholar]
  18. Soto-Adames, F.N.; Barra, J.A.; Christiansen, K.; Jordana, R. Suprageneric classification of Collembola Entomobryomorpha. Ann. Entomol. Soc. Am. 2008, 101, 501–513. [Google Scholar] [CrossRef]
  19. Cassagnau, P. Des hexapodes vieux de 400 millions d’années: Les collemboles. 1.—Biologie et Évolution. Ann. Biol. 1990, 29, 1–69. [Google Scholar]
  20. Deharveng, L. Recent advances in Collembola systematics. Pedobiologia 2004, 48, 415–433. [Google Scholar] [CrossRef]
  21. Fjellberg, A. Redescription of Mackenziella psocoides Hammer, 1953 and discussion of its systematic position (Collembola, Mackenziellidae). Proc. 2nd Intern. Sem. Apterygota 1989, 1, 93–105. [Google Scholar]
  22. Fanciulli, P.P.; Melegari, D.; Carapelli, A.; Frati, F.; Dallai, R. Population structure, gene flow and evolutionary relationships in four species of the genera Tomocerus and Pogonognathellus (Collembola, Tomoceridae). Biol. J. Linn. Soc. 2000, 70, 221–238. [Google Scholar] [CrossRef] [Green Version]
  23. Fanciulli, P.P.; Summa, D.; Dallai, R.; Frati, F. High levels of genetic variability and population differentiation in Gressittacantha terranova (Collembola, Hexapoda) from Victoria Land, Antarctica. Antarct. Sci. 2001, 13, 246–254. [Google Scholar] [CrossRef]
  24. Timmermans, M.J.T.N.; Ellers, J.; Mariën, J.; Verhoef, S.C.; Ferwerda, E.B.; Van Straalen, N.M. Genetic structure in Orchesella cincta (Collembola): Strong subdivision of European populations inferred from mtDNA and AFLP markers. Mol. Ecol. 2005, 14, 2017–2024. [Google Scholar] [CrossRef] [PubMed]
  25. Stevens, M.I.; Frati, F.; McGaughran, A.; Spinsanti, G.; Hogg, I.D. Phylogeographic structure suggests multiple glacial refugia in northern Victoria Land for the endemic Antarctic springtail Desoria klovstadi (Collembola, Isotomidae). Zool. Scr. 2007, 36, 201–212. [Google Scholar] [CrossRef]
  26. Torricelli, G.; Frati, F.; Convey, P.; Telford, M.; Carapelli, A. Population structure of Friesea grisea (Collembola, Neanuridae) in the Antarctic Peninsula and Victoria Land: Evidence for local genetic differentiation of pre-Pleistocene origin. Antarct. Sci. 2010, 22, 757–765. [Google Scholar] [CrossRef] [Green Version]
  27. Schneider, C.; Cruaud, C.; D’Haese, C.A. Unexpected diversity in Neelipleona revealed by molecular phylogeny approach (Hexapoda, Collembola). Soil Org. 2011, 83, 383–398. [Google Scholar]
  28. Schneider, C.; Zon, S.D.; D’Haese, C.A. Megalothorax laevis (Neelipleona, Neelidae): Account of a neglected springtail widely distributed in the intertropical zone. Int. J. Trop. Insect Sci. 2018, 38, 168–191. [Google Scholar] [CrossRef]
  29. Katz, A.D.; Giordano, R.; Soto-Adames, F.N. Operational criteria for cryptic species delimitation when evidence is limited, as exemplified by North American Entomobrya (Collembola: Entomobryidae). Zool. J. Linn. Soc. 2015, 173, 818–840. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, F.; Jantarit, S.; Nilsai, A.; Stevens, M.I.; Ding, Y.; Satasook, C. Species delimitation in the morphologically conserved Coecobrya (Collembola: Entomobryidae): A case study integrating morphology and molecular traits to advance current taxonomy. Zool. Scr. 2018, 47, 342–356. [Google Scholar] [CrossRef]
  31. Yu, D.; Qin, C.; Ding, Y.; Hu, F.; Zhang, F.; Liu, M. Revealing species diversity of Tomocerus ocreatus complex (Collembola: Tomoceridae): Integrative species delimitation and evaluation of taxonomic characters. Arthropod Syst. Phylogeny 2018, 76, 147–172. [Google Scholar]
  32. Mateos, E.; Winkler, D.; Riutort, M.; Álvarez-Presas, M. New morphological and molecular data reveal an important underestimation of species diversity and indicate evolutionary patterns in European Lepidocyrtus (Collembola: Entomobryidae). Invertebr. Syst. 2021, 35, 471–492. [Google Scholar] [CrossRef]
  33. Dukes, C.D.; Janssens, F.; Recuero, E.; Caterino, M.S. Specific and Intraspecific Diversity of Symphypleona and Neelipleona (Hexapoda: Collembola) in Southern High Appalachia (USA). Diversity 2022, 14, 847. [Google Scholar] [CrossRef]
  34. Soto-Adames, F.N. Molecular phylogeny of the Puerto Rican Lepidocyrtus and Pseudosinella (Hexapoda: Collembola), a validation of Yoshii’s “color pattern species”. Mol. Phylogenet. Evol. 2002, 25, 27–42. [Google Scholar] [CrossRef] [PubMed]
  35. Schneider, C.; D’Haese, C.A. Morphological and molecular insights on Megalothorax: The largest Neelipleona genus revisited (Collembola). Invertebr. Syst. 2013, 27, 317–364. [Google Scholar] [CrossRef]
  36. Mateos, E.; Escuer, P.; Busmachiu, G.; Riutort, M.; Álvarez-Presas, M. Untangling Lepidocyrtus (Collembola, Entomobryidae): New molecular data shed light on the relationships of the European groups. Invertebr. Syst. 2018, 32, 639–651. [Google Scholar] [CrossRef]
  37. Winkler, D.; Mateos, E.; Traser, G.; Lakatos, F.; Tóth, V. New Insight into the Systematics of European Lepidocyrtus (Collembola: Entomobryidae) Using Molecular and Morphological Data. Insects 2020, 11, 302. [Google Scholar] [CrossRef]
  38. Zhang, F.; Pan, Z.; Wu, J.; Ding, Y.; Yu, D.; Wang, B. Dental scales could occur in all scaled subfamilies of Entomobryidae (Collembola): New definition of Entomobryinae with description of a new genus and three new species. Invertebr. Syst. 2016, 30, 598–615. [Google Scholar] [CrossRef]
  39. Zhang, F.; Ma, Y.; Greenslade, P. New Australian Paronellidae (Collembola) reveal anomalies in existing tribal diagnoses. Invertebr. Syst. 2017, 31, 375–393. [Google Scholar] [CrossRef]
  40. Nunes, R.C.; Godeiro, N.N.; Pacheco, G.; Liu, S.; Gilbert, M.T.P.; Alvarez-Valin, F.; Zhang, F.; Bellini, B.C. The discovery of Neotropical Lepidosira (Collembola, Entomobryidae) and its systematic position. Zool. Scr. 2019, 48, 783–800. [Google Scholar] [CrossRef]
  41. Godeiro, N.N.; Pacheco, G.; Liu, S.; Cipola, N.G.; Berbel-Filho, W.M.; Zhang, F.; Gilbert, M.T.; Bellini, B.C. Phylogeny of Neotropical Seirinae (Collembola, Entomobryidae) based on mitochondrial genomes. Zool. Scr. 2020, 49, 329–339. [Google Scholar] [CrossRef]
  42. Zhang, F.; Chen, Z.; Dong, R.R.; Deharveng, L.; Stevens, M.I.; Huang, Y.H.; Zhu, C.D. Molecular phylogeny reveals independent origins of body scales in Entomobryidae (Hexapoda: Collembola). Mol. Phylogenet. Evol. 2014, 70, 231–239. [Google Scholar] [CrossRef]
  43. Zhang, F.; Sun, D.; Yu, D.; Wang, B. Molecular phylogeny supports S-chaetae as a key character better than jumping organs and body scales in classification of Entomobryoidea (Collembola). Sci. Rep. 2015, 5, 12471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Zhang, F.; Bellini, B.C.; Soto-Adames, F.N. New insights into the systematics of Entomobryoidea (Collembola: Entomobryomorpha): First instar chaetotaxy, homology and classification. Zool. Syst. 2019, 44, 249–278. [Google Scholar] [CrossRef]
  45. Zhang, F.; Deharveng, L. Systematic revision of Entomobryidae (Collembola) by integrating molecular and new morphological evidence. Zool. Scr. 2015, 44, 298–311. [Google Scholar] [CrossRef]
  46. Yu, D.; Zhang, F.; Stevens, M.I.; Yan, Q.; Liu, M.; Hu, F. New insight into the systematics of Tomoceridae (Hexapoda, Collembola) by integrating molecular and morphological evidence. Zool. Scr. 2016, 45, 286–299. [Google Scholar] [CrossRef]
  47. Yu, D.; Ding, Y.; Tihelka, E.; Cai, C.; Hu, F.; Liu, M.; Zhang, F. Phylogenomics of Elongate-Bodied Springtails Reveals Independent Transitions from Aboveground to Belowground Habitats in Deep Time. Syst. Biol. 2022, 71, 1023–1031. [Google Scholar] [CrossRef] [PubMed]
  48. D’Haese, C.A. Were the first springtails semi-aquatic? A phylogenetic approach by means of 28S rDNA and optimization alignment. Proc. R. Soc. Lond. 2002, 269, 1143–1151. [Google Scholar] [CrossRef] [Green Version]
  49. Luan, Y.-X.; Mallatt, J.M.; Xie, R.D.; Yang, Y.M.; Yin, W.Y. The phylogenetic positions of three basal-hexapod groups (Protura, Diplura, and Collembola) based on ribosomal RNA gene sequences. Mol. Biol. Evol. 2005, 22, 1579–1592. [Google Scholar] [CrossRef] [Green Version]
  50. Gao, Y.; Bu, Y.; Luan, Y.-X. Phylogenetic Relationships of Basal Hexapods Reconstructed from Nearly Complete 18S and 28S rRNA Gene Sequences. Zool. Sci. 2008, 25, 1139–1145. [Google Scholar] [CrossRef]
  51. Xiong, Y.; Gao, Y.; Yin, W.Y.; Luan, Y.X. Molecular phylogeny of Collembola inferred from ribosomal RNA genes. Mol. Phylogenet. Evol. 2008, 49, 728–735. [Google Scholar] [CrossRef]
  52. Carapelli, A.; Convey, P.; Nardi, F.; Frati, F. The mitochondrial genome of the antarctic springtail Folsomotoma octooculata (Hexapoda; Collembola), and an update on the phylogeny of collembolan lineages based on mitogenomic data. Entomologia 2014, 2, 46–55. [Google Scholar] [CrossRef]
  53. Sun, X.; Yu, D.; Xie, Z.; Dong, J.; Ding, Y.; Yao, H.; Greenslade, P. Phylomitogenomic analyses on collembolan higher taxa with enhanced taxon sampling and discussion on method selection. PLoS ONE 2020, 15, e0230827. [Google Scholar] [CrossRef] [PubMed]
  54. Sun, X.; Ding, Y.; Orr, M.C.; Zhang, F. Streamlining universal single-copy orthologue and ultraconserved element design: A case study in Collembola. Mol. Ecol. Resour. 2020, 20, 706–717. [Google Scholar] [CrossRef] [PubMed]
  55. Leo, C.; Carapelli, A.; Cicconardi, F.; Frati, F.; Nardi, F. Mitochondrial genome diversity in Collembola: Phylogeny, dating and gene order. Diversity 2019, 11, 169. [Google Scholar] [CrossRef] [Green Version]
  56. Cucini, C.; Fanciulli, P.P.; Frati, F.; Convey, P.; Nardi, F.; Carapelli, A. Re-evaluating the internal phylogenetic relationships of Collembola by means of Mitogenome Data. Genes 2021, 12, 44. [Google Scholar] [CrossRef]
  57. Ma, Y.; Huang, C.-W.; Luan, Y.-X.; Chen, W.J. The mitochondrial genome of a minute springtail species Megalothorax incertus (Collembola: Neelipleona: Neelidae). Mitochondrial DNA B Resour. 2021, 6, 2430–2432. [Google Scholar] [CrossRef]
  58. Song, F.; Li, H.; Jiang, P.; Zhou, X.G.; Liu, J.P.; Sun, C.H.; Vogler, A.P.; Cai, W.Z. Capturing the phylogeny of Holometabola with mitochondrial genome data and Bayesian site-heterogeneous mixture models. Genome Biol. Evol. 2016, 8, 1411–1426. [Google Scholar] [CrossRef]
  59. Wang, Y.; Cao, J.J.; Li, N.; Ma, G.Y.; Li, W.H. The first mitochondrial genome from Scopuridae (Insecta: Plecoptera) reveals structural features and phylogenetic implications. Int. J. Biol. Macromol. 2019, 122, 893–902. [Google Scholar] [CrossRef]
  60. Godeiro, N.N.; Bellini, B.C.; Ding, N.; Xu, C.; Ding, Y.; Zhang, F. A mitogenomic phylogeny of the Entomobryoidea (Collembola): A comparative perspective. Zool. Scr. 2021, 50, 658–666. [Google Scholar] [CrossRef]
  61. Godeiro, N.N.; Zhang, F.; Cipola, N.G. First partial mitogenome of a new Seira Lubbock species (Collembola, Entomobryidae, Seirinae) from Cambodia reveals a possible separate lineage from the Neotropical Seirinae. Zootaxa 2020, 4890, 451–472. [Google Scholar] [CrossRef]
  62. Godeiro, N.N.; Zhang, F. First record of Seira dowlingi (Wray, 1953) (Collembola, Entomobryidae, Seirinae) from China and mitogenome comparison with the New World specimens. Zootaxa 2021, 5020, 191–196. [Google Scholar] [CrossRef]
  63. Godeiro, N.N.; Palacios-Vargas, J.G.; Gao, Y.; Bu, Y. Complete mitochondrial genome of the Mexican marine littoral hygrophilous Spinactaletes boneti (Collembola: Actaletidae) and its phylogenetic placement. Mitochondrial DNA B Resour. 2022, 7, 755–757. [Google Scholar] [CrossRef] [PubMed]
  64. Nardi, F.; Cucini, C.; Leo, C.; Frati, F.; Fanciulli, P.P.; Carapelli, A. The complete mitochondrial genome of the springtail Allacma fusca, the internal phylogenetic relationships and gene order of Symphypleona. Mitochondrial DNA B Resour. 2020, 5, 3103–3105. [Google Scholar] [CrossRef] [PubMed]
  65. Beutel, R.G.; Yavorskaya, M.I.; Mashimo, Y.; Fukui, M.; Meusemann, K. The Phylogeny of Hexapoda (Arthropoda) and the Evolution of Megadiversity. Proc. Arthropod. Embryol. Soc. Jpn. 2017, 51, 1–15. [Google Scholar]
  66. Sasaki, G.; Ishiwata, K.; Machida, R.; Miyata, T.; Su, Z.H. Molecular phylogenetic analyses support the monophyly of Hexapoda and suggest the paraphyly of Entognatha. BMC Evol. Biol. 2013, 13, 236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Misof, B.; Liu, S.; Meusemann, K.; Peters, R.S.; Donath, A.; Mayer, C.; Frandsen, P.B.; Ware, J.; Flouri, T.; Beutel, R.G.; et al. Phylogenomics resolves the timing and pattern of insect evolution. Science 2014, 346, 763–767. [Google Scholar] [CrossRef]
  68. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 2013, 8, 1494–1512. [Google Scholar] [CrossRef]
  69. Smirnov, V.; Warnow, T. MAGUS: Multiple sequence alignment using graph clustering. Bioinformatics 2021, 37, 1666–1672. [Google Scholar] [CrossRef]
  70. 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]
  71. Criscuolo, A.; Gribaldo, S. BMGE (Block Mapping and Gathering with Entropy): A new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol. Biol. 2010, 10, 2–21. [Google Scholar] [CrossRef] [Green Version]
  72. Steenwyk, J.L.; Buida, T.J.; Labella, A.L.; Li, Y.; Shen, X.X.; Rokas, A. PhyKIT: A broadly applicable UNIX shell toolkit for processing and analyzing phylogenomic data. Bioinformatics 2021, 37, 2325–2331. [Google Scholar] [CrossRef]
  73. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Hoang, D.T.; Chernomor, O.; von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the Ultrafast Bootstrap Approximation. Mol. Biol. Evol. 2018, 35, 518–522. [Google Scholar] [CrossRef] [PubMed]
  75. Kalyaanamoorthy, S.; Minh, B.; Wong, T.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Wang, H.; Susko, E.; Minh, B.; Roger, A.J. Modeling site heterogeneity with posterior mean site frequency profiles accelerates accurate phylogenomic estimation. Syst. Biol. 2018, 67, 216–235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Lartillot, N.; Rodrigue, N.; Stubbs, D.; Richer, J. PhyloBayes MPI: Phylogenetic reconstruction with infinite mixtures of profiles in a parallel environment. Syst. Biol. 2013, 62, 611–615. [Google Scholar] [CrossRef]
  78. Hoang, D.T.; Vinh, L.S.; Flouri, T.; Stamatakis, A.; von Haeseler, A.; Minh, B.Q. MPBoot: Fast phylogenetic maximum parsimony tree inference and bootstrap approximation. BMC Evol. Biol. 2018, 18, 11. [Google Scholar] [CrossRef] [Green Version]
  79. Rambaut, A.; Institute of Evolutionary Biology, University of Edinburgh, Edinburgh. FigTree V1.3.1. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 16 November 2022).
  80. Soto-Adames, F.N. Revision de la familia Actaletidae Borner, 1902 (Insects: Collembola). Caribb. J. Sci. 1988, 24, 161–196. [Google Scholar]
  81. Bellini, B.C. Colêmbolos: Uma riqueza microscópica no Semiárido. In Conhecendo os artrópodes do Semiárido, 1st ed.; Bravo, F., Calor, A.R., Eds.; Métis Produção Editorial: São Paulo, Brazil, 2016; Volume 1, pp. 43–55. [Google Scholar]
  82. Sánchez-García, A.; Engel, M.S. Long-term stasis in a diverse fauna of Early Cretaceous springtails (Collembola: Symphypleona). J. Syst. Palaeontol. 2016, 15, 513–537. [Google Scholar] [CrossRef]
  83. Bretfeld, G. Synopses on Palaeartic Collembola, Volume 2. Symphypleona. Abh. Ber. Naturkundemus. Gorlitz 1999, 71, 1–318. [Google Scholar]
  84. D’Haese, C.A. Homology and morphology in Poduromorpha (Hexapoda, Collembola). Eur. J. Entomol. 2003, 101, 385–407. [Google Scholar] [CrossRef]
  85. Hopkin, S.P. Biology of the Springtails (Insecta: Collembola), 1st ed.; Oxford University Press: New York, NY, USA, 1997; pp. 1–330. [Google Scholar]
  86. Greenslade, P.; Stevens, M.I.; Torricelli, G.; D’Haese, C.A. An ancient Antarctic endemic genus restored: Morphological and molecular support for Gomphiocephalus hodgsoni (Collembola: Hypogastruridae). Syst. Entomol. 2011, 36, 223–240. [Google Scholar] [CrossRef]
  87. D’Haese, C.A. Morphological appraisal of Collembola phylogeny with special emphasis on Poduromorpha and a test of the aquatic origin hypothesis. Zool. Scr. 2003, 32, 563–586. [Google Scholar] [CrossRef]
  88. Zhang, F.; Deharveng, L. First instar tibiotarsal chaetotaxy supports the Entomobryidae and Symphypleona (Collembola) forming a cluster in a phylogenetic tree. Zootaxa 2015, 3955, 487–504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Bonet, F. Monografia de la familia Neelidae (Collembola). Rev. Soc. Mex. Hist. Nat. 1947, 8, 131–192. [Google Scholar]
  90. Betsch, J.M.; Waller, A. Chaetotaxic nomenclature of the head, thorax and abdomen in Symphypleona (Insecta, Collembola). Acta Zool. Fenn. 1994, 195, 5–12. [Google Scholar]
  91. Schneider, C. Morphological review of the order Neelipleona (Collembola) through the redescription of the type species of Acanthoneelidus, Neelides and Neelus. Zootaxa 2017, 4308, 1–94. [Google Scholar] [CrossRef]
  92. Betsch, J.M. An ontogenetically focused chaetotaxial scheme in Symphypleona (collembolan): The 6th abdominal segment. Pedobiologia 1997, 41, 13–18. [Google Scholar]
  93. Papáč, V.; Palacios-Vargas, J.G. A new genus of Neelidae (Collembola) from Mexican caves. ZooKeys 2016, 569, 37–51. [Google Scholar] [CrossRef] [Green Version]
  94. Yosii, R. Studies on the Collembolan Genus Hypogastrura II. Neartic Forms collected by Prof. F. Bonet. Contrib. Biol. Lab. Kyoto Univ. 1962, 13, 1–25. [Google Scholar]
  95. Nayrolles, P. Chetotaxie tibiotarsale des collemboles symphyplèones. Trav. Lab. Ecobiol. Arthropodes Edaphiques 1988, 5, 1–19. [Google Scholar]
  96. Bretfeld, G.; Griegel, A. Description of a new Neelidae genus and species and of new specimens of Sminthurides annulicornis Axelson 1905 from Poland. Senckenberg. Biol. 1999, 79, 211–223. [Google Scholar]
  97. Havird, J.C.; Sloan, D.B. The roles of mutation, selection, and expression in determining relative rates of evolution in mitochondrial versus nuclear genomes. Mol. Biol. Evol. 2016, 33, 3042–3053. [Google Scholar] [CrossRef] [PubMed]
  98. Zhang, F.; Cipola, N.G.; Pan, Z.-X.; Ding, Y. New insight into the systematics of Heteromurini (Collembola: Entomobryidae: Heteromurinae) with special reference to Alloscopus and Sinodicranocentrus gen.n. Arth. Syst. Phylogeny 2020, 78, 1–16. [Google Scholar] [CrossRef]
  99. Ding, Y.-H.; Yu, D.-Y.; Guo, W.B.; Li, J.-N.; Zhang, F. Molecular phylogeny of Entomobrya (Collembola: Entomobryidae) from China: Color pattern groups and multiple origins. Insect Sci. 2018, 26, 587–597. [Google Scholar] [CrossRef] [PubMed]
  100. Soto-Adames, F.N. Postembryonic development of the dorsal chaetotaxy in Seira dowlingi (Collembola, Entomobryidae); with an analysis of the diagnostic and phylogenetic significance of primary chaetotaxy in Seira. Zootaxa 2008, 1683, 1–31. [Google Scholar] [CrossRef]
  101. Potapov, M. Isotomidae. In Synopses on Palearctic Collembola, 1st ed.; Dunger, W., Ed.; Abdhandlungen und Berichte des Naturkundemuseums Görlitz: Görlitz, Germany, 2001; Volume 3, pp. 1–603. [Google Scholar]
  102. Nunes, R.C.; Bellini, B.C. A new species of Nothobrya Arlé, 1961 (Collembola: Entomobryidae) from Brazil and notes on key characters for Nothobryinae taxonomy, with an identification key to the species of the subfamily. Zootaxa 2019, 4615, 375–391. [Google Scholar] [CrossRef] [PubMed]
  103. Nunes, R.C.; Santos-Costas, R.C.; Bellini, B.C. The first Neotropical Capbrya Barra, 1999 (Collembola: Orchesellidae: Nothobryinae) and the reinterpretation of Nothobryinae systematics. Zool. Anz. 2020, 288, 24e42. [Google Scholar] [CrossRef]
  104. Massoud, Z. Monographie des Neanuridae, Collemboles Poduromorphes apiéces buccales modifiées. In Biologie de l’Amerique Australe, 1st ed.; Delamare Deboutteville, C., Rapoport, E.H., Eds.; Éditions du CNRS: Paris, France, 1967; Volume 3, pp. 7–399. [Google Scholar]
  105. Börner, C. Die Familien der Collembolen. Zool. Anz. 1913, 61, 315–322. [Google Scholar]
  106. Hammer, M. Investigations on the Microfauna of Northern Canada, Part II, Collembola. Acta Arct. 1953, 6, 1–108. [Google Scholar]
  107. Huther, W. Die systematische Stellung von Mackenziella psocoides Hammer (Collembola). Zool. Anz. 1964, 173, 119–126. [Google Scholar]
  108. Bretfeld, G. Sturmius epiphytus n. gen. n. spec. from Colombia, a taxon of the Symphypleona (Insecta, Colembola) with an unexpected character combination. Description and position in non-Linnean and Linnean classifications of the Symphypleona. Zool. Syst. Evol. Res. 1994, 32, 264–281. [Google Scholar] [CrossRef]
  109. Souza, P.G.C.; Medeiros, G.S.; Ferreira, R.L.; Souza-Silva, M.; Bellini, B.C. A Highly Troglomorphic New Genus of Sminthuridae (Collembola, Symphypleona) from the Brazilian Semiarid Region. Insects 2022, 13, 650. [Google Scholar] [CrossRef] [PubMed]
  110. Bernard, C.E.; Wynne, J.J. Disparrhopalites naasaveqw n. sp. from caves at Wupatki National Monument, Arizona, synonymy of Dietersminthurus Palacios-Vargas, Cuéllar & Vázquez, 1998 with Disparrhopalites Stach, 1956 and composition of Songhaicinae (Collembola: Sminthuridae). Zootaxa 2017, 4319, 77–90. [Google Scholar] [CrossRef]
Figure 1. Phylogeny of Collembola inferred from Maximum Likelihood unpartitioned dataset substitution model mtZOA + F + R10 (ML_1). Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively; colors represent the following scores: black = 100; yellow = 95–99.9; grey = 90–94.9; white < 90. Branches of each Collembola order are represented in different colors.
Figure 1. Phylogeny of Collembola inferred from Maximum Likelihood unpartitioned dataset substitution model mtZOA + F + R10 (ML_1). Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively; colors represent the following scores: black = 100; yellow = 95–99.9; grey = 90–94.9; white < 90. Branches of each Collembola order are represented in different colors.
Diversity 15 00007 g001
Figure 2. Summarized phylogenetic relationships of the Collembola orders obtained in this study: (A) topology of ML_1, ML_2 and ML_4 trees (nodes depict only ML_1 support levels); (B) topology of ML_3 tree; (C) topology of BI tree; (D) topology of MP tree. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML trees, posterior probability for BI tree, and maximum parsimony bootstrap for MP tree. Colors represent the following scores: black = 100 (1 for BI); yellow = 95–99.9 (0.95–0.99 for BI); grey = 90–94.9 (0.9–0.94 for BI); white < 90 (<0.9 for BI).
Figure 2. Summarized phylogenetic relationships of the Collembola orders obtained in this study: (A) topology of ML_1, ML_2 and ML_4 trees (nodes depict only ML_1 support levels); (B) topology of ML_3 tree; (C) topology of BI tree; (D) topology of MP tree. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML trees, posterior probability for BI tree, and maximum parsimony bootstrap for MP tree. Colors represent the following scores: black = 100 (1 for BI); yellow = 95–99.9 (0.95–0.99 for BI); grey = 90–94.9 (0.9–0.94 for BI); white < 90 (<0.9 for BI).
Diversity 15 00007 g002
Figure 3. Summarized phylogenetic relationships of the higher Entomobryomorpha: (A) topology of ML_1 and BI trees (nodes depict only ML_1 support levels); (B) topology of ML_2 tree; (C) topology of ML_3 tree; (D) topology of ML_4 tree; (E) topology of MP tree (Actaletidae omitted since it was recovered outside of the main Entomobryomorpha branch). Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML trees, and maximum parsimony bootstrap for MP tree. Colors represent the following scores: black = 100; yellow = 95–99.9; grey = 90–94.9; white < 90.
Figure 3. Summarized phylogenetic relationships of the higher Entomobryomorpha: (A) topology of ML_1 and BI trees (nodes depict only ML_1 support levels); (B) topology of ML_2 tree; (C) topology of ML_3 tree; (D) topology of ML_4 tree; (E) topology of MP tree (Actaletidae omitted since it was recovered outside of the main Entomobryomorpha branch). Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML trees, and maximum parsimony bootstrap for MP tree. Colors represent the following scores: black = 100; yellow = 95–99.9; grey = 90–94.9; white < 90.
Diversity 15 00007 g003
Figure 4. Summarized phylogenetic relationships of the higher Poduromorpha: (A) topology of ML_1, ML_2, ML_3 and ML_4 trees (nodes depict only ML_1 support levels); (B) topology of BI tree; (C) topology of MP tree. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML tree, posterior probability for BI tree, and maximum parsimony bootstrap for MP tree. Colors represent the following scores: black = 100 (1 for BI); yellow = 95–99.9 (0.95–0.99 for BI); grey = 90–94.9 (0.9–0.94 for BI); white < 90 (<0.9 for BI).
Figure 4. Summarized phylogenetic relationships of the higher Poduromorpha: (A) topology of ML_1, ML_2, ML_3 and ML_4 trees (nodes depict only ML_1 support levels); (B) topology of BI tree; (C) topology of MP tree. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML tree, posterior probability for BI tree, and maximum parsimony bootstrap for MP tree. Colors represent the following scores: black = 100 (1 for BI); yellow = 95–99.9 (0.95–0.99 for BI); grey = 90–94.9 (0.9–0.94 for BI); white < 90 (<0.9 for BI).
Diversity 15 00007 g004
Figure 5. Summarized phylogenetic relationships of the higher Symphypleona: (A) topology of ML_1, ML_2, ML_3 and ML_4 trees (nodes depict only ML_1 support levels); (B) topology of BI tree; (C) topology of MP tree. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML tree, posterior probability for BI tree, and maximum parsimony bootstrap for MP tree. Colors represent the following scores: black = 100 (1 for BI); yellow = 95–99.9 (0.95–0.99 for BI); grey = 90–94.9 (0.9–0.94 for BI); white < 90 (<0.9 for BI).
Figure 5. Summarized phylogenetic relationships of the higher Symphypleona: (A) topology of ML_1, ML_2, ML_3 and ML_4 trees (nodes depict only ML_1 support levels); (B) topology of BI tree; (C) topology of MP tree. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML tree, posterior probability for BI tree, and maximum parsimony bootstrap for MP tree. Colors represent the following scores: black = 100 (1 for BI); yellow = 95–99.9 (0.95–0.99 for BI); grey = 90–94.9 (0.9–0.94 for BI); white < 90 (<0.9 for BI).
Diversity 15 00007 g005
Figure 6. Summarized phylogenetic relationships of the Neelipleona genera obtained by the ML_1, ML_2, ML_3, ML_4, BI, and MP trees. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML_1 tree. Colors represent the following scores: black = 100; yellow = 95–99.9; grey = 90–94.9; white < 90.
Figure 6. Summarized phylogenetic relationships of the Neelipleona genera obtained by the ML_1, ML_2, ML_3, ML_4, BI, and MP trees. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively, for ML_1 tree. Colors represent the following scores: black = 100; yellow = 95–99.9; grey = 90–94.9; white < 90.
Diversity 15 00007 g006
Figure 7. Summarized phylogenetic relationships of the sampled higher Collembola based on ML_1 tree. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively. Colors represent the following scores: black = 100; yellow = 95–99.9; grey = 90–94.9; white < 90.
Figure 7. Summarized phylogenetic relationships of the sampled higher Collembola based on ML_1 tree. Node circles represent SH-aLRT support (left side) and bootstrap values (right side), respectively. Colors represent the following scores: black = 100; yellow = 95–99.9; grey = 90–94.9; white < 90.
Diversity 15 00007 g007
Table 1. Different types of analyses, partition schemes (for ML analyses), models, and respective tree codes used in this study and their respective resulting figures.
Table 1. Different types of analyses, partition schemes (for ML analyses), models, and respective tree codes used in this study and their respective resulting figures.
Type of AnalysisPartition SchemeModelsTree CodeFigure
Maximum LikelihoodUnpartitionedmtZOA + F + R10ML_1 *Figure 1 and Figure S6
Maximum LikelihoodPartitionedEX-EHOML_2Figure S1
Maximum LikelihoodUnpartitionedmtART + C60 + FO + RML_3Figure S2
Maximum LikelihoodPartitionedTable S2ML_4Figure S3
Bayesian Inference-CAT + GTRBIFigure S4
Maximum Parsimony--MPFigure S5
Legends: ‘*’ = main phylogeny of this study; ‘-’ = not applicable.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bellini, B.C.; Zhang, F.; de Souza, P.G.C.; dos Santos-Costa, R.C.; Medeiros, G.d.S.; Godeiro, N.N. The Evolution of Collembola Higher Taxa (Arthropoda, Hexapoda) Based on Mitogenome Data. Diversity 2023, 15, 7. https://doi.org/10.3390/d15010007

AMA Style

Bellini BC, Zhang F, de Souza PGC, dos Santos-Costa RC, Medeiros GdS, Godeiro NN. The Evolution of Collembola Higher Taxa (Arthropoda, Hexapoda) Based on Mitogenome Data. Diversity. 2023; 15(1):7. https://doi.org/10.3390/d15010007

Chicago/Turabian Style

Bellini, Bruno Cavalcante, Feng Zhang, Paolla Gabryelle Cavalcante de Souza, Renata Clicia dos Santos-Costa, Gleyce da Silva Medeiros, and Nerivânia Nunes Godeiro. 2023. "The Evolution of Collembola Higher Taxa (Arthropoda, Hexapoda) Based on Mitogenome Data" Diversity 15, no. 1: 7. https://doi.org/10.3390/d15010007

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