Draft genome and reference transcriptomic resources for the urticating pine defoliator Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

Large-scale genomic data reveal the phylogeny and evolution of owlet moths (Noctuoidea)

The owlet moths (Noctuoidea; ~43-45K described species) are one of the most ecologically diverse and speciose superfamilies of animals. Moreover, they comprise some of the world's most notorious pests of agriculture and forestry. Despite their contributions to terrestrial biodiversity and impacts on ecosystems and economies, the evolutionary history of Noctuoidea remains unclear because the superfamily lacks a statistically robust phylogenetic and temporal framework. We reconstructed the phylogeny of Noctuoidea using data from 1234 genes (946.4 kb nucleotides) obtained from the genome and transcriptome sequences of 76 species. The relationships among the six families of Noctuoidea were well resolved and consistently recovered based on both concatenation and gene coalescence approaches, supporting the following relationships: Oenosandridae + (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))). A Yule tree prior with three unlinked molecular clocks was identified as the preferred BEAST analysis using marginal-likelihood estimations. The crown age of Noctuoidea was estimated at 74.5 Ma, with most families originating before the end of the Paleogene (23 Ma). Our study provides the first statistically robust phylogenetic and temporal framework for Noctuoidea, including all families of owlet moths, based on large-scale genomic data.

Sequence, assembly and count datasets of viruses associated to the pine processionary moth Thaumetopoea pityocampa (Denis & Schiffermüller) (Lepidoptera, Notodontidae) identified from transcriptomic high-throughput sequencing

The pine processionary moth Thaumetopoea pityocampa is a Lepidopteran pest species occurring in the Western Mediterranean. It causes heavy pine defoliations and it is a public and animal health concern because of its urticating caterpillars. Very little is known about the viruses associated to this species, as only two viruses were described so far. We here present a dataset corresponding to 34 viral transcripts, among which 27 could be confidently assigned to 9 RNA and DNA viral families (Iflaviridae, Reoviridae, Partitiviridae, Permutotetraviridae, Flaviviridae, Rhabdoviridae, Parvoviridae, Baculoviridae and PolyDNAviridae). These transcripts were identified from an original transcriptome assembled for the insect host, using both blast search and phylogenetic approaches. The data were acquired from 2 populations in Portugal and 2 populations in Italy. The transcripts were de novo assembled and used to identify viral sequences by homology searches. We also provide information about the populations and life stages in which each virus was identified. The data produced will allow to enrich the virus taxonomy in Lepidopteran hosts, and to develop PCR-based diagnostic tools to screen colonies across the range and determine the distribution and prevalence of the identified viral species.

SNP4OrphanSpecies: A bioinformatics pipeline to isolate molecular markers for studying genetic diversity of orphan species

Background

For several decades, an increase in disease or pest emergences due to anthropogenic introduction or environmental changes has been recorded. This increase leads to serious threats to the genetic and species diversity of numerous ecosystems. Many of these events involve species with poor or no genomic resources (called here "orphan species"). This lack of resources is a serious limitation to our understanding of the origin of emergent populations, their ability to adapt to new environments and to predict future consequences to biodiversity. Analyses of genetic diversity are an efficient method to obtain this information rapidly, but require available polymorphic genetic markers.

New information

We developed a generic bioinformatics pipeline to rapidly isolate such markers with the goal for the pipeline to be applied in studies of invasive taxa from different taxonomic groups, with a special focus on forest fungal pathogens and insect pests. This pipeline is based on: 1) an automated de novo genome assembly obtained from shotgun whole genome sequencing using paired-end Illumina technology; 2) the isolation of single-copy genes conserved in species related to the studied emergent organisms; 3) primer development for multiplexed short sequences obtained from these conserved genes. Previous studies have shown that intronic regions of these conserved genes generally contain several single nucleotide polymorphisms within species. The pipeline's functionality was evaluated with sequenced genomes of five invasive or expanding pathogen and pest species in Europe (Armillariaostoyae (Romagn.) Herink 1973, Bursaphelenchusxylophilus Steiner & Buhrer 1934, Sphaeropsissapinea (fr.) Dicko & B. Sutton 1980, Erysiphealphitoides (Griffon & Maubl.) U. Braun & S. Takam. 2000, Thaumetopoeapityocampa Denis & Schiffermüller, 1775). We successfully isolated several pools of one hundred short gene regions for each assembled genome, which can be amplified in multiplex. The bioinformatics pipeline is user-friendly and requires little computational resources. This easy-to-set-up and run method for genetic marker identification will be useful for numerous laboratories studying biological invasions, but with limited resources and expertise in bioinformatics.

Complete mitogenome data from a summer population specimen of the urticating pine defoliator Thaumetopoea pityocampa (Lepidoptera, Notodontidae, Thaumetopoeinae, Thaumetopoea)

