Comparative analysis of two phenologically divergent populations of the pine processionary moth (Thaumetopoea pityocampa) by de novo transcriptome sequencing

Transcriptome profiling for developmental stages Protaetia brevitarsis seulensis with focus on wing development and metamorphosis

A white-spotted flower chafer Protaetia brevitarsis seulensis widely distributed in Asian countries is traditionally used in oriental medicine. This study explored gene expression abundance with respect to wing development and metamorphosis in P. b. seulensis based on the large-scale RNA-seq data. The transcriptome assembly consists of 23,551 high-quality transcripts which are approximately 96.7% covered. We found 265 wing development genes, 19 metamorphosis genes, and 1,314 candidates. Of the 1,598 genes, 1,594 are included exclusively in cluster 4 with similar gene co-expression patterns. The network centrality analyses showed that wing development- and metamorphosis-related genes have a high degree of betweenness centrality and are expressed most highly in eggs, moderately in pupa and adults, and lowest in larva. This study provides some meaningful clues for elucidating the genetic modulation mechanism of wing development and metamorphosis in P. b. seulensis.

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.

Functional characterization of the Hyles euphorbiae hawkmoth transcriptome reveals strong expression of phorbol ester detoxification and seasonal cold hardiness genes

Background

The European spurge hawkmoth, Hyles euphorbiae (Lepidoptera, Sphingidae), has been intensively studied as a model organism for insect chemical ecology, cold hardiness and evolution of species delineation. To understand species isolation mechanisms at a molecular level, this study aims at determining genetic factors underlying two adaptive ecological trait candidates, phorbol ester (TPA) detoxification and seasonal cold acclimation.

Method

A draft transcriptome of H. euphorbiae was generated using Illumina sequencing, providing the first genomic resource for the hawkmoth subfamily Macroglossinae. RNA expression levels in tissues of experimental TPA feeding larvae and cooled pupae was compared to levels in control larvae and pupae using 26 bp RNA sequence tag libraries (DeepSuperSAGE). Differential gene expression was assessed by homology searches of the tags in the transcriptome.

Results

In total, 389 and 605 differentially expressed transcripts for detoxification and cold hardiness, respectively, could be identified and annotated with proteins. The majority (22 of 28) of differentially expressed detox transcripts of the four 'drug metabolism' enzyme groups (cytochrome P450 (CYP), carboxylesterases (CES), glutathione S-transferases (GST) and lipases) are up-regulated. Triacylglycerol lipase was significantly over proportionally annotated among up-regulated detox transcripts. We record several up-regulated lipases, GSTe2, two CESs, CYP9A21, CYP6BD6 and CYP9A17 as candidate genes for further H. euphorbiae TPA detoxification analyses. Differential gene expression of the cold acclimation treatment is marked by metabolic depression with enriched Gene Ontology terms among down-regulated transcripts almost exclusively comprising metabolism, aerobic respiration and dissimilative functions. Down-regulated transcripts include energy expensive respiratory proteins like NADH dehydrogenase, cytochrome oxidase and ATP synthase. Gene expression patterns show shifts in carbohydrate metabolism towards cryoprotectant production. The Glycolysis enzymes, G1Pase, A1e, Gpi and an Akr isoform are up-regulated. Glycerol, an osmolyte which lowers the body liquid supercooling point, appears to be the predominant polyol cryoprotectant in H. euphorbiae diapause pupae. Several protein candidates involved in glucose, glycerol, myo-inositol and potentially sorbitol and trehalose synthesis were identified.

Conclusions

A majority of differently expressed transcripts unique for either detoxification or cold hardiness indicates highly specialized functional adaptation which may have evolved from general cell metabolism and stress response.The transcriptome and extracted candidate biomarkers provide a basis for further gene expression studies of physiological processes and adaptive traits in H. euphorbiae.

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

