Environmental stresses induce the expression of putative glycine-rich insect cuticular protein genes in adult Leptinotarsa decemlineata (Say)

Genome-wide annotation and comparative analysis revealed conserved cuticular protein evolution among non-biting midges with varied environmental adaptability

Chironomidae, non-biting midges, a diverse and abundant insect group in global aquatic ecosystems, represent an exceptional model for investigating genetic adaptability mechanisms in aquatic insects due to their extensive species diversity and resilience to various environmental conditions. The cuticle in insects acts as the primary defense against ecological pressures. Cuticular Proteins (CPs) determine cuticle characteristics, facilitating adaptation to diverse challenges. However, systematic annotation of CP genes has only been conducted for one Chironomidae species, Propsilocerus akamusi, by our team. In this study, we expanded this annotation by identifying CP genes in eight additional Chironomidae species, covering all Chironomidae species with available genome data. We identified a total of 889 CP genes, neatly categorized into nine CP families: 215 CPR RR1 genes, 272 CPR RR2 genes, 23 CPR RR3 genes, 21 CPF genes, 16 CPLCA genes, 19 CPLCG genes, 28 CPLCP genes, 77 CPAP genes, and 37 Tweedle genes. Subsequently, we conducted a comprehensive phylogenetic analysis of CPs within the Chironomidae family. This expanded annotation of CP genes across diverse Chironomidae species significantly contributes to our understanding of their remarkable adaptability.

Identification of insect cuticular protein genes LCP17 and SgAbd5 from Helicoverpa armigera and evaluation their roles in fenvalerate resistance

Insect cuticular protein (ICP) plays an important role in insect growth and development. However, research on the role of ICP in insecticide resistance is very limited. In this study, insect cuticular protein genes LCP17 and SgAbd5 were cloned and characterized in Helicoverpa armigera based on previous transcriptome data. The functions of LCP17 and SgAbd5 genes in fenvalerate resistance were assessed by RNA interference (RNAi), and their response to fenvalerate was further detected. The results showed that LCP17 and SgAbd5 were overexpressed in the fenvalerate-resistant strain comparing with a susceptible strain. The open reading frames of LCP17 and SgAbd5 genes were 423 bp and 369 bp, encoding 141 and 123 amino acids, respectively. LCP17 and SgAbd5 genes were highly expressed in the larval stage, but less expressed in the adult and pupal stages. The expression level of LCP17 and SgAbd5 genes increased significantly after fenvalerate treatment at 24 h. When the cotton bollworms larvae were exposed to fenvalerate at LD50 level, RNAi-mediated silencing of LCP17 and SgAbd5 genes increased the mortality from 50.68% to 68.67% and 63.89%, respectively; the mortality increased to even higher level, which was 73.61%, when these two genes were co-silenced. Moreover, silencing of these two genes caused the cuticle lamellar structure to become loose, which led to increased penetration of fenvalerate into the larvae. The results suggested that LCP17 and SgAbd5 may be involved in the resistance of cotton bollworm to fenvalerate, and LCP17 and SgAbd5 could serve as potential targets for H. armigera control.

Assessing genome-wide adaptations associated with range expansion in the pink rice borer, Sesamia inferens

Understanding the genetic basis of adaptive evolution following habitat expansion can have important implications for pest management. The pink rice borer (PRB), Sesamia inferens (Walker), is a destructive pest of rice that was historically restricted to regions south of 34° N latitude in China. However, with changes in global climate and farming practices, the distribution of this moth has progressively expanded, encompassing most regions in North China. Here, 3 highly differentiated subpopulations were discovered using high-quality single-nucleotide polymorphism and structural variant datasets across China, corresponding to northern, southern China regions, and the Yunnan-Guizhou Plateau, with significant patterns of isolation by geographic and environmental distances. Our estimates of evolutionary history indicate asymmetric migration with varying population sizes across the 3 subpopulations. Selective sweep analyses estimated strong selection at insect cuticle glycine-rich cuticular protein genes which are associated with enhanced desiccation adaptability in the northern group, and at the histone-lysine-N-methyltransferase gene associated with range expansion and local adaptation in the Shandong population. Our findings have significant implications for the development of effective strategies to control this pest.

GARP: A family of glycine and alanine-rich proteins that helps spider mites feed on plants

Spider mites (Tetranychidae) are destructive agricultural pests which have evolved strategies to overcome plant defenses, such as the ability to puncture the leaves of their hosts to feed. The expression of many genes with unknown functions is altered during feeding, but little is known about the role of these genes in plant-mite interactions. Here, we identified 3 novel gene families through analysis of genomic and transcriptomic data from 3 spider mite species. These GARP family genes encode glycine and alanine-rich proteins; they are present in mites (Acariformes) but absent in ticks (Parasitiformes) in the subclass Acari, indicating that these genes have undergone a significant expansion in spider mites and thus play important adaptive roles. Transcriptomic analysis revealed that the expression of GARP genes is strongly correlated with feeding and the transfer to new hosts. We used RNA interference to silence GARP1d in the two-spotted spider mite Tetranychus urticae, which inhibited feeding and egg laying and significantly increased mortality when the mites were transferred to soybean shoots; a similar effect was observed after TuVATPase was silenced. However, no changes in mite mortality were observed after TuGARP1d-silenced mites were placed on an artificial diet, which was different from the effect of TuVATPase silencing. Our results indicate that GARP family members play important roles in mite-plant interactions. Additional studies are needed to clarify the mechanisms underlying these interactions.

Identifying potential insecticide resistance markers through genomic-level comparison of Bemisia tabaci (Gennadius) lines

The invasive whitefly (Bemisia tabaci) MED is one of the most economically damaging plant pests. The extensive use of insecticide over decades has led to that the invasive B. tabaci MED has developed resistance to a wide range of insecticide classes, but little is known about the genetic background associated with resistance. To this end, we conducted a comparative genome-wide analysis of single-base nucleotide polymorphisms between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line collected in 1976. First, low-coverage genome sequencings were conducted on DNA isolated from individual whiteflies. The sequencing results were evaluated using an available B. tabaci MED genome as a reference. Significant genetic differences were discovered between MED whitefly lines collected from fields that were recently infested and an insecticide-susceptible MED whitefly line based on the principal component analyses. Top GO categories and KEGG pathways that might be involved in insecticide resistance development were identified, and several of them have not been previously associated with resistance. Additionally, we identified several genetic loci with novel variations including Cytochrome P450 monooxygenases (P450s), UDP-glucuronosyltransferases (UGTs), Glutathione S-transferases (GSTs), esterase, carboxyl-esterases (COE), ABC transporters, fatty acyl-CoA reductase, voltage-gated sodium channels, GABA receptor, and cuticle proteins (CPs) that were previously reported to have close associations with pesticide resistance in well-studied insect groups that provide an essential resource for the design of insecticide resistance-linked loci arrays insecticide. Our results was obtained solely on resequencing genome data sets, more pesticide bio-assays combined with omics datasets should be further used to verify the markers identified here.

Honey Bee Larval Hemolymph as a Source of Key Nutrients and Proteins Offers a Promising Medium for Varroa destructor Artificial Rearing

Varroa destructor, a major ectoparasite of the Western honey bee Apis mellifera, is a widespread pest that damages colonies in the Northern Hemisphere. Throughout their lifecycle, V. destructor females feed on almost every developmental stage of their host, from the last larval instar to the adult. The parasite is thought to feed on hemolymph and fat body, although its exact diet and nutritional requirements are poorly known. Using artificial Parafilm™ dummies, we explored the nutrition of V. destructor females and assessed their survival when fed on hemolymph from bee larvae, pupae, or adults. We compared the results with mites fed on synthetic solutions or filtered larval hemolymph. The results showed that the parasites could survive for several days or weeks on different diets. Bee larval hemolymph yielded the highest survival rates, and filtered larval plasma was sufficient to maintain the mites for 14 days or more. This cell-free solution therefore theoretically contains all the necessary nutrients for mite survival. Because some bee proteins are known to be hijacked without being digested by the parasite, we decided to run a proteomic analysis of larval honey bee plasma to highlight the most common proteins in our samples. A list of 54 proteins was compiled, including several energy metabolism proteins such as Vitellogenin, Hexamerin, or Transferrins. These molecules represent key nutrient candidates that could be crucial for V. destructor survival.

