Genome assembly, sex-biased gene expression and dosage compensation in the damselfly Ischnura elegans

Behavioral and Transcriptomic Analyses in the Indoxacarb Response of a Non-Target Damselfly Species

Ischnura senegalensis, which widely spreads in paddy fields, has the potential to be used as a natural predator of insect pests. However, the application of insecticides in the field could pose a threat to the survival of I. senegalensis. Among these pesticides, indoxacarb, an oxadiazine insecticide, is renowned for its broad-spectrum efficacy against numerous insect pests. In this study, we examined the toxicity of indoxacarb towards the larvae of I. senegalensis. Behavioral experiments and transcriptome analyses were conducted under indoxacarb treatments. Results revealed that indoxacarb induced abnormal body gestures and significant locomotory impairments, which could ultimately reduce the survival rate of the larvae in their natural habitat. Moreover, transcriptome analyses indicated that genes related to muscle function were significantly affected. Interestingly, at lower concentrations of indoxacarb (0.004 mg/L), the larvae seem to detoxify the indoxacarb with the aid of the cytochrome P450 gene. However, under higher concentrations (0.4 mg/L), the sensory abilities of the larvae were significantly diminished, and they were unable to degrade the toxicity of indoxacarb. Our study underscores the importance of carefully evaluating the impact of insecticides on non-target predatory insects before their widespread application.

The chromosome-level genome assembly of the giant dobsonfly Acanthacorydalis orientalis (McLachlan, 1899)

Acanthacorydalis orientalis (McLachlan, 1899) (Megaloptera: Corydalidae) is an important freshwater-benthic invertebrate species that serves as an indicator for water-quality biomonitoring and is valuable for conservation from East Asia. Here, a high-quality reference genome for A. orientalis was constructed using Oxford Nanopore sequencing and High throughput Chromosome Conformation Capture (Hi-C) technology. The final genome size is 547.98 Mb, with the N50 values of contig and scaffold being 7.77 Mb and 50.53 Mb, respectively. The longest contig and scaffold are 20.57 Mb and 62.26 Mb in length, respectively. There are 99.75% contigs anchored onto 13 pseudo-chromosomes. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis showed that the completeness of the genome assembly is 99.01%. There are 10,977 protein-coding genes identified, of which 84.00% are functionally annotated. The genome contains 44.86% repeat sequences. This high-quality genome provides substantial data for future studies on population genetics, aquatic adaptation, and evolution of Megaloptera and other related insect groups.

The Scorpionfly (Panorpa cognata) Genome Highlights Conserved and Derived Features of the Peculiar Dipteran X Chromosome

Many insects carry an ancient X chromosome-the Drosophila Muller element F-that likely predates their origin. Interestingly, the X has undergone turnover in multiple fly species (Diptera) after being conserved for more than 450 My. The long evolutionary distance between Diptera and other sequenced insect clades makes it difficult to infer what could have contributed to this sudden increase in rate of turnover. Here, we produce the first genome and transcriptome of a long overlooked sister-order to Diptera: Mecoptera. We compare the scorpionfly Panorpa cognata X-chromosome gene content, expression, and structure to that of several dipteran species as well as more distantly related insect orders (Orthoptera and Blattodea). We find high conservation of gene content between the mecopteran X and the dipteran Muller F element, as well as several shared biological features, such as the presence of dosage compensation and a low amount of genetic diversity, consistent with a low recombination rate. However, the 2 homologous X chromosomes differ strikingly in their size and number of genes they carry. Our results therefore support a common ancestry of the mecopteran and ancestral dipteran X chromosomes, and suggest that Muller element F shrank in size and gene content after the split of Diptera and Mecoptera, which may have contributed to its turnover in dipteran insects.

Exploring chromosome evolution in 250 million year old groups of dragonflies and damselflies (Insecta:Odonata)

Using recently published chromosome-length genome assemblies of two damselfly species, Ischnura elegans and Platycnemis pennipes, and two dragonfly species, Pantala flavescens and Tanypteryx hageni, we demonstrate that the autosomes of Odonata have undergone few fission, fusion, or inversion events, despite 250 million years of separation. In the four genomes discussed here, our results show that all autosomes have a clear ortholog in the ancestral karyotype. Despite this clear chromosomal orthology, we demonstrate that different factors, including concentration of repeat dynamics, GC content, relative position on the chromosome, and the relative proportion of coding sequence all influence the density of syntenic blocks across chromosomes. However, these factors do not interact to influence synteny the same way in any two pairs of species, nor is any one factor retained in all four species. Furthermore, it was previously unknown whether the micro-chromosomes in Odonata are descended from one ancestral chromosome. Despite structural rearrangements, our evidence suggests that the micro-chromosomes in the sampled Odonata do indeed descend from an ancestral chromosome, and that the micro-chromosome in P. flavescens was lost through fusion with autosomes.

Chromosome-level Assembly, Dosage Compensation and Sex-biased Gene Expression in the Small Brown Planthopper, Laodelphax striatellus

