Worldwide revision of synanthropic silverfish (Insecta: Zygentoma: Lepismatidae) combining morphological and molecular data.

Synanthropic silverfish are the best-known and most widely distributed insects of the order Zygentoma. However, there is a great gap in the knowledge and confusion about the geographic distribution and the diagnostic characteristics that allow their identification. In this work, we provide an exhaustive and deep analysis of the most common 9 synanthropic silverfish of the world, combining previously published and newly derived morphological and molecular data. Updated descriptions of Ctenolepisma calvum (Ritter, 1910) and Ctenolepisma (Sceletolepisma) villosum (Fabricius, 1775) are included, and morphological remarks, illustrations, and photographs of the remaining synanthropic species are provided to clarify their diagnosis and differentiation among them and from other free-living species. In addition, Ctenolepisma targionii (Grassi and Rovelli, 1889) is synonymized with C. villosum. A molecular phylogeny is presented based on the COI sequences of all the synanthropic species deposited in BOLD and GenBank, with 15 new sequences provided by this study. This has allowed us to detect and correct a series of identification errors based on the lack of morphological knowledge of several species. Moreover, 2 different lineages of Ctenolepisma longicaudatumEscherich, 1905 have also been detected. To help future studies, we also provide a taxonomic interpretation guide for the most important diagnostic characters of the order Zygentoma, as well as an identification key for all the Synanthropic studied species. Finally, an approximation of the global distribution of synanthropic silverfish is discussed. Several new records indicate that the expansion of these species, generally associated with the transport of goods and people, is still far from over.

A new suite of reporter vectors and a novel landing site survey system to study cis-regulatory elements in diverse insect species.

Comparative analyses between traditional model organisms, such as the fruit fly Drosophila melanogaster, and more recent model organisms, such as the red flour beetle Tribolium castaneum, have provided a wealth of insight into conserved and diverged aspects of gene regulation. While the study of trans-regulatory components is relatively straightforward, the study of cis-regulatory elements (CREs, or enhancers) remains challenging outside of Drosophila. A central component of this challenge has been finding a core promoter suitable for enhancer-reporter assays in diverse insect species. Previously, we demonstrated that a Drosophila Synthetic Core Promoter (DSCP) functions in a cross-species manner in Drosophila and Tribolium. Given the over 300 million years of divergence between the Diptera and Coleoptera, we reasoned that DSCP-based reporter constructs will be useful when studying cis-regulation in a variety of insect models across the holometabola and possibly beyond. To this end, we sought to create a suite of new DSCP-based reporter vectors, leveraging dual compatibility with piggyBac and PhiC31-integration, the 3xP3 universal eye marker, GATEWAY cloning, different colors of reporters and markers, as well as Gal4-UAS binary expression. While all constructs functioned properly with a Tc-nub enhancer in Drosophila, complications arose with tissue-specific Gal4-UAS binary expression in Tribolium. Nevertheless, the functionality of these constructs across multiple holometabolous orders suggests a high potential compatibility with a variety of other insects. In addition, we present the piggyLANDR (piggyBac-LoxP AttP Neutralizable Destination Reporter) platform for the establishment of proper PhiC31 landing sites free from position effects. As a proof-of-principle, we demonstrated the workflow for piggyLANDR in Drosophila. The potential utility of these tools ranges from molecular biology research to pest and disease-vector management, and will help advance the study of gene regulation beyond traditional insect models.

Comprehensive analysis of silk proteins and gland compartments in Limnephilus lunatus, a case-making trichopteran.

Caddisfly larvae produce silk containing heavy and light fibroins, similar to the silk of Lepidoptera, for the construction of underwater structures. We analyzed the silk of Limnephilus lunatus belonging to the case-forming suborder Integripalpia. We analyzed the transcriptome, mapped the transcripts to a reference genome and identified over 80 proteins using proteomic methods, and checked the specificity of their expression. For comparison, we also analyzed the transcriptome and silk proteome of Limnephilus flavicornis. Our results show that fibroins and adhesives are produced together in the middle and posterior parts of the silk glands, while the anterior part produces enzymes and an unknown protein AT24. The number of silk proteins of L. lunatus far exceeds that of the web-spinning Plectrocnemia conspersa, a previously described species from the suborder Annulipalpia. Our results support the idea of increasing the structural complexity of silk in rigid case builders compared to trap web builders.

Evolution repeats itself in replicate long-term studies in the wild.