Transcriptomic and Proteomic Analyses of Myzus persicae Carrying Brassica Yellows Virus

Viruses in the genus Polerovirus infect a wide range of crop plants and cause severe economic crop losses. BrYV belongs to the genus Polerovirus and is transmitted by Myzus persicae. However, the changes in transcriptome and proteome profiles of M. persicae during viral infection are unclear. Here, RNA-Seq and TMT-based quantitative proteomic analysis were performed to compare the differences between viruliferous and nonviruliferous aphids. In total, 1266 DEGs were identified at the level of transcription with 980 DEGs being upregulated and 286 downregulated in viruliferous aphids. At the protein level, among the 18 DEPs identified, the number of upregulated proteins in viruliferous aphids was twice that of the downregulated DEPs. Enrichment analysis indicated that these DEGs and DEPs were mainly involved in epidermal protein synthesis, phosphorylation, and various metabolic processes. Interestingly, the expressions of a number of cuticle proteins and tubulins were upregulated in viruliferous aphids. Taken together, our study revealed the complex regulatory network between BrYV and its vector M. persicae from the perspective of omics. These findings should be of great benefit to screening key factors involved in the process of virus circulation in aphids and provide new insights for BrYV prevention via vector control in the field.

Leveraging RNA Interference to Impact Insecticide Resistance in the Colorado Potato Beetle, Leptinotarsa decemlineata

The Colorado potato beetle, Leptinotarsa decemlineata Say, is a potato pest that can cause important economic losses to the potato industry worldwide. Diverse strategies have been deployed to target this insect such as biological control, crop rotation, and a variety of insecticides. Regarding the latter, this pest has demonstrated impressive abilities to develop resistance against the compounds used to regulate its spread. Substantial work has been conducted to better characterize the molecular signatures underlying this resistance, with the overarching objective of leveraging this information for the development of novel approaches, including RNAi-based techniques, to limit the damage associated with this insect. This review first describes the various strategies utilized to control L. decemlineata and highlights different examples of reported cases of resistances against insecticides for this insect. The molecular leads identified as potential players modulating insecticide resistance as well as the growing interest towards the use of RNAi aimed at these leads as part of novel means to control the impact of L. decemlineata are described subsequently. Finally, select advantages and limitations of RNAi are addressed to better assess the potential of this technology in the broader context of insecticide resistance for pest management.

Comparative transcriptome analysis of Callosobruchus chinensis (L.) (Coleoptera: Chrysomelidae-Bruchinae) after heat and cold stress exposure

Callosobruchus chinensis is regarded as one of the cosmopolitan pests of legume crops and can cause tremendous losses to a variety of beans. This study focused on comparative transcriptome analyses of C. chinensis exposed to 45 °C (heat stress), 27 °C (ambient temperature) and -3 °C (cold stress) for 3 h to investigate the gene differences and underlying molecular mechanisms. There were 402 and 111 differentially expressed genes (DEGs) identified in the heat and cold stress treatments, respectively. "cell process", "cell" and "binding" were the main enriched functions and biological processes revealed by gene ontology (GO) analysis. The clusters of orthologous genes (COG) showed that DEGs were assigned to the categories: "posttranslational modification, protein turnover, chaperones", "lipid transport and metabolism", and "general function prediction only". With respect to the Kyoto Encyclopedia of Genes and Genomes (KEGG), the "longevity regulating pathway-multiple species", "carbon metabolism", "peroxisome", "protein processing in endoplasmic", "glyoxylate and dicarboxylate metabolism" pathways were significantly enriched. The annotation and enrichment analysis revealed that genes encoding heat shock proteins (Hsps) and cuticular proteins were significantly upregulated under high and low-temperature stresses, respectively. In addition, some DEGs encoding "Protein lethal essential for life", "Reverse transcriptase", "DnaJ domain", "Cytochrome" and "Zinc finger protein" were also upregulated to varying degrees. Transcriptomic data were validated using qRT‒PCR, which confirmed that they were consistent. In this paper, the temperature tolerance of C. chinensis adults was evaluated and the results showed that female adults were more sensitive to heat and cold stress than males, and the upregulation of heat shock protein and epidermal protein was the largest in DEGs after heat and cold stress, respectively. These findings provide a reference for further understanding the biological characteristics of C. chinensis adults and the molecular mechanisms underlying the response to low and high temperatures.

Genome-wide annotation of cuticular protein genes in non-biting midge Propsilocerus akamusi and transcriptome analysis of their response to heavy metal pollution

The insect cuticle is a sophisticated chitin-protein extracellular structure for mutable functions. The cuticles varied their structures and properties in different species, and the same species but in different regions or at different stages, to fill the requirements of different functions. The alteration of cuticle structures may also be induced due to challenges by some environmental crises, such as pollution exposures. The physical properties of the cuticle were determined by the cuticle proteins (CPs) they contain. The cuticle proteins are large protein groups in all insects, which are commonly divided into different families according to their conserved protein sequence motifs. Although Chironomidae is an abundant and universal insect in global aquatic ecosystems and a popular model for aquatic toxicology, no systematic annotation of CPs was done for any species in Chironomidae before. In this work, we annotated the CP genes of Propsilocerus akamusi, the most abundant Chironomidae species in Asia. A total of 160 CP genes were identified, and 97 of them could be well classified into eight CP families: 76 CPR genes can be subdivided into three groups (further divided into three subgroups: 36 RR1 genes, 37 RR2 genes, and 3 RR3 genes), 2 CPF genes, 3 CPLCA genes, 1 CPLCG gene, 8 CPAP genes, and 3 Tweedle genes. Additionally, we analyzed the response of P. akamusi CP genes at expression level to Cu exposure, which is related to the high heavy metal tolerance and the earlier onset of pupariation in heavy metal polluted water.

Comparative transcriptome analysis reveals differences in gene expression in whitefly following individual or combined applications of Akanthomyces attenuatus (Zare & Gams) and matrine

BACKGROUND

Bemisia tabaci Gennadius (Hemiptera: Aleyrodidae) is a serious pest of crops in different regions of the world. Our recent studies on the joint application of Akanthomyces attenuatus (a pathogenic insect fungus) and matrine (a botanical insecticide) against B. tabaci have shown promising results. Using RNA sequencing (RNA-Seq), we identified differentially expressed genes involved in whitefly responses to single or mixed applications of A. attenuatus and matrine.

METHODS

In this study, we compared the transcriptome profiles of B. tabaci treated with individual and combined treatments of A. attenuatus and matrine to determine variations in gene expression among whiteflies in response to different treatments.

RESULTS

Transcriptomic data analysis showed differential expression of 71, 1194, and 51 genes in response to A. attenuatus (BtA), matrine (BtM), and A. attenuatus + matrine (BtAM) treatment, respectively. A total of 65 common differentially expressed genes (DEGs) were identified between whiteflies treated with A. attenuatus (BtA) and matrine (BtM). A comparison of DEGs across the three treatments (BtA, BtM, and BtAM) revealed two common DEGs. The results also revealed that AMPK signaling, apoptosis, and drug metabolism pathways are likely involved in whitefly defense responses against A. attenuatus and matrine infection. Furthermore, a notable suppression of general metabolism and immune response genes was observed in whiteflies treated with A. attenuatus + matrine (BtAM) compared to whiteflies treated with individual A. attenuatus (BtA) or matrine (BtM) treatments.

CONCLUSION

Dynamic changes in the number of differentially expressed genes were observed in B. tabaci subjected to different treatments (BtA, BtM, and BtAM). To the best of our knowledge, this is the first report on the molecular interactions between whitefly and individual or combined treatments of A. attenuatus and matrine. These results will further improve our knowledge of the infection mechanism and complex biochemical processes involved in the synergistic action of A. attenuatus and matrine against B. tabaci.

Characterization of cuticular proteins in CPR family in the wolf spider, Pardosa pseudoannulata, and the response of one subfamily genes to environmental stresses

