DNA methylation differences between stick insect ecotypes

Epigenetic mechanisms, such as DNA methylation, can influence gene regulation and affect phenotypic variation, raising the possibility that they contribute to ecological adaptation. Beginning to address this issue requires high-resolution sequencing studies of natural populations to pinpoint epigenetic regions of potential ecological and evolutionary significance. However, such studies are still relatively uncommon, especially in insects, and are mainly restricted to a few model organisms. Here, we characterize patterns of DNA methylation for natural populations of Timema cristinae adapted to two host plant species (i.e. ecotypes). By integrating results from sequencing of whole transcriptomes, genomes and methylomes, we investigate whether environmental, host and genetic differences of these stick insects are associated with methylation levels of cytosine nucleotides in the CpG context. We report an overall genome-wide methylation level for T. cristinae of ~14%, with methylation being enriched in gene bodies and impoverished in repetitive elements. Genome-wide DNA methylation variation was strongly positively correlated with genetic distance (relatedness), but also exhibited significant host-plant effects. Using methylome-environment association analysis, we pinpointed specific genomic regions that are differentially methylated between ecotypes, with these regions being enriched for genes with functions in membrane processes. The observed association between methylation variation and genetic relatedness, and with the ecologically important variable of host plant, suggests a potential role for epigenetic modification in T. cristinae adaptation. To substantiate such adaptive significance, future studies could test whether methylation can be transmitted across generations and the extent to which it responds to experimental manipulation in field and laboratory studies.

Comparative transcriptome analysis to unveil genes affecting the host cuticle destruction in Metarhizium rileyi

Insect pathogenic fungi, also known as entomopathogenic fungi, are one of the largest insect pathogenic microorganism communities, represented by Beauveria spp. and Metarhizium spp. Entomopathogenic fungi have been proved to be a great substitute for chemical pesticide in agriculture. In fact, a lot of functional genes were also already characterized in entomopathogenic fungi, but more depth of exploration is still needed to reveal their complicated pathogenic mechanism to insects. Metarhizium rileyi (Nomuraea rileyi) is a great potential biocontrol fungus that can parasitize more than 40 distinct species (mainly Lepidoptera: Noctuidae) to cause large-scale infectious diseases within insect population. In this study, a comparative analysis of transcriptome profile was performed with topical inoculation and hemolymph injection to character the infectious pattern of M. rileyi. Appressorium and multiple hydrolases are indispensable constituents to break the insect host primary cuticle defense in entomopathogenic fungi. Within our transcriptome data, numerous transcripts related to destruction of insect cuticle rather growth regulations were obtained. Most importantly, some unreported ribosomal protein genes and novel unannotated protein (hypothetical protein) genes were proved to participate in the course of pathogenic regulation. Our current data provide a higher efficiency gene library for virulence factors screen in M. rileyi, and this library may be also useful for furnishing valuable information on entomopathogenic fungal pathogenic mechanisms to host.

The loci of insect phenotypic evolution

Insects are important elements of terrestrial ecosystems because they pollinate plants, destroy crops, transmit diseases to livestock and humans, and are important components of food chains. Here, I used Gephebase, a manually curated database of genetic variants associated with natural and domesticated trait variation, to explore current knowledge about the genes and the mutations known to contribute to natural phenotypic variation in insects. Analysis of over 600 mutations reveals that data are concentrated toward certain species and traits and that experimental approaches have changed over time. The distribution of coding and cis-regulatory changes varies with traits, experimental approaches, and identified gene loci. Recent studies highlight the important role of standing variation, repeated mutations in hotspot genes, recombination, inversions, and introgression.

Complete metamorphosis and microbiota turnover in insects

