Search

Background: Gene duplication is a dominant mechanism for the evolution of genomes and plays a key role in genome expansion. Gene duplication via retroposition produces RNA-mediated intron-less copies called retrocopies, that may gain regulatory sequence and biological function to generate retrogenes. Retrocopies dynamics have been reported in several model insect species, but there is still a huge knowledge gap about retrocopies dynamics in most insects, and their role in adaptation., Results: In this study, we reported retrocopy dynamics in 40 species of insect pests of plants belonging to six insect orders. We identified a total of 9,930 retrocopies, which is so far the largest set of retrocopies identified in insects. The identified retrocopies were further grouped into 2,599 Retrogenes, 4,578 Chimeras, 1,241 Intact retrocopies, and 1,512 Pseudogene. We also analyzed all the identified retrogenes that were annotated into 506 gene families. The highest number of retrogenes annotated belong to the heat shock proteins superfamily and are present across all the 40 species from the six orders. We found a significant expansion of the heat shock protein superfamily in the studied species. Almost all the retrogenes, including those belonging to heat shock proteins, are under purifying selection. In summary, we report the retrocopies and retrogenes dynamics in a large set of insect pests of plants and the expansion of the heat shock protein family due to retroposition., Conclusion: This study unveils retrocopy dynamics in the insect pests of plants and highlights the evolution of new genes due to retroposition, and their role in important gene families' expansion., Competing Interests: Declarations. Ethics approval and consent to participate: Not applicable. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Background: Adenosine-to-inosine (A-to-I) RNA editing is a co-/post-transcriptional modification introducing A-to-G variations in RNAs. There is extensive discussion on whether the flexibility of RNA editing exerts a proteomic diversification role, or it just acts like hardwired mutations to correct the genomic allele. Eusocial insects evolved the ability to generate phenotypically differentiated individuals with the same genome, indicating the involvement of epigenetic/transcriptomic regulation., Methods: We obtained the genomes of 104 Hymenoptera insects and the transcriptomes of representative species. Comparative genomic analysis was performed to parse the evolutionary trajectory of a regulatory Ile > Met auto-recoding site in Adar gene., Results: At genome level, the pre-editing Ile codon is conserved across a node containing all eusocial hymenopterans. At RNA level, the editing events are confirmed in representative species and shows considerable condition-specificity. Compared to random expectation, the editable Ile codon avoids genomic substitutions to Met or to uneditable Ile codons, but does not avoid mutations to other unrelated amino acids., Conclusions: The flexibility of Adar auto-recoding site in Hymenoptera is selectively maintained, supporting the flexible RNA editing hypothesis. We proposed a new angle to view the adaptation of RNA editing, providing another layer to explain the great phenotypical plasticity of eusocial insects., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

Background: The Argentine stem weevil (ASW, Listronotus bonariensis) is a significant pasture pest in Aotearoa New Zealand, primarily controlled by the parasitoid biocontrol agent Microctonus hyperodae. Despite providing effective control of ASW soon after release, M. hyperodae parasitism rates have since declined significantly, with ASW hypothesised to have evolved resistance to its biocontrol agent. While the parasitism arsenal of M. hyperodae has previously been investigated, revealing many venom components and an exogenous novel DNA virus Microctonus hyperodae filamentous virus (MhFV), the effects of said arsenal on gene expression in ASW during parasitism have not been examined. In this study, we performed a multi-species transcriptomic analysis to investigate the biology of ASW parasitism by M. hyperodae, as well as the decline in efficacy of this biocontrol system., Results: The transcriptomic response of ASW to parasitism by M. hyperodae involves modulation of the weevil's innate immune system, flight muscle components, and lipid and glucose metabolism. The multispecies approach also revealed continued expression of venom components in parasitised ASW, as well as the transmission of MhFV to weevils during parasitism and some interrupted parasitism attempts. Transcriptomics did not detect a clear indication of parasitoid avoidance or other mechanisms to explain biocontrol decline., Conclusions: This study has expanded our understanding of interactions between M. hyperodae and ASW in a biocontrol system of critical importance to Aotearoa-New Zealand's agricultural economy. Transmission of MhFV to ASW during successful and interrupted parasitism attempts may link to a premature mortality phenomenon in ASW, hypothesised to be a result of a toxin-antitoxin system. Further research into MhFV and its potential role in ASW premature mortality is required to explore whether manipulation of this viral infection has the potential to increase biocontrol efficacy in future., (© 2024. The Author(s).)

Background: Euholognatha is a monophyletic group within stoneflies comprised by a superfamily Nemouroidea and a family Scopuridae. Based on morphological data, the family-level phylogenetic relationships within Euholognatha are widely accepted, but there is still controversy among different molecular studies. To better understand the phylogeny of all six extant euholognathan families, we sequenced and analyzed seven euholognathan mitogenomes., Results: The sequence heterogeneity analysis observed a low degree of compositional heterogeneity in euholognathan mitogenomes. Meanwhile, leuctrid mitogenomes were more heterogeneous than other euholognathan families, which may affect the phylogenetic reconstruction. Phylogenetic analyses with various datasets generated three topologies. The Leuctridae was recovered as the earliest branching lineage, and the sister relationship of Capniidae and Taeniopterygidae was supported by most tree topologies and FcLM analyses. When separately excluding sparsely sampled Scopuridae or high heterogeneity leuctrid taxa, phylogenetic analyses under the same methods generated more stable and consistent tree topologies. Finally, based on the results of this study, we reconstructed the relationships within Euholognatha as: Leuctridae + (Scopuridae + ((Taeniopterygidae + Capniidae) + (Nemouridae + Notonemouridae)))., Conclusion: Our research shows the potential of data optimizing strategies in reconstructing phylogeny within Euholognatha and provides new insight into the phylogeny of this group., (© 2024. The Author(s).)

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., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: Freshwater ecosystems, such as streams, are facing increasing pressures from agricultural land use and recent literature stresses the importance of robust biomonitoring to detect trends in insect decline globally. Aquatic insects and other macroinvertebrates are often used as indicators of ecological condition in freshwater biomonitoring programs; however, these diverse groups can present challenges to morphological identification and coarse-level taxonomic resolution can mask patterns in community composition. Here, we incorporate molecular identification (DNA metabarcoding) into a stream biomonitoring sampling design to explore the diversity and variability of aquatic macroinvertebrate communities at small spatial scales. While individual stream reaches can be very heterogenous, most community ecology studies focus on larger, landscape-level patterns of community composition. A high degree of community variability at the local scale has important implications for both biomonitoring and ecological research, and the incorporation of DNA metabarcoding into local biodiversity assessments will inform future sampling protocols., Results: We sampled twenty streams in southern Ontario, Canada, for aquatic macroinvertebrates across multiple time points and assessed local community variability by comparing field replicates taken ten meters apart within the same stream. Using bulk-tissue DNA metabarcoding, we revealed that aquatic macroinvertebrate communities are highly diverse at small spatial scales with unprecedented levels of local taxonomic turnover. We detected over 1600 Operational Taxonomic Units (OTUs) from 149 families, and a single insect family, the Chironomidae, contained over one third of the total number of OTUs detected in our study. Benthic communities were largely comprised of rare taxa detected only once per stream despite multiple biological replicates (24-94% rare taxa per site). In addition to numerous rare taxa, our species pool estimates indicated that there was a large proportion of taxa that remained undetected by our sampling regime (14-94% per site). Our sites were located across a gradient of agricultural activity, and while we predicted that increased land use would homogenize benthic communities, this was not supported as within-stream dissimilarity was unrelated to land use. Within-stream dissimilarity estimates were consistently high for all levels of taxonomic resolution (invertebrate families, invertebrate OTUs, chironomid OTUs), indicating stream communities are very dissimilar at small spatial scales., (© 2023. The Author(s).)

Background: The western flower thrips Frankliniella occidentalis is an insect pest that damages various crops, including hot peppers. It is a vector of a plant pathogen, tomato spotted wilt virus. To control this pest, chemical insecticides have been used in the past, but the control efficacy is unsatisfactory owing to rapid resistance development by F. occidentalis., Methodology: This study reports a novel control technology against this insect pest using RNA interference (RNAi) of the vacuolar-type ATPase (vATPase) expression. Eight subunit genes (vATPase-A ∼ vATPase-H) of vATPase were obtained from the F. occidentalis genome and confirmed for their expressions at all developmental stages., Results: Double-stranded RNAs (dsRNAs) specific to the eight subunit genes were fed to larvae and adults, which significantly suppressed the corresponding gene expressions after 24-h feeding treatment. These RNAi treatments resulted in significant mortalities, in which the dsRNA treatments at ∼2,000 ppm specific to vATPase-A or vATPase-B allowed complete control efficacy near 100% mortality in 7 days after treatment. To prevent dsRNA degradation by the digestive proteases during oral feeding, dsRNAs were formulated in a liposome and led to an enhanced mortality of the larvae and adults of F. occidentalis. The dsRNAs were then sprayed at 2,000 ppm on F. occidentalis infesting hot peppers in a greenhouse, which resulted in 53.5-55.9% control efficacy in 7 days after treatment. Even though the vATPases are conserved in different organisms, the dsRNA treatment was relatively safe for non-target insects owing to the presence of mismatch sequences compared to the dsRNA region of F. occidentalis., Conclusion: These results demonstrate the practical feasibility of spraying dsRNA to control F. occidentalis infesting crops., (© 2023. The Author(s).)