Cuticular protein (CP) plays an essential role in the construction and function of exoskeleton in arthropods. CPR family, CP with Rebers and Riddiford (R&R) Consensus, is the largest CP family in insects, but it lacks systematic research in non-insect arthropods. In this study, we explored CPRs in the wolf spider, Pardosa pseudoannulata, a predator to many insect pests. We totally identified 152 CPRs in P. pseudoannulata genome, which were divided into two subgroups based on R&R Consensus sequences, with 12 CPRs in RR-1 and 140 in RR-2. All RR-2 members presented a novel Consensus with 34 amino acids, G-x(8)-G-x(6)-Y-x-A-x(3)-G-x(7)-N-E-x-G, which was a common characteristic for RR-2 CPRs in chelicerates. Transcriptome data was used to document the expression patterns of CPR genes in different tissues and ecdysis processes. The specific expressions were found for part CPR genes, such as five RR-2 genes that were specifically expressed in male genital bulbs and eleven RR-1 genes that were highly expressed in the integument. Due to the limited number and integument-specific expression of RR-1 genes, we further analyzed their responses to different environmental stresses at the transcriptional level. Except for PapsCPR11, ten RR-1 genes responded to at least one environmental stress, among with the expression of PapsCPR12 was significantly changed by three stresses (dryness, low temperature and imidacloprid treatments). Silencing PapsCPR12 increased the tolerance of P. pseudoannulata to imidacloprid. Overall, the results presented novel Consensus characteristics of CPRs in P. pseudoannulata, which was helpful for the identification and evolution analysis of CPRs in non-insect arthropods.

Cuticular protein genes involve heat acclimation of insect larvae under global warming

Cuticular proteins (CPs) play important roles in insect growth and development. However, it is unknown whether CPs are related to heat tolerance. Cnaphalocrocis medinalis, a serious pest of rice, occurs in summer and exhibits strong adaptability to high temperature, but the underlying mechanism is unclear. Here, the role of CP genes in heat acclimation was studied. Heat tolerance of the heat-acclimated larvae was significantly stronger than the unacclimated larvae. The cuticular protein content in the heat-acclimated larvae was higher than that of the unacclimated larvae. 191 presumed CP genes of C. medinalis (CmCPs) were identified. Expression patterns of 14 CmCPs were different between the heat acclimated (S39) and unacclimated (S27) larvae under heat stress. CmCPs were specifically expressed in epidermis and the head except CmCPR20 mainly expressed in Malpighian tubules. CmCPR20 was upregulated in S39 while downregulated in S27, but CmTweedle1 and CmCPG1 were upregulated in S27 and downregulated in S39. RNAi CmTweedle1 or CmCPG1 remarkably decreased heat tolerance and cuticular protein content of the heat-acclimated larvae but not the unacclimated larvae. RNAi CmCPR20 decreased heat tolerance and cuticular protein content of the unacclimated larvae but not the heat-acclimated larvae. CmTweedle1 and CmCPG1 genes involve heat acclimation of C. medinalis.

Convergent adaptation of ootheca formation as a reproductive strategy in Polyneoptera

Insects have evolved numerous adaptations and colonized diverse terrestrial environments. Several polyneopterans, including dictyopterans (cockroaches and mantids) and locusts, have developed oothecae, but little is known about the molecular mechanism, physiological function, and evolutionary significance of ootheca formation. Here, we demonstrate that the cockroach asymmetric colleterial glands produce vitellogenins, proline-rich protein, and glycine-rich protein as major ootheca structural proteins (OSPs) that undergo sclerotization and melanization for ootheca formation through the cooperative protocatechuic acid pathway and dopachrome and dopaminechrome subpathway. Functionally, OSP sclerotization and melanization prevent eggs from losing water at warm and dry conditions, and thus effectively maintain embryo viability. Dictyopterans and locusts convergently evolved vitellogenins, apolipoprotein D, and laminins as OSPs, whereas within Dictyoptera, cockroaches and mantids independently developed glycine-rich protein and fibroins as OSPs. Highlighting the ecological-evolutionary importance, convergent ootheca formation represents a successful reproductive strategy in Polyneoptera that promoted the radiation and establishment of cockroaches, mantids, and locusts.

Genomic and Transcriptomic Analysis Reveals Cuticular Protein Genes Responding to Different Insecticides in Fall Armyworm Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is a serious pest of crucial crops causing great threats to the food security of the world. It has evolved resistance to various insecticides, while the underlying molecular mechanisms remain largely unknown. Cuticular proteins (CPs), as primary components in cuticle, play an important role in insects' protection against environmental stresses. Few of them have been documented as participating in insecticide resistance in several insect species. In order to explore whether CP genes of the FAW exhibit a functional role in responding to insecticides stress, a total of 206 CPs, classified into eight families, were identified from the genome of the FAW through a homology-based approach coupled with manual efforts. The temporal expression profiles of all identified CP genes across developmental stages and their responses to 23 different insecticides were analyzed using the RNA-seq data. Expression profiling indicated that most of the CP genes displayed stage-specific expression patterns. It was found that the expression of 51 CP genes significantly changed after 48 h exposure to 17 different insecticides. The expression of eight CP genes responding to four insecticides were confirmed by RT-PCR analysis. The results showed that their overall expression profiles were consistent with RNA-seq analysis. The findings provide a basis for further functional investigation of CPs implied in insecticide stress in FAW.

Immune-Related Genes of Megalurothrips usitatus (Bagrall) Against Beauveria brongniartii and Akanthomyces attenuatus Identified Using RNA Sequencing

Megalurothrips usitatus (Bagrall) is an important pest of legumes worldwide, causing great economic loss every year. Beauveria brongniartii and Akanthomyces attenuatus have shown considerable pathogenicity against M. usitatus in our previous studies. The medial lethal concentration (LC₅₀) and the sublethal lethal concentration (LC₂₅) of B. brongniartii isolate SB010 against M. usitatus were 8.38 × 10⁵ and 1.73 × 10⁵ conidia mL⁻¹, respectively, whereas those of A. attenuatus isolate SCAUDCL-53 against M. usitatus were 4.37 × 10⁵ and 2.97 × 10⁴ conidia mL⁻¹, respectively. This study reports the transcriptome-based explanation of the stress responses of M. usitatus following the application of B. brongniartii and A. attenuatus. The analysis of the transcriptomic data revealed the expression of 254, 207, 195, and 234 immunity-related unigenes by M. usitatus in response to B. brongniartii LC₅₀ (SB1), B. brongniartii LC₂₅ (SB2), A. attenuatus LC₅₀ (V1), and A. attenuatus LC₂₅ (V2), respectively. The biological function and metabolic pathway analyses showed that these unigenes were mainly related to pattern recognition receptors, information transduction factors, and reaction factors, such as scavenger receptor, cytochrome b5, cuticle protein, lysozyme, and serine protease.

Comparison of gene expression in the red imported fire ant (Solenopsis invicta) under different temperature conditions

The red imported fire ant (RIFA), Solenopsis invicta Buren is native to South America and is known as a global problematic invasive species. This study focused on the molecular response of RIFA by comparing gene expression profiles after exposing ants to low (10 °C) and high (40 °C) temperature stress and comparing them to untreated controls (30 °C). A total of 99,085 unigenes (the clustered non-redundant transcripts that are filtered from the longest assembled contigs) were obtained, of which 19,154 were annotated with gene descriptions, gene ontology terms, and metabolic pathways. 86 gene ontology (GO) functional sub-groups and 23 EggNOG terms resulted. Differentially expressed genes (DEGs) with log₂FC ≥ 10 were screened and were compared at different temperatures. We found 203, 48, and 66 specific DEGs co-regulated at 10, 20, and 40 °C. Comparing transcriptome profiles for differential gene expression resulted in various DE genes, including cytochrome P450, NADH dehydrogenase subunit 1, cuticle protein and heat shock protein (HSP), which have previously been reported to be involved in cold and high temperature resistance. GO analysis revealed that antioxidant activity is up-regulated under high temperature stress. We verified the RNA-seq data by qPCR on 20 up- and down-regulated DEGs. These findings provide a basis for future understanding of the adaptation mechanisms of RIFA and the molecular mechanisms underlying the response to low and high temperatures.

Inotocin, a potential modulator of reproductive behaviours in a biparental beetle, Lethrus apterus

