The functional evolution of collembolan Ubx on the regulation of abdominal appendage formation.

Folsomia candida is a tiny soil-living arthropod belonging to the Collembola, which is an outgroup to Insecta. It resembles insects as having a pair of antennae and three pairs of thorax legs, while it also possesses three abdominal appendages: a ventral tube located in the first abdominal segment (A1), a retinaculum in A3, and a furca in A4. Collembolan Ubx and AbdA specify abdominal appendages, but they are unable to repress appendage marker gene Dll. The genetic basis of collembolan appendage formation and the mechanisms by which Ubx and AbdA regulate Dll transcription and appendage development remains unknown. In this study, we analysed the developmental transcriptomes of F. candida and identified candidate appendage formation genes, including Ubx (FcUbx). The expression data revealed the dominance of Dll over Ubx during the embryonic 3.5 and 4.5 days, suggesting that Ubx is deficient in suppressing Dll at early appendage formation stages. Furthermore, via electrophoretic mobility shift assays and dual luciferase assays, we found that the binding and repression capacity of FcUbx on Drosophila Dll resembles those of the longest isoform of Drosophila Ubx (DmUbx_Ib), while the regulatory mechanism of the C-terminus of FcUbx on Dll repression is similar to that of the crustacean Artemia franciscana Ubx (AfUbx), demonstrating that the function of collembolan Ubx is intermediate between that of Insecta and Crustacea. In summary, our study provides novel insights into collembolan appendage formation and sheds light on the functional evolution of Ubx. Additionally, we propose a model that collembolan Ubx regulates abdominal segments in a context-specific manner.

E93 is indispensable for reproduction in ametabolous and hemimetabolous insects.

Ecdysone-induced protein 93 (E93), known as the 'adult-specifier' transcription factor in insects, triggers metamorphosis in both hemimetabolous and holometabolous insects. Although E93 is conserved in ametabolous insects, its spatiotemporal expression and physiological function remain poorly understood. In this study, we first discover that, in the ametabolous firebrat Thermobia domestica, the previtellogenic ovary exhibits cyclically high E93 expression, and E93 mRNA is broadly distributed in previtellogenic ovarioles. E93 homozygous mutant females of T. domestica exhibit severe fecundity deficiency due to impaired previtellogenic development of the ovarian follicles, likely because E93 induces the expression of genes involved in ECM (extracellular matrix)-receptor interactions during previtellogenesis. Moreover, we reveal that in the hemimetabolous cockroach Blattella germanica, E93 similarly promotes previtellogenic ovarian development. In addition, E93 is also essential for vitellogenesis that is necessary to guarantee ovarian maturation and promotes the vitellogenesis-previtellogenesis switch in the fat body of adult female cockroaches. Our findings deepen the understanding of the roles of E93 in controlling reproduction in insects, and of E93 expression and functional evolution, which are proposed to have made crucial contributions to the origin of insect metamorphosis.

300 Million years of coral treaders (Insecta: Heteroptera: Hermatobatidae) back to the ocean in the phylogenetic context of Arthropoda.

Among hundreds of insect families, Hermatobatidae (commonly known as coral treaders) is one of the most unique. They are small, wingless predaceous bugs in the suborder Heteroptera. Adults are almost black in colour, measuring about 5 mm in body length and 3 mm in width. Thirteen species are known from tropical coral reefs or rocky shores, but their origin and evolutionary adaptation to their unusual marine habitat were unexplored. We report here the genome and metagenome of Hermatobates lingyangjiaoensis, hitherto known only from its type locality in the South China Sea. We further reconstructed the evolutionary history and origin of these marine bugs in the broader context of Arthropoda. The dated phylogeny indicates that Hexapoda diverged from their marine sister groups approximately 498 Ma and that Hermatobatidae originated 192 Ma, indicating that they returned to an oceanic life some 300 Myr after their ancestors became terrestrial. Their origin is consistent with the recovery of tropical reef ecosystems after the end-Triassic mass extinction, which might have provided new and open niches for them to occupy and thrive. Our analyses also revealed that both the genome changes and the symbiotic bacteria might have contributed to adaptations necessary for life in the sea.

IFT52 plays an essential role in sensory cilia formation and neuronal sensory function in Drosophila.