We present a de novo mitogenome assembly obtained from specimens sampled in the so-called summer population (SP) of Thaumetopoea pityocampa (Denis and Schiffermüller, 1775) in Portugal. Contrary to the typical larval development occurring in winter in this species, the larvae of this unique population develop during summer. The sequencing data used were obtained from genomic libraries originally generated to assemble the nuclear genome of T. pityocampa[1]. We also provide a complete annotation and a phylogenetic representation which positions the Portuguese summer population of T. pityocampa and an Italian typical individual of the same species among the Notodontidae family and more distant Noctuoidea species. This data represents a valuable new resource for an expanding and urticating insect pest.

Expanding the Menu: Are Polyphagy and Gene Family Expansions Linked across Lepidoptera?

Evolutionary expansions and contractions of gene families are often correlated with key innovations and/or ecological characteristics. In butterflies and moths (Lepidoptera), expansions of gene families involved in detoxification of plant specialized metabolites are hypothesized to facilitate a polyphagous feeding style. However, analyses supporting this hypothesis are mostly based on a limited number of lepidopteran species. We applied a phylogenomics approach, using 37 lepidopteran genomes, to analyze if gene family evolution (gene gain and loss) is associated with the evolution of polyphagy. Specifically, we compared gene counts and evolutionary gene gain and loss rates of gene families involved in adaptations with plant feeding. We correlated gene evolution to host plant family range (phylogenetic diversity) and specialized metabolite content of plant families (functional metabolite diversity). We found a higher rate for gene loss than gene gain in Lepidoptera, a potential consequence of genomic rearrangements and deletions after (potentially small-scale) duplication events. Gene family expansions and contractions varied across lepidopteran families, and were associated to host plant use and specialization levels. Within the family Noctuidae, a higher expansion rate for gene families involved in detoxification can be related to the large number of polyphagous species. However, gene family expansions are observed in both polyphagous and monophagous lepidopteran species and thus seem to be species-specific in the taxa sampled. Nevertheless, a significant positive correlation of gene counts of the carboxyl- and choline esterase and glutathione-S-transferase detoxification gene families with the level of polyphagy was identified across Lepidoptera.

Coding-Complete Genome Sequence of a Partitivirus Isolated from Pine Processionary Moth Eggs

Two coding-complete nucleotide sequences of a partitivirus (family Partitiviridae) were discovered in transcriptomic data sets obtained from eggs of the Lepidoptera Thaumetopoea pityocampa. Each segment encodes a single open reading frame, and these two segments are predicted to encode an RNA-dependent RNA polymerase and a coat protein, respectively.

Two coding-complete nucleotide sequences of a partitivirus (family Partitiviridae) were discovered in transcriptomic data sets obtained from eggs of the Lepidoptera Thaumetopoea pityocampa. Each segment encodes a single open reading frame, and these two segments are predicted to encode an RNA-dependent RNA polymerase and a coat protein, respectively.

Coding-Complete Genome Sequences of an Iteradensovirus and an Alphapermutotetra-Like Virus Identified from the Pine Processionary Moth (Thaumetopoea pityocampa) in Portugal

The coding-complete nucleotide sequences of an iteradensovirus (family Parvoviridae) and an alphapermutotetra-like virus (family Permutotetraviridae) were discovered from transcriptomic data sets obtained from Thaumetopoea pityocampa larvae collected in Portugal. Each of the coding-complete genome sequences of these viruses contain three main open reading frames (ORFs).

The coding-complete genome sequences of an iteradensovirus (family Parvoviridae) and an alphapermutotetra-like virus (family Permutotetraviridae) were discovered from transcriptomic data sets obtained from Thaumetopoea pityocampa larvae collected in Portugal. Each of the coding-complete genome sequences of these viruses contains three main open reading frames (ORFs).

Joyce A, Pereira de Souza A, Gabriela Murúa M, J. Clough S, Imaculada Zucchi M. A first draft genome of the Sugarcane borer, Diatraea saccharalis

Background: The sugarcane borer (Diatraea saccharalis), a widely distributed moth throughout the Americas, is a pest that affects economically important crops such as sugarcane, sorghum, wheat, maize and rice. Given its significant impact on yield reduction, whole-genome information of the species is needed. Here, we report the first draft assembly of the D. saccharalis genome. Methods: The genomic sequences were obtained using the Illumina HiSeq 2500 whole-genome sequencing of a single adult male specimen. We assembled the short-reads using the SPAdes software and predicted protein-coding genes using MAKER. Genome assembly completeness was assessed through BUSCO and the repetitive content by RepeatMasker. Results: The 453 Mb assembled sequences contain 1,445 BUSCO gene orthologs and 1,161 predicted gene models identified based on homology evidence to the domestic silk moth, Bombyx mori. The repeat content composes 41.18% of the genomic sequences which is in the range of other lepidopteran species. Conclusions: Functional annotation reveals that predicted gene models are involved in important cellular mechanisms such as metabolic pathways and protein synthesis. Thus, the data generated in this study expands our knowledge on the genomic characteristics of this devastating pest and provides essential resources for future genetic studies of the species.