Several members of the highly conserved oxytocin/vasopressin neuropeptide family are involved in the regulation of reproductive and affiliative behaviours in numerous vertebrate and invertebrate species. Here we investigate gene expression patterns of inotocin, the insect ortholog of this peptide family, and its receptor to decipher their possible role in the control of reproductive behaviour in a beetle, Lethrus apterus, with biparental care. In an experiment performed on individuals of a wild population, we found that inotocin is not related to the control of water balance in this species because expression patterns did not change as a response to drought exposure. The expression levels of inotocin and its receptor, however, increased over the reproductive season i.e., when behaviour shifts from pair formation to parental care, suggesting that inotocin might be involved in the regulation of parental care in this insect. No difference was, however, found between sexes; a finding which might indicate that inotocin plays a similar role in both parents.

Transcriptome analysis of Apis mellifera under benomyl stress to discriminate the gene expression in response to development and immune systems

The health and safety of the honeybees are seriously threatened due to the abuse of chemical pesticides in modern agriculture and apiculture. In this study, the RNA Seq approach was used to assess the effects of the honeybees treated with benomyl. The results showed that there were a total of 11,902 differentially expressed genes (DEGs). Among them, 5,759 DEGs were up-regulated and involved in the functions of immunity, detoxification, biological metabolism, and regulation. The DEGs were clustered in the GO terms of epidermal structure and response to external stimuli, and most of the DEGs were enriched in 15 pathways, such as light conduction, MAPK, calcium ion pathway, and so on. Moreover, the pathway of the toll signal transduction was activated. The data investigated that the expression of functional genes involved in the growth, development, foraging, and immunity of honeybees were significantly affected by benomyl stress, which would seriously threaten the health of the honeybees. This study provided a theoretical basis for revealing the response mechanism of honeybees to pesticides stress.

De novo transcriptome assembly of the deep-sea hydrothermal vent, shrimp Rimicaris exoculate (Crustacea: Decapoda), from the south Mid-Atlantic Ridge

The deep-sea hydrothermal vent is a special ecosystem, which is different from terrestrial or coastal ecosystems. Rimicaris exoculata, which adapts well to several deep-sea hydrothermal vent environments, is the ideal model for studying hydrothermal vent fauna. In the present study, we obtained R. exoculata from a newly found hydrothermal vent in the south Mid-Atlantic Ridge, and the Illumina next-generation sequencing and de novo assembly were performed by Beijing Genomics Institution. A total of 17,258 annotated Unigenes were obtained. Several Unigenes associated with sulfide metabolism, which might contribute to well adaptation to high concentration of sulfide for R. exoculata, were annotated. This study is the first report on the high-throughput sequencing of R. exoculata. Our data can allow for further studies on the ability of R. exoculata adaptation to harsh conditions and provide abundant gene resources for research and development.

Comparative Proteomic Analysis of Rhipicephalus sanguineus sensu lato (Acari: Ixodidae) Tropical and Temperate Lineages: Uncovering Differences During Ehrlichia canis Infection

The tick vector Rhipicephalus sanguineus is established as a complex of closely related species with high veterinary-medical significance, in which the presence of different genetic, morphological, and biological traits has resulted in the recognition of different lineages within taxa. One of the most striking differences in the "temperate" and "tropical" lineages of R. sanguineus (s.l.) is the vector competence to Ehrlichia canis, suggesting that these ticks tolerate and react differently to pathogen infection. The present study addresses the SG and MG proteome of the R. sanguineus tropical and temperate lineages and compares their proteomic profile during E. canis infection. Batches of nymphs from the two lineages were allowed to feed on naïve and experimentally E. canis infected dogs and after molting, adults were dissected, and salivary glands and midgut tissues separated. Samples were screened for the presence of E. canis before proteomic analyses. The representation of the proteins identified in infected and non-infected tissues of each lineage was compared and gene ontology used for protein classification. Results highlight important differences in those proteomic profiles that added to previous reported genetic, biological, behavioral, and morphological differences, strengthening the hypothesis of the existence of two different species. Comparing infected and non-infected tissues, the results show that, while in midgut tissues the response to E. canis infection is similar in the salivary glands, the two lineages show a different pattern of protein representation. Focusing on the proteins found only in the infected condition, the data suggests that the cement cone produced during tick feeding may be implicated in pathogen infection. This study adds useful information to the debate on the controversial R. sanguineus systematic status, to the discussion related with the different vectorial competence occurring between the two lineages and identifies potential targets for efficient tick and tick-borne disease control.

Predator-induced stress responses in insects: A review

Predators can induce extreme stress and profound physiological responses in prey. Insects are the most dominant animal group on Earth and serve as prey for many different predators. Although insects have an extraordinary diversity of anti-predator behavioral and physiological responses, predator-induced stress has not been studied extensively in insects, especially at the molecular level. Here, we review the existing literature on physiological predator-induced stress responses in insects and compare what is known about insect stress to vertebrate stress systems. We conclude that many unrelated insects share a baseline pathway of predator-induced stress responses that we refer to as the octopamine-adipokinetic hormone (OAH) axis. We also present best practices for studying predator-induced stress responses in prey insects. We encourage investigators to compare neurophysiological responses to predator-related stress at the organismal, neurohormonal, tissue, and cellular levels within and across taxonomic groups. Studying stress-response variation between ecological contexts and across taxonomic levels will enable the field to build a holistic understanding of, and distinction between, taxon- and stimulus-specific responses relative to universal stress responses.

Effects of pesticide exposure and predation risk on nutrient cycling and primary production

Understanding how pesticides and natural stressors shape ecosystem functions remains a major challenge. A largely overlooked way how stressors may affect nutrient cycling and primary production is through effects on body stoichiometry and the egestion of elements. We investigated how exposure to the pesticide chlorpyrifos and to predation risk, an abundant natural stressor in aquatic systems, altered the stoichiometry of the bodies and the egested faecal pellets of Enallagma cyathigerum damselfly larvae and how this further cascaded into effects on primary production (algae growth). Chlorpyrifos exposure reduced egestion rates while predation risk had no effect. Chlorpyrifos exposure and predation risk affected both elemental composition of bodies and faecal pellets, and this in an additive way. Chlorpyrifos exposure increased body C(carbon), N(nitrogen), and P(phosphorous) contents, and increased the C content of the faecal pellets. Predation risk induced an increase of the N content, resulting in a decreased C:N ratio, of both the bodies and faecal pellets. The changes in the composition of the faecal pellets caused by predation risk but not by chlorpyrifos exposure increased algae growth under control conditions. This indicated that algae growth was N limited. Our results provide an important proof-of-principle how a stressor may shape nutrient cycling and subsequently primary productivity.

TMT-Based Quantitative Proteomic Profiling of Overwintering Lissorhoptrus oryzophilus

Adaptations to low temperature play a critical role in restricting the geographical distribution of insects. Decreasing day lengths and temperatures trigger seasonal cold adaptations in insects. These adaptions include changes in expression at the miRNA, mRNA and protein levels. The rice water weevil (RWW), Lissorhoptrus oryzophilus, introduced from the Mississippi River, is a globally invasive pest of wetland rice that can survive at the northern border of China. To investigate the changes in expression at the protein level in overwintering female RWW adults, 6-plex tandem mass tags (TMTs) were used in overwintering and summer adults. By using a proteome database available for Curculionidae, 1077 proteins were quantified, 183 of which differed significantly between the overwintering and summer samples. To further understand these differentially expressed proteins (DEPs), bioinformatics analyses such as gene ontology (GO) enrichment analyses were performed. DEPs associated with the terms binding, structural molecule activity, catalytic activity, multicellular organismal process, extracellular region, chitin binding, metabolic process, intracellular part and organic cyclic compound binding were altered by selection during winter. The changes in the expression of these proteins suggest that the proteins are important for RWW survival in winter.

Repurposing of Glycine-Rich Proteins in Abiotic and Biotic Stresses in the Lone-Star Tick (Amblyomma americanum)