Cilia are microtubule-based, hair-like organelles involved in sensory function or motility, playing critical roles in many physiological processes such as reproduction, organ development, and sensory perception. In insects, cilia are restricted to certain sensory neurons and sperms, being important for chemical and mechanical sensing, and fertility. Although great progress has been made regarding the mechanism of cilia assembly, the formation of insect cilia remains poorly understand, even in the insect model organism Drosophila. Intraflagellar transport (IFT) is a cilia-specific complex that traffics protein cargos bidirectionally along the ciliary axoneme and is essential for most cilia. Here we investigated the role of IFT52, a core component of IFT-B, in cilia/flagellar formation in Drosophila. We show that Drosophila IFT52 is distributed along the sensory neuronal cilia, and is essential for sensory cilia formation. Deletion of Ift52 results in severe defects in cilia-related sensory behaviors. It should be noted that IFT52 is not detected in spermatocyte cilia or sperm flagella of Drosophila. Accordingly, ift52 mutants can produce sperms with normal motility, supporting a dispensable role of IFT in Drosophila sperm flagella formation. Altogether, IFT52 is a conserved protein essential for sensory cilia formation and sensory neuronal function in insects.

High-quality genomes reveal significant genetic divergence and cryptic speciation in the model organism Folsomia candida (collembola).

The collembolan Folsomia candida Willem, 1902, is widely distributed throughout the world and has been frequently used as a test organism in soil ecology and ecotoxicology studies. However, it is questioned as an ideal "standard" because of differences in reproductive modes and cryptic genetic diversity between strains from various geographical origins. In this study, we obtained two high-quality chromosome-level genomes of F. candida, for a parthenogenetic strain (named FCDK, 219.08 Mb, 25,139 protein-coding genes) and a sexual strain (named FCSH, 153.09 Mb, 21,609 protein-coding genes), reannotated the genome of the parthenogenetic strain reported by Faddeeva-Vakhrusheva et al. in 2017 (named FCBL, 221.7 Mb, 25,980 protein-coding genes) and conducted comparative genomic analyses of the three strains. High genome similarities between FCDK and FCBL based on synteny, genome architecture, mitochondrial and nuclear gene sequences suggest that they are conspecific. The seven chromosomes of FCDK are each 25%-54% larger than the corresponding chromosomes of FCSH, showing obvious repetitive element expansions and large-scale inversions and translocations but no whole-genome duplication. The strain-specific genes, expanded gene families and genes in nonsyntenic chromosomal regions identified in FCDK are highly related to the broader environmental adaptation of parthenogenetic strains. In addition, FCDK has fewer strain-specific microRNAs than FCSH, and their mitochondrial and nuclear genes have diverged greatly. In conclusion, FCDK/FCBL and FCSH have accumulated independent genetic changes and evolved into distinct species after 10 million years ago. Our work provides important genomic resources for studying the mechanisms of rapidly cryptic speciation and soil arthropod adaptation to soil ecosystems.

Convergent Adaptation of Ootheca Formation as a Reproductive Strategy in Polyneoptera.

Insects have evolved numerous adaptations and colonized diverse terrestrial environments. Several polyneopterans, including dictyopterans (cockroaches and mantids) and locusts, have developed oothecae, but little is known about the molecular mechanism, physiological function, and evolutionary significance of ootheca formation. Here, we demonstrate that the cockroach asymmetric colleterial glands produce vitellogenins, proline-rich protein, and glycine-rich protein as major ootheca structural proteins (OSPs) that undergo sclerotization and melanization for ootheca formation through the cooperative protocatechuic acid pathway and dopachrome and dopaminechrome subpathway. Functionally, OSP sclerotization and melanization prevent eggs from losing water at warm and dry conditions, and thus effectively maintain embryo viability. Dictyopterans and locusts convergently evolved vitellogenins, apolipoprotein D, and laminins as OSPs, whereas within Dictyoptera, cockroaches and mantids independently developed glycine-rich protein and fibroins as OSPs. Highlighting the ecological-evolutionary importance, convergent ootheca formation represents a successful reproductive strategy in Polyneoptera that promoted the radiation and establishment of cockroaches, mantids, and locusts.

The mitochondrial genome of a minute springtail species Megalothorax incertus (Collembola: Neelipleona: Neelidae).