The extent to which evolution is repeatable remains debated. Here, we study changes over time in the frequency of cryptic color-pattern morphs in 10 replicate long-term field studies of a stick insect, each spanning at least a decade (across 30 years of total data). We find predictable "up-and-down" fluctuations in stripe frequency in all populations, representing repeatable evolutionary dynamics based on standing genetic variation. A field experiment demonstrates that these fluctuations involve negative frequency-dependent natural selection (NFDS). These fluctuations rely on demographic and selective variability that pushes populations away from equilibrium, such that they can reliably move back toward it via NFDS. Last, we show that the origin of new cryptic forms is associated with multiple structural genomic variants such that which mutations arise affects evolution at larger temporal scales. Thus, evolution from existing variation is predictable and repeatable, but mutation adds complexity even for traits evolving deterministically under natural selection.

Measuring the state of aquatic environments using eDNA-upscaling spatial resolution of biotic indices.

Aquatic macroinvertebrates, including many aquatic insect orders, are a diverse and ecologically relevant organismal group yet they are strongly affected by anthropogenic activities. As many of these taxa are highly sensitive to environmental change, they offer a particularly good early warning system for human-induced change, thus leading to their intense monitoring. In aquatic ecosystems there is a plethora of biotic monitoring or biomonitoring approaches, with more than 300 assessment methods reported for freshwater taxa alone. Ultimately, monitoring of aquatic macroinvertebrates is used to calculate ecological indices describing the state of aquatic systems. Many of the methods and indices used are not only hard to compare, but especially difficult to scale in time and space. Novel DNA-based approaches to measure the state and change of aquatic environments now offer unprecedented opportunities, also for possible integration towards commonly applicable indices. Here, we first give a perspective on DNA-based approaches in the monitoring of aquatic organisms, with a focus on aquatic insects, and how to move beyond traditional point-based biotic indices. Second, we demonstrate a proof-of-concept for spatially upscaling ecological indices based on environmental DNA, demonstrating how integration of these novel molecular approaches with hydrological models allows an accurate evaluation at the catchment scale. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Towards holistic insect monitoring: species discovery, description, identification and traits for all insects.

Holistic insect monitoring needs scalable techniques to overcome taxon biases, determine species abundances, and gather functional traits for all species. This requires that we address taxonomic impediments and the paucity of data on abundance, biomass and functional traits. We here outline how these data deficiencies could be addressed at scale. The workflow starts with large-scale barcoding (megabarcoding) of all specimens from mass samples obtained at biomonitoring sites. The barcodes are then used to group the specimens into molecular operational taxonomic units that are subsequently tested/validated as species with a second data source (e.g. morphology). New species are described using barcodes, images and short diagnoses, and abundance data are collected for both new and described species. The specimen images used for species discovery then become the raw material for training artificial intelligence identification algorithms and collecting trait data such as body size, biomass and feeding modes. Additional trait data can be obtained from vouchers by using genomic tools developed by molecular ecologists. Applying this pipeline to a few samples per site will lead to greatly improved insect monitoring regardless of whether the species composition of a sample is determined with images, metabarcoding or megabarcoding. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Three steps towards comparability and standardization among molecular methods for characterizing insect communities.

Molecular methods are currently some of the best-suited technologies for implementation in insect monitoring. However, the field is developing rapidly and lacks agreement on methodology or community standards. To apply DNA-based methods in large-scale monitoring, and to gain insight across commensurate data, we need easy-to-implement standards that improve data comparability. Here, we provide three recommendations for how to improve and harmonize efforts in biodiversity assessment and monitoring via metabarcoding: (i) we should adopt the use of synthetic spike-ins, which will act as positive controls and internal standards; (ii) we should consider using several markers through a multiplex polymerase chain reaction (PCR) approach; and (iii) we should commit to the publication and transparency of all protocol-associated metadata in a standardized fashion. For (i), we provide a ready-to-use recipe for synthetic cytochrome c oxidase spike-ins, which enable between-sample comparisons. For (ii), we propose two gene regions for the implementation of multiplex PCR approaches, thereby achieving a more comprehensive community description. For (iii), we offer guidelines for transparent and unified reporting of field, wet-laboratory and dry-laboratory procedures, as a key to making comparisons between studies. Together, we feel that these three advances will result in joint quality and calibration standards rather than the current laboratory-specific proof of concepts. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

DNA Barcoding and Metabarcoding Protocols for Species Identification.