Background: Generating the most contiguous, accurate genome assemblies given available sequencing technologies is a long-standing challenge in genome science. With the rise of long-read sequencing, assembly challenges have shifted from merely increasing contiguity to correctly assembling complex, repetitive regions of interest, ideally in a phased manner. At present, researchers largely choose between two types of long read data: longer, but less accurate sequences, or highly accurate, but shorter reads (i.e., >Q20 or 99% accurate). To better understand how these types of long-read data as well as scale of data (i.e., mean length and sequencing depth) influence genome assembly outcomes, we compared genome assemblies for a caddisfly, Hesperophylax magnus, generated with longer, but less accurate, Oxford Nanopore (ONT) R9.4.1 and highly accurate PacBio HiFi (HiFi) data. Next, we expanded this comparison to consider the influence of highly accurate long-read sequence data on genome assemblies across 6750 plant and animal genomes. For this broader comparison, we used HiFi data as a surrogate for highly accurate long-reads broadly as we could identify when they were used from GenBank metadata., Results: HiFi reads outperformed ONT reads in all assembly metrics tested for the caddisfly data set and allowed for accurate assembly of the repetitive ~ 20 Kb H-fibroin gene. Across plants and animals, genome assemblies that incorporated HiFi reads were also more contiguous. For plants, the average HiFi assembly was 501% more contiguous (mean contig N50 = 20.5 Mb) than those generated with any other long-read data (mean contig N50 = 4.1 Mb). For animals, HiFi assemblies were 226% more contiguous (mean contig N50 = 20.9 Mb) versus other long-read assemblies (mean contig N50 = 9.3 Mb). In plants, we also found limited evidence that HiFi may offer a unique solution for overcoming genomic complexity that scales with assembly size., Conclusions: Highly accurate long-reads generated with HiFi or analogous technologies represent a key tool for maximizing genome assembly quality for a wide swath of plants and animals. This finding is particularly important when resources only allow for one type of sequencing data to be generated. Ultimately, to realize the promise of biodiversity genomics, we call for greater uptake of highly accurate long-reads in future studies., (© 2023. The Author(s).)

Background: The nose of most animals comprises multiple sensory subsystems, which are defined by the expression of different olfactory receptor families. Drosophila melanogaster antennae contain two morphologically and functionally distinct subsystems that express odorant receptors (Ors) or ionotropic receptors (Irs). Although these receptors have been thoroughly characterized in this species, the subsystem-specific expression and roles of other genes are much less well-understood., Results: Here we generate subsystem-specific transcriptomic datasets to identify hundreds of genes, encoding diverse protein classes, that are selectively enriched in either Or or Ir subsystems. Using single-cell antennal transcriptomic data and RNA in situ hybridization, we find that most neuronal genes-other than sensory receptor genes-are broadly expressed within the subsystems. By contrast, we identify many non-neuronal genes that exhibit highly selective expression, revealing substantial molecular heterogeneity in the non-neuronal cellular components of the olfactory subsystems. We characterize one Or subsystem-specific non-neuronal molecule, Osiris 8 (Osi8), a conserved member of a large, insect-specific family of transmembrane proteins. Osi8 is expressed in the membranes of tormogen support cells of pheromone-sensing trichoid sensilla. Loss of Osi8 does not have obvious impact on trichoid sensillar development or basal neuronal activity, but abolishes high sensitivity responses to pheromone ligands., Conclusions: This work identifies a new protein required for insect pheromone detection, emphasizes the importance of support cells in neuronal sensory functions, and provides a resource for future characterization of other olfactory subsystem-specific genes., (© 2022. The Author(s).)

Background: Hoverflies (Diptera: Syrphidae) including Eupeodes corollae are important insects worldwide that provide dual ecosystem services including pest control and pollination. The larvae are dominant predators of aphids and can be used as biological control agents, and the adults are efficient pollinators. The different feeding habits of larvae and adults make hoverflies a valuable genetic resource for understanding the mechanisms underlying the evolution and adaptation to predation and pollination in insects., Results: Here, we present a 595-Mb high-quality reference genome of the hoverfly E. corollae, which is typical of an aphid predator and a pollinator. Comparative genomic analyses of E. corollae and Coccinellidae (ladybugs, aphid predators) shed light on takeout genes (3), which are involved in circadian rhythms and feeding behavior and might regulate the feeding behavior of E. corollae in a circadian manner. Genes for sugar symporter (12) and lipid transport (7) related to energy production in E. corollae had homologs in pollinator honeybees and were absent in predatory ladybugs. A number of classical cytochrome P450 detoxification genes, mainly CYP6 subfamily members, were greatly expanded in E. corollae. Notably, comparative genomic analyses of E. corollae and other aphidophagous hoverflies highlighted three homologous trypsins (Ecor12299, Ecor12301, Ecor2966). Transcriptome analysis showed that nine trypsins, including Ecor12299, Ecor12301, and Ecor2966, are strongly expressed at the larval stage, and 10 opsin genes, which are involved in visual perception, are significantly upregulated at the adult stage of E. corollae., Conclusions: The high-quality genome assembly provided new insights into the genetic basis of predation and pollination by E. corollae and is a valuable resource for advancing studies on genetic adaptations and evolution of hoverflies and other natural enemies., (© 2022. The Author(s).)

Background: The re-evolution of complex characters is generally considered impossible, yet, studies of recent years have provided several examples of phenotypic reversals shown to violate Dollo's law. Along these lines, the regain of wings in stick and leaf insects (Phasmatodea) was hypothesised to have occurred several times independently after an ancestral loss, a scenario controversially discussed among evolutionary biologists due to overestimation of the potential for trait reacquisition as well as to the lack of taxonomic data., Results: We revisited the recovery of wings by reconstructing a phylogeny based on a comprehensive taxon sample of over 500 representative phasmatodean species to infer the evolutionary history of wings. We additionally explored the presence of ocelli, the photoreceptive organs used for flight stabilisation in winged insects, which might provide further information for interpreting flight evolution. Our findings support an ancestral loss of wings and that the ancestors of most major lineages were wingless. While the evolution of ocelli was estimated to be dependent on the presence of (fully-developed) wings, ocelli are nevertheless absent in the majority of all examined winged species and only appear in the members of few subordinate clades, albeit winged and volant taxa are found in every euphasmatodean lineage., Conclusion: In this study, we explored the evolutionary history of wings in Phasmatodea and demonstrate that the disjunct distribution of ocelli substantiates the hypothesis on their regain and thus on trait reacquisition in general. Evidence from the fossil record as well as future studies focussing on the underlying genetic mechanisms are needed to validate our findings and to further assess the evolutionary process of phenotypic reversals., (© 2022. The Author(s).)

Background: Insects have exceptionally fast smelling capabilities, and some can track the temporal structure of odour plumes at rates above 100 Hz. It has been hypothesized that this fast smelling capability is an adaptation for flying. We test this hypothesis by comparing the olfactory acuity of sympatric flighted versus flightless lineages within a wing-polymorphic stonefly species., Results: Our analyses of olfactory receptor neuron responses reveal that recently-evolved flightless lineages have reduced olfactory acuity. By comparing flighted versus flightless ecotypes with similar genetic backgrounds, we eliminate other confounding factors that might have affected the evolution of their olfactory reception mechanisms. Our detection of different patterns of reduced olfactory response strength and speed in independently wing-reduced lineages suggests parallel evolution of reduced olfactory acuity., Conclusions: These reductions in olfactory acuity echo the rapid reduction of wings themselves, and represent an olfactory parallel to the convergent phenotypic shifts seen under selective gradients in other sensory systems (e.g. parallel loss of vision in cave fauna). Our study provides evidence for the hypothesis that flight poses a selective pressure on the speed and strength of olfactory receptor neuron responses and emphasizes the energetic costs of rapid olfaction., (© 2022. The Author(s).)