The insects constitute the majority of animal diversity. Most insects are holometabolous: during complete metamorphosis their bodies are radically reorganized. This reorganization poses a significant challenge to the gut microbiota, as the gut is replaced during pupation, a process that does not occur in hemimetabolous insects. In holometabolous hosts, it offers the opportunity to decouple the gut microbiota between the larval and adult life stages resulting in high beta diversity whilst limiting alpha diversity. Here, we studied 18 different herbivorous insect species from five orders of holometabolous and three orders of hemimetabolous insects. Comparing larval and adult specimens, we find a much higher beta-diversity and hence microbiota turnover in holometabolous insects compared to hemimetabolous insects. Alpha diversity did not differ between holo- and hemimetabolous insects nor between developmental stages within these groups. Our results support the idea that pupation offers the opportunity to change the gut microbiota and hence might facilitate ecological niche shifts. This possible effect of niche shift facilitation could explain a selective advantage of the evolution of complete metamorphosis, which is a defining trait of the most speciose insect taxon, the holometabola.

Efficient nanoparticle-based CRISPR-Cas13d induced mRNA disruption of an eye pigmentation gene in the white-backed planthopper, Sogatella furcifera

The discovery of the clustered regularly interspaced short palindromic repeat (CRISPR) system has driven gene manipulation technology to a new era with applications reported in organisms that span the tree of life. The utility of CRISPR-mediated editing was further expanded to mRNA following identification of the RNA-targeting Cas13 family of smaller endonuclease proteins. Application of this family to insect research, however, has been more limited. In this study, the smallest Cas13 family member, Cas13d, and guide RNAs (gRNAs) were complexed with a versatile nanomaterial (star polycation, SPc) to generate a proof-of-concept RNA-editing platform capable of disrupting mRNA expression of the eye pigmentation gene tryptophan 2,3-dioxygenase (SfTO) in white-backed planthoppers (WBPHs). The resulting red-eye phenotype was present in 19.76% (with SPc) and 22.99% (without SPc) of the treatment groups and was comparable to the red-eye phenotype generated following conventional RNA interference knockdown (22.22%). Furthermore, the Cas13/gRNA phenotype manifested more quickly than RNA interference. Consistent with the expected Cas13d mechanism, SfTO transcript levels were significantly reduced. Taken together, the results indicate that the SPc-CRISPR-Cas13d/gRNA complex negatively impacted expression of the target gene. These findings confirm the utility of this novel mRNA disruption system in insects and lay the foundation for further development of these tools in the implementation of green agricultural pest management tactics.

Pollen DNA metabarcoding reveals cryptic diversity and high spatial turnover in alpine plant-pollinator networks

Alpine plant-pollinator communities play an important role in the functioning of alpine ecosystems, which are highly threatened by climate change. However, we still have a poor understanding of how environmental factors and spatiotemporal variability shape these communities. Here, we investigate what drives structure and beta diversity in a plant-pollinator metacommunity from the Australian alpine region using two approaches: pollen DNA metabarcoding (MB) and observations. Individual pollinators often carry pollen from multiple plant species, and therefore we expected MB to reveal a more diverse and complex network structure. We used two gene regions (ITS2 and trnL) to identify plant species present in the pollen loads of 154 insect pollinator specimens from three alpine habitats and construct MB networks, and compared them to networks based on observations alone. We compared species and interaction turnover across space for both types of networks, and evaluated their differences for plant phylogenetic diversity and beta diversity. We found significant structural differences between the two types of networks; notably, MB networks were much less specialized but more diverse than observation networks, with MB detecting many cryptic plant species. Both approaches revealed that alpine pollination networks are very generalized, but we estimated a high spatial turnover of plant species (0.79) and interaction rewiring (0.6) as well as high plant phylogenetic diversity (0.68) driven by habitat differences based on the larger diversity of plant species and species interactions detected with MB. Overall, our findings show that habitat and microclimatic heterogeneity drives diversity and fine-scale spatial turnover of alpine plant-pollinator networks.

Opsin Gene Duplication in Lepidoptera: Retrotransposition, Sex Linkage, and Gene Expression