The article presents the several steps to be performed on a plant, fungal, insect, or soil sample to obtain DNA sequences for DNA barcode analysis. The chapter begins with a description of sample preparation including procedures for cleaning and proceeds to DNA extraction with methods adapted for the specific type of sample. Next, DNA quantification is described so the proper amount is used for the amplification of the selected barcode regions. Information is provided for reaction mixes and amplification conditions for several referenced barcode primer pairs tuned for the individual sample of interest. This is followed by a description of procedures to access the success of amplification, cleanup, and quantification of the product ready for either Sanger sequencing or library preparation for massive parallel sequencing (MPS). Finally, procedures are provided for Sanger sequencing, library preparation, and MPS sequencing. The chapter provides several references of barcode regions for different sample types.

Looking across the gap: Understanding the evolution of eyes and vision among insects.

The compound eyes of insects exhibit stunning variation in size, structure, and function, which has allowed these animals to use their vision to adapt to a huge range of different environments and lifestyles, and evolve complex behaviors. Much of our knowledge of eye development has been learned from Drosophila, while visual adaptations and behaviors are often more striking and better understood from studies of other insects. However, recent studies in Drosophila and other insects, including bees, beetles, and butterflies, have begun to address this gap by revealing the genetic and developmental bases of differences in eye morphology and key new aspects of compound eye structure and function. Furthermore, technical advances have facilitated the generation of high-resolution connectomic data from different insect species that enhances our understanding of visual information processing, and the impact of changes in these processes on the evolution of vision and behavior. Here, we review these recent breakthroughs and propose that future integrated research from the development to function of visual systems within and among insect species represents a great opportunity to understand the remarkable diversification of insect eyes and vision.

Parthenogenetic Stick Insects Exhibit Signatures of Preservation in the Molecular Architecture of Male Reproduction.

After the loss of a trait, theory predicts that the molecular machinery underlying its phenotypic expression should decay. Yet, empirical evidence is contrasting. Here, we test the hypotheses that (i) the molecular ground plan of a lost trait could persist due to pleiotropic effects on other traits and (ii) that gene co-expression network architecture could constrain individual gene expression. Our testing ground has been the Bacillus stick insect species complex, which contains close relatives that are either bisexual or parthenogenetic. After the identification of genes expressed in male reproductive tissues in a bisexual species, we investigated their gene co-expression network structure in two parthenogenetic species. We found that gene co-expression within the male gonads was partially preserved in parthenogens. Furthermore, parthenogens did not show relaxed selection on genes upregulated in male gonads in the bisexual species. As these genes were mostly expressed in female gonads, this preservation could be driven by pleiotropic interactions and an ongoing role in female reproduction. Connectivity within the network also played a key role, with highly connected-and more pleiotropic-genes within male gonad also having a gonad-biased expression in parthenogens. Our findings provide novel insight into the mechanisms which could underlie the production of rare males in parthenogenetic lineages; more generally, they provide an example of the cryptic persistence of a lost trait molecular architecture, driven by gene pleiotropy on other traits and within their co-expression network.

Exploring gene regulation and biological processes in insects: Insights from omics data using gene regulatory network models.

Gene regulatory network (GRN) comprises complicated yet intertwined gene-regulator relationships. Understanding the GRN dynamics will unravel the complexity behind the observed gene expressions. Insect gene regulation is often complicated due to their complex life cycles and diverse ecological adaptations. The main interest of this review is to have an update on the current mathematical modelling methods of GRNs to explain insect science. Several popular GRN architecture models are discussed, together with examples of applications in insect science. In the last part of this review, each model is compared from different aspects, including network scalability, computation complexity, robustness to noise and biological relevancy.

Four phospholipase A2 genes encoded in the western flower thrips genome and their functional differentiation in mediating development and immunity.

Eicosanoids are synthesized from phospholipids by the catalytic activity of phospholipase A2 (PLA2). Even though several PLA2s are encoded in the genome of different insect species, their physiological functions are not clearly discriminated. This study identified four PLA2 genes encoded in the western flower thrips, Frankliniella occidentalis. Two PLA2s (Fo-PLA2C and Fo-PLA2D) are predicted to be secretory while the other two PLA2s (Fo-PLA2A and Fo-PLA2B) are intracellular. All four PLA2 genes were expressed in all developmental stages, of which Fo-PLA2B and Fo-PLA2C were highly expressed in larvae while Fo-PLA2A and Fo-PLA2D were highly expressed in adults. Their expressions in different tissues were also detected by fluorescence in situ hybridization. All four PLA2s were detected in the larval and adult intestines and the ovary. Feeding double-stranded RNAs specific to the PLA2 genes specifically suppressed the target transcript levels. Individual RNA interference (RNAi) treatments led to significant developmental retardation, especially in the treatments specific to Fo-PLA2B and Fo-PLA2D. The RNAi treatments also showed that Fo-PLA2B and Fo-PLA2C expressions were required for the induction of immune-associated genes, while Fo-PLA2A and Fo-PLA2D expressions were required for ovary development. These results suggest that four PLA2s are associated with different physiological processes by their unique catalytic activities and expression patterns.