Tick feeding requires the secretion of a huge number of pharmacologically dynamic proteins and other molecules which are vital for the formation of the cement cone, the establishment of the blood pool and to counter against the host immune response. Glycine-rich proteins (GRP) are found in many organisms and can function in a variety of cellular processes and structures. The functional characterization of the GRPs in the tick salivary glands has not been elucidated. GRPs have been found to play a role in the formation of the cement cone; however, new evidence suggests repurposing of GRPs in the tick physiology. In this study, an RNA interference approach was utilized to silence two glycine-rich protein genes expressed in early phase of tick feeding to determine their functional role in tick hematophagy, cement cone structure, and microbial homeostasis within the tick host. Additionally, the transcriptional regulation of GRPs was determined after exposure to biotic and abiotic stresses including cold and hot temperature, injury, and oxidative stress. This caused a significant up-regulation of AamerSigP-34358, Aam-40766, AamerSigP-39259, and Aam-36909. Our results suggest ticks repurpose these proteins and further functional characterization of GRPs may help to design novel molecular strategies to disrupt the homeostasis and the pathogen transmission.

Applying microarray-based technique to study and analyze silkworm (Bombyx mori) transcriptomic response to long-term high iron diet

To investigate the biological processes affected by long-term iron supplementation, newly hatched silkworms were exposed to high iron mulberry diet (10 and 100 ppm) and its effect on silkworm transcriptom was determined. The results showed that the silkworm was responsive to iron by increasing iron concentration and ferritin levels in the hemolymph and by regulating the expression of many other genes. A total of 523 and 326 differentially expressed genes were identified in 10 and 100 ppm Fe group compared to the control, respectively. Of these genes, 249 were shared between in both the 10 ppm and 100 ppm Fe group, including 152 up-regulated and 97 down-regulated genes. These shared genes included 19 known Fe regulated, 24 immune-related, 12 serine proteases and serine proteases homologs, 41 cuticular and cuticle genes. Ten genes (carboxypeptidases A, serine protease homologs 85, fibrohexamerin/P25, transferrin, sex-specific storage-protein 2, fungal protease inhibitor F, insect intestinal mucin, peptidoglycan recognition protein B, cuticle protein CPH45, unknown gene) were involved in the regulation of iron overload responses.

Cuticular protein defective Bamboo mutant of Bombyx mori is sensitive to environmental stresses

Insect cuticle acts as a primary protective barrier against environment stresses that may directly impact the insect body. Here, we report the mechanical defense function of a structural cuticular protein, BmorCPH24, to environmental stresses using a silkworm Bamboo (Bo) mutant with this gene mutation. Ultraviolet (UV) irradiation and topical application of an acetone insecticide were used as environmental stresses to determine the differences in susceptibility between Bo and wild-type larvae. UV irradiation resulted in a sunburn phenotype in the Bo strains earlier than the wild-type indicating the sensitivity of Bo. Higher malondialdehyde (MDA) content and a lower survival ratio were also observed in the Bo strains. Treatment with deltamethrin revealed that Bo larvae were more sensitive to insecticides than the wild-type. Furthermore, cuticle analysis by microsection revealed thinner cuticle and a significant decrease in the endocuticle layer (∼64.0%) in Bo. These results suggest that BmorCPH24 mutation can lead to deficiency in resources required to construct the cuticle in Bo resulting in thin cuticle and reduced resistance to UV and insecticides. These results provide us new insight into the role of structural cuticular proteins in insect cuticle against environment stresses.

Genome-wide annotation of cuticular proteins in the oriental fruit fly (Bactrocera dorsalis), changes during pupariation and expression analysis of CPAP3 protein genes in response to environmental stresses

Cuticular proteins (CPs) are essential components of the insect cuticle as they create a structural and protective shield and may have a role in insect development. In this paper, we studied the CPs in the oriental fruit fly (Bactrocera dorsalis), one of the most economically important pests in the Tephritidae family around the world. The availability of a complete genome sequence (NCBI Assembly: ASM78921v2) allowed the identification of 164 CP genes in B. dorsalis. Comparative analysis of the CPs in B. dorsalis with those in the model insect Drosophila melanogaster and the closely related Ceratitis capitata, and CPs from mosquitoes, Lepidoptera, Hymenoptera and Coleoptera identified Diptera-specific genes and cuticle development patterns. Analysis of their evolutionary relationship revealed that some CP families had evolved according to the phylogeny of the different insect species, while others shared a closer relationship based on domain architecture. Subsequently, transcriptome analysis showed that while most of the CPs (60-100% of the family members) are expressed in the epidermis, some were also present in internal organs such as the fat body and the reproductive organs. Furthermore, the study of the expression profiles throughout development revealed a profound change in the expression of CPs during the formation of the puparium (pupariation). Further analysis of the expression profiles of the CPAP3 genes under various environmental stresses revealed them to be involved in the response to pesticides and arid and extreme temperatures conditions. In conclusion, the data provide a particular overview of CPs and their evolutionary and transcriptional dynamics, and in turn they lay a molecular foundation to explore their roles in the unique developmental process of insect metamorphosis and stress responses.

ASGDB: a specialised genomic resource for interpreting Anopheles sinensis insecticide resistance

BACKGROUND

Anopheles sinensis is an important malaria vector in Southeast Asia. The widespread emergence of insecticide resistance in this mosquito species poses a serious threat to the efficacy of malaria control measures, particularly in China. Recently, the whole-genome sequencing and de novo assembly of An. sinensis (China strain) has been finished. A series of insecticide-resistant studies in An. sinensis have also been reported. There is a growing need to integrate these valuable data to provide a comprehensive database for further studies on insecticide-resistant management of An. sinensis.

RESULTS

A bioinformatics database named An. sinensis genome database (ASGDB) was built. In addition to being a searchable database of published An. sinensis genome sequences and annotation, ASGDB provides in-depth analytical platforms for further understanding of the genomic and genetic data, including visualization of genomic data, orthologous relationship analysis, GO analysis, pathway analysis, expression analysis and resistance-related gene analysis. Moreover, ASGDB provides a panoramic view of insecticide resistance studies in An. sinensis in China. In total, 551 insecticide-resistant phenotypic and genotypic reports on An. sinensis distributed in Chinese malaria-endemic areas since the mid-1980s have been collected, manually edited in the same format and integrated into OpenLayers map-based interface, which allows the international community to assess and exploit the high volume of scattered data much easier. The database has been given the URL: http://www.asgdb.org /.

CONCLUSIONS

ASGDB was built to help users mine data from the genome sequence of An. sinensis easily and effectively, especially with its advantages in insecticide resistance surveillance and control.

Microscopic cuticle structure comparison of pupal melanic and wild strain of Spodoptera exigua and their gene expression profiles in three time points

The beet armyworm, Spodoptera exigua (Hubner), is one of the major crop pests and is a target for current pest control approaches using insecticides. S. exigua melanic mutants (SEM) spontaneously occurred in the S. exigua wild type (SEW) strain and have been maintained under laboratory conditions on an artificial diet. Scanning electron microscopy showed that the inner cuticle of the SEM had a denser and less orderly structure. We investigated the cuticle protein genes using RNA-seq at three different developmental stages of both SEM and SEW. Comparison of cDNA libraries showed that 7257 CPs were significantly up-regulated and 664 genes were significantly downregulated in SEM at the developmental stage of 46-h in the fifth instar. In addition, 460 genes were significantly up-regulated and 439 genes were significantly down-regulated in the SEM at the development stage of 4-h before pupation. Moreover, 162 genes were significantly up-regulated and 293 genes were significantly downregulated in the SEM, just after pupation. Two genes CPR63 and CPR97 were identified from RNA sequences to verify the differentially expressed gene (DEG) results through quantitative real-time PCR (qRT-PCR). The results show that expression of both CPR63 and CPR97 structural cuticular proteins were significantly different between SEM and SEW. This functional analysis may help in understanding the role that these genes play in the cuticle pattern of the SEM.

Characterization of a glycine-rich protein from Rhipicephalus microplus: tissue expression, gene silencing and immune recognition

Salivary molecules, as glycine-rich proteins (GRPs), are essential to tick attachment and feeding on the host and are suggested to be involved in the host's immune system evasion, therefore representing natural candidates in the search for protective vaccine antigens. This work shows the molecular characterization of a GRP from Rhipicephalus microplus (RmGRP). The cDNA and putative amino acid sequences were analysed, as well as the transcription level in tick tissues/developmental stages, showing the highest levels of gene expression in 1-day-old larvae and salivary glands of fully engorged females. RmGRP gene silencing resulted in a lower hatching rate of larvae from treated females. In addition, recombinant RmGRP (rRmGRP) was recognized by sera from naturally and experimentally infested bovines, displaying considerable differences among the individuals tested. rRmGRP was recognized by anti-saliva and anti-salivary glands sera, while anti-rRmGRP serum recognized RmGRP in saliva and salivary glands, indicating its secretion into the host. The data collected indicate that RmGRP may present roles other than in the tick–host relationship, especially in embryo development. In addition, the high expression in adult females, antigenicity and presence of shared characteristics with other tick protective GRPs turns RmGRP a potential candidate to compose an anti-tick vaccine cocktail.