Background: Visual opsins are expressed in the compound eyes and ocelli of insects and enable light detection. Three distinct phylogenetic groups of visual opsins are found in insects, named long (LW), short (SW) and ultraviolet (UV) wavelength sensitive opsins. Recently, the LW group was found to be duplicated into the LW2b and the LW2a opsins. The expression of LW2b opsins is ocelli specific in some insects (e.g., bees, cricket, scorpion flies), but the gene was not found in other orders possessing three or less ocelli (e.g., dragonflies, beetles, moths, bugs). In flies, two LW2b homologs have been characterised, with one expressed in the ocelli and the other in the compound eyes. To date, it remains unclear which evolutionary forces have driven gains and losses of LW opsins in insects. Here we take advantage of the recent rapid increase in available sequence data (i.e., from insect genomes, targeted PCR amplification, RNAseq) to characterize the phylogenetic relationships of 1000 opsin sequences in 18 orders of Insects. The resulting phylogeny discriminates between four main groups of opsins, and onto this phylogeny we mapped relevant morphological and life history traits., Results: Our results demonstrate a conserved LW2b opsin only present in insects with three ocelli. Only two groups (Brachycera and Odonata) possess more than one LW2b opsin, likely linked to their life history. In flies, we hypothesize that the duplication of the LW2b opsin occurred after the transition from aquatic to terrestrial larvae. During this transition, higher flies (Brachycera) lost a copy of the LW2a opsin, still expressed and duplicated in the compound eyes of lower flies (Nematocera). In higher flies, the LW2b opsin has been duplicated and expressed in the compound eyes while the ocelli and the LW2b opsin were lost in lower flies. In dragonflies, specialisation of flight capabilities likely drove the diversification of the LW2b visual opsins., Conclusion: The presence of the LW2b opsin in insects possessing three ocelli suggests a role in specific flight capabilities (e.g., stationary flight). This study provides the most complete view of the evolution of visual opsin genes in insects yet, and provides new insight into the influence of ocelli and life history traits on opsin evolution in insects., (© 2022. The Author(s).)

Background: Conserved syntenic gene complexes are rare in Arthropods and likely only retained due to functional constraint. Numerous sHSPs have been identified in the genomes of insects, some of which are located clustered in close proximity. Previous phylogenetic analyses of these clustered sHSP have been limited to a small number of holometabolous insect species and have not determined the pattern of evolution of the clustered sHSP genes (sHSP-C) in insect or Arthropod lineages., Results: Using eight genomes from representative insect orders and three non-insect arthropod genomes we have identified that a syntenic cluster of sHSPs (sHSP-C) is a hallmark of most Arthropod genomes. Using 11 genomes from Hymenopteran species our phylogenetic analyses have refined the evolution of the sHSP-C in Hymenoptera and found that the sHSP-C is order-specific with evidence of birth-and-death evolution in the hymenopteran lineage. Finally we have shown that the honeybee sHSP-C is co-ordinately expressed and is marked by genomic features, including H3K27me3 histone marks consistent with coordinate regulation, during honeybee ovary activation., Conclusions: The syntenic sHSP-C is present in most insect genomes, and its conserved coordinate expression and regulation implies that it is an integral genomic component of environmental response in arthropods., (© 2021. The Author(s).)

Background: Fragmented mitochondrial (mt) genomes and extensive mt gene rearrangements have been frequently reported from parasitic lice (Insecta: Phthiraptera). However, relatively little is known about the mt genomes from the family Philopteridae, the most species-rich family within the suborder Ischnocera., Methods: Herein, we use next-generation sequencing to decode the mt genome of Falcolipeurus suturalis and compare it with the mt genome of F. quadripustulatus. Phylogenetic relationships within the family Philopteridae were inferred from the concatenated 13 protein-coding genes of the two Falcolipeurus lice and members of the family Philopteridae using Bayesian inference (BI) and maximum likelihood (ML) methods., Results: The complete mt genome of F. suturalis is a circular, double-stranded DNA molecule 16,659 bp in size that contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and three non-coding regions. The gene order of the F. suturalis mt genome is rearranged relative to that of F. quadripustulatus, and is radically different from both other louse species and the putative ancestral insect. Phylogenetic analyses revealed clear genetic distinctiveness between F. suturalis and F. quadripustulatus (Bayesian posterior probabilities = 1.0 and bootstrapping frequencies = 100), and that the genus Falcolipeurus is sister to the genus Ibidoecus (Bayesian posterior probabilities = 1.0 and bootstrapping frequencies = 100)., Conclusions: These datasets help to better understand gene rearrangements in lice and the phylogenetic position of Falcolipeurus and provide useful genetic markers for systematic studies of bird lice.

Objective: Insects are the most evolutionarily successful groups of organisms, and this success is largely due to their flight ability. Interestingly, some stick insects have lost their flight ability despite having wings. To elucidate the shift from wingless to flying forms during insect evolution, we compared the nutritional metabolism system among flight-winged, flightless-winged, and flightless-wingless stick insect groups., Results: Here, we report RNA sequencing of midgut transcriptome of Entoria okinawaensis, a prominent Japanese flightless-wingless stick insect, and the comparative analysis of its transcriptome in publicly available midgut transcriptomes obtained from seven stick insect species. A gene enrichment analysis for differentially expressed genes, including those obtained from winged vs wingless and flight vs flightless genes comparisons, revealed that carbohydrate metabolic process-related genes were highly expressed in the winged stick insect group. We also found that the expression of the mitochondrial enolase superfamily member 1 transcript was significantly higher in the winged stick insect group than in the wingless stick insect group. Our findings could indicate that carbohydrate metabolic processes are related to the evolutionary process through which stick insects gain the ability of flight.

Background: The blood-feeding behavior evolved multiple times in Insecta lineages and it represents an excellent opportunity to study patterns of convergent molecular evolution regarding this habit. In insects the expansion of some gene families is linked with blood-feeding behavior, but a wide study comparing the evolution of these gene families among different lineages is still missing. Here we gathered genomic data from six independently-evolved hematophagous lineages, aiming to identify convergent expansions and/or contractions of gene families in hematophagous lineages of insects., Results: We found four rapidly evolving gene families shared by at least two hematophagous independently-evolved lineages, including a heat-shock and a chemosensory protein. On the expression of these four rapidly evolving gene families we found more genes expressed in mated individuals compared with virgin individuals in rapidly-expanded families and more genes expressed in non-blood-feeding individuals compared with blood-feeding individuals in rapidly-contracted families., Conclusion: Our results reveal a new set of candidate genes to be explored in further analysis to help the development of new strategies to deal with blood-feeding vectors and also presents a new perspective to study the evolution of hematophagy identifying convergent molecular patterns.

Background: Gene duplication events play an important role in the evolution and adaptation of organisms. Duplicated genes can arise through different mechanisms, including whole-genome duplications (WGDs). Recently, WGD was suggested to be an important driver of evolution, also in hexapod animals., Results: Here, we analyzed 20 high-quality hexapod genomes using whole-paranome distributions of estimated synonymous distances (K S ), patterns of within-genome co-linearity, and phylogenomic gene tree-species tree reconciliation methods. We observe an abundance of gene duplicates in the majority of these hexapod genomes, yet we find little evidence for WGD. The majority of gene duplicates seem to have originated through small-scale gene duplication processes. We did detect segmental duplications in six genomes, but these lacked the within-genome co-linearity signature typically associated with WGD, and the age of these duplications did not coincide with particular peaks in K S distributions. Furthermore, statistical gene tree-species tree reconciliation failed to support all but one of the previously hypothesized WGDs., Conclusions: Our analyses therefore provide very limited evidence for WGD having played a significant role in the evolution of hexapods and suggest that alternative mechanisms drive gene duplication events in this group of animals. For instance, we propose that, along with small-scale gene duplication events, episodes of increased transposable element activity could have been an important source for gene duplicates in hexapods.

Background: Transposable elements (TEs) are a major component of metazoan genomes and are associated with a variety of mechanisms that shape genome architecture and evolution. Despite the ever-growing number of insect genomes sequenced to date, our understanding of the diversity and evolution of insect TEs remains poor., Results: Here, we present a standardized characterization and an order-level comparison of arthropod TE repertoires, encompassing 62 insect and 11 outgroup species. The insect TE repertoire contains TEs of almost every class previously described, and in some cases even TEs previously reported only from vertebrates and plants. Additionally, we identified a large fraction of unclassifiable TEs. We found high variation in TE content, ranging from less than 6% in the antarctic midge (Diptera), the honey bee and the turnip sawfly (Hymenoptera) to more than 58% in the malaria mosquito (Diptera) and the migratory locust (Orthoptera), and a possible relationship between the content and diversity of TEs and the genome size., Conclusion: While most insect orders exhibit a characteristic TE composition, we also observed intraordinal differences, e.g., in Diptera, Hymenoptera, and Hemiptera. Our findings shed light on common patterns and reveal lineage-specific differences in content and evolution of TEs in insects. We anticipate our study to provide the basis for future comparative research on the insect TE repertoire.