RNAi efficiency is enhanced through knockdown of double-stranded RNA-degrading enzymes in butterfly Papilio xuthus.

The efficiency of RNA interference (RNAi) has always limited the research on the phenotype innovation of Lepidoptera insects. Previous studies have found that double-stranded RNA-degrading enzyme (dsRNase) is an important factor in RNAi efficiency, but there have been no relevant reports in butterflies (Papilionoidea). Papilio xuthus is one of the important models in butterflies with an extensive experimental application value. To explore the effect of dsRNase in the RNAi efficiency on butterflies, six dsRNase genes (PxdsRNase 1-6) were identified in P. xuthus genome, and their dsRNA-degrading activities were subsequently detected by ex vivo assays. The result shows that the dsRNA-degrading ability of gut content (<1 h) was higher than hemolymph content (>12 h). We then investigated the expression patterns of these PxdsRNase genes during different tissues and developmental stages, and related RNAi experiments were carried out. Our results show that different PxdsRNase genes had different expression levels at different developmental stages and tissues. The expression of PxdsRNase2, PxdsRNase3, and PxdsRNase6 were upregulated significantly through dsGFP injection, and PxdsRNase genes can be silenced effectively by injecting their corresponding dsRNA. RNAi-of-RNAi studies with PxEbony, which acts as a reporter gene, observed that silencing PxdsRNase genes can increase RNAi efficiency significantly. These results confirm that silencing dsRNase genes can improve RNAi efficiency in P. xuthus significantly, providing a reference for the functional study of insects such as butterflies with low RNAi efficiency.

Aberrant splicing of a nicotinic acetylcholine receptor alpha 6 subunit is associated with spinosad tolerance in the thrips predator Orius laevigatus.

Susceptibility to insecticides is one of the limiting factors preventing wider adoption of natural enemies to control insect pest populations. Identification and selective breeding of insecticide tolerant strains of commercially used biological control agents (BCAs) is one of the approaches to overcome this constraint. Although a number of beneficial insects have been selected for increased tolerance to insecticides the molecular mechanisms underpinning these shifts in tolerance are not well characterised. Here we investigated the molecular mechanisms of enhanced tolerance of a lab selected strain of Orius laevigatus (Fieber) to the commonly used biopesticide spinosad. Transcriptomic analysis showed that spinosad tolerance is not a result of overexpressed detoxification genes. Molecular analysis of the target site for spinosyns, the nicotinic acetylcholine receptor (nAChR), revealed increased expression of truncated transcripts of the nAChR α6 subunit in the spinosad selected strain, a mechanism of resistance which was described previously in insect pest species. Collectively, our results demonstrate the mechanisms by which some beneficial biological control agents can evolve insecticide tolerance and will inform the development and deployment of insecticide-tolerant natural enemies in integrated pest management strategies.

Mixed insect pest populations of Diaspididae species under control of oligonucleotide insecticides: 3'-end nucleotide matters.