Color vision in insects is determined by signaling cascades, central to which are opsin proteins, resulting in sensitivity to light at different wavelengths. In certain insect groups, lineage-specific evolution of opsin genes, in terms of copy number, shifts in expression patterns, and functional amino acid substitutions, has resulted in changes in color vision with subsequent behavioral and niche adaptations. Lepidoptera are a fascinating model to address whether evolutionary change in opsin content and sequence evolution are associated with changes in vision phenotype. Until recently, the lack of high-quality genome data representing broad sampling across the lepidopteran phylogeny has greatly limited our ability to accurately address this question. Here, we annotate opsin genes in 219 lepidopteran genomes representing 33 families, reconstruct their evolutionary history, and analyze shifts in selective pressures and expression between genes and species. We discover 44 duplication events in opsin genes across ∼300 million years of lepidopteran evolution. While many duplication events are species or family specific, we find retention of an ancient long-wavelength-sensitive (LW) opsin duplication derived by retrotransposition within the speciose superfamily Noctuoidea (in the families Nolidae, Erebidae, and Noctuidae). This conserved LW retrogene shows life stage-specific expression suggesting visual sensitivities or other sensory functions specific to the early larval stage. This study provides a comprehensive order-wide view of opsin evolution across Lepidoptera, showcasing high rates of opsin duplications and changes in expression patterns.

The X chromosome of insects likely predates the origin of class Insecta

Sex chromosomes have evolved independently multiple times, but why some are conserved for more than 100 million years whereas others turnover rapidly remains an open question. Here, we examine the homology of sex chromosomes across nine orders of insects, plus the outgroup springtails. We find that the X chromosome is likely homologous across insects and springtails; the only exception is in the Lepidoptera, which has lost the X and now has a ZZ/ZW sex-chromosome system. These results suggest the ancestral insect X chromosome has persisted for more than 450 million years-the oldest known sex chromosome to date. Further, we propose that the shrinking of gene content the dipteran X chromosome has allowed for a burst of sex-chromosome turnover that is absent from other speciose insect orders.

Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex

As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.

Knockdown of double-stranded RNases (dsRNases) enhances oral RNA interference (RNAi) in the corn leafhopper, Dalbulus maidis

The leafhopper Dalbulus maidis is a harmful pest that causes severe damage to corn crops. Conventional chemical pesticides have negative environmental impacts, emphasizing the need for alternative solutions. RNA interference (RNAi) is a more specific and environmentally friendly method for controlling pests and reducing the negative impacts of current pest management practices. Previous studies have shown that orally administered double-stranded RNA (dsRNA) is less effective than injection protocols in silencing genes. This study focuses on identifying and understanding the role of double-stranded ribonucleases (dsRNases) in limiting the efficiency of oral RNAi in D. maidis. Three dsRNases were identified and characterized, with Dmai-dsRNase-2 being highly expressed in the midgut and salivary glands. An ex vivo degradation assay revealed significant nuclease activity, resulting in high instability of dsRNA when exposed to tissue homogenates. Silencing Dmai-dsRNase-2 improved the insects' response to the dsRNA targeting the gene of interest, providing evidence of dsRNases involvement in oral RNAi efficiency. Therefore, administering both dsRNase-specific and target gene-specific-dsRNAs simultaneously is a promising approach to increase the efficiency of oral RNAi and should be considered in future control strategies.

Untangling an insect's virome from its endogenous viral elements

BACKGROUND

Insects are an important reservoir of viral biodiversity, but the vast majority of viruses associated with insects have not been discovered. Recent studies have employed high-throughput RNA sequencing, which has led to rapid advances in our understanding of insect viral diversity. However, insect genomes frequently contain transcribed endogenous viral elements (EVEs) with significant homology to exogenous viruses, complicating the use of RNAseq for viral discovery.

METHODS

In this study, we used a multi-pronged sequencing approach to study the virome of an important agricultural pest and prolific vector of plant pathogens, the potato aphid Macrosiphum euphorbiae. We first used rRNA-depleted RNAseq to characterize the microbes found in individual insects. We then used PCR screening to measure the frequency of two heritable viruses in a local aphid population. Lastly, we generated a quality draft genome assembly for M. euphorbiae using Illumina-corrected Nanopore sequencing to identify transcriptionally active EVEs in the host genome.