Caterpillar Venom: A Health Hazard of the 21st Century

Caterpillar envenomation is a global health threat in the 21st century. Every direct or indirect contact with the urticating hairs of a caterpillar results in clinical manifestations ranging from local dermatitis symptoms to potentially life-threatening systemic effects. This is mainly due to the action of bioactive components in the venom that interfere with targets in the human body. The problem is that doctors are limited to relieve symptoms, since an effective treatment is still lacking. Only for Lonomia species an effective antivenom does exist. The health and economical damage are an underestimated problem and will be even more of a concern in the future. For some caterpillar species, the venom composition has been the subject of investigation, while for many others it remains unknown. Moreover, the targets involved in the pathophysiology are poorly understood. This review aims to give an overview of the knowledge we have today on the venom composition of different caterpillar species along with their pharmacological targets. Epidemiology, mode of action, clinical time course and treatments are also addressed. Finally, we briefly discuss the future perspectives that may open the doors for future research in the world of caterpillar toxins to find an adequate treatment.

From refugia to contact: Pine processionary moth hybrid zone in a complex biogeographic setting

Contact zones occur at the crossroad between specific dispersal routes and are facilitated by biogeographic discontinuities. Here, we focused on two Lepidoptera sister species that come in contact near the Turkish Straits System (TSS). We aimed to infer their phylogeographic histories in the Eastern Mediterranean and finely analyze their co-occurrence and hybridization patterns in this biogeographic context. We used molecular mitochondrial and nuclear markers to study 224 individuals from 42 localities. We used discordances between markers and complementary assignment methods to identify and map hybrids and parental individuals. We confirmed the parapatric distribution of Thaumetopoea pityocampa (Lepidoptera: Notodontidae) in the west and Thaumetopoea wilkinsoni in the east and identified a narrow contact zone. We identified several glacial refugia of T. wilkinsoni in southern Turkey with a strong east-west differentiation in this species. Unexpectedly, T. pityocampa crossed the TSS and occur in northern Aegean Turkey and some eastern Greek islands. We found robust evidence of introgression between the two species in a restricted zone in northwestern Turkey, but we did not identify any F1 individuals. The identified hybrid zone was mostly bimodal. The distributions and genetic patterns of the studied species were strongly influenced both by the Quaternary climatic oscillations and the complex geological history of the Aegean region. T. pityocampa and T. wilkinsoni survived the last glacial maximum in disjoint refugia and met in western Turkey at the edge of the recolonization routes. Expanding population of T. wilkinsoni constrained T. pityocampa to the western Turkish shore. Additionally, we found evidence of recurrent introgression by T. wilkinsoni males in several T. pityocampa populations. Our results suggest that some prezygotic isolation mechanisms, such as differences in timing of the adult emergences, might be a driver of the isolation between the sister species.

Candidate detoxification-related genes in brown planthopper, Nilaparvata lugens, in response to β-asarone based on transcriptomic analysis

Nilaparvata lugens(Stål) is a serious pest of rice and has evolved different levels of resistance against most chemical pesticides. β-asarone is the main bioactive insecticidal compound of Acorus calamus L. that shows strong insecticidal activity against pests. In this study, we conducted a bioassay experiment to determine the contact toxicity of β-asarone to N. lugens nymphs. The LD₃₀ sublethal dose was 0.106 μg per nymph, with 95% confidence limits of 0.070-0.140 μg. We applied the LD₃₀ concentration of β-asarone to nymphs for 24 h or 72 h and then performed a transcriptome sequence analysis by referencing the N. lugens genome to characterize the variation. The transcriptomic analysis showed that several GO terms and KEGG pathways presented significant changes. Individually, 126 differentially expressed genes (DEGs), including 72 upregulated and 54 downregulated genes, were identified at 24 h, and 1771 DEGs, including 882 upregulated and 889 downregulated genes, were identified at 72 h. From the DEGs, we identified a total of 40 detoxification-related genes, including eighteen Cytochrome P450 monooxygenase genes (P450s), three Glutathione S-transferase genes, one Carboxylesterase gene, twelve UDP-glucosyltransferases and six ATP-binding cassette genes. We selected the eighteen P450s for subsequent verification by quantitative PCR. These findings indicated that β-asarone presented strong contact toxicity to N. lugens nymphs and induced obvious variation of detoxification-related genes that may be involved in the response to β-asarone.