In this study, the complete mitochondrial genome for the Wolbachia infected parthenogenetic collembola Megalothorax incertus Börner, 1903 was determined. It represents the first report of a complete mitochondrial genome from Neelipleona, one of the four orders of Collembola. The circularized 14,994 bp mitochondrial genome sequence consists of canonical 37 mito-genes, including 13 protein-coding genes (PCGs), 22 tRNA genes, and two rRNA genes. The base composition of the majority strand (same direction as most of the mitochondrial genes) is 32.0% for A, 24.1% for C, 11.9% for G, and 32.1% for T. The phylogenetic trees inferred from 13 PCGs using the Maximum-likelihood (ML) and Bayesian inference (BI) methods suggested that Neelidae is basal to the remaining springtails. This mitochondrial genome provides new insights to decipher the phylogeny of Collembola.

Grainy head signaling regulates epithelium development and ecdysis in Blattella germanica.

The transcription factor grainy head (Grh) functions in the protection of the epithelium against the external environment by generating strongly adhesive layers, and this function is conserved in vertebrates and invertebrates. In Drosophila, the top model for holometabolous insects, Grh is necessary during embryonic development, epidermal differentiation, central nervous system specification and epithelial repair. However, the function of this gene in hemimetabolous insect epithelia remains unknown. To examine the function of Grh signaling in regulating epithelium development in Hemimetabola, we focused on the Blattella germanica epidermal layer using a gene knockdown strategy. The spatiotemporal expression pattern of BgGrh was detected, and knockdown of BgGrh and BgCad96ca, which provide positive feedback to BgGrh, caused severe defects in new epithelium development and impeded the molting process required to discard the old integument. Knockdown of the expression of BgGrh and BgCad96ca caused increased expression of chitin synthase gene (BgCHS1) and chitinase gene (BgCht5), the upregulations of which should be mediated by the higher level of hormone receptor 3 (BgHr3) gene. In conclusion, epithelium development is regulated by Grh signaling, which might represent a potential target for the control of urban pest cockroaches.

Selection of Reference Genes for Normalization of Gene Expression in Thermobia domestica (Insecta: Zygentoma: Lepismatidae).

Zygentoma occupies a key evolutionary position for understanding the evolution of insect metamorphosis but has received little attention in terms of genetic analysis. To develop functional genomic studies in this insect, we evaluated five candidate internal reference genes for quantitative RT-PCR (qPCR) studies from Thermobia domestica, a representative species of Zygentoma, including Actin 5C (Actin5C), Elongation factor-1 alpha (EF1A), Ribosome protein S26 (RPS26), Ribosome protein L32 (RPL32), and Superoxide dismutase 2 (SOD2), at different developmental stages, in various body parts, and under dsRNA microinjection and starvation stresses, using four algorithms (delta Ct, geNorm, NormFinder and BestKeeper) and a comparative algorithm (RefFinder). Specific suitable reference genes were recommended across specific experimental conditions, and the combination of RPS26 and RPL32 was appropriate for all tested samples. Employing our selected reference gene combination, we investigated the gene expression pattern of Myoglianin (Myo), a crucial gene-regulating insect metamorphosis, in ametabolous T. domestica, and demonstrated the efficiency of RNA interference (RNAi) in firebrat nymphs. This study provides a basis for reliable quantitative studies of genes and greatly benefits evolutionary and functional genomics studies in Zygentoma.

New and little known Isotomidae (Collembola) from the shore of Lake Baikal and saline lakes of continental Asia.

Collembola of the family Isotomidae from the shores of Lake Baikal and from six saline lake catenas of the Buryat Republic (Russia) and Inner Mongolia Province (China) were studied. Pseudanurophorus barathrum Potapov & Gulgenova, sp. nov. and Parisotoma baicalica Potapov & Gulgenova, sp. nov. from Baikal and Ephemerotoma buryatica Potapov, Huang & Gulgenova, sp. nov. and Folsomia mongolica Huang & Potapov, sp. nov. from saline lakes are described here. A morphological description of epitokous males of Scutisotoma acorrelata Potapov, Babenko & Fjellberg, 2006 is given. A list of 23 species of the family Isotomidae found in the shores of studied lakes is provided based on literature sources and newly collected material.

Protura from Hainan Island, China: new species, checklist and distribution.

More than 1500 proturan specimens from Hainan Island are systematically studied. An annotated list of all species of Protura from Hainan Island is provided and their geographical distribution is discussed. The genus Paracondeellum is reported from Hainan Island for the first time, and Paracondeellum paradisum sp. nov. is described. The type species Paracondeellum dukouense (Tang & Yin, 1988) is redescribed based on syntype, and the lectotype and paralectotype are designated. The characters of the genus Paracondeellum are redefined, and the two known species are compared in detail. The Protura fauna of Hainan Island is mainly composed of species from the Oriental region, with 91% of the species belonging to the families Berberentulidae and Eosentomidae.