RESULTS

We found reads from two insect-specific viruses (a Flavivirus and an Ambidensovirus) in our RNAseq data, as well as a parasitoid virus (Bracovirus), a plant pathogenic virus (Tombusvirus), and two phages (Acinetobacter and APSE). However, our genome assembly showed that part of the 'virome' of this insect can be attributed to EVEs in the host genome.

CONCLUSION

Our work shows that EVEs have led to the misidentification of aphid viruses from RNAseq data, and we argue that this is a widespread challenge for the study of viral diversity in insects.

Neofunctionalization of Toll Signaling in Insects: From Immunity to Dorsoventral Patterning

Toll signaling plays a crucial role in pathogen defense throughout the animal kingdom. It was discovered, however, for its function in dorsoventral (DV) axis formation in Drosophila. In all other insects studied so far, but not outside the insects, Toll is also required for DV patterning. However, in insects more distantly related to Drosophila, Toll's patterning role is frequently reduced and substituted by an expanded influence of BMP signaling, the pathway implicated in DV axis formation in all major metazoan lineages. This suggests that Toll was integrated into an ancestral BMP-based patterning system at the base of the insects or during insect evolution. The observation that Toll signaling has an immune function in the extraembryonic serosa, an early differentiating tissue of most insect embryos, suggests a scenario of how Toll was co-opted from an ancestral immune function for its new role in axis formation.

Phylogeographical evidence for historical long-distance dispersal in the flightless stick insect Ramulus mikado

Exploring how organisms overcome geographical barriers to dispersal is a fundamental question in biology. Passive long-distance dispersal events, although infrequent and unpredictable, have a considerable impact on species range expansions. Despite limited active dispersal capabilities, many stick insect species have vast geographical ranges, indicating that passive long-distance dispersal is vital for their distribution. A potential mode of passive dispersal in stick insects is via the egg stage within avian digestive tracts, as suggested by experimental evidence. However, detecting such events under natural conditions is challenging due to their rarity. Therefore, to indirectly assess the potential of historical avian-mediated dispersal, we examined the population genetic structure of the flightless stick insect Ramulus mikado across Japan, based on a multifaceted molecular approach [cytochrome oxidase subunit I (COI) haplotypes, nuclear simple sequence repeat markers and genome-wide single nucleotide polymorphisms]. Subsequently, we identified unique phylogeographic patterns, including the discovery of identical COI genotypes spanning considerable distances, which substantiates the notion of passive long-distance genotypic dispersal. Overall, all the molecular data revealed the low and mostly non-significant genetic differentiation among populations, with identical or very similar genotypes across distant populations. We propose that long-distance dispersal facilitated by birds is the plausible explanation for the unique phylogeographic pattern observed in this flightless stick insect.

ALKBH8 as a potential N6 -methyladenosine (m6 A) eraser in insects

The N6 -methyladenosine (m6 A) machinery functions through three groups of proteins in eukaryotic cells, including m6 A writers, erasers and readers. The m6 A cellular machinery has mostly been characterised in mammalian species, and the relevant literature on insects is currently scant. While homologues of m6 A writers and readers have been reported from insects, no erasers have been described so far. Here, using BLAST search, we searched for potential erasers in insects. While we found homologues of human m6 A eraser ALKBH5 in termites, beetles and true bugs, they could not be found in representative dipteran and lepidopteran species. However, a potential m6 A eraser, ALKBH8, was identified and experimentally investigated. Our results showed that ALKBH8 can reduce the m6 A levels of Aedes aegypti and Drosophila melanogaster RNAs, suggesting that AeALKBH8 could be a candidate m6 A eraser in insects.

Launching insectphylo.org; a new hub facilitating construction and use of synthesis molecular phylogenies of insects