Background: Members of the transient receptor potential (TRP) superfamily are proteins that are critical for insects to detect changes in environmental stimuli and also play key roles in their sensory physiology. Moreover, this family provides potential targets for the design of insecticides. In contrast to a large number of studies conducted on Drosophila melanogaster, molecular studies to characterize TRP channels in agricultural pests are lacking., Results: In this study, we identified 15 TRP channel genes in the genome of a notorious agricultural pest, the oriental fruit fly (Bactrocera dorsalis). Comparative analysis of the TRP channels (TRPs) in B. dorsalis with those in D. melanogaster, Glossina morsitans, Musca domestica and the closely related Ceratitis capitata, and TRPs from mosquitoes, Hymenoptera, Lepidoptera, Coleoptera and Hemiptera reveals that members of TRPA and TRPP subfamily are most diverse among insects. The results also suggest that Tephritidae family have two TRP-Polycystin 2 members even though most insects either possess just one or none. The highest expression levels of these two genes are in the testes of B. dorsalis, implying a role in regulating sperm function. We analyzed the expression profiles of the TRP channels identified in this study at different life stages using quantitative real time PCR. The results of this study demonstrate that all TRP channels are mainly expressed in adults, especially at mature stages. The one exception to this trend is BdTRPM, which is more highly expressed in the eggs of B. dorsalis, implying an important role in early development. We also detected the spatial expression of TRP channels in mature adult fruit flies by investigating expression levels within various tissues including those involved in sensory function, such as antennae, compound eyes, mouthparts, legs, and wings, as well as tissues critical for homeostasis and physiology (i.e., Malpighian tubules, the brain and gut as well as fat bodies, ovaries, and testes)., Conclusion: The results of this study establish a solid foundation for future functional characterization of B. dorsalis TRP channels as well as those of other insects and will help future insecticide design targeting these channels.

Background: The ever increasing availability of genomes makes it possible to investigate and compare not only the genomic complements of genes and proteins, but also of RNAs. One class of RNAs, the long noncoding RNAs (lncRNAs) and, in particular, their subclass of long intergenic noncoding RNAs (lincRNAs) have recently gained much attention because of their roles in regulation of important biological processes such as immune response or cell differentiation and as possible evolutionary precursors for protein coding genes. lincRNAs seem to be poorly conserved at the sequence level but at least some lincRNAs have conserved structural elements and syntenic genomic positions. Previous studies showed that transposable elements are a main contribution to the evolution of lincRNAs in mammals. In contrast, plant lincRNA emergence and evolution has been linked with local duplication events. However, little is known about their evolutionary dynamics in general and in insect genomes in particular., Results: Here we compared lincRNAs between seven insect genomes and investigated possible evolutionary changes and functional roles. We find very low sequence conservation between different species and that similarities within a species are mostly due to their association with transposable elements (TE) and simple repeats. Furthermore, we find that TEs are less frequent in lincRNA exons than in their introns, indicating that TEs may have been removed by selection. When we analysed the predicted thermodynamic stabilities of lincRNAs we found that they are more stable than their randomized controls which might indicate some selection pressure to maintain certain structural elements. We list several of the most stable lincRNAs which could serve as prime candidates for future functional studies. We also discuss the possibility of de novo protein coding genes emerging from lincRNAs. This is because lincRNAs with high GC content and potentially with longer open reading frames (ORF) are candidate loci where de novo gene emergence might occur., Conclusion: The processes responsible for the emergence and diversification of lincRNAs in insects remain unclear. Both duplication and transposable elements may be important for the creation of new lincRNAs in insects.

Background: Taxonomy offers precise species identification and delimitation and thus provides basic information for biological research, e.g. through assessment of species richness. The importance of molecular taxonomy, i.e., the identification and delimitation of taxa based on molecular markers, has increased in the past decade. Recently developed exploratory tools now allow estimating species-level diversity in multi-locus molecular datasets., Results: Here we use molecular species delimitation tools that either quantify differences in intra- and interspecific variability of loci, or divergence times within and between species, or perform coalescent species tree inference to estimate species-level entities in molecular genetic datasets. We benchmark results from these methods against 14 morphologically readily differentiable species of a well-defined subgroup of the diverse Drusinae subfamily (Trichoptera, Limnephilidae). Using a 3798 bp (6 loci) molecular data set we aim to corroborate a geographically isolated new species by integrating comparative morphological studies and molecular taxonomy., Conclusions: Our results indicate that only multi-locus species delimitation provides taxonomically relevant information. The data further corroborate the new species Drusus zivici sp. nov. We provide differential diagnostic characters and describe the male, female and larva of this new species and discuss diversity patterns of Drusinae in the Balkans. We further discuss potential and significance of molecular species delimitation. Finally we argue that enhancing collaborative integrative taxonomy will accelerate assessment of global diversity and completion of reference libraries for applied fields, e.g., conservation and biomonitoring.

Background: The entomopathogenic nematode Steinernema carpocapsae has been used worldwide as a biocontrol agent for insect pests, making it an interesting model for understanding parasite-host interactions. Two models propose that these interactions are co-evolutionary processes in such a way that equilibrium is never reached. In one model, known as "arms race", new alleles in relevant genes are fixed in both host and pathogens by directional positive selection, producing recurrent and alternating selective sweeps. In the other model, known as"trench warfare", persistent dynamic fluctuations in allele frequencies are sustained by balancing selection. There are some examples of genes evolving according to both models, however, it is not clear to what extent these interactions might alter genome-level evolutionary patterns and intraspecific diversity. Here we investigate some of these aspects by studying genomic variation in S. carpocapsae and other pathogenic and free-living nematodes from phylogenetic clades IV and V., Results: To look for signatures of an arms-race dynamic, we conducted massive scans to detect directional positive selection in interspecific data. In free-living nematodes, we detected a significantly higher proportion of genes with sites under positive selection than in parasitic nematodes. However, in these genes, we found more enriched Gene Ontology terms in parasites. To detect possible effects of dynamic polymorphisms interactions we looked for signatures of balancing selection in intraspecific genomic data. The observed distribution of Tajima's D values in S. carpocapsae was more skewed to positive values and significantly different from the observed distribution in the free-living Caenorhabditis briggsae. Also, the proportion of significant positive values of Tajima's D was elevated in genes that were differentially expressed after induction with insect tissues as compared to both non-differentially expressed genes and the global scan., Conclusions: Our study provides a first portrait of the effects that lifestyle might have in shaping the patterns of selection at the genomic level. An arms-race between hosts and pathogens seems to be affecting specific genetic functions but not necessarily increasing the number of positively selected genes. Trench warfare dynamics seem to be acting more generally in the genome, likely focusing on genes responding to the interaction, rather than targeting specific genetic functions.

Background: The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), has become the most damaging insect pest of corn in Asia. However, the lack of genome or transcriptome information heavily hinders our further understanding of ACB in every aspect at a molecular level and on a genome-wide scale. Here, we used the Ion Torrent Personal Genome Machine (PGM) Sequencer to explore the ACB transcriptome and to identify relevant genes in response to flubendiamide, showing high selective activity against ACB., Results: We obtained 35,430 unigenes, with an average length of 716 bp, representing a dramatic expansion of existing cDNA sequences available for ACB. These sequences were annotated with Non-redundant Protein (Nr), Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to better understand their functions. A total of 31 cytochrome P450 monooxygenases (P450s), 27 carboxyl/cholinesterases (CCEs) and 19 glutathione S-transferases (GSTs) were manually curated to construct phylogenetic trees, and 25 unigenes encoding target proteins (acetylcholinesterase, nicotinic acetylcholine receptor, gamma-aminobutyric acid receptor, glutamate-gated chloride channel, voltage-gated sodium channel and ryanodine receptor) were identified. In addition, we compared and validated the differentially expressed unigenes upon flubendiamide treatment, revealing that the genes for detoxification enzymes (P450s and esterase), calcium signaling pathways and muscle control pathways (twitchin and tropomyosin), immunoglobulin (hemolin), chemosensory protein and heat shock protein 70 were significantly overexpressed in response to flubendiamide, while the genes for cuticular protein, protease and oxidoreductase showed much lower expression levels., Conclusion: The obtained transcriptome information provides large genomic resources available for further studies of ACB. The differentially expressed gene data will elucidate the molecular mechanisms of ACB in response to the novel diamide insecticide, flubendiamide. In particular, these findings will facilitate the identification of the genes involved in insecticide resistance and the development of new compounds to control the ACB.