Differential transcriptome analysis reveals genes related to cold tolerance in seabuckthorn carpenter moth, Eogystia hippophaecolus

Seabuckthorn carpenter moth, Eogystia hippophaecolus (Lepidoptera: Cossidae), is an important pest of sea buckthorn (Hippophae rhamnoides), which is a shrub that has significant ecological and economic value in China. E. hippophaecolus is highly cold tolerant, but limited studies have been conducted to elucidate the molecular mechanisms underlying its cold resistance. Here we sequenced the E. hippophaecolus transcriptome using RNA-Seq technology and performed de novo assembly from the short paired-end reads. We investigated the larval response to cold stress by comparing gene expression profiles between treatments. We obtained 118,034 unigenes, of which 22,161 were annotated with gene descriptions, conserved domains, gene ontology terms, and metabolic pathways. These resulted in 57 GO terms and 193 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. By comparing transcriptome profiles for differential gene expression, we identified many differentially expressed proteins and genes, including heat shock proteins and cuticular proteins which have previously been reported to be involved in cold resistance of insects. This study provides a global transcriptome analysis and an assessment of differential gene expression in E. hippophaecolus under cold stress. We found seven differential expressed genes in common between developmental stages, which were verified with qPCR. Our findings facilitate future genomic studies aimed at improving our understanding of the molecular mechanisms underlying the response of insects to low temperatures.

Osiris9a is a major component of silk fiber in lepidopteran insects

In a previous high-throughput proteomics study, it was found that the silkworm cocoon contains hundreds of complex proteins, many of which have unknown functions, in addition to fibroins, sericins, and some protease inhibitors. Osiris was one of the proteins with no known function. In this study, we identified the Osiris gene family members and constructed a phylogenetic tree based on the sequences from different species. Our results indicate that the Osiris9 gene subfamily contains six members; it is specifically expressed in lepidopteran insects and has evolved by gene duplication. An Osiris gene family member from Bombyx mori was designated as BmOsiris9a (BmOsi9a) on the basis of its homology to Drosophila melanogaster Osiris9. The expression pattern of BmOsi9a showed that it was highly expressed only in the middle silk gland of silkworm larvae, similar to Sericin1 (Ser1). BmOsi9a was visualized as two bands in western blot analysis, implying that it probably undergoes post-translational modifications. Immunohistochemistry analysis revealed that BmOsi9a was synthesized and secreted into the lumen of the middle silk gland, and was localized in the sericin layer in the silk fiber. BmOsi9a was found in the silk fibers of not only three Bombycidae species, viz. B. mori, B. mandarina, and B. huttoni, but also in the fibers collected from Saturniidae species, including Antheraea assama, Antheraea mylitta, and Samia cynthia. Although the exact biological function of Osi9a in the silk fibers is unknown, our results are important because they demonstrate that Osi9a is a common structural component of silk fiber and is expressed widely among the silk-producing Bombycidae and Saturniidae insects. Our results should help in understanding the role of Osi9a in silk fibers.

Comparative transcriptome analysis of Rimicaris sp. reveals novel molecular features associated with survival in deep-sea hydrothermal vent

Shrimp of the family Alvinocarididae are the predominant megafauna of deep-sea hydrothermal vents. However, genome information on this family is currently unavailable. In the present study, by employing Illumina sequencing, we performed the first de novo transcriptome analysis of the gills of the shrimp Rimicaris sp. from the hydrothermal vent in Desmos, Manus Basin. The analysis was conducted in a comparative manner with the shrimp taken directly from the vent (GR samples) and the shrimp that had been maintained for ten days under normal laboratory condition (mGR samples). Among the 128,938 unigenes identified, a large number of differentially expressed genes (DEGs) between the GR and mGR samples were detected, including 2365 and 1607 genes significantly upregulated and downregulated, respectively, in GR. The DEGs covered diverse functional categories. Most of the DEGs associated with immunity were downregulated in GR, while most of the DEGs associated with sulfur metabolism and detoxification were upregulated in GR. These results provide the first comprehensive transcriptomic resource for hydrothermal vent Rimicaris and revealed varied categories of genes likely involved in deep-sea survival.

Implicating ABC Transporters in Insecticide Resistance: Research Strategies and a Decision Framework

Pest insects damage crops, transmit diseases, and are household nuisances. Historically, they have been controlled with insecticides, but overuse often leads to resistance to one or more of these chemicals. Insects gain resistance to insecticides through behavioral, metabolic, genetic, and physical mechanisms. One frequently overlooked strategy is through the use of ATP-binding cassette (ABC) transporters. ABC transporters, present in all domains of life, perform natural excretory functions, thus the exploitation of these transporters to excrete insecticides and contribute to resistance is highly plausible. Previous work has implicated ABC transporters in some cases of insecticide resistance. Proposed herein is a framework meant as a formal guide for more easily incorporating the analysis of ABC transporters into existing resistance monitoring using suggested simple research methods. This framework functions as a simple decision tree and its utility is demonstrated using case examples. Determining a role for ABC transporters in insecticide resistance would help to shape future resistance management plans and guide the design of new insecticides.

DIRProt: a computational approach for discriminating insecticide resistant proteins from non-resistant proteins

BACKGROUND

Insecticide resistance is a major challenge for the control program of insect pests in the fields of crop protection, human and animal health etc. Resistance to different insecticides is conferred by the proteins encoded from certain class of genes of the insects. To distinguish the insecticide resistant proteins from non-resistant proteins, no computational tool is available till date. Thus, development of such a computational tool will be helpful in predicting the insecticide resistant proteins, which can be targeted for developing appropriate insecticides.

RESULTS

Five different sets of feature viz., amino acid composition (AAC), di-peptide composition (DPC), pseudo amino acid composition (PAAC), composition-transition-distribution (CTD) and auto-correlation function (ACF) were used to map the protein sequences into numeric feature vectors. The encoded numeric vectors were then used as input in support vector machine (SVM) for classification of insecticide resistant and non-resistant proteins. Higher accuracies were obtained under RBF kernel than that of other kernels. Further, accuracies were observed to be higher for DPC feature set as compared to others. The proposed approach achieved an overall accuracy of >90% in discriminating resistant from non-resistant proteins. Further, the two classes of resistant proteins i.e., detoxification-based and target-based were discriminated from non-resistant proteins with >95% accuracy. Besides, >95% accuracy was also observed for discrimination of proteins involved in detoxification- and target-based resistance mechanisms. The proposed approach not only outperformed Blastp, PSI-Blast and Delta-Blast algorithms, but also achieved >92% accuracy while assessed using an independent dataset of 75 insecticide resistant proteins.

CONCLUSIONS

This paper presents the first computational approach for discriminating the insecticide resistant proteins from non-resistant proteins. Based on the proposed approach, an online prediction server DIRProt has also been developed for computational prediction of insecticide resistant proteins, which is accessible at http://cabgrid.res.in:8080/dirprot/ . The proposed approach is believed to supplement the efforts needed to develop dynamic insecticides in wet-lab by targeting the insecticide resistant proteins.

Temporal patterns of imidacloprid resistance throughout a growing season in Leptinotarsa decemlineata populations

BACKGROUND

The Colorado potato beetle, Leptinotarsa decemlineata (Say), is a major agricultural pest of commercial potatoes. Pest managers use a combination of control tactics to limit populations, including multiple insecticides. Finding a window of insecticide susceptibility and understanding genetic responses to insecticide exposure during a growing season may provide novel management recommendations for L. decemlineata.

RESULTS

We examined temporal changes (during one growing season) in phenotypic response between a susceptible population and an imidacloprid-resistant population. Beetles remained more susceptible to imidacloprid in the susceptible population throughout the growing season. Estimated mean LC₅₀ values varied throughout the growing season in the resistant population, with increased susceptibility among overwintered and recently emerged adult beetles compared with a heightened level of resistance in the second generation. RNA transcript abundance was compared among multiple time points through the growing season, showing that cuticular proteins and cytochrome p450s were highly upregulated during peaks of measured resistance.