A multi-generational risk assessment of Cry1F on the non-target soil organism Folsomia candida (Collembola) based on whole transcriptome profiling.

The Bacillus thuringiensis toxin Cry1F has been used to develop insect-resistant genetically engineered crops. There has been great interest in evaluating its potential risk to non-target organisms (NTOs). However, the majority of previous risk assessments only examined one generation of NTOs using several physiological indicators, which cannot comprehensively detect some potential sub-lethal effects at the molecular level. In this study, we conducted a laboratory-based, multi-generational risk assessment of Cry1F for the collembolan Folsomia candida, an important representative of soil arthropods in terms of survival, reproduction, and differentially expressed genes (DEGs) identified from whole transcriptome profiles. Our results demonstrated that Cry1F was continuously ingested by collembolans over three consecutive generations, but it did not affect the survival or reproduction of F. candida. There were no significant differences in the global gene expression between F. candida-fed diets with and without Cry1F, and no consistent co-expressed DEGs over three generations. In addition, Cry1F did not obviously alter the expression profiles of seven sensitive biological markers. Our composite data indicates that Cry1F had no long-term harmful effects on collembolan F. candida.

Going Deeper into High and Low Phylogenetic Relationships of Protura.

Proturans are small, wingless, soil-dwelling arthropods, generally associated with the early diversification of Hexapoda. Their bizarre morphology, together with conflicting results of molecular studies, has nevertheless made their classification ambiguous. Furthermore, their limited dispersal capability (due to the primarily absence of wings) and their euedaphic lifestyle have greatly complicated species-level identification. Mitochondrial and nuclear markers have been applied herein to investigate and summarize proturan systematics at different hierarchical levels. Two new mitochondrial genomes are described and included in a phylum-level phylogenetic analysis, but the position of Protura could not be resolved with confidence due to an accelerated rate of substitution and extensive gene rearrangements. Mitochondrial and nuclear loci were also applied in order to revise the intra-class systematics, recovering three proturan orders and most of the families/subfamilies included as monophyletic, with the exception of the subfamily Acerentominae. At the species level, most morphologically described species were confirmed using molecular markers, with some exceptions, and the advantages of including nuclear, as well as mitochondrial, markers and morphology are discussed. At all levels, an enlarged taxon sampling and the integration of data from different sources may be of significant help in solving open questions that still persist on the evolutionary history of Protura.

DNA barcoding and an updated key to the genus Hesperentomon (Protura: Acerentomata: Hesperentomidae), with a new species from Northwest China.

Hesperentomon bolense sp. n., described from Northwest China, is distinguished from other Hesperentomon spp. by having foretarsal sensillum b clearly shorter than c, a robust sensillum c broader than the other foretarsal sensilla, 18 posterior setae on the mesonotum and 16 posterior setae on the metanotum, 12 posterior setae on urotergites II-VI (P1a and P2a absent), 8 posterior setae on urosternites IV-VI (Pc absent), and absence of seta P2a on urotergite VII. The morphological characters of all described species of Hesperentomon are compared, and an updated identification key to species is given. The DNA barcodes of H. bolense sp. n. were compared to similar species of the genus; the K2P genetic divergences were 0-5.1% between individuals of H. bolense sp. n., but 21.4%-27.3% between H. bolense sp. n. and other congeneric species. The molecular data further confirmed H. bolense as a distinct species.

A new species and a new record of Sinella Brook (Collembola: Entomobryidae) from Guizhou Province, China.

Two species of Sinella Brook were discovered in an expedition to Maolan National Nature Reserve, Guizhou Province, China: S. maolanensis sp. nov. and S. colorata Zhang, Qu Deharveng. The new species is characterized by the postlabial chaetae X, X2 and X4 ciliate and X3 smooth and the 6+6 central mac on the fourth abdominal segment. It differs from the closely related species by the combination of postlabial chaetae and abdominal chaetotaxy. Sinella colorata is a new record from the province. A key to the eyeless Sinella species from China is provided.

Study on the Pauropoda from Tibet, China. Part I. The genera Decapauropus and Hemipauropus (Myriapoda).

Three new species of family Pauropodidae: Decapauropus biconjugarus Qian & Bu, sp. n., D. tibeticus Qian & Bu, sp. n. and Hemipauropus quadrangulus Qian & Bu, sp. n. are described and illustrated from southeastern Tibet, China. The genus Hemipauropus is recorded for the first time from China. This is the second report of pauropods from Tibet.