The Holy Grail of an Insect Tree of Life can only be 'discovered' through extensive collaboration among taxon specialists, phylogeneticists and centralized frameworks such as Open Tree of Life, but insufficient effort from stakeholders has so far hampered this promising approach. The resultant unavailability of synthesis phylogenies is an unfortunate situation given the numerous practical usages of phylogenies in the near term and against the backdrop of the ongoing biodiversity crisis. To resolve this issue, we establish a new online hub that centralizes the collation of relevant phylogenetic data and provides the resultant synthesis molecular phylogenies. This is achieved through key developments in a proposed pipeline for the construction of a species-level insect phylogeny. The functionality of the framework is demonstrated through the construction of a highly supported, species-comprehensive phylogeny of Diptera, built from integrated omics data, COI DNA barcodes, and a compiled database of over 100 standardized, published Diptera phylogenies. Machine-readable forms of the phylogeny (and subsets thereof) are publicly available at insectphylo.org, a new public repository for species-comprehensive phylogenies for biological research.

ICTV Virus Taxonomy Profile: Qinviridae 2023

Qinviridae is a family of negative-sense RNA viruses with genomes of 7.3-8.2 kb that have been associated with crustaceans, insects, gastropods, and nematodes. The qinvirid genome consists of two segments, each with at least one open reading frame (ORF). The large (L) segment ORF encodes a large protein containing an RNA-directed RNA polymerase (RdRP) domain. The small (S) segment ORF encodes a nucleocapsid protein. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Qinviridae, which is available at ictv.global/report/qinviridae.

Tiny insects, big troubles: a review of BOLD's COI database for Thysanoptera (Insecta)

DNA Barcoding is an important tool for disciplines such as taxonomy, phylogenetics and phylogeography, with Barcode of Life Data System (BOLD) being the largest database of partial cytochrome c oxidase subunit I (COI) sequences. We provide the first extensive revision of the information available in this database for the insect order Thysanoptera, to assess: how many COI sequences are available; how representative these sequences are for the order; and the current potential of BOLD as a reference library for specimen identification and species delimitation. The COI database at BOLD currently represents only about 5% of the over 6400 valid thrips species, with a heavy bias towards a few species of economic importance. Clear Barcode gaps were observed for 24 out of 33 genera evaluated, but many outliers were also observed. We suggest that the COI sequences available in BOLD as a reference would not allow for accurate identifications in about 30% of Thysanoptera species in this database, which rises to 40% of taxa within Thripidae, the most sampled family within the order. Thus, we call for caution and a critical evaluation in using BOLD as a reference library for thrips Barcodes, and future efforts should focus on improving the data quality of this database.

Inefficient uptake of small interfering RNAs is responsible for their inability to trigger RNA interference in Colorado potato beetle cells

There has been limited success in the usage of exogenous small interference RNA (siRNA) or small hairpin RNA (shRNA) to trigger RNA interference (RNAi) in insects. Instead, long double-stranded RNAs (dsRNA) are used to induce knockdown of target genes in insects. Here, we compared the potency of si/sh RNAs and dsRNA in Colorado potato beetle (CPB) cells. CPB cells showed highly efficient RNAi response to dsRNA. However, si/sh RNAs were inefficient in triggering RNAi in CPB cells. Confocal microscopy observations of Cy3 labeled-si/sh RNA cellular uptake revealed reduced si/sh RNA uptake compared to dsRNA. si/sh RNAs were stable in the conditioned media of CPB cells. Although in a small amount, when internalized by CPB cells, the si/sh RNAs were processed by the Dicer enzyme. Lipid-mediated transfection and chimeric dsRNA approaches were used to improve the delivery of si/sh RNAs. Our results suggest that the uptake of si/sh RNAs is inefficient in CPB cells, resulting in ineffective RNAi response. However, with the help of effective delivery methods, si/sh RNA could be a useful option for developing target-specific RNAi-mediated biopesticides.

Reduced-representation libraries in insect genetics