New Data on the Range Expansion of the Thaumetopoea pityocampa (Lepidoptera: Notodontidae) 'ENA clade' in Greece: The Role of Bacterial Endosymbionts

The pine processionary moth, Thaumetopoea pityocampa (Denis and Schiffermüller), is an important insect in the Mediterranean region, as it defoliates pines and its urticating hairs can cause allergic reactions in humans and animals. Moreover, this species exhibits an interesting genetic structure as recently a distinct East-North African mtDNA lineage ('ENA clade') has been described. This clade has been recently detected in Greek populations where it has currently expanded its range by replacing the 'endemic' T. pityocampa lineages. Here, we report new data on the rapid spread of 'ENA clade' in the Greek island Evoia in only a few years. As the underlying mechanisms of the 'ENA clade' range expansion has not been studied so far, we screened T. pityocampa for an infection with the heritable bacterial endosymbionts Wolbachia (Bacteria: Anaplasmataceae), Cardinium (Bacteria: Bacteroidaceae), Rickettsia (Bacteria: Rickettsiaceae) and Spiroplasma (Bacteria: Spiroplasmataceae). These bacteria can manipulate the reproduction of infected hosts, something that could potentially explain the rapid spread of 'ENA clade' lineage. Therefore, we screened 28 individuals that exhibited T. pityocampa 'ENA clade' and 'endemic' T. pityocampa haplotypes from nine populations scattered all over Greece. None of them was infected with any of the four endosymbionts, suggesting that these bacteria do not cause reproductive manipulations in T. pityocampa lineages and, thus, other factors should be explored in future research efforts.

Expansion of LINEs and species-specific DNA repeats drives genome expansion in Asian Gypsy Moths

Two subspecies of Asian gypsy moth (AGM), Lymantria dispar asiatica and L. dispar japonica, pose a serious alien invasive threat to North American forests. Despite decades of research on the ecology and biology of this pest, limited AGM-specific genomic resources are currently available. Here, we report on the genome sequences and functional content of these AGM subspecies. The genomes of L.d. asiatica and L.d. japonica are the largest lepidopteran genomes sequenced to date, totaling 921 and 999 megabases, respectively. Large genome size in these subspecies is driven by the accumulation of specific classes of repeats. Genome-wide metabolic pathway reconstructions suggest strong genomic signatures of energy-related pathways in both subspecies, dominated by metabolic functions related to thermogenesis. The genome sequences reported here will provide tools for probing the molecular mechanisms underlying phenotypic traits that are thought to enhance AGM invasiveness.

Identification and analysis of the complete mitochondrial genome of Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

The mitochondrial genome (mitogenome) provides important information for phylogenetic analysis and understanding evolutionary origins. Thaumetopoea pityocampa is a forest pest that harms nearly all cedar and pine species. In this study, the T. pityocampa mitochondrial genome was sequenced, assembled, and annotated. The sequence length of the genome was found to be 15,737 bp, containing 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and an A + T-rich region compared with the genomes of other lepidopterans. The overall nucleotide composition is: 37.3% T, 40.5% A, 14.6% C, and 7.6% G, demonstrating an AT bias (A + T: 77.8%). Our phylogenetic tree analysis results showed that T. pityocampa and Ochrogaster lunifer were the most similar species, with the closest evolutionary distance. The mitogenome sequence determined in this study will contribute to improved understanding of Notodontidae evolution.

The Genomic Basis of Color Pattern Polymorphism in the Harlequin Ladybird

Many animal species comprise discrete phenotypic forms. A common example in natural populations of insects is the occurrence of different color patterns, which has motivated a rich body of ecological and genetic research [1-6]. The occurrence of dark, i.e., melanic, forms displaying discrete color patterns is found across multiple taxa, but the underlying genomic basis remains poorly characterized. In numerous ladybird species (Coccinellidae), the spatial arrangement of black and red patches on adult elytra varies wildly within species, forming strikingly different complex color patterns [7, 8]. In the harlequin ladybird, Harmonia axyridis, more than 200 distinct color forms have been described, which classic genetic studies suggest result from allelic variation at a single, unknown, locus [9, 10]. Here, we combined whole-genome sequencing, population-based genome-wide association studies, gene expression, and functional analyses to establish that the transcription factor Pannier controls melanic pattern polymorphism in H. axyridis. We show that pannier is necessary for the formation of melanic elements on the elytra. Allelic variation in pannier leads to protein expression in distinct domains on the elytra and thus determines the distinct color patterns in H. axyridis. Recombination between pannier alleles may be reduced by a highly divergent sequence of ∼170 kb in the cis-regulatory regions of pannier, with a 50 kb inversion between color forms. This most likely helps maintain the distinct alleles found in natural populations. Thus, we propose that highly variable discrete color forms can arise in natural populations through cis-regulatory allelic variation of a single gene.