A new species of Lobellina and first record of Vietnura from China (Collembola: Neanuridae: Neanurinae).

A new species of Lobellina Yosii, 1956 and a key to all species of the genus is provided. It is distinguished from all known members of the genus by its unique set of morphological characters: mandible with six teeth, cephalic chaeta O present, and free from tubercle Fr, cephalic tubercle Oc with three chaetae, cephalic tubercle Di separate, and tubercle Dl with four (sometimes three) chaetae, Ant. I with eight chaetae, and claw with an inner tooth. Vietnuracaerulea Deharveng & Bedos, 2000 is recorded from China for the first time. New localities of Rambutanurahunanensis Jiang & Dong, 2018 and Vitronuraparaacuta Wang, Wang & Jiang, 2016 from southwest China are also provided.

Revisiting the phylogeny of Wolbachia in Collembola.

The endosymbiont Wolbachia has been detected in a few parthenogenetic collembolans sampled in Europe and America, including three species of Poduromorpha, two species of Entomobryomorpha, and two species of Neelipleona. Based on 16S rRNA and ftsZ gene sequences, most of the Wolbachia infecting parthenogenetic collembolans were characterized as members of supergroup E and showed concordant phylogeny with their hosts. However, the two neelipleonan symbionts form another unique group, indicating that Wolbachia has infected parthenogenetic collembolans multiple times. In this study, five parthenogenetic collembolan species were identified as hosts of Wolbachia, and four new Wolbachia strains were reported for four collembolan species sampled in China, respectively, including a neelipleonan strain from Megalothorax incertus (wMinc). Our results demonstrated that the Wolbachia multilocus sequence typing (MLST) system is superior to the 16S rRNA + ftsZ approach for phylogenetic analyses of collembolan Wolbachia. The MLST system assigned these Wolbachia of parthenogenetic collembolans to supergroup E as a unique clade, which included wMinc, supporting the monophyletic origin of Wolbachia in parthenogenetic collembolan species. Moreover, our data suggested supergroup E as one of the most divergent lineages in Wolbachia and revealed the discrepancy between the phylogenies of Wolbachia from parthenogenetic collembolans and their hosts, which may result from the high level of genetic divergence between collembolan Wolbachia, in association with the geographic differentiation of their hosts, or the possible horizontal transmission of Wolbachia between different collembolan species.

Paracerella Imadaté in China: the description of a new species and the analysis of genetic differences between populations (Protura, Acerentomata, Nipponentomidae).

The genus Paracerella Imadaté, 1980 is recorded from China for the first time, with the description of a new species, Paracerella sinensis sp. n. Paracerella sinensis is characterized by four pairs of A-setae on tergite I, the presence of setae Pc and P3a on tergite VII, eight A-setae on tergite VIII, the presence of seta Pc on both sternites VI and VII, and 4/2 setae on sternite VIII, which are different from all other members of the genus. The key to the four species of the genus is updated. In addition, DNA barcodes of four populations are sequenced and their genetic differences are analyzed.

Pax6 in Collembola: Adaptive Evolution of Eye Regression.

Unlike the compound eyes in insects, collembolan eyes are comparatively simple: some species have eyes with different numbers of ocelli (1 + 1 to 8 + 8), and some species have no apparent eye structures. Pax6 is a universal master control gene for eye morphogenesis. In this study, full-length Pax6 cDNAs, Fc-Pax6 and Cd-Pax6, were cloned from an eyeless collembolan (Folsomia candida, soil-dwelling) and an eyed one (Ceratophysella denticulata, surface-dwelling), respectively. Their phylogenetic positions are between the two Pax6 paralogs in insects, eyeless (ey) and twin of eyeless (toy), and their protein sequences are more similar to Ey than to Toy. Both Fc-Pax6 and Cd-Pax6 could induce ectopic eyes in Drosophila, while Fc-Pax6 exhibited much weaker transactivation ability than Cd-Pax6. The C-terminus of collembolan Pax6 is indispensable for its transactivation ability, and determines the differences of transactivation ability between Fc-Pax6 and Cd-Pax6. One of the possible reasons is that Fc-Pax6 accumulated more mutations at some key functional sites of C-terminus under a lower selection pressure on eye development due to the dark habitats of F. candida. The composite data provide a first molecular evidence for the monophyletic origin of collembolan eyes, and indicate the eye degeneration of collembolans is caused by adaptive evolution.