Diaspididae are one of the most serious small herbivorous insects with piercing-sucking mouth parts and are major economic pests as they attack and destroy perennial ornamentals and food crops. Chemical control is the primary management approach for armored scale infestation. However, chemical insecticides do not possess selectivity in action and not always effective enough for the control of armored scale insects. Our previous work showed that green oligonucleotide insecticides (olinscides) are highly effective against armored and soft scale insects. Moreover, olinscides possess affordability, selectivity in action, fast biodegradability, and a low carbon footprint. Insect pest populations undergo microevolution and olinscides should take into account the problem of insecticide resistance. Using sequencing results, it was found that in the mixed populations of insect pests Dynaspidiotus britannicus Newstead and Aonidia lauri Bouche, predominates the population of A. lauri. Individuals of A. lauri comprised for 80% of individuals with the sequence 3'-ATC-GTT-GGC-AT-5' in the 28S rRNA site, and 20% of the population comprised D. britannicus individuals with the sequence 3'-ATC-GTC-GGT-AT-5'. We created olinscides Diasp80-11 (5'-ATG-CCA-ACG-AT-3') and Diasp20-11 (5'-ATA-CCG-ACG-AT-3') with perfect complementarity to each of the sequences. Mortality of insects on the 14th day comprised 98.19 ± 3.12% in Diasp80-11 group, 64.66 ± 0.67% in Diasp20-11 group (p < 0.05), and 3.77 ± 0.94% in the control group. Results indicate that for maximum insecticidal effect it is necessary to use an oligonucleotide insecticide that corresponds to the dominant species. Mortality in Diasp80-11 group was accompanied with significant decrease in target 28S rRNA concentration and was 8.44 ± 0.14 and 1.72 ± 0.36 times lower in comparison with control (p < 0.05) on the 10th and 14th days, respectively. We decided to make single nucleotide substitutions in Diasp20-11 olinscide to understand which nucleotide will play the most important role in insecticidal effect. We created three sequences with single nucleotide transversion substitutions at the 5'-end - Diasp20(5')-11 (A to T), 3'-end - Diasp20(3')-11 (T to A), and in the middle of the sequence - Diasp20(6)-11 (6th nitrogenous base of the sequence; G to C), respectively. As a result, mortality of mixed population of the field experiment decreased and comprised 53.89 ± 7.25% in Diasp20(5')-11 group, 40.68 ± 4.33% in Diasp20(6)-11 group, 35.74 ± 5.51% in Diasp20(3')-11 group, and 3.77 ± 0.94% in the control group on the 14th day. Thus, complementarity of the 3'-end nucleotide to target 28S rRNA was the most important for pronounced insecticidal effect (significance of complementarity of nucleotides for insecticidal effect: 5' nt < 6 nt < 3' nt). As was found in our previous research works, the most important rule to obtain maximum insecticidal effect is complete complementarity to the target rRNA sequence and maximum coverage of target sequence in insect pest populations. However, in this article we also show that the complementarity of 3'-end is a second important factor for insecticidal potential of olinscides. Also in this article we propose 2-step DNA containment mechanism of action of olinscides, recruiting RNase H. The data obtained indicate the selectivity of olinscides and at the same time provide a simple and flexible platform for the creation of effective plant protection products, based on antisense DNA oligonucleotides.

CRISPR-Cas9 mediated dsRNase knockout improves RNAi efficiency in the fall armyworm.

Lepidopteran insects are refractory to RNA interference (RNAi) response, especially to orally delivered double-stranded RNA (dsRNA). High nuclease activity in the midgut lumen is proposed as one of the major reasons for RNAi insensitivity. We identified three dsRNase genes highly expressed in the midgut of fall armyworm (FAW), Spodoptera frugiperda. The genomic region harboring those three dsRNase genes was deleted using the CRISPR-Cas9-mediated genome editing method. A homozygous line with deletion of three dsRNase genes was produced. dsRNA degradation by midgut lumen contents of mutant larvae was lower than in wild-type larvae. Feeding dsRNA targeting the inhibitor of apoptosis (IAP) gene increased knockdown of the target gene and mortality in mutants compared to wild-type larvae. These results suggest that dsRNases in the midgut contribute to RNAi inefficiency in FAW. Formulations that protect dsRNA from dsRNase degradation may improve RNAi efficiency in FAW and other lepidopteran insects.

Molecular identification of Hymenopteran insects collected by using Malaise traps from Hazarganji Chiltan National Park Quetta, Pakistan.