In insects, sex chromosome differentiation often results in unequal gene dosages between sexes. Dosage compensation mechanisms evolve to balance gene expression, but the degree and mechanism of regulation often vary by insect species. In hemipteran species, the small brown planthopper (SBPH), Laodelphax striatellus, is an injurious crop pest, with a sex chromosome type XX in females and XO in males. This species offers the opportunity to study dosage compensation and sex-biased gene expression. In this study, we generated a chromosome-level genome of SBPH using Oxford Nanopore Technologies and high-throughput chromatin conformation capture (Hi-C) technology. We also sequenced RNA-seq data from 16 tissue samples to annotate the genome and analyze gene dosage compensation. We finally obtained a 510.2 megabases (Mb) genome with 99.12% of the scaffolds anchored on 15 chromosomes (14 autosomes and 1 X chromosome) and annotated 16,160 protein-coding genes based on full-length cDNA sequencing data. Furthermore, we found complete dosage compensation in all L. striatellus somatic tissues, but lack of dosage compensation in gonad tissue testis. We also found that female-biased genes were significantly enriched on the X chromosome in all tissues, whereas male-biased genes in gonad tissues were enriched on autosomes. This study not only provides a high-quality genome assembly but also lays a foundation for a better understanding of the sexual regulatory network in hemipteran insects.

X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

Restricted X chromosome introgression and support for Haldane's rule in hybridizing damselflies

Contemporary hybrid zones act as natural laboratories for the investigation of species boundaries and may shed light on the little understood roles of sex chromosomes in species divergence. Sex chromosomes are considered to function as a hotspot of genetic divergence between species; indicated by less genomic introgression compared to autosomes during hybridization. Moreover, they are thought to contribute to Haldane's rule, which states that hybrids of the heterogametic sex are more likely to be inviable or sterile. To test these hypotheses, we used contemporary hybrid zones of Ischnura elegans, a damselfly species that has been expanding its range into the northern and western regions of Spain, leading to chronic hybridization with its sister species Ischnura graellsii. We analysed genome-wide SNPs in the Spanish I. elegans and I. graellsii hybrid zone and found (i) that the X chromosome shows less genomic introgression compared to autosomes, and (ii) that males are underrepresented among admixed individuals, as predicted by Haldane's rule. This is the first study in Odonata that suggests a role of the X chromosome in reproductive isolation. Moreover, our data add to the few studies on species with X0 sex determination system and contradict the hypothesis that the absence of a Y chromosome causes exceptions to Haldane's rule.

A high-quality genome of the dobsonfly Neoneuromus ignobilis reveals molecular convergences in aquatic insects

Neoneuromus ignobilis is an archaic holometabolous aquatic predatory insect. However, a lack of genomic resources hinders the use of whole genome sequencing to explore their genetic basis and molecular mechanisms for adaptive evolution. Here, we provided a high-contiguity, chromosome-level genome assembly of N. ignobilis using high coverage Nanopore and PacBio reads with the Hi-C technique. The final assembly is 480.67 MB in size, containing 12 telomere-ended pseudochromosomes with only 17 gaps. We compared 42 hexapod species genomes including six independent lineages comprising 11 aquatic insects, and found convergent expansions of long wavelength-sensitive and blue-sensitive opsins, thermal stress response TRP channels, and sulfotransferases in aquatic insects, which may be related to their aquatic adaptation. We also detected strong nonrandom signals of convergent amino acid substitutions in aquatic insects. Collectively, our comparative genomic analysis revealed the evidence of molecular convergences in aquatic insects during both gene family evolution and convergent amino acid substitutions.

Chromosome-level genome of the globe skimmer dragonfly (Pantala flavescens)

Background

The globe skimmer dragonfly (Pantala flavescens) is a notable Odonata insect distributed in nature fields and farmlands worldwide, and it is commonly recognized as a natural enemy because it preys on agricultural pests and health pests. As one of the sister groups of winged insects, odonatan species are key to understanding the evolution of insect wings.

Findings

We present a high-quality reference genome of P. flavescens, which is the first chromosome-level genome in the Palaeoptera (Odonata and Ephemeroptera). The assembled genome size was 662 Mb, with a contig N50 of 16.2 Mb. Via Hi-C scaffolding, 648 Mb (97.9%) of contig sequences were clustered, ordered, and assembled into 12 large scaffolds, each corresponding to a natural chromosome. The X chromosome was identified by sequence coverage depth. The repetitive sequences and gene density of the X chromosome are similar to those of autosomal sequences, but the X chromosome shows a much lower degree of heterozygosity. Our analysis shows that the effective population size experienced 3 declining events, which may have been caused by climate change and environmental pollution.

Conclusions

The genome of P. flavescens provides more information on the biology and evolution of insects and will help for the use of this species in pest control.

Molecular mechanisms underlying metamorphosis in the most-ancestral winged insect