Genotyping-by-sequencing of reduced-representation libraries has ushered in an era where genome-wide data can be gotten for any species. Here, I review research on this topic during the last two years, report meta-analysis of the results, and discuss analysis methods and issues. Scanning the literature from 2021 to 2022 identified 21 papers, the majority of which were on population differences, including local adaptation and migration, but several papers were on genetic maps and their use in assembly scaffolding or analysis of quantitative trait loci, on the origin of incursions of pest insects, or on infection rates of a pathogen in a disease vector. The research reviewed includes 33 species from 25 families and 11 orders. Meta-analysis showed that less than 16%, and most often, less than 1% of the genome was implicated in local adaptation and that the number of adaptive loci correlated with genetic divergence among populations.

Autorecoding A-to-I RNA editing sites in the Adar gene underwent compensatory gains and losses in major insect clades

As one of the most prevalent RNA modifications in animals, adenosine-to-inosine (A-to-I) RNA editing facilitates the environmental adaptation of organisms by diversifying the proteome in a temporal-spatial manner. In flies and bees, the editing enzyme Adar has independently gained two different autorecoding sites that form an autofeedback loop, stabilizing the overall editing efficiency. This ensures cellular homeostasis by keeping the normal function of target genes. However, in a broader range of insects, the evolutionary dynamics and significance of this Adar autoregulatory mechanism are unclear. We retrieved the genomes of 377 arthropod species covering the five major insect orders (Hemiptera, Hymenoptera, Coleoptera, Diptera, and Lepidoptera) and aligned the Adar autorecoding sites across all genomes. We found that the two autorecoding sites underwent compensatory gains and losses during the evolution of two orders with the most sequenced species (Diptera and Hymenoptera), and that the two editing sites were mutually exclusive among them: One editable site is significantly linked to another uneditable site. This autorecoding mechanism of Adar could flexibly diversify the proteome and stabilize global editing activity. Many insects independently selected different autorecoding sites to achieve a feedback loop and regulate the global RNA editome, revealing an interesting phenomenon during evolution. Our study reveals the evolutionary force acting on accurate regulation of RNA editing activity in insects and thus deepens our understanding of the functional importance of RNA editing in environmental adaptation and evolution.

Diversity of transgenes in sustainable management of insect pests

Insecticidal transgenes, when incorporated and expressed in plants, confer resistance against insects by producing several products having insecticidal properties. Protease inhibitors, lectins, amylase inhibitors, and chitinase genes are associated with the natural defenses developed by plants to counter insect attacks. Several toxin genes are also derived from spiders and scorpions for protection against insects. Bacillus thuringiensis Berliner is a microbial source of insecticidal toxins. Several methods have facilitated the large-scale production of transgenic plants. Bt-derived cry, cyt, vip, and sip genes, plant-derived genes such as lectins, protease inhibitors, and alpha-amylase inhibitors, insect cell wall-degrading enzymes like chitinase and some proteins like arcelins, plant defensins, and ribosome-inactivating proteins have been successfully utilized to impart resistance to insects. Besides, transgenic plants expressing double-stranded RNA have been developed with enhanced resistance. However, the long-term effects of transgenes on insect resistance, the environment, and human health must be thoroughly investigated before they are made available for commercial planting. In this chapter, the present status, prospects, and future scope of transgenes for insect pest management have been summarized and discussed.

Characterization of the fragmented mitochondrial genome of domestic pig louse Haematopinus suis (Insecta: Haematopinidae) from China

The domestic pig louse Haematopinus suis (Linnaeus, 1758) (Phthiraptera: Anoplura) is a common ectoparasite of domestic pigs, which can act as a vector of various infectious disease agents. Despite its significance, the molecular genetics, biology and systematics of H. suis from China have not been studied in detail. In the present study, the entire mitochondrial (mt) genome of H. suis isolate from China was sequenced and compared with that of H. suis isolate from Australia. We identified 37 mt genes located on nine circular mt minichromosomes, 2.9 kb-4.2 kb in size, each containing 2-8 genes and one large non-coding region (NCR) (1,957 bp-2,226 bp). The number of minichromosomes, gene content, and gene order in H. suis isolates from China and Australia are identical. Total sequence identity across coding regions was 96.3% between H. suis isolates from China and Australia. For the 13 protein-coding genes, sequence differences ranged from 2.8%-6.5% consistent nucleotides with amino acids. Our result is H. suis isolates from China and Australia being the same H. suis species. The present study determined the entire mt genome of H. suis from China, providing additional genetic markers for studying the molecular genetics, biology and systematics of domestic pig louse.

