Abstract:Background A robust and stable phylogenetic framework is a fundamental goal of evolutionary biology. As the third largest insect order following by Diptera and Coleoptera in the world, lepidoptera (butterflies and moths) play a central role in almost every terrestrial ecosystem as the indicators of environmental change and serve as important models for biologists exploring questions related to ecology and evolutionary biology. However, for such charismatic insect group, the higher-level phylogenetic relationships among its superfamilies are still poorly unresolved. Results we increased taxon sampling among Lepidoptera (40 superfamilies and 76 families contained 286 taxa) and filtered the unqualified samples, then acquired a series of large amino-acid datasets from 69,680 to 400,330 for phylogenomic reconstructions. Using these datasets, we explored the effect of different taxon sampling on tree topology by considering a series of systematic errors using ML and BI methods. Moreover, we also tested the effectiveness in topology robustness among the three ML-based models. The results showed that taxon sampling is an important determinant in tree robustness of accurate lepidopteran phylogenetic estimation. Long-branch attraction (LBA) caused by site-wise heterogeneity is a significant source of bias given rise to topologies divergence of ditrysia in phylogenomic reconstruction. Phylogenetic inference showed a most comprehensive framework by far to reveal the relationships among lepidopteran superfamilies, but limited by taxon sampling, it could only represent the current understanding of the lepidopteran tree of life. The relationships within the species-rich and relatively rapid radiation Ditrysia and especially Apoditrysia remain poorly unresolved, which need to increase taxon sampling and adopt lineage-specific genes for further phylogenomic reconstruction. Conclusions The present study further expands the taxon sampling of lepidopteran phylogeny and provides a potential phylogenomic foundation for further understanding its current higher-level relationships. ### Competing Interest Statement The authors have declared no competing interest. * BUSCO : Benchmarking universal single-copy orthologue LBA : Long-branch attraction NGS : next-generation sequencing LCWGS : low-coverage whole-genome sequencing AHE : anchored hybrid enrichment UCE : ultraconserved element enrichment aTRAM : automated target restricted assembly method aUBF : Average ultrafast bootstraps GHOST : General Heterogeneous evolution On a Single Topology PMSF : posterior mean site frequency MSC : multi-species coalescent gCF : gene concordance factor sCF : site concordance factor ML : maximum-likelihood BI : Bayesian inference SH-aLRT : Shimodaira-Hasegawa-like approximate likelihood ratio test UFBS : ultrafast bootstraps ILS : incomplete lineage sorting EBP : Earth BioGenome Project HGT : horizontal gene transfer NCBI : National center for biotechnology information BOLD : Barcode of life data system COI : cytochrome c oxidase subunit I LPP : local posterior probabilities AU : approximately the unbiased WKH : weighted Kishino-Hasegawa WSH : weighted Shimodaira-Hasegawa SRH : stationary, reversible and homogeneous PISs : parsimony informative sites.

Phylogeny of the superfamily Gelechioidea (Lepidoptera: Obtectomera), with an exploratory application on geometric morphometrics

Comparative mitochondrial genome analysis and phylogenetic relationship among lepidopteran species

Molecular Phylogeny of Hesperiidae (lepidoptera) with an Emphasis on Asian and African Genera

Phylogenetic Analysis Within Lepidoptera Based on the Concatenated 11 Protein-coding Gene Sequences of Mitochondrial Genomes

Phylogenomic data exploration with increased sampling provides new insights into the higher-level relationships of butterflies and moths (Lepidoptera)

Phylogenomics reveals the relationships of butterflies and moths (Lepidoptera): providing the potential landscape using universal single copy orthologues

Enhanced Resolution of Evolution and Phylogeny of the Moths Inferred from Nineteen Mitochondrial Genomes.

Phylogenomics and systematics of Entomobryoidea (Collembola): marker design, phylogeny and classification

The Genus Leucophenga (diptera, Drosophilidae), Part II: the Ornata Species Group from East Asia, with Morphological and Molecular Evidence (I).

Phylogeny of Membracoidea (hemiptera: Auchenorrhyncha) Based on Transcriptome Data

Phylogeny of the thaphium butterflies inferred from nuclear 28 S rDNA and mitochondrial ND 5 gene sequences

Characterization of Seventeen Complete Mitochondrial Genomes: Structural Features and Phylogenetic Implications of the Lepidopteran Insects

The Complete Mitochondrial Genome of Brachmia Macroscopa (lepidoptera: Gelechiidae) and Its Related Phylogenetic Analysis.

Evolution of Green Lacewings (neuroptera: Chrysopidae): a Molecular Supermatrix Approach

Evolution of lacewings and allied orders using anchored phylogenomics ( <scp>N</scp> europtera, <scp>M</scp> egaloptera, <scp>R</scp> aphidioptera)

Phylogenomics of the Superfamily Dytiscoidea (coleoptera: Adephaga) with an Evaluation of Phylogenetic Conflict and Systematic Error

Mitogenomic phylogenetic analyses provide novel insights into the taxonomic problems of several hesperiid taxa (Lepidoptera: Hesperiidae)

Phylogeny and Classification of the Leafhopper Subfamily Eurymelinae (hemiptera: Cicadellidae) Inferred from Molecules and Morphology

The Complete Mitochondrial Genome of the Golden Birdwing Butterfly Troides Aeacus: Insights into Phylogeny and Divergence Times of the Superfamily Papilionoidea

A Molecular Phylogeny of the Hawkmoth Genus Hyles (lepidoptera: Sphingidae, Macroglossinae)

Structural Features and Phylogenetic Implications of Cicadellidae Subfamily and Two New Mitogenomes Leafhoppers.