Background: The last time the phylogenetic relationships among members of the family Hemerobiidae were studied quantitatively was over 12 years ago and based exclusively on morphology. Our study builds upon this morphological evidence by adding sequence data from three gene loci to provide a total evidence phylogeny of brown lacewings (Neuroptera: Hemerobiidae). Thirty-seven species representing nineteen Hemerobiidae genera were compared with outgroups from the families Ithonidae, Psychopsidae and Chrysopidae in Bayesian and parsimony analyses using a single nuclear gene (CAD) and two mitochondrial (16S rDNA and Cytochrome Oxidase I) genes. We compare divergence time estimates of Hemerobiidae cladogenesis under the two most commonly used relaxed clock models and discuss the evolution of wing venation in the family., Results: We recovered a phylogeny largely incongruent with previously published morphological studies, although all but two subfamilies (i.e., Notiobiellinae and Drepanacrinae) were recovered as monophyletic. We found the subfamily Drepanacrinae paraphyletic with respect to Psychobiellinae, and Notiobiellinae to be polyphyletic. We thus offer a revised concept of Notiobiellinae, comprising only Notiobiella Banks, and erect a new subfamily Zachobiellinae including the remaining genera previously placed in Notiobiellinae. Psychobiellinae is synonymized with Drepanacrinae. Unlike the previous hypothesis that proposed a remarkably laddered topology, our tree suggests that hemerobiids diverged as three main clades. Moreover, in contrast to the vein proliferation hypothesis, we found that hemerobiids have instead undergone multiple reductions in the number of radial veins, this scenario questions the relevance of this character as diagnostic of various subfamilies, Conclusions: Our phylogenetic hypothesis and divergence times analysis suggest that extant hemerobiids originated around the end of the Triassic and evolved as three distinct clades that diverged from one another during the Late Jurassic to Early Cretaceous. Contrary to earlier phylogenetic hypotheses, Carobius Banks (Carobiinae) is sister to the previously unplaced genus Notherobius New in a clade more closely related to Sympherobiinae, Megalominae and Zachobiellinae subfam. nov. The addition of taxa which are not available for DNA sequencing should be the focus of future studies, especially Adelphohemerobius Oswald, which is particularly important to test our inferences regarding the evolution of wing venation in Hemerobiidae.

Background: Grapevine phylloxera, an insect related to true aphids, is a major historic pest of viticulture only controlled through the selection of resistant rootstocks or through quarantine regulations where grapevine is cultivated own-rooted. Transcriptomic data could help understand the bases of its original life-traits, including a striking case of polyphenism, with forms feeding on roots and forms feeding in leaf-galls. Comparisons with true aphids (for which complete genomes have been sequenced) should also allow to link differences in life-traits of the two groups with changes in gene repertoires or shifts in patterns of expression., Results: We sequenced transcriptomes of the grapevine phylloxera (Illumina technology), choosing three life-stages (adults on roots or on leaf galls, and eggs) to cover a large catalogue of transcripts, and performed a de novo assembly. This resulted in 105,697 contigs, which were annotated: most contigs had a best blastx hit to the pea aphid (phylogenetically closest complete genome), while very few bacterial hits were recorded (except for Probionibacterium acnes). Coding sequences were predicted from this data set (17,372 sequences), revealing an extremely high AT-bias (at the third codon position). Differential expression (DE) analysis among root-feeding and gall-feeding showed that i) the root-feeding form displayed a much larger number of differentially expressed transcripts ii) root-feeding biased genes were enriched in some categories, for example cuticular proteins and genes associated with cell-cell signaling iii) leaf-galling-biased genes were enriched in genes associated with the nucleus and DNA-replication, suggesting a metabolism more oriented towards fast and active multiplication. We also identified a gene family with a very high expression level (copies totaling nearly 10% of the reads) in the grapevine phylloxera (both in root and leaf galling forms), but usually expressed at very low levels in true aphids (except in sexual oviparous females). These transcripts thus appear to be associated with oviparity., Conclusions: Our study illustrated major intraspecific changes in transcriptome profiles, related with different life-styles (and the feeding on roots versus in leaf-galls). At a different scale, we could also illustrate one major shift in expression levels associated with changes in life-traits that occurred along evolution and that respectively characterize (strictly oviparous) grapevine phylloxera and (mostly viviparous) true aphids.

Background: Neuropeptides are key players in information transfer and act as important regulators of development, growth, metabolism, and reproduction within multi-cellular animal organisms (Metazoa). These short protein-like substances show a high degree of structural variability and are recognized as the most diverse group of messenger molecules. We used transcriptome sequences from the 1KITE (1K Insect Transcriptome Evolution) project to search for neuropeptide coding sequences in 24 species from the non-pterygote hexapod lineages Protura (coneheads), Collembola (springtails), Diplura (two-pronged bristletails), Archaeognatha (jumping bristletails), and Zygentoma (silverfish and firebrats), which are often referred to as "basal" hexapods. Phylogenetically, Protura, Collembola, Diplura, and Archaeognatha are currently placed between Remipedia and Pterygota (winged insects); Zygentoma is the sistergroup of Pterygota. The Remipedia are assumed to be among the closest relatives of all hexapods and belong to the crustaceans., Results: We identified neuropeptide precursor sequences within whole-body transcriptome data from these five hexapod groups and complemented this dataset with homologous sequences from three crustaceans (including Daphnia pulex), three myriapods, and the fruit fly Drosophila melanogaster. Our results indicate that the reported loss of several neuropeptide genes in a number of winged insects, particularly holometabolous insects, is a trend that has occurred within Pterygota. The neuropeptide precursor sequences of the non-pterygote hexapods show numerous amino acid substitutions, gene duplications, variants following alternative splicing, and numbers of paracopies. Nevertheless, most of these features fall within the range of variation known from pterygote insects. However, the capa/pyrokinin genes of non-pterygote hexapods provide an interesting example of rapid evolution, including duplication of a neuropeptide gene encoding different ligands., Conclusions: Our findings delineate a basic pattern of neuropeptide sequences that existed before lineage-specific developments occurred during the evolution of pterygote insects.

Background: The Mesopsychidae is an extinct family of Mecoptera, comprising eleven described genera from Upper Permian to Lower Cretaceous deposits. In 2009, several well-preserved mesopsychids with long proboscides were reported from the mid Mesozoic of Northeastern China, suggesting the presence of pollination mutualisms with gymnosperm plants and highlighting their elevated genus-level diversity. Since that time, additional mesopsychid taxa have been described. However, the phylogeny of genera within Mesopsychidae has not been studied formally, attributable to the limited number of well-preserved fossils., Results: Here, we describe two new species, Lichnomesopsyche prochorista sp. nov. and Vitimopsyche pristina sp. nov. and revise the diagnosis of Lichnomesopsyche daohugouensis Ren, Labandeira and Shih, 2010, based on ten specimens from the latest Middle Jurassic Jiulongshan Formation of Inner Mongolia, China. After compiling data from these new fossil species and previously reported representative taxa, we conducted phylogenetic analyses and geometric morphometric studies that now shed light on the taxonomy and phylogeny of Mesopsychidae. We also evaluate the recurring origin of the siphonate proboscis in the Mecoptera and propose an evolutionary developmental model for its multiple origins., Conclusions: Phylogenetic and geometric morphometric results confirm the establishment of two new species, each to Lichnomesopsyche and Vitimopsyche. Vitimopsyche pristina sp. nov. extends the existence of the genus Vitimopsyche Novokshonov and Sukacheva, 2001, from the mid Lower Cretaceous to the latest Middle Jurassic. Two methods of analyses indicate an affiliation of Mesopsyche dobrokhotovae Novokshonov, 1997 with Permopsyche Bashkuev, 2011. A phylogenetic analysis of the Mesopsychidae supports: 1), Mesopsychidae as a monophyletic group; 2), Mesopsyche as a paraphyletic group, to be revised pending future examination of additional material; and 3), the independent origin of the proboscis in the Pseudopolycentropodidae, its subsequent loss in earliest Mesopsychidae such as Epicharmesopsyche, its re-origination in the common ancestor (or perhaps independently) in the Vitimopsyche and Lichnomesopsyche clades of the Mesopsychidae. The third conclusion indicates that the proboscis originated four or five times within early Mecoptera, whose origin is explained by an evolutionary developmental model.

Background: Phylogeographic studies of aquatic insects provide valuable insights into mechanisms that shape the genetic structure of communities, yet studies that include broad geographic areas are uncommon for this group. We conducted a broad scale phylogeographic analysis of the least salmonfly Pteronarcella badia (Plecoptera) across western North America. We tested hypotheses related to mode of dispersal and the influence of historic climate oscillations on population genetic structure. In order to generate a larger mitochondrial data set, we used 454 sequencing to reconstruct the complete mitochondrial genome in the early stages of the project., Results: Our analysis revealed high levels of population structure with several deeply divergent clades present across the sample area. Evidence from five mitochondrial genes and one nuclear locus identified a potentially cryptic lineage in the Pacific Northwest. Gene flow estimates and geographic clade distributions suggest that overland flight during the winged adult stage is an important dispersal mechanism for this taxon. We found evidence of multiple glacial refugia across the species distribution and signs of secondary contact within and among major clades., Conclusions: This study provides a basis for future studies of aquatic insect phylogeography at the inter-basin scale in western North America. Our findings add to an understanding of the role of historical climate isolations in shaping assemblages of aquatic insects in this region. We identified several geographic areas that may have historical importance for other aquatic organisms with similar distributions and dispersal strategies as P. badia. This work adds to the ever-growing list of studies that highlight the potential of next-generation DNA sequencing in a phylogenetic context to improve molecular data sets from understudied groups.

