Research Article |
Corresponding author: Xing Wang ( xingwanghjt@163.com ) Corresponding author: Guo-Hua Huang ( ghhuang@hunau.edu.cn ) Academic editor: Anna Hundsdoerfer
© 2023 Min Deng, Andreas Zwick, Qi Chen, Cheng-Qing Liao, Wei Wang, Xing Wang, Guo-Hua Huang.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
|
The small, relict-like moth family Endromidae is well-established within the superfamily Bombycoidea, but relationships within the family have remained vague for the last decade, primarily due to very limited taxon sampling. This resulted in the explicit removal of all internal suprageneric classification by
Endromidae, mitochondrial genome, phylogenetic analysis, revision
The moth family Endromidae Boisduval, 1828 is relatively species poor (72 species) and occurs primarily in Asia, with just a single species extending into Europe. This species, the very distinct and widespread Palearctic species Endromis versicolora (Linnaeus, 1758), is placed in the monobasic genus Endromis Ochsenheimer, 1810 and, for over a century, its own family Endromidae. It was regarded as an isolated taxon within the “Bombyces”, until
The use of DNA sequence data has greatly contributed to the phylogenetic hypotheses and consequential changes in the classification of Bombycoidea. Based on five protein-coding nuclear genes, phylogenetic analyses of the ‘bombycoid complex’ grouped Endromidae, Mirinidae and the bombycid subfamilies Oberthueriinae Kuznetsov & Stekolnikov, 1985 and Prismostictinae Forbes, 1955 into a single clade (
With further studies on the phylogeny of Bombycoidea (
Several endromid genera (e.g., Dalailama Staudinger, 1896, Sesquiluna Forbes, 1955 and Falcogona) are rarely collected and underrepresented in collections, making it difficult to obtain comprehensive taxon sampling for molecular phylogenetic studies. DNA sequencing of old collection specimens helps to improve taxon sampling of rarely collected species (
The mt-genomes of 26 taxa of Endromidae, belonging to 13 genera in four subfamilies (Table S1), were newly sequenced for this study. Following collection, the three right legs or thorax of each specimen were preserved in absolute ethanol and then stored at –20°C. The remainders of all specimens are deposited as vouchers in the Insect Museum of Hunan Agricultural University, Changsha City, Hunan Province, China.
Total genomic DNA was extracted from the legs or thoracic tissue of each specimen using a TaKaRa MiniBEST Universal Genomic DNA Extraction Kit Ver.5.0 (Shiga Prefecture, Kusatsu City, Japan). Purified DNA was preserved at –20°C prior to sequencing.
Illumina TruSeq libraries with 350 bp insert size were prepared for each species, and these Whole Genome Shotgun libraries were sequenced by Novogene (Beijing, China) on the Illumina Hiseq platform with 150 bp paired-end reads. For each library, 6 Gbp of clean data were obtained after removing reads containing adaptor contamination, poly-Ns (>15 bp Ns), or >75 bp bases with quality scores ≤ 20. Cleaned reads were assembled into contigs and scaffolds using IDBA (v1.1.3;
Twelve publicly available mt-genomes were obtained from NCBI GenBank (http://www.ncbi.nlm.nih.gov), including nine ingroup species and three lasiocampid species used as outgroups. All mt-genome sequences were imported and standardized in Geneious R8. All PCGs were exported from Geneious R8. The 13 PCGs were aligned with the TranslatorX server (http://www.translatorx.co.uk) (Castresana 2000), with the “ALL”-parameter. Concatenation of single gene alignments was performed in Geneious R8, resulting in 2 datasets: 1) 13 protein-coding genes (13PCGs); 2) 13 PCGs as amino acids (13PCGs-AA). Partitionfinder 2.1.1 was used to search the optimal partitioning scheme and models for each data set. IQ-Tree (v1.5.5;
Our study provides the first mt-genomes for six endromid genera, i.e., Endromis, Mirina, Pseudandraca Miyata, 1970, Smerkata Zolotuhin, 2007, Dalailama and Promustilia Zolotuhin, 2007. Most mt-genomes, except Andraca lawa_21, Mustilizans dierli and Prismosticta tiretta_32, used in this study comprised a total of 37 genes (13 PCGs, 22 tRNAs and 2 rRNAs), and the total length of all sequences ranges from an incomplete 6,350 bp (Andraca lawa_21) to 15,880 bp (Andraca olivacea−GD) (Figure S1). As is usual for the mt-genomes of Lepidoptera (
Based on the 13PCGs dataset, the two phylogenetic trees estimated with BI (Figure S4) and ML analyses (Figure S5) are almost identical. Only the relationships between the three clades of Oberthuriinae and the relationships between species of Primosticta differ. And the Bayesian tree of the 13-PCGs data set is used to label the values of other trees. Statistical support (PP and BP) is strong for 34 of the 40 nodes, with weak support restricted to backbone nodes and within a clade of Oberthueriinae stat. rev. Within the limits of taxon sampling, the results strongly support the monophyly of the family Endromidae, as well as a division into four subfamilies, with Prismostictinae stat. rev. sister to all other taxa (PP = 1, BP = 100%). The phylogenetic relationship between Endrominae and Mirininae stat. rev., which are nested between the two other subfamilies, is strongly supported in the BI tree (PP = 0.992). Within Prismostictinae stat. rev., the genera Prismosticta Butler, 1880 and Prismostictoides are sister to each other (PP = 1, BP = 99.5%). A distinct monophyletic group, subfamily Oberthueriinae stat. rev. comprises three major, well-supported clades (PP = 1, BP = 100%). Clade 1 includes the genera Pseudandraca Miyata, 1970 and Andraca Walker, 1865. In addition, Andraca gongshanensis and Pseudandraca flavamaculata are shown as sisters with strong support (PP = 1, BP = 99.8%). The clades (Andraca apodecta + Andraca melli) and (Andraca trilochoides + (Andraca draco + Andraca lawa) are grouped together (PP = 1, BP = 100%). Clade 2 includes four genera, with Comparmustilia Wang, X. & Zolotuhin, 2015 sister to the other genera. Clade 3 is divided into two major groups comprising the monophyletic genus Oberthueria Kirby, 1892 (PP = 1, BP = 100%) and the genus Mustilizans Yang, 1995, which is a paraphyletic relative to the species Promustilia yajiangensis Wang, X. & Zolotuhin, 2015. The relationship of these three clades is shown as (Clade 1 + (Clade 2 + Clade 3)) (PP = 0.771, BP = 49.8). Although both methods of analysis of the different datasets resulted in largely congruent topologies, there are still obvious differences compared to the analysis results of the 13PCGs-AA dataset (BI in Figure S6, ML in Figure S7, and combined in Figure
All the newly sequenced Endromidae mitochondrial genomes have the same gene order as in the other known Lepidoptera (Cao et al. 2012;
Mitochondrial genomes are widely used for studying population genetics, comparative and evolutionary genomics, the reconstruction of phylogenetic relationships, and evolutionary biology (e.g.,
Our results are likewise in agreement with morphologically recognized group. Our analyses recovered both the family Endromidae and its four major lineages as monophyletic and strongly supported. These four lineages correspond to the morphologically recognized (sub)families that were synonymized with Endromidae (
Endrominae, which is just consisting of one species, has typical characteristics different from other endromid moths as follows: forewing with three triangular white spots, thorax and abdomen bright-colored (Figure
Phylogeny of Endromidae inferred from different data sets (13PCGs-AA, 13PCGs) using Bayesian inference and maximum likelihood analyses. Numbers above branches are posterior probabilities (BI PP), beneath which are bootstrap percentages (ML BP) for 1000 replicates; nodes with maximum support values are marked with a black dot instead. Dashed arrows (two in total) identify alternative topologies (relative to the topology shown) that receive at least 70% bootstrap support by one or more of the approaches. The asterisks indicate newly sequenced mitochondrial genomes.
Mirininae stat. rev. is consisting of only one genus, Mirina, which was been controversial. Some scholars considered that it should be a separate family (
The subfamily Prismostictinae stat. rev., which is sister to a clade comprising all the other subfamilies, consists of only two morphologically similar genera, Prismosticta and Prismostictoides. The monobasic genus Prismostictoides, with the type species Prismosticta unihyala Chu & Wang, 1993, was separated from Prismosticta based on a broad postmedial line and broad, dark yellow submarginal band on the forewing upperside, uncus with a long uncuslike projiection arising from the base of uncus, and valva asymmetrical. Otherwise, the two genera are rather similar and share a characteristic transparent spot near the apex of the forewing (
Oberthueriinae stat. rev., the largest subfamily of Endromidae, is divided into three major clades that are strongly supported as monophyletic groups. Based on our results, we treat these clades as three tribes: Andracini tribe nov. (Clade 1), Mustiliini tribe nov. (Clade 2) and Oberthueriini Kuznetzov & Stekolnikov, 1985 stat. rev. (Clade 3). The sister relationship between Andraca and Pseudandraca in Andracini tribe nov. is consistent with findings of previous studies (
Previously, Mustilia Walker, 1865 had been split into six separate genera (
Our analyses place three genera in the tribe Oberthueriini, i.e., Oberthueria, Mustilizans and Promustilia, of which the latter two were previously included in Mustilia. The genus Oberthueria is monophyletic, with strong statistical support for the three species included in this study, and the six currently recognized species are morphologically very similar and show only moderate differences in COX1 barcode sequences (
Unfortunately, specimens of the remaining three endromid genera, Falcogona, Sesquiluna and Theophoba Fletcher & Nye, 1982 were unavailable to us and we lack mt-genomes for them. Based on similarities in adult morphology, such as eyes surrounded by setae, completely pectinate antennae and similar size, we tentatively include Sesquiluna and Theophoba in subfamily Prismostictinae stat. rev., as postulated by
Oberthueria Kirby, 1892
Members of this tribe share the following characters: 1) forewings long and narrow; 2) labial palpi of moderate length, about 2/3 of the vertical diameter of the compound eye; 3) distal half of antenna devoid of well-developed rami (Figure
Oberthueria Kirby, 1892, Syn. Cat. Lepid. Het., 1: 720. Type species: Euphranor caeca Oberthür, 1880, by monotypy
Oberthueria Staudinger, 1892, in Romanoff, Mém. Lépid.: 337 Type species: Euphranor caeca Oberthür, 1880, by monotypy (a junior homonym and junior objective synonym of Oberthueria Kirby, 1892)
Oberthüria : Staudinger, 1892, in Romanoff, Mémoires sur les lepidoptères (Mém. Lépid.) 6: 337. (incorrect original spelling)
Euphraor: Kirby, 1892, Syn. Cat. Lepid. Het. 1: 720 (incorrect subsequent spelling)
Euphranor Oberthür, 1880, Etudes d’Entomologie (Étud. ent.) 5: 40. Type species: Euphranor caeca Oberthür, 1880, by monotypy (a junior homonym of Euphranor Herrich-Schäffer, 1855 (Lepidoptera, Saturniidae))
Mustilizans Yang, 1995, Insects of Baishanzu Mountain, Eastern China: 355. Type species: Mustilizans drepaniformis Yang, 1995, by original designation
Promustilia Zolotuhin, 2007, Neue ent. Nachr. 60: 199. Originally erected as a subgenus of Mustilizans Yang, 1993. Type species: Mustilizans (Promustilia) andracoides Zolotuhin, 2007, by original designation
Mustilia Walker, 1865
Mustiliini tribe nov. is similar to Oberthueriini stat. rev. in having narrow forewings and the distal half of the antenna with underdeveloped rami, but it can be easily distinguished by particularly short or completely reduced labial palpi. Larvae of this tribe also possess a hairless body, but the thoracic tergites are laterally strongly expanded, and the anal horn is relatively long and stout (Figure
Based on the above morphological characteristics, we here establish the new tribe Mustiliini tribe nov., which is also supported by molecular data. Although we lack molecular data for Falcogona, it is included in this tribe because of its similarity in habitus and male genital structures to Smerkata.
Comparmustilia Wang, X. & Zolotuhin, 2015, Zootaxa, 3989: 79. Type species: Mustilia sphingiformis Moore, 1879, by present designation
Smerkata Zolotuhin, 2007, Neue ent. Nachr. 60: 193. Originally proposed as a subgenus of Mustilia Walker, 1865. Type species: Mustilia phaeopera Hampson, 1910, by original designation
Dalailama Staudinger, 1896, Dt. ent. Z. Iris 8 (2): 303. Type species: Dalailama bifurca Staudinger, 1896, by monotypy
Dailalama: Staudinger, 1901, Cat. Lepid. palaearct. Faunengeb. (1): 128. Incorrect subsequent spelling
Deilelamia: Pagenstecher, 1909, Geschichte eur. Schmett.: 433. Incorrect subsequent spelling
Mustilia Walker, 1865, List Specimens lepid. Insects Colln Br. Mus. 32: 580. Type species: Mustilia falcipennis Walker, 1865, by monotypy
Falcogona Zolotuhin, 2007, Neue ent. Nachr. 60: 199. Type species: Falcogona gryphea Zolotuhin, 2007, by original designation
Andraca Walker, 1865
Morphological synapomorphies supporting the monophyly of Andracini tribe nov. are the relatively broad forewings, the very long labial palpi (longer than the vertical diameter of the compound eye), and underdeveloped rami over the distal 1/3 of the antenna (Figure
Pseudandraca was established by
Andraca Walker, 1865, List Specimens lepid. Insects Colln Br. Mus. 32: 581. Type species: Andraca bipunctata Walker, 1865, by monotypy
Pseudoeupterote Shiraki, 1911, Catalogue Insectorum Noxiorum Formosarum: 48. Type species: Oreta theae Matsumura, 1909, by monotypy
Pseudandraca Miyata, 1970, Tinea. 8: 190. Type species: Andraca gracilis Butler, 1885, by original designation
The complete sequences were uploaded to the NCBI (GenBank accession numbers: OQ472264–OQ472285). The incomplete sequences provided in this article can be accessed from Zenodo, DOI: https://doi.org/10.5281/zenodo.7655269
The authors have no conflicts of interests to declare.
M.D. was responsible for drafting the manuscript, as well as the acquisition, analysis and interpretation of data. A.Z. provided part molecular sequences and contributed to the conception and design of the current study. Q.C. analyzed and interpreted the data. W.W. made suggestions and revised the manuscript. X.W. confirmed the insect species and revised the manuscript. G.-H.H. helped perform the analysis with constructive discussions and provided financial support. All authors read and approved the final manuscript.
The authors thank Prof. Min Wang, Dr. Hou-Shuai Wang (South China Agricultural University, China), Ms. Zhuang-Mei Chen, Meng-Yue Chen and Si-Jia Yi, Mr. Bin Chen and Lu Chen (Hunan Agricultural University) for their kind help in collecting the samples. This study was supported by the National Natural Science Foundation of China (31970450, 32111540167), and China Agriculture Research System (CARS-23-C08), and Guangdong Provincial Academy of Sciences Special Project of Science and Technology Development (2020GDASYL-20200102021).
Table S1
Data type: .docx
Explanation note: List of taxa (genera and species) examined for the study and sources of information.
Table S2
Data type: .docx
Explanation note: The characteristics of the mitochondrial genomes of Endromidae.
Figure S1
Data type: .pdf
Explanation note: Characters of thirty-five newly sequenced endromid species’ mitochondrial genomes. Gene names are annotated using standard abbreviations; single letters are IUPAC amino acid abbreviations.
Figure S2
Data type: .pdf
Explanation note: Boxplot showing P−distances between all 47 samples for each of the 13 genes analyzed. Outlier values are depicted as points outside of the box.
Figure S3
Data type: .pdf
Explanation note: Boxplot showing the GC-ratios of all 47 samples for each of the 13 genes analyzed. Outlier values are depicted as points outside of the box.
Figure S4
Data type: .pdf
Explanation note: Bayesian inference phylogram constructed with the 13PCGs data set of Endromidae.
Figure S5
Data type: .pdf
Explanation note: Maximum likelihood phylogram constructed with the 13PCGs data set of Endromidae.
Figure S6
Data type: .pdf
Explanation note: Bayesian inference phylogram constructed with the 13PCGs-AA data set of Endromidae.
Figure S7
Data type: .pdf
Explanation note: Maximum likelihood phylogram constructed with the 13PCGs-AA data set of Endromidae.