CONCLUSION

Temporal variation in imidacloprid susceptibility of L. decemlineata was observed, which included early time points of susceptibility and later peaks in resistance. Heightened resistance occurred during the second generation and correlated to increased transcript abundance of multiple mechanisms of resistance, including multiple cuticular protein and cytochrome p450 transcripts. © 2016 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Under-Expression of Chemosensory Genes in Domiciliary Bugs of the Chagas Disease Vector Triatoma brasiliensis

BACKGROUND

In Latin America, the bloodsucking bugs Triatominae are vectors of Trypanosoma cruzi, the parasite that causes Chagas disease. Chemical elimination programs have been launched to control Chagas disease vectors. However, the disease persists because native vectors from sylvatic habitats are able to (re)colonize houses-a process called domiciliation. Triatoma brasiliensis is one example. Because the chemosensory system allows insects to interact with their environment and plays a key role in insect adaption, we conducted a descriptive and comparative study of the chemosensory transcriptome of T. brasiliensis samples from different ecotopes.

METHODOLOGY/PRINCIPAL FINDING

In a reference transcriptome built using de novo assembly, we found transcripts encoding 27 odorant-binding proteins (OBPs), 17 chemosensory proteins (CSPs), 3 odorant receptors (ORs), 5 transient receptor potential channel (TRPs), 1 sensory neuron membrane protein (SNMPs), 25 takeout proteins, 72 cytochrome P450s, 5 gluthatione S-transferases, and 49 cuticular proteins. Using protein phylogenies, we showed that most of the OBPs and CSPs for T. brasiliensis had well supported orthologs in the kissing bug Rhodnius prolixus. We also showed a higher number of these genes within the bloodsucking bugs and more generally within all Hemipterans compared to the other species in the super-order Paraneoptera. Using both DESeq2 and EdgeR software, we performed differential expression analyses between samples of T. brasiliensis, taking into account their environment (sylvatic, peridomiciliary and domiciliary) and sex. We also searched clusters of co-expressed contigs using HTSCluster. Among differentially expressed (DE) contigs, most were under-expressed in the chemosensory organs of the domiciliary bugs compared to the other samples and in females compared to males. We clearly identified DE genes that play a role in the chemosensory system.

CONCLUSION/SIGNIFICANCE

Chemosensory genes could be good candidates for genes that contribute to adaptation or plastic rearrangement to an anthropogenic system. The domiciliary environment probably includes less diversity of xenobiotics and probably has more stable abiotic parameters than do sylvatic and peridomiciliary environments. This could explain why both detoxification and cuticle protein genes are less expressed in domiciliary bugs. Understanding the molecular basis for how vectors adapt to human dwellings may reveal new tools to control disease vectors; for example, by disrupting chemical communication.

Characterization of the transcriptomes and cuticular protein gene expression of alate adult, brachypterous neotenic and adultoid reproductives of Reticulitermes labralis

The separation of primary reproductive and secondary reproductive roles based on the differentiation of alate adults and neotenic reproductives is the most prominent characteristic of termites. To clarify the mechanism underlying this differentiation, we sequenced the transcriptomes of alate adults (ARs), brachypterous neotenics (BNs) and adultoid reproductives (ANs) from the last instar nymphs of Reticulitermes labralis. A total of 404,152,188 clean sequencing reads was obtained and 61,953 unigenes were assembled. Of the 54 identified cuticular protein (CP) genes of the reproductives, 22 were classified into the CPR family and 7 were classified into the CPG family. qRT-PCR analyses of the 6 CP genes revealed that the CP genes involved in exocuticle sclerotization were highly expressed in the ARs and RR-1 involved in soft endocuticle was highly expressed in the ARs and ANs. These results suggest that the alate adults might increase cuticular component deposition to adapt to new or changing environments and that the development of reproductive individuals into primary or secondary reproductives is controlled by the expression of cuticular protein genes involved in the hardening of the exocuticle. In addition, the AN caste is a transitional type between the BN and AR castes in the process of evolution.

Identification of Immunity-Related Genes in Dialeurodes citri against Entomopathogenic Fungus Lecanicillium attenuatum by RNA-Seq Analysis

Dialeurodes citri is a major pest in citrus producing areas, and large-scale outbreaks have occurred increasingly often in recent years. Lecanicillium attenuatum is an important entomopathogenic fungus that can parasitize and kill D. citri. We separated the fungus from corpses of D. citri larvae. However, the sound immune defense system of pests makes infection by an entomopathogenic fungus difficult. Here we used RNA sequencing technology (RNA-Seq) to build a transcriptome database for D. citri and performed digital gene expression profiling to screen genes that act in the immune defense of D. citri larvae infected with a pathogenic fungus. De novo assembly generated 84,733 unigenes with mean length of 772 nt. All unigenes were searched against GO, Nr, Swiss-Prot, COG, and KEGG databases and a total of 28,190 (33.3%) unigenes were annotated. We identified 129 immunity-related unigenes in transcriptome database that were related to pattern recognition receptors, information transduction factors and response factors. From the digital gene expression profile, we identified 441 unigenes that were differentially expressed in D. citri infected with L. attenuatum. Through calculated Log₂Ratio values, we identified genes for which fold changes in expression were obvious, including cuticle protein, vitellogenin, cathepsin, prophenoloxidase, clip-domain serine protease, lysozyme, and others. Subsequent quantitative real-time polymerase chain reaction analysis verified the results. The identified genes may serve as target genes for microbial control of D. citri.

Comparative proteomic analysis in Aphis glycines Mutsumura under lambda-cyhalothrin insecticide stress

Lambda-cyhalothrin is now widely used in China to control the soybean aphid Aphis glycines. To dissect the resistance mechanism, a laboratory-selected resistant soybean aphid strain (CRR) was established with a 43.42-fold resistance ratio to λ-cyhalothrin than the susceptible strain (CSS) in adult aphids. In this study, a comparative proteomic analysis between the CRR and CSS strains revealed important differences between the susceptible and resistant strains of soybean aphids for λ-cyhalothrin. Approximately 493 protein spots were detected in two-dimensional polyacrylamide gel electrophoresis (2-DE). Thirty-six protein spots displayed differential expression of >2-fold in the CRR strain compared to the CSS strain. Out of these 36 protein spots, 21 had elevated and 15 had decreased expression. Twenty-four differentially expressed proteins were identified by MALDI TOF MS/MS and categorized into the functional groups cytoskeleton-related protein, carbohydrate and energy metabolism, protein folding, antioxidant system, and nucleotide and amino acid metabolism. Function analysis showed that cytoskeleton-related proteins and energy metabolism proteins have been associated with the λ-cyhalothrin resistance of A. glycines. The differential expression of λ-cyhalothrin responsive proteins reflected the overall change in cellular structure and metabolism after insecticide treatment in aphids. In summary, our studies improve understanding of the molecular mechanism resistance of soybean aphid to lambda-cyhalothrin, which will facilitate the development of rational approaches to improve the management of this pest and to improve the yield of soybean.

QTLs for Woolly Poplar Aphid (Phloeomyzus passerinii L.) Resistance Detected in an Inter-Specific Populus deltoides x P. nigra Mapping Population

The genus Populus represents one of the most economically important groups of forest trees. It is composed by approximately 30 species used for wood and non-wood products, phytoremediation and biomass. Poplar is subjected to several biological and environmental threats although, compared to annual crops, we know far less about the genetic bases of biotic stress resistance. Woolly poplar aphid (Phloeomyzus passerinii) is considered a main pest of cultivated poplars in European and American countries. In this work we present two high density linkage maps in poplar obtained by a genotyping by sequencing (GBS) approach and the identification of QTLs involved in Ph. passerinii resistance. A total of 5,667 polymorphic markers (5,606 SNPs and 61 SSRs) identified on expressed sequences have been used to genotype 131 plants of an F1 population P ×canadensis obtained by an interspecific mate between Populus deltoides (resistant to woolly poplar aphid) and Populus nigra (susceptible to woolly poplar aphid). The two linkage maps, obtained following the two-way pseudo-testcross mapping strategy, have been used to investigate the genetic bases of woolly poplar aphid resistance. One major QTL and two QTLs with minor effects (mapped on LGV, LGXVI and LG XIX) explaining the 65.8% of the genetic variance observed in the progeny in response to Ph. passerinii attack were found. The high density coverage of functional markers allowed the identification of three genes belonging to disease resistance pathway as putative candidates for P. deltoides resistance to woolly poplar aphid. This work is the first report on genetic of woolly poplar aphid genetic resistance and the resistant loci associated markers identified represent a valuable tool in resistance poplar breeding programs.