Background: Trachipleistophora hominis was isolated from an HIV/AIDS patient and is a member of a highly successful group of obligate intracellular parasites., Methods: Here we have investigated the evolution of the parasite and the interplay between host and parasite gene expression using transcriptomics of T. hominis-infected rabbit kidney cells., Results: T. hominis has about 30% more genes than small-genome microsporidians. Highly expressed genes include those involved in growth, replication, defence against oxidative stress, and a large fraction of uncharacterised genes. Chaperones are also highly expressed and may buffer the deleterious effects of the large number of non-synonymous mutations observed in essential T. hominis genes. Host expression suggests a general cellular shutdown upon infection, but ATP, amino sugar and nucleotide sugar production appear enhanced, potentially providing the parasite with substrates it cannot make itself. Expression divergence of duplicated genes, including transporters used to acquire host metabolites, demonstrates ongoing functional diversification during microsporidian evolution. We identified overlapping transcription at more than 100 loci in the sparse T. hominis genome, demonstrating that this feature is not caused by genome compaction. The detection of additional transposons of insect origin strongly suggests that the natural host for T. hominis is an insect., Conclusions: Our results reveal that the evolution of contemporary microsporidian genomes is highly dynamic and innovative. Moreover, highly expressed T. hominis genes of unknown function include a cohort that are shared among all microsporidians, indicating that some strongly conserved features of the biology of these enormously successful parasites remain uncharacterised.

Background: Identification of biological specimens is a requirement for a range of applications. Reference-free methods analyse unprocessed sequencing data without relying on prior knowledge, but generally do not scale to arbitrarily large genomes and arbitrarily large phylogenetic distances., Results: We present Cnidaria, a practical tool for clustering genomic and transcriptomic data with no limitation on genome size or phylogenetic distances. We successfully simultaneously clustered 169 genomic and transcriptomic datasets from 4 kingdoms, achieving 100% identification accuracy at supra-species level and 78% accuracy at the species level., Conclusion: CNIDARIA allows for fast, resource-efficient comparison and identification of both raw and assembled genome and transcriptome data. This can help answer both fundamental (e.g. in phylogeny, ecological diversity analysis) and practical questions (e.g. sequencing quality control, primer design).

Background: Phylogenetic footprinting is a comparative method based on the principle that functional sequence elements will acquire fewer mutations over time than non-functional sequences. Successful comparisons of distantly related species will thus yield highly important sequence elements likely to serve fundamental biological roles. RNA regulatory elements are less well understood than those in DNA. In this study we use the emerging model organism Nasonia vitripennis, a parasitic wasp, in a comparative analysis against 12 insect genomes to identify deeply conserved non-coding elements (CNEs) conserved in large groups of insects, with a focus on 5' UTRs and promoter sequences., Results: We report the identification of 322 CNEs conserved across a broad range of insect orders. The identified regions are associated with regulatory and developmental genes, and contain short footprints revealing aspects of their likely function in translational regulation. The most ancient regions identified in our analysis were all found to overlap transcribed regions of genes, reflecting stronger conservation of translational regulatory elements than transcriptional elements. Further expanding sequence analyses to non-insect species we also report the discovery of, to our knowledge, the two oldest and most ubiquitous CNE's yet described in the animal kingdom (700 MYA). These ancient conserved non-coding elements are associated with the two ribosomal stalk genes, RPLP1 and RPLP2, and were very likely functional in some of the earliest animals., Conclusions: We report the identification of the most deeply conserved CNE's found to date, and several other deeply conserved elements which are without exception, part of 5' untranslated regions of transcripts, and occur in a number of key translational regulatory genes, highlighting translational regulation of translational regulators as a conserved feature of insect genomes.

Background: Holometabolous insects are the most diverse, speciose and ubiquitous group of multicellular organisms in terrestrial and freshwater ecosystems. The enormous evolutionary and ecological success of Holometabola has been attributed to their unique postembryonic life phases in which nonreproductive and wingless larvae differ significantly in morphology and life habits from their reproductive and mostly winged adults, separated by a resting stage, the pupa. Little is known of the evolutionary developmental mechanisms that produced the holometabolous larval condition and their Paleozoic origin based on fossils and phylogeny., Results: We provide a detailed anatomic description of a 311 million-year-old specimen, the oldest known holometabolous larva, from the Mazon Creek deposits of Illinois, U.S.A. The head is ovoidal, downwardly oriented, broadly attached to the anterior thorax, and bears possible simple eyes and antennae with insertions encircled by molting sutures; other sutures are present but often indistinct. Mouthparts are generalized, consisting of five recognizable segments: a clypeo-labral complex, mandibles, possible hypopharynx, a maxilla bearing indistinct palp-like appendages, and labium. Distinctive mandibles are robust, triangular, and dicondylic. The thorax is delineated into three, nonoverlapping regions of distinctive surface texture, each with legs of seven elements, the terminal-most bearing paired claws. The abdomen has ten segments deployed in register with overlapping tergites; the penultimate segment bears a paired, cercus-like structure. The anterior eight segments bear clawless leglets more diminutive than the thoracic legs in length and cross-sectional diameter, and inserted more ventrolaterally than ventrally on the abdominal sidewall., Conclusions: Srokalarva berthei occurred in an evolutionary developmental context likely responsible for the early macroevolutionary success of holometabolous insects. Srokalarva berthei bore head and prothoracic structures, leglet series on successive abdominal segments - in addition to comparable features on a second taxon eight million-years-younger - that indicates Hox-gene regulation of segmental and appendage patterning among earliest Holometabola. Srokalarva berthei body features suggest a caterpillar-like body plan and head structures indicating herbivory consistent with known, contemporaneous insect feeding damage on seed plants. Taxonomic resolution places Srokalarva berthei as an extinct lineage, apparently possessing features closer to neuropteroid than other holometabolous lineages.

Background: The study of islands as model systems plays a key role in understanding many evolutionary processes. Knowledge of the historical events leading to present-day island communities is pivotal for exploring fundamental mechanisms of speciation and adaptation. The remote Mascarene archipelago (Mauritius, Réunion, Rodrigues), considered to be the product of an age-progressive trend of north-to-south volcanic activity in the Indian Ocean, hosts a remarkably diverse, endemic and threatened concentration of flora and fauna that has traditionally been considered to be biogeographically related to Madagascar and Africa. To explore the evolutionary diversity of the Mascarene stick insects (Phasmatodea), we constructed a global phylogeny from approximately 2.4 kb of mitochondrial and nuclear sequence data of more than 120 species representing all major phasmatodean lineages., Results: Based on the obtained time-calibrated molecular tree we demonstrate that the current phasmid community of the Mascarene archipelago, which consists of members of four presumably unrelated traditional subfamilies, is the result of a single ancient dispersal event from Australasia and started radiating between 16-29 million years ago, significantly predating the age of Mauritius (8-10 million years)., Conclusions: We propose that the Mascarene stick insects diversified on landmasses now eroded away, presumably to the north of Mauritius. In consequence, ancient islands have probably persisted in the Indian Ocean until the emergence of Mauritius and not only served as stepping stones for colonisation events during sea-level lowstands, but as long-lasting cradles of evolution. These ancient landmasses most likely allowed for adaptive speciation and served as significant sources of diversity that contributed to the biomes of the Mascarene archipelago and the megadiverse Madagascar.

Background: Insects belong to a class that accounts for the majority of animals on earth. With over one million identified species, insects display a huge diversity and occupy extreme environments. At present, there are dozens of fully sequenced insect genomes that cover a range of habitats, social behavior and morphologies. In view of such diverse collection of genomes, revealing evolutionary trends and charting functional relationships of proteins remain challenging., Results: We analyzed the relatedness of 17 complete proteomes representative of proteomes from insects including louse, bee, beetle, ants, flies and mosquitoes, as well as an out-group from the crustaceans. The analyzed proteomes mostly represented the orders of Hymenoptera and Diptera. The 287,405 protein sequences from the 18 proteomes were automatically clustered into 20,933 families, including 799 singletons. A comprehensive analysis based on statistical considerations identified the families that were significantly expanded or reduced in any of the studied organisms. Among all the tested species, ants are characterized by an exceptionally high rate of family gain and loss. By assigning annotations to hundreds of species-specific families, the functional diversity among species and between the major clades (Diptera and Hymenoptera) is revealed. We found that many species-specific families are associated with receptor signaling, stress-related functions and proteases. The highest variability among insects associates with the function of transposition and nucleic acids processes (collectively coined TNAP). Specifically, the wasp and ants have an order of magnitude more TNAP families and proteins relative to species that belong to Diptera (mosquitoes and flies)., Conclusions: An unsupervised clustering methodology combined with a comparative functional analysis unveiled proteomic signatures in the major clades of winged insects. We propose that the expansion of TNAP families in Hymenoptera potentially contributes to the accelerated genome dynamics that characterize the wasp and ants.