The order Hymenoptera holds great significance for humans, particularly in tropical and subtropical regions, due to its role as a pollinator of wild and cultivated flowering plants, parasites of destructive insects and honey producers. Despite this importance, limited attention has been given to the genetic diversity and molecular identification of Hymenopteran insects in most protected areas. This study provides insights into the first DNA barcode of Hymenopteran insects collected from Hazarganji Chiltan National Park (HCNP) and contributes to the global reference library of DNA barcodes. A total of 784 insect specimens were collected using Malaise traps, out of which 538 (68.62%) specimens were morphologically identified as Hymenopteran insects. The highest abundance of species of Hymenoptera (133/538, 24.72%) was observed during August and least in November (16/538, 2.97%). Genomic DNA extraction was performed individually from 90/538 (16.73%) morphologically identified specimens using the standard phenol-chloroform method, which were subjected separately to the PCR for their molecular confirmation via the amplification of cytochrome c oxidase subunit 1 (cox1) gene. The BLAST analyses of obtained sequences showed 91.64% to 100% identities with related sequences and clustered phylogenetically with their corresponding sequences that were reported from Australia, Bulgaria, Canada, Finland, Germany, India, Israel, and Pakistan. Additionally, total of 13 barcode index numbers (BINs) were assigned by Barcode of Life Data Systems (BOLD), out of which 12 were un-unique and one was unique (BOLD: AEU1239) which was assigned for Anthidium punctatum. This indicates the potential geographical variation of Hymenopteran population in HCNP. Further comprehensive studies are needed to molecularly confirm the existing insect species in HCNP and evaluate their impacts on the environment, both as beneficial (for example, pollination, honey producers and natural enemies) and detrimental (for example, venomous stings, crop damage, and pathogens transmission).

Mitochondrial genomes of Nemourinae species (Plecoptera: Nemouridae) and the phylogenetic implications.

Currently, the classification system of 2 subfamilies within Nemouridae has been widely accepted. However, monophyly of 2 subfamilies has not been well supported by molecular evidence. To date, only mitogenomes from genus Nemoura of the subfamily Nemourinae were used in previous phylogenetic studies and produced conflicting results with morphological studies. Herein, we analyzed mitogenomes of 3 Nemourinae species to reveal their mitogenomic characteristics and to examine genus-level classification among Nemouridae. In this study, the genome organization of 3 mitogenomes is highly conserved in gene order, nucleotide composition, codon usage, and amino acid composition. In 3 Nemourinae species, there is a high variation in nucleotide diversity among the 13 protein-coding genes (PCGs). The Ka/Ks values for all PCGs were far lower than 1, indicating that these genes were evolving under purifying selection. The phylogenetic analyses highly support Nemurella as the sister group to Ostrocerca. Meanwhile, Nemoura is recovered as the sister group of Malenka; they are grouped with other Amphinemurinae and emerged from a paraphyletic Nemourinae. More molecular data from different taxonomic groups are needed to understand stoneflies phylogeny and evolution.

The first A-to-I RNA editome of hemipteran species Coridius chinensis reveals overrepresented recoding and prevalent intron editing in early-diverging insects.

BACKGROUND: Metazoan adenosine-to-inosine (A-to-I) RNA editing resembles A-to-G mutation and increases proteomic diversity in a temporal-spatial manner, allowing organisms adapting to changeable environment. The RNA editomes in many major animal clades remain unexplored, hampering the understanding on the evolution and adaptation of this essential post-transcriptional modification. METHODS: We assembled the chromosome-level genome of Coridius chinensis belonging to Hemiptera, the fifth largest insect order where RNA editing has not been studied yet. We generated ten head RNA-Seq libraries with DNA-Seq from the matched individuals. RESULTS: We identified thousands of high-confidence RNA editing sites in C. chinensis. Overrepresentation of nonsynonymous editing was observed, but conserved recoding across different orders was very rare. Under cold stress, the global editing efficiency was down-regulated and the general transcriptional processes were shut down. Nevertheless, we found an interesting site with "conserved editing but non-conserved recoding" in potassium channel Shab which was significantly up-regulated in cold, serving as a candidate functional site in response to temperature stress. CONCLUSIONS: RNA editing in C. chinensis largely recodes the proteome. The first RNA editome in Hemiptera indicates independent origin of beneficial recoding during insect evolution, which advances our understanding on the evolution, conservation, and adaptation of RNA editing.

Divergent dynamics of sexual and habitat isolation at the transition between stick insect populations and species.

Speciation is often viewed as a continuum along which populations diverge until they become reproductively-isolated species. However, such divergence may be heterogeneous, proceeding in fits and bursts, rather than being uniform and gradual. We show in Timema stick insects that one component of reproductive isolation evolves non-uniformly across this continuum, whereas another does not. Specifically, we use thousands of host-preference and mating trials to study habitat and sexual isolation among 42 pairs of taxa spanning a range of genomic differentiation and divergence time. We find that habitat isolation is uncoupled from genomic differentiation within species, but accumulates linearly with it between species. In contrast, sexual isolation accumulates linearly across the speciation continuum, and thus exhibits similar dynamics to morphological traits not implicated in reproductive isolation. The results show different evolutionary dynamics for different components of reproductive isolation and highlight a special relevance for species status in the process of speciation.