Conifer flavonoid compounds inhibit detoxification enzymes and synergize insecticides

Detoxification by glutathione S-transferases (GSTs) and esterases are important mechanisms associated with insecticide resistance. Discovery of novel GST and esterase inhibitors from phytochemicals could provide potential new insecticide synergists. Conifer tree species contain flavonoids, such as taxifolin, that inhibit in vitro GST activity. The objectives were to test the relative effectiveness of taxifolin as an enzyme inhibitor and as an insecticide synergist in combination with the organophosphorous insecticide, Guthion (50% azinphos-methyl), and the botanical insecticide, pyrethrum, using an insecticide-resistant Colorado potato beetle (CPB) Leptinotarsa decemlineata (Say) strain. Both taxifolin and its isomer, quercetin, increased the mortality of 1(st) instar CPB larvae after 48h when combined with Guthion, but not pyrethrum. Taxifolin had greater in vitro esterase inhibition compared with the commonly used esterase inhibitor, S, S, S-tributyl phosphorotrithioate (DEF). An in vivo esterase and GST inhibition effect after ingestion of taxifolin was measured, however DEF caused a greater suppression of esterase activity. This study demonstrated that flavonoid compounds have both in vitro and in vivo esterase inhibition, which is likely responsible for the insecticide synergism observed in insecticide-resistant CPB.

Starvation-responsive glycine-rich protein gene in the silkworm Bombyx mori

Four glycine-rich protein (GRP) genes were identified from expressed sequence tags of the maxillary galea of the silkworm. All four genes were expressed in the maxillary pulp, antenna, labrum, and labium, but none of the genes were expressed in most internal organs. Expression of one of the genes, termed bmSIGRP, was further increased approximately fivefold in the mouth region (including the maxilla, antenna, labrum, labium, and mandible) after 24 h of starvation. bmSIGRP expression peaked at 24 h and gradually declined during the subsequent 2 days. When a synthetic diet not containing proteins was fed, bmSIGRP expression increased significantly in the mouth region to levels similar to that observed in starved larvae. Synthetic diets that lacked vitamins or salts but contained amino acids did not significantly affect bmSIGRP expression. These results suggest that amino acid depletion increases bmSIGRP expression.

Gene expression profile of Bombyx mori hemocyte under the stress of destruxin A

Destruxin A (DA) is a cyclo-peptidic mycotoxin from the entomopathogenic fungus Metarhizium anisopliae. To uncover potential genes associated with its molecular mechanisms, a digital gene expression (DGE) profiling analysis was used to compare differentially expressed genes in the hemocytes of silkworm larvae treated with DA. Ten DGE libraries were constructed, sequenced, and assembled, and the unigenes with least 2.0-fold difference were further analyzed. The numbers of up-regulated genes were 10, 20, 18, 74 and 8, as well as the numbers of down-regulated genes were 0, 1, 8, 13 and 3 at 1, 4, 8, 12 and 24 h post treatment, respectively. Totally, the expression of 132 genes were significantly changed, among them, 1, 3 and 12 genes were continually up-regulated at 4, 3 and 2 different time points, respectively, while 1 gene was either up or down-regulated continually at 2 different time points. Furthermore, 68 genes were assigned to one or multiple gene ontology (GO) terms and 89 genes were assigned to specific Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology. In-depth analysis identified that these genes putatively involved in insecticide resistance, cell apoptosis, and innate immune defense. Finally, twenty differentially expressed genes were randomly chosen and validated by quantitative real-time PCR (qRT-PCR). Our studies provide insights into the toxic effect of this microbial insecticide on silkworm's hemocytes, and are helpful to better understanding of the molecular mechanisms of DA as a biological insecticide.

Cuticular differences associated with aridity acclimation in African malaria vectors carrying alternative arrangements of inversion 2La

BACKGROUND

Principal malaria vectors in Africa, An. gambiae and An. coluzzii, share an inversion polymorphism on the left arm of chromosome 2 (2La/2L+a) that is distributed non-randomly in the environment. Genomic sequencing studies support the role of strong natural selection in maintaining steep clines in 2La inversion frequency along environmental gradients of aridity, and physiological studies have directly implicated 2La in heat and desiccation tolerance, but the precise genetic basis and the underlying behavioral and physiological mechanisms remain unknown. As the insect cuticle is the primary barrier to water loss, differences in cuticle thickness and/or epicuticular waterproofing associated with alternative 2La arrangements might help explain differences in desiccation resistance.

METHODS

To test that hypothesis, two subcolonies of both An. gambiae and An. coluzzii were established that were fixed for alternative 2La arrangements (2La or 2L+a) on an otherwise homosequential and shared genetic background. Adult mosquitoes reared under controlled environmental conditions (benign or arid) for eight days post-eclosion were collected and analyzed. Measurements of cuticle thickness were made based on scanning electron microscopy, and cuticular hydrocarbon (CHC) composition was evaluated by gas chromatography-mass spectrometry.

RESULTS

After removing the allometric effects of body weight, differences in mean cuticle thickness were found between alternative 2La karyotypes, but not between alternative environments. Moreover, the thicker cuticle of the An. coluzzii 2La karyotype was contrary to the known higher rate of water loss of this karyotype relative to 2L+a. On the other hand, quantitative differences in individual CHCs and overall CHC profiles between alternative karyotypes and environmental conditions were consistent with expectation based on previous physiological studies.

CONCLUSIONS

Our results suggest that alternative arrangements of the 2La inversion are associated with differences in cuticle thickness and CHC composition, but that only CHC composition appears to be relevant for desiccation resistance. Differences in the CHC composition were consistent with previous findings of a lower rate of water loss for the 2L+a karyotype at eight days post-eclosion, suggesting that CHC composition is an important strategy for maintaining water balance in this genetic background, but not for 2La. Despite a higher rate of water loss at eight days, higher body water content of the 2La karyotype confers a level of desiccation resistance equivalent to that of the 2L+a karyotype.

Extreme temperatures in the adult stage shape delayed effects of larval pesticide stress: a comparison between latitudes

Global warming and pesticide pollution are major threats for aquatic biodiversity. Yet, how pesticide effects are influenced by the increased frequency of extreme temperatures under global warming and how local thermal adaptation may mitigate these effects is unknown. We therefore investigated the combined impact of larval chlorpyrifos exposure, larval food stress and adult heat exposure on a set of fitness-related traits in replicated low- and high-latitude populations of the damselfly Ischnura elegans. Larval pesticide exposure resulted in lighter adults with a higher water content, lower fat content, higher Hsp70 levels and a lower immune function (PO activity). Heat exposure reduced water content, mass, fat content and flying ability. Importantly, both stressors interacted across metamorphosis: adult heat exposure lowered the reduction of fat content, and generated a stronger decrease in PO activity in pesticide-exposed animals. Larval pesticide exposure and larval food stress also reduced the defense response to the adult heat stress in terms of increased Hsp70 levels. In line with strong life history differences in the unstressed control situation, high-latitude animals were less sensitive to food stress (body mass and water content), but more sensitive to pesticide stress (development time and PO activity) and heat exposure (PO activity and Hsp70 levels). While low-latitude adults could better withstand the extreme temperature as suggested by the weaker increase in Hsp70, heat exposure similarly affected the delayed effects of larval pesticide exposure at both latitudes. Our study highlighted two key findings relevant for ecological risk assessment under global warming. Firstly, the delayed effects of larval pesticide exposure on adult damselflies depended upon subsequent adult heat exposure, indicating that larval pesticide stress and adult heat stress interacted across metamorphosis. Secondly, low- and high-latitude animals responded differently to the imposed stressors, highlighting that intraspecific evolution along natural thermal gradients may shape sensitivity to pesticides.