We previously characterized the β-actin gene promoter of Indian domesticated rohu carp (Labeo rohita) and made a reporter construct via fusion to green fluorescence protein (GFP) cDNA. In this study, the same construct was used to breed transgenic rohu fish. About 20% of the transgenic offspring showed ubiquitous expression of the reporter GFP gene. In a few of the transgenic fish, we documented massive epithelial and/or muscular expression with visible green color under normal light. The expression of GFP mRNA was higher in the muscle tissue of transgenic fish than in that of non-transgenic fish. A highly efficient nucleofection protocol was optimized to transfect proliferating spermatogonial cells of rohu using this reporter construct. The β-actin promoter also drove expressions in HEK293 (derived from human embryonic kidney cells), K562 (human leukemic cells) and SF21 (insect ovarian cells) lines. These findings imply conserved regulatory mechanisms of β-actin gene expression across eukaryotes. Furthermore, the isolated β-actin promoter with consensus regulatory elements has the potential to be used in generating transgenic carp with genes of interest and in basic biology research.

Background: Several methods exist for the prediction of precursor miRNAs (pre-miRNAs) in genomic or sRNA-seq (small RNA sequences) data produced by NGS (Next Generation Sequencing). One key information used for this task is the characteristic hairpin structure adopted by pre-miRNAs, that in general are identified using RNA folders whose complexity is cubic in the size of the input. The vast majority of pre-miRNA predictors then rely on further information learned from previously validated miRNAs from the same or a closely related genome for the final prediction of new miRNAs. With this paper, we wished to address three main issues. The first was methodological and aimed at obtaining a more time-efficient predictor, however without losing in accuracy which represented a second issue. We indeed aimed at better predicting miRNAs at a genome scale, but also from sRNAseq data where in some cases, notably of plants, the current folding methods often infer the wrong structure. The third issue is related to the fact that it is important to rely as little as possible on previously recorded examples of miRNAs. We therefore also sought a method that is less dependent on previous miRNA records., Results: As concerns the first and second issues, we present a novel alternative to a classical folder based on a thermodynamic Nearest-Neighbour (NN) model for computing the free energy and predicting the classical hairpin structure of a pre-miRNA. We show that the free energies thus computed correlate well with those of RNAFOLD. This novel method, called MIRINHO, has quadratic instead of cubic complexity and is much more efficient also in practice. When applied to sRNAseq data of plants, it gives in general better results than classical folders. On the third issue, we show that MIRINHO, which uses as only knowledge the length of the loops and stem-arms and the free energy of the pre-miRNA hairpin, compares well with algorithms that require more information. The results, obtained with different datasets, are indeed similar to those of other approaches with which such a comparison was possible. These needed to be publicly available softwares that could be used on a large input. In some cases, MIRINHO is even better in terms of sensitivity or precision., Conclusion: We provide a simpler and much faster method with very reasonable sensitivity and precision, which can be applied without special adaptation to the prediction of both animal and plant pre-miRNAs, using as input either genomic sequences or sRNA-seq data.

Background: Social animals have the unique capability of mounting social defenses against pathogens. Over the last decades, social immunity has been extensively studied in species with obligatory and permanent forms of social life. However, its occurrence in less derived social systems and thus its role in the early evolution of group-living remains unclear. Here, we investigated whether lining nests with feces is a form of social immunity against microbial growth in the European earwig Forficula auricularia, an insect with temporary family life and facultative maternal care., Results: Using a total of 415 inhibition zone assays, we showed that earwig feces inhibit the growth of two GRAM+ bacteria, two fungi, but not of a GRAM- bacteria. These inhibitions did not result from the consumed food or the nesting environment. We then demonstrated that the antimicrobial activity against fungus was higher in offspring than maternal feces, but that this difference was absent against bacteria. Finally, we showed that family interactions inhibited the antibacterial activity of maternal feces against one of the two GRAM+ bacteria, whereas it had no effect on the one of nymphal feces. By contrast, antifungal activities of the feces were independent of mother-offspring interactions., Conclusion: These results demonstrate that social immunity occurs in a species with simple and facultative social life, and thus shed light on the general importance of this process in the evolution of group-living. These results also emphasize that defecation can be under selection for other life-history traits than simple waste disposal.

Background: Ophiocordyceps sinensis, a worm and fungus combined mixture which Hirsutella sinensis is parasitic on the caterpillar body, has been used as a traditional medicine or healthy food in China for thousands of years. H. sinensis is reported as the only correct anamorph of O. sinensis and its main active ingredients are similar to the natural O. sinensis., Results: H. sinensis L0106, asexual strain of O. sinensis, was isolated and identified in this study. Three transcriptomes of H. sinensis at different cultivation periods (growth period 3d, pre-stable period 6d and stable period 9d) were sequenced for the first time by RNA-Seq method, and 25,511 unigenes (3d), 25,214 unigenes (6d) and 16,245 unigenes (9d) were assembled and obtained, respectively. These unigenes of the three samples were further assembled into 20,822 unigenes (All), and 62.3 percent of unigenes (All) could be annotated based on protein databases. Subsequently, the genes and enzymes involved in the biosynthesis of the active ingredients according to the sequencing and annotation results were predicted. Based on the predictions, we further investigated the interaction of different pathway networks and the corresponding enzymes. Furthermore, the differentially expressed genes (DEGs) of H. sinensis grown during different developmental stages (3d-VS-6d, 3d-VS-9d and 6d-VS-9d) were globally detected and analyzed based on the data from RNA-Seq, and 764 DEGs between 3d and 6d, 1,869 DEGs between 3d and 9d, and 770 DEGs between 6d and 9d were found, respectively., Conclusions: This work presented here would aid in understanding and carrying out future studies on the genetic basis of H. sinensis and contribute to the further artificial production and application of this organism. This study provided a substantial contribution and basis to further characterize the gene expression profiles of H. sinensis in the metabolic pathways of active ingredients.

Background: Zoraptera, generally regarded as a member of Polyneoptera, represents one of the most enigmatic insect orders. Although phylogenetic analyses based on a wide array of morphological and/or nuclear data have been performed, the position of Zoraptera is still under debate. Mitochondrial genome (mitogenome) information is commonly considered to be preferable to reconstruct phylogenetic relationships, but no efforts have been made to incorporate it in Zorapteran phylogeny. To characterize Zoraptera mitogenome features and provide insights into its phylogenetic placement, here we sequenced, for the first time, one complete mitogenome of Zoraptera and reconstructed the phylogeny of Polyneoptera., Results: The mitogenome of Zorotypus medoensis with an A+T content of 72.50% is composed of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a noncoding A+T-rich region. The gene content and arrangement are identical to those considered ancestral for insects. This mitogenome shows a number of very unusual features. First, it is very compact, comprising 14,572 bp, and is the smallest among all known polyneopteran mitogenomes. Second, both noncoding sequences and coding genes exhibit a significant decrease in size compared with those of other polyneopterans. Third, Z. medoensis mitogenome has experienced an accelerated substitution rate. Fourth, truncated secondary structures of tRNA genes occur with loss of dihydrouridine (DHU) arm in trnC, trnR, and trnS(AGN) and loss of TΨC arm in trnH and trnT. The phylogenetic analyses based on the mitogenome sequence information indicate that Zoraptera, represented by Z. medoensis, is recovered as sister to Embioptera. However, both Zoraptera and Embioptera exhibit very long branches in phylogenetic trees., Conclusions: Characterization of Z. medoensis mitogenome contributes to our understanding of the enigmatic Zoraptera. Mitogenome data demonstrate an overall strong resolution of deep-level phylogenies of Polyneoptera but not Insecta. It is preferable to expand taxon sampling of Zoraptera and other poorly represented orders in future to break up long branches.