As caterpillars metamorphose to butterflies, insects change their appearance dramatically through metamorphosis. Some insects have an immobile pupal stage for morphological remodeling (homometaboly). Other insects, such as cockroaches, have no pupal stage, and the juveniles and adults are morphologically similar (hemimetaboly). Notably, among the most-ancestral hemimetabolous insects, dragonflies drastically alter their appearance from aquatic nymphs to aerial adults. In dragonflies, we showed that transcription factors Kr-h1 and E93 are essential for regulating metamorphosis as in other insects, while broad, the master gene for pupation in holometabolous insects, regulates a number of both nymph-specific genes and adult-specific genes, providing insight into what evolutionary trajectory the key transcription factor broad has experienced before ending up with governing pupation and holometaboly.

Insects comprise over half of the described species, and the acquisition of metamorphosis must have contributed to their diversity and prosperity. The order Odonata (dragonflies and damselflies) is among the most-ancestral insects with drastic morphological changes upon metamorphosis, in which understanding of the molecular mechanisms will provide insight into the evolution of incomplete and complete metamorphosis in insects. In order to identify metamorphosis-related genes in Odonata, we performed comprehensive RNA-sequencing of the blue-tailed damselfly Ischnura senegalensis at different developmental stages. Comparative RNA-sequencing analyses between nymphs and adults identified eight nymph-specific and seven adult-specific transcripts. RNA interference (RNAi) of these candidate genes demonstrated that three transcription factors, Krüppel homolog 1 (Kr-h1), broad, and E93 play important roles in metamorphosis of both I. senegalensis and a phylogenetically distant dragonfly, Pseudothemis zonata. E93 is essential for adult morphogenesis, and RNAi of Kr-h1 induced precocious metamorphosis in epidermis via up-regulation of E93. Precocious metamorphosis was also induced by RNAi of the juvenile hormone receptor Methoprene-tolerant (Met), confirming that the regulation of metamorphosis by the MEKRE93 (Met-Kr-h1-E93) pathway is conserved across diverse insects including the basal insect lineage Odonata. Notably, RNAi of broad produced unique grayish pigmentation on the nymphal abdominal epidermis. Survey of downstream genes for Kr-h1, broad, and E93 uncovered that unlike other insects, broad regulates a substantial number of nymph-specific and adult-specific genes independently of Kr-h1 and E93. These findings highlight the importance of functional changes and rewiring of the transcription factors Kr-h1, broad, and E93 in the evolution of insect metamorphosis.

The genome sequence of the blue-tailed damselfly, Ischnura elegans (Vander Linden, 1820)

We present a genome assembly from an individual female Ischnura elegans (the blue-tailed damselfly; Arthropoda; Insecta; Odonata; Coenagrionidae). The genome sequence is 1,723 megabases in span. The majority of the assembly (99.55%) is scaffolded into 14 chromosomal pseudomolecules, with the X sex chromosome assembled.

Wolbachia-driven selective sweep in a range expanding insect species

Background

Evolutionary processes can cause strong spatial genetic signatures, such as local loss of genetic diversity, or conflicting histories from mitochondrial versus nuclear markers. Investigating these genetic patterns is important, as they may reveal obscured processes and players. The maternally inherited bacterium Wolbachia is among the most widespread symbionts in insects. Wolbachia typically spreads within host species by conferring direct fitness benefits, and/or by manipulating its host reproduction to favour infected over uninfected females. Under sufficient selective advantage, the mitochondrial haplotype associated with the favoured maternally-inherited symbiotic strains will spread (i.e. hitchhike), resulting in low mitochondrial genetic variation across the host species range.

Method

The common bluetail damselfly (Ischnura elegans: van der Linden, 1820) has recently emerged as a model organism for genetics and genomic signatures of range expansion during climate change. Although there is accumulating data on the consequences of such expansion on the genetics of I. elegans, no study has screened for Wolbachia in the damselfly genus Ischnura. Here, we present the biogeographic variation in Wolbachia prevalence and penetrance across Europe and Japan (including samples from 17 populations), and from close relatives in the Mediterranean area (i.e. I. genei: Rambur, 1842; and I. saharensis: Aguesse, 1958).

Results

Our data reveal (a) multiple Wolbachia-strains, (b) potential transfer of the symbiont through hybridization, (c) higher infection rates at higher latitudes, and (d) reduced mitochondrial diversity in the north-west populations, indicative of hitchhiking associated with the selective sweep of the most common strain. We found low mitochondrial haplotype diversity in the Wolbachia-infected north-western European populations (Sweden, Scotland, the Netherlands, Belgium, France and Italy) of I. elegans, and, conversely, higher mitochondrial diversity in populations with low penetrance of Wolbachia (Ukraine, Greece, Montenegro and Cyprus). The timing of the selective sweep associated with infected lineages was estimated between 20,000 and 44,000 years before present, which is consistent with the end of the last glacial period about 20,000 years.

Conclusions

Our findings provide an example of how endosymbiont infections can shape spatial variation in their host evolutionary genetics during postglacial expansion. These results also challenge population genetic studies that do not consider the prevalence of symbionts in many insects, which we show can impact geographic patterns of mitochondrial genetic diversity.