Evidence for cryptic sex in parthenogenetic stick insects of the genus Timema

Obligately parthenogenetic species are expected to be short lived since the lack of sex and recombination should translate into a slower adaptation rate and increased accumulation of deleterious alleles. Some, however, are thought to have been reproducing without males for millions of years. It is not clear how these old parthenogens can escape the predicted long-term costs of parthenogenesis, but an obvious explanation is cryptic sex. In this study, we screen for signatures of cryptic sex in eight populations of four parthenogenetic species of Timema stick insects, some estimated to be older than 1 Myr. Low genotype diversity, homozygosity of individuals and high linkage disequilibrium (LD) unaffected by marker distances support exclusively parthenogenetic reproduction in six populations. However, in two populations (namely, of the species Timema douglasi and T. monikensis) we find strong evidence for cryptic sex, most likely mediated by rare males. These populations had comparatively high genotype diversities, lower LD, and a clear LD decay with genetic distance. Rare sex in species that are otherwise largely parthenogenetic could help explain the unusual success of parthenogenesis in the Timema genus and raises the question whether episodes of rare sex are in fact the simplest explanation for the persistence of many old parthenogens in nature.

Evolution of an Alternative Genetic Code in the Providencia Symbiont of the Hematophagous Leech Haementeria acuecueyetzin

Strict blood-feeding animals are confronted with a strong B-vitamin deficiency. Blood-feeding leeches from the Glossiphoniidae family, similarly to hematophagous insects, have evolved specialized organs called bacteriomes to harbor symbiotic bacteria. Leeches of the Haementeria genus have two pairs of globular bacteriomes attached to the esophagus which house intracellular "Candidatus Providencia siddallii" bacteria. Previous work analyzing a draft genome of the Providencia symbiont of the Mexican leech Haementeria officinalis showed that, in this species, the bacteria hold a reduced genome capable of synthesizing B vitamins. In this work, we aimed to expand our knowledge on the diversity and evolution of Providencia symbionts of Haementeria. For this purpose, we sequenced the symbiont genomes of three selected leech species. We found that all genomes are highly syntenic and have kept a stable genetic repertoire, mirroring ancient insect endosymbionts. Additionally, we found B-vitamin pathways to be conserved among these symbionts, pointing to a conserved symbiotic role. Lastly and most notably, we found that the symbiont of H. acuecueyetzin has evolved an alternative genetic code, affecting a portion of its proteome and showing evidence of a lineage-specific and likely intermediate stage of genetic code reassignment.

Convergent resistance to GABA receptor neurotoxins through plant-insect coevolution

The molecular mechanisms of coevolution between plants and insects remain elusive. GABA receptors are targets of many neurotoxic terpenoids, which represent the most diverse array of natural products known. Over deep evolutionary time, as plant terpene synthases diversified in plants, so did plant terpenoid defence repertoires. Here we show that herbivorous insects and their predators evolved convergent amino acid changing substitutions in duplicated copies of the Resistance to dieldrin (Rdl) gene that encodes the GABA receptor, and that the evolution of duplicated Rdl and terpenoid-resistant GABA receptors is associated with the diversification of moths and butterflies. These same substitutions also evolved in pests exposed to synthetic insecticides that target the GABA receptor. We used in vivo genome editing in Drosophila melanogaster to evaluate the fitness effects of each putative resistance mutation and found that pleiotropy both facilitates and constrains the evolution of GABA receptor resistance. The same genetic changes that confer resistance to terpenoids across 300 Myr of insect evolution have re-evolved in response to synthetic analogues over one human lifespan.