Background: Ecological studies routinely show genotype-genotype interactions between insects and their parasites. The mechanisms behind these interactions are not clearly understood. Using the bumblebee Bombus terrestris/trypanosome Crithidia bombi model system (two bumblebee colonies by two Crithidia strains), we have carried out a transcriptome-wide analysis of gene expression and alternative splicing in bees during C. bombi infection. We have performed four analyses, 1) comparing gene expression in infected and non-infected bees 24 hours after infection by Crithidia bombi, 2) comparing expression at 24 and 48 hours after C. bombi infection, 3) determining the differential gene expression associated with the bumblebee-Crithidia genotype-genotype interaction at 24 hours after infection and 4) determining the alternative splicing associated with the bumblebee-Crithidia genotype-genotype interaction at 24 hours post infection., Results: We found a large number of genes differentially regulated related to numerous canonical immune pathways. These genes include receptors, signaling pathways and effectors. We discovered a possible interaction between the peritrophic membrane and the insect immune system in defense against Crithidia. Most interestingly, we found differential expression and alternative splicing of immunoglobulin related genes (Dscam and Twitchin) are associated with the genotype-genotype interactions of the given bumblebee colony and Crithidia strain., Conclusions: In this paper we have shown that the expression and alternative splicing of immune genes is associated with specific interactions between different host and parasite genotypes in this bumblebee/trypanosome model.

Background: As attested by the fossil record, Cretaceous environmental changes have significantly impacted the diversification dynamics of several groups of organisms. A major biome turnover that occurred during this period was the rise of angiosperms starting ca. 125 million years ago. Though there is evidence that the latter promoted the diversification of phytophagous insects, the response of other insect groups to Cretaceous environmental changes is still largely unknown. To gain novel insights on this issue, we assess the diversification dynamics of a hyperdiverse family of detritivorous beetles (Tenebrionidae) using molecular dating and diversification analyses., Results: Age estimates reveal an origin after the Triassic-Jurassic mass extinction (older than previously thought), followed by the diversification of major lineages during Pangaean and Gondwanan breakups. Dating analyses indicate that arid-adapted species diversified early, while most of the lineages that are adapted to more humid conditions diversified much later. Contrary to other insect groups, we found no support for a positive shift in diversification rates during the Cretaceous; instead there is evidence for an 8.5-fold increase in extinction rates that was not compensated by a joint increase in speciation rates., Conclusions: We hypothesize that this pattern is better explained by the concomitant reduction of arid environments starting in the mid-Cretaceous, which likely negatively impacted the diversification of arid-adapted species that were predominant at that time.

Background: Stick and leaf insects (Phasmatodea) are an exclusively leaf-feeding order of insects with no record of omnivory, unlike other "herbivorous" Polyneoptera. They represent an ideal system for investigating the adaptations necessary for obligate folivory, including plant cell wall degrading enzymes (PCWDEs). However, their physiology and internal anatomy is poorly understood, with limited genomic resources available., Results: We de novo assembled transcriptomes for the anterior and posterior midguts of six diverse Phasmatodea species, with RNA-Seq on one exemplar species, Peruphasma schultei. The latter's assembly yielded >100,000 transcripts, with over 4000 transcripts uniquely or more highly expressed in specific midgut sections. Two to three dozen PCWDE encoding gene families, including cellulases and pectinases, were differentially expressed in the anterior midgut. These genes were also found in genomic DNA from phasmid brain tissue, suggesting endogenous production. Sequence alignments revealed catalytic sites on most PCWDE transcripts. While most phasmid PCWDE genes showed homology with those of other insects, the pectinases were homologous to bacterial genes., Conclusions: We identified a large and diverse PCWDE repertoire endogenous to the phasmids. If these expressed genes are translated into active enzymes, then phasmids can theoretically break plant cell walls into their monomer components independently of microbial symbionts. The differential gene expression between the two midgut sections provides the first molecular hints as to their function in living phasmids. Our work expands the resources available for industrial applications of animal-derived PCWDEs, and facilitates evolutionary analysis of lower Polyneopteran digestive enzymes, including the pectinases whose origin in Phasmatodea may have been a horizontal transfer event from bacteria.

Background: The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes., Results: To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes., Conclusions: Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

Background: The Kalligrammatidae are distinctive, large, conspicuous, lacewings found in Eurasia from the Middle Jurassic to mid Early Cretaceous. Because of incomplete and often inadequate fossil preservation, an absence of detailed morphology, unclear relationships, and unknown evolutionary trends, the Kalligrammatidae are poorly understood., Results: We describe three new subfamilies, four new genera, twelve new species and four unassigned species from the late Middle Jurassic Jiulongshan and mid Early Cretaceous Yixian Formations of China. These kalligrammatid taxa exhibit diverse morphological characters, such as mandibulate mouthparts in one major clade and siphonate mouthparts in the remaining four major clades, the presence or absence of a variety of distinctive wing markings such as stripes, wing spots and eyespots, as well as multiple major wing shapes. Based on phylogenetic analyses, the Kalligrammatidae are divided into five principal clades: Kalligrammatinae Handlirsch, 1906, Kallihemerobiinae Ren & Engel, 2008, Meioneurinae subfam. nov., Oregrammatinae subfam. nov. and Sophogrammatinae subfam. nov., each of which is accorded subfamily-level status. Our results show significant morphological and evolutionary differentiation of the Kalligrammatidae family during a 40 million-year-interval of the mid Mesozoic., Conclusion: A new phylogeny and classification of five subfamilies and their constituent genera is proposed for the Kalligrammatidae. These diverse, yet highly specialized taxa from northeastern China suggest that eastern Eurasia likely was an important diversification center for the Kalligrammatidae. Kalligrammatids possess an extraordinary morphological breadth and panoply of adaptations during the mid-Mesozoic that highlight our conclusion that their evolutionary biology is much more complex than heretofore realized.

Background: Raphidioptera (snakeflies) is a holometabolous order of the superorder Neuropterida characterized by the narrowly elongate adult prothorax and the long female ovipositor. Mesozoic snakeflies were markedly more diverse than the modern ones are. However, the evolutionary history of Raphidioptera is largely unexplored, as a result of the poorly studied phylogeny among fossil and extant lineages within the order., Results: In this paper, we report a new snakefly family, Juroraphidiidae fam. nov., based on exquisitely preserved fossils, attributed to a new species Juroraphidia longicollumgen. et sp. nov., from the Jiulongshan Formation (Middle Jurassic) in Inner Mongolia, China. The new family is characterized by an unexpected combination of plesiomorphic and apomorphic characters of Raphidioptera. Based on our phylogenetic analysis, Juroraphidiidae fam. nov. together with Raphidiomorpha form a monophyletic clade, which is the sister to Priscaenigmatomorpha. The snakefly affinity of Priscaenigmatomorpha is confirmed and another new family, Chrysoraphidiidae fam. nov., is erected in this suborder., Conclusions: Juroraphidiidae fam. nov. is determined to be a transitional lineage between Priscaenigmatomorpha and Raphidiomorpha. Diversification of higher snakefly taxa had occurred by the Early Jurassic, suggesting that these insects had already had a long but undocumented history by this time.

Background: Despite considerable progress in systematics, a comprehensive scenario of the evolution of phenotypic characters in the mega-diverse Holometabola based on a solid phylogenetic hypothesis was still missing. We addressed this issue by de novo sequencing transcriptome libraries of representatives of all orders of holometabolan insects (13 species in total) and by using a previously published extensive morphological dataset. We tested competing phylogenetic hypotheses by analyzing various specifically designed sets of amino acid sequence data, using maximum likelihood (ML) based tree inference and Four-cluster Likelihood Mapping (FcLM). By maximum parsimony-based mapping of the morphological data on the phylogenetic relationships we traced evolutionary transformations at the phenotypic level and reconstructed the groundplan of Holometabola and of selected subgroups., Results: In our analysis of the amino acid sequence data of 1,343 single-copy orthologous genes, Hymenoptera are placed as sister group to all remaining holometabolan orders, i.e., to a clade Aparaglossata, comprising two monophyletic subunits Mecopterida (Amphiesmenoptera + Antliophora) and Neuropteroidea (Neuropterida + Coleopterida). The monophyly of Coleopterida (Coleoptera and Strepsiptera) remains ambiguous in the analyses of the transcriptome data, but appears likely based on the morphological data. Highly supported relationships within Neuropterida and Antliophora are Raphidioptera + (Neuroptera + monophyletic Megaloptera), and Diptera + (Siphonaptera + Mecoptera). ML tree inference and FcLM yielded largely congruent results. However, FcLM, which was applied here for the first time to large phylogenomic supermatrices, displayed additional signal in the datasets that was not identified in the ML trees., Conclusions: Our phylogenetic results imply that an orthognathous larva belongs to the groundplan of Holometabola, with compound eyes and well-developed thoracic legs, externally feeding on plants or fungi. Ancestral larvae of Aparaglossata were prognathous, equipped with single larval eyes (stemmata), and possibly agile and predacious. Ancestral holometabolan adults likely resembled in their morphology the groundplan of adult neopteran insects. Within Aparaglossata, the adult's flight apparatus and ovipositor underwent strong modifications. We show that the combination of well-resolved phylogenies obtained by phylogenomic analyses and well-documented extensive morphological datasets is an appropriate basis for reconstructing complex morphological transformations and for the inference of evolutionary histories.