The pine processionary moth Thaumetopoea pityocampa (Lepidoptera: Notodontidae) is the main pine defoliator in the Mediterranean region. Its urticating larvae cause severe human and animal health concerns in the invaded areas. This species shows a high phenotypic variability for various traits, such as phenology, fecundity and tolerance to extreme temperatures. This study presents the construction and analysis of extensive genomic and transcriptomic resources, which are an obligate prerequisite to understand their underlying genetic architecture. Using a well-studied population from Portugal with peculiar phenological characteristics, the karyotype was first determined and a first draft genome of 537 Mb total length was assembled into 68,292 scaffolds (N50 = 164 kb). From this genome assembly, 29,415 coding genes were predicted. To circumvent some limitations for fine-scale physical mapping of genomic regions of interest, a 3X coverage BAC library was also developed. In particular, 11 BACs from this library were individually sequenced to assess the assembly quality. Additionally, de novo transcriptomic resources were generated from various developmental stages sequenced with HiSeq and MiSeq Illumina technologies. The reads were de novo assembled into 62,376 and 63,175 transcripts, respectively. Then, a robust subset of the genome-predicted coding genes, the de novo transcriptome assemblies and previously published 454/Sanger data were clustered to obtain a high-quality and comprehensive reference transcriptome consisting of 29,701 bona fide unigenes. These sequences covered 99% of the cegma and 88% of the busco highly conserved eukaryotic genes and 84% of the busco arthropod gene set. Moreover, 90% of these transcripts could be localized on the draft genome. The described information is available via a genome annotation portal (http://bipaa.genouest.org/sp/thaumetopoea_pityocampa/).

Proteome Analysis of Urticating Setae From Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

Thaumetopoea pityocampa (Denis & Schiffermüller) (Lepidoptera: Notodontidae) is harmful to conifer trees because of defoliation and to public health because of the release of urticating setae from the caterpillars. Contact with setae by humans and domestic animals induces dermatitis, usually localized to the exposed areas. Recent studies demonstrated the presence of a complex urticating mechanism where proteins present in the setae may play a role as activators of immune responses. Yet, limited information is available at present about the proteins occurring in the setae of T. pityocampa. Using a refined method for protein extraction from the setae, and a combination of liquid chromatography tandem-mass spectrometry (LC-MS/MS), de novo assembly of transcriptomic data, and sequence similarity searches, an extensive data set of 353 proteins was obtained. These were further categorized by molecular function, biological process, and cellular location. All the 353 proteins identified were found to match through BLAST search with at least one Lepidoptera sequence available in databases. We found the previously known allergens Tha p 1 and Tha p 2 described from T. pityocampa, as well as enzymes involved in chitin biosynthesis, one of the principal components of the setae, and serine proteases that were responsible for inflammatory and allergic reactions in other urticating Lepidoptera. This new proteomic database may allow for a better understanding of the complexity of allergenic reactions due to T. pityocampa and to other Lepidoptera sharing similar defense systems.

Effect of heat waves on embryo mortality in the pine processionary moth

Extreme climate events such as heat waves are predicted to become more frequent with climate change, representing a challenge for many organisms. The pine processionary moth Thaumetopoea pityocampa is a Mediterranean pine defoliator, which typically lays eggs during the summer. We evaluated the effects of heat waves on egg mortality of three populations with different phenologies: a Portuguese population with a classical life cycle (eggs laid in summer), an allochronic Portuguese population reproducing in spring, and a Tunisian population from the extreme southern limit of T. pityocampa distribution range, in which eggs are laid in fall. We tested the influence of three consecutive hot days on egg survival and development time, using either constant (CT) or daily cycling temperatures (DT) with equivalent mean temperatures. Maximum temperatures (T max) used in the experiment ranged from 36 to 48°C for DT and from 30 to 42°C for CT. Heat waves had a severe negative effect on egg survival when T max reached 42°C for all populations. No embryo survived above this threshold. At high mean temperatures (40°C), significant differences were observed between populations and between DT and CT regimes. Heat waves further increased embryo development time. The knowledge we gained about the upper lethal temperature to embryos of this species will permit better prediction of the potential expansion of this insect under different climate warming scenarios.

In search of pathogens: transcriptome-based identification of viral sequences from the pine processionary moth (Thaumetopoea pityocampa)

Thaumetopoea pityocampa (pine processionary moth) is one of the most important pine pests in the forests of Mediterranean countries, Central Europe, the Middle East and North Africa. Apart from causing significant damage to pinewoods, T. pityocampa occurrence is also an issue for public and animal health, as it is responsible for dermatological reactions in humans and animals by contact with its irritating hairs. High throughput sequencing technologies have allowed the fast and cost-effective generation of genetic information of interest to understand different biological aspects of non-model organisms as well as the identification of potential pathogens. Using these technologies, we have obtained and characterized the transcriptome of T. pityocampa larvae collected in 12 different geographical locations in Turkey. cDNA libraries for Illumina sequencing were prepared from four larval tissues, head, gut, fat body and integument. By pooling the sequences from Illumina platform with those previously published using the Roche 454-FLX and Sanger methods we generated the largest reference transcriptome of T. pityocampa. In addition, this study has also allowed identification of possible viral pathogens with potential application in future biocontrol strategies.