Transcriptome analysis of the Chinese white wax scale Ericerus pela with focus on genes involved in wax biosynthesis

Variation of gene expression of fatty acid acyl CoA reductase associated with wax secretion of a scale insect, Ericerus pela, and identification of its regulation factors through the accessible chromatin analyses and yeast one-hybrid

The Chinese white wax scale insect (CWWSI), Ericerus pela, can secret an amount of wax equivalent to their body weight. Previous studies demonstrated the fatty acyl-CoA reductase (far3) plays a pivotal role in wax secretion of CWWSI. The high expression of far3 is crucial for the massive wax secretion. However, the transcription regulation of far3 was not clear. To identify regulatory factors that control the expression of far3, the assay for transposase-accessible chromatin (ATAC) and yeast one-hybrid (Y1H) were carried out in this study. The ATAC sequencing of the CWWSI at the early wax-secretion stage ATAC-seq resulted in 22.75 GB raw data, generated 75,827,225 clean reads and revealed 142,771 peaks. There was one significant peak in the 3 kb upstream regulation regions. The peak sequence is located between -1000 and -670 bp upstream of the far3 transcription start site, spanning a length of 331 bp. This peak sequence served as bait for creating the pAbAi-peak recombinant vector, used in Y1H screenings to identify proteins interacting with far3 gene. The results indicate a successful CWWSI cDNA library construction with a capacity of 1.2 × 107 colony forming unit, a 95.8% recombination rate, and insert sizes between 1,000 and 2,000 bp. Self-activation tests established that 100 ng/mL of AbA effectively inhibited bait vector self-activation. Finally, a total of 88 positive clones were selected. After sequencing and removal of duplication, 63 unique clones were obtained from these screened colonies. By aligning the clone sequences with full-length transcriptome and genome of CWWSI, the full-length coding sequences of these clones were obtained. BlastX analysis identified a transcription factor, nuclear transcription factor Y beta, and two co-activators, cAMP-response-element-binding-protein-binding protein and WW domain binding protein 2. Reverse transcription quantitative polymerase chain reaction analysis confirmed that their expression patterns were consistent with the developmental stages preceding wax secretion and matched the wax secretion characteristics during ovulation periods. These results are beneficial for further research into the regulatory mechanisms of wax secretion of CWWSI.

The transcriptome of Icerya aegyptiaca (Hemiptera: Monophlebidae) and comparison with neococcoids reveal genetic clues of evolution in the scale insects

BACKGROUND

Scale insects are worldwide sap-sucking parasites, which can be distinguished into neococcoids and non-neococcoids. Neococcoids are monophyletic with a peculiar reproductive system, paternal genome elimination (PGE). Different with neococcoids, Iceryini, a tribe in non-neococcoids including several damaging pests, has abdominal spiracles, compound eyes in males, relatively abundant wax, unique hermaphrodite system, and specific symbionts. However, the current studies on the gene resources and genomic mechanism of scale insects are mainly limited in the neococcoids, and lacked of comparison in an evolution frame.

RESULT

We sequenced and de novo assembled a transcriptome of Icerya aegyptiaca (Douglas), a worldwide pest of Iceryini, and used it as representative of non-neococcoids to compare with the genomes or transcriptomes of other six species from different families of neococcoids. We found that the genes under positive selection or negative selection intensification (simplified as "selected genes" below) in I. aegyptiaca included those related to neurogenesis and development, especially eye development. Some genes related to fatty acid biosynthesis were unique in its transcriptome with relatively high expression and not detected in neococcoids. These results may indicate a potential link to the unique structures and abundant wax of I. aegyptiaca compared with neococcoids. Meanwhile, genes related to DNA repair, mitosis, spindle, cytokinesis and oogenesis, were included in the selected genes in I. aegyptiaca, which is possibly associated with cell division and germ cell formation of the hermaphrodite system. Chromatin-related process were enriched from selected genes in neococcoids, along with some mitosis-related genes also detected, which may be related to their unique PGE system. Moreover, in neococcoid species, male-biased genes tend to undergo negative selection relaxation under the PGE system. We also found that the candidate horizontally transferred genes (HTGs) in the scale insects mainly derived from bacteria and fungi. bioD and bioB, the two biotin-synthesizing HTGs were exclusively found in the scale insects and neococcoids, respectively, which possibly show potential demand changes in the symbiotic relationships.

CONCLUSION

Our study reports the first I. aegyptiaca transcriptome and provides preliminary insights for the genetic change of structures, reproductive systems and symbiont relationships at an evolutionary aspect. This will provide a basis for further research and control of scale insects.

The Complete Mitochondrial Genome of the Chinese White Wax Scale Insect, Ericerus pela Chavannes (Hemiptera: Coccidae), with Novel Gene Arrangement and Truncated tRNA Genes

The Chinese white wax scale insect, Ericerus pela Chavannes (Hemiptera: Coccidae), is one of the scale insects with great economic value and has been dispersed and reared in China for over one thousand years. Its mitochondrial genome provides essential information for the molecular identification and genetic study of this species. We assembled the complete mitochondrial genome of E. pela based on PacBio sequencing and analyzed its genomic features. The genome was 17,766 bp in length with 13 protein-coding genes, 22 tRNAs, and two rRNA genes. The analysis results showed E. pela had significant gene rearrangements involving tRNAs compared with other Coccoidea species. Furthermore, E. pela's nine tRNAs were identified to have obvious truncated structures. The phylogenetic tree compiled of the species showed a long branch of the Coccoidea lineage, which indicated the high evolutionary rate in this group. Our study revealed the mitochondrial characteristics of E. pela and enriched the mitochondrial genetic information on Coccoidea species. It also determined the occurrence of gene rearrangement for the species in this superfamily.

Differential gene expression underlying the biosynthesis of Dufour's gland signals in Bombus impatiens

Pheromones regulating social behavior are one of the most explored phenomena in social insects. However, compound identity, biosynthesis and their genetic basis are known in only a handful of species. Here we examined the gene expression associated with pheromone biosynthesis of two main chemical classes: esters and terpenes, using the social bee Bombus impatiens. We conducted chemical and RNA-seq analyses of the Dufour's gland, an exocrine gland producing a plethora of pheromones regulating social behavior in hymenopteran species. The Dufour's gland contains mostly long-chained hydrocarbons, terpenes and esters that signal reproductive and social status in several bee species. In bumble bees, the Dufour's gland contains queen- and worker-specific esters, in addition to terpenes and terpene-esters only found in gynes and queens. These compounds are assumed to be synthesized de novo in the gland, however, their genetic basis is unknown. A whole transcriptome gene expression analysis of the gland in queens, gynes, queenless and queenright workers showed distinct transcriptomic profiles, with thousands of differentially expressed genes between the groups. Workers and queens express genes associated with key enzymes in the biosynthesis of wax esters, while queens and gynes preferentially express key genes in terpene biosynthesis. Overall, our data demonstrate gland-specific regulation of chemical signals associated with social behavior and identifies candidate genes and pathways regulating caste-specific chemical signals in social insects.

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.

Identification of the Key Pathways and Genes Involved in the Wax Biosynthesis of the Chinese White Wax Scale Insect (Ericerus pela Chavannes) by Integrated Weighted Gene Coexpression Network Analysis

The white wax secreted by the male insects of the Chinese white wax scale (CWWS) is a natural high-molecular-weight compound with important economic value. However, its regulatory mechanism of wax biosynthesis is still unclear. In this study, a weighted gene coexpression network analysis (WGCNA) was used to analyze transcriptome data of first- and second-instar females, early and late female adults, and first- and second-instar males. A total of 19 partitioned modules with different topological overlaps were obtained, and three modules were identified as highly significant for wax secretion (p < 0.05). A total of 30 hub genes were obtained through screening, among which elongation of very-long-chain fatty acids protein (ELOVL) and fatty acyl-CoA reductase (FAR) are important catalytic enzymes of fatty acid metabolism. Furthermore, their metabolic catalytic products are involved in the synthesis of wax biosynthesis. The results demonstrate that WGCNA can be used for insect transcriptome analysis and effectively screen out the key genes related to wax biosynthesis.

Two Cladosporium Fungi with Opposite Functions to the Chinese White Wax Scale Insect Have Different Genome Characters

Insects encounter infection of microorganisms, and they also harbor endosymbiosis to participate in nutrition providing and act as a defender against pathogens. We previously found the Chinese white wax scale insect, Ericerus pela, was infected and killed by Cladosporium sp. (pathogen). We also found it harbored Cladosporium sp. (endogensis). In this study, we cultured these two Cladosporium fungi and sequenced their genome. The results showed Cladosporium sp. (endogensis) has a larger genome size and more genes than Cladosporium sp. (pathogen). Pan-genome analysis showed Cladosporium sp. (endogensis)-specific genes enriched in pathways related to nutrition production, such as amino acid metabolism, carbohydrate metabolism, and energy metabolism. These pathways were absent in that of Cladosporium sp. (pathogen). Gene Ontology analysis showed Cladosporium sp. (pathogen)-specific genes enriched in the biosynthesis of asperfuranone, emericellamide, and fumagillin. These terms were not found in that of Cladosporium sp. (endogensis). Pathogen Host Interactions analysis found Cladosporium sp. (endogensis) had more genes related to loss of pathogenicity and reduced virulence than Cladosporium sp. (pathogen). Cytotoxicity assay indicated Cladosporium sp. (pathogen) had cytotoxicity, while Cladosporium sp. (endogensis) had no cytotoxicity. These characters reflect the adaptation of endosymbiosis to host-restricted lifestyle and the invader of the entomopathogen to the host.

Genomic Analyses of the Fungus Paraconiothyrium sp. Isolated from the Chinese White Wax Scale Insect Reveals Its Symbiotic Character

The Chinese white wax scale, Ericerus pela, is an insect native to China. It harbors a variety of microbes. The Paraconiothyrium fungus was isolated from E. pela and genome sequenced in this study. A fungal cytotoxicity assay was performed on the Aedes albopictus cell line C6/36. The assembled Paraconiothyrium sp. genome was 39.55 Mb and consisted of 14,174 genes. The coding sequences accounted for 50.75% of the entire genome. Functional pathway analyses showed that Paraconiothyrium sp. possesses complete pathways for the biosynthesis of 20 amino acids, 10 of which E. pela lacks. It also had complementary genes in the vitamin B groups synthesis pathways. Secondary metabolism prediction showed many gene clusters that produce polyketide. Additionally, a large number of genes associated with ‘reduced virulence’ in the genome were annotated with the Pathogen–Host Interaction database. A total of 651 genes encoding carbohydrate-active enzymes were predicted to be mostly involved in plant polysaccharide degradation. Pan-specific genomic analyses showed that genes unique to Paraconiothyrium sp. were enriched in the pathways related to amino acid metabolism and secondary metabolism. GO annotation analysis yielded similar results. The top COG categories were ‘carbohydrate transport and metabolism’, ‘lipid transport and metabolism’, and ‘secondary metabolite biosynthesis, transport and catabolism’. Phylogenetic analyses based on gene family and pan genes showed that Paraconiothyrium sp is clustered together with species from the Didymosphaeriaceae family. A multi-locus sequence analysis showed that it converged with the same branch as P. brasiliense and they formed one group with fungi from the Paraconiothyrium genus. To validate the in vitro toxicity of Paraconiothyrium sp., a cytotoxicity assay was performed. The results showed that medium-cultured Paraconiothyrium sp. had no harmful effect on cell viability. No toxins were secreted by the fungus during growth. Our results imply that Paraconiothyrium sp. may establish a symbiotic relationship with the host to supply complementary nutrition to E. pela.

Reference Gene Selection for Analyzing the Transcription Patterns of Two Fatty Acyl-CoA Reductase Genes From Paracoccus marginatus (Hemiptera: Pseudococcidae)

Paracoccus marginatus (Hemiptera: Pseudococcidae), known as the papaya mealybug, could cause considerable yield loss of several plants. To date, there is no molecular-based study of P. marginatus. Fatty acyl-CoA reductases (FARs) are key enzymes involved in wax synthesis. In the present study, we cloned and characterized coding sequences (CDS) of two FAR genes from P. marginatus. The results showed that PmFAR1 and PmFAR2 CDS were 1,590 and 1,497 bp in length, respectively, and sequence analysis indicated that these two genes both had the conservative motifs belonging to FAR_C superfamily. Furthermore, seven candidate reference genes were analyzed for their expression stability by using common algorithms including comparative ΔCq method, geNorm, NormFinder, BestKeeper, and RefFinder. Eventually, β-actin and GAPDH were the best reference genes in evaluating the expression of those two FAR genes. We found that PmFAR1 and PmFAR2 showed distinct expression patterns in different life stages. Moreover, the transcription of PmFAR1 and PmFAR2 in P. marginatus fed on resistant cassava cultivars was significantly lower compared with those fed on susceptible ones, indicating the potential function of FAR genes in cassava resistance to P. marginatus. The present study might help in better understanding the molecular mechanism of cassava resistance to mealybug.

Transcriptome analyses and weighted gene coexpression network analysis reveal key pathways and genes involved in the rapid cold resistance of the Chinese white wax scale insect

The Chinese white wax scale insect, Ericerus pela, is an important resource insect in China. The rapid response of E. pela to decreasing temperatures plays key roles in the population distribution. In this study, we analyzed the gene expression of E. pela treated with low temperature using transcriptome analyses and weighted gene coexpression network analysis (WGCNA). The results showed that the cold resistance of E. pela involved changes in the expression of many genes. The genes were mainly involved in alcohol formation activity, lipid metabolism, membrane and structure, and oxidoreductase activity. According to the WGCNA results, some pathways related to cold resistance were found in the genes in the modules, such as cytoskeleton proteins, cytoskeleton protein pathway, biosynthesis of unsaturated fatty acids, glycerophospholipid metabolism, ether lipid metabolism, and thermogenesis. Some of the hub genes were nonspecific lipid-transfer proteins, DnaJ homolog subfamily C member 13, paramyosin, tropomodulin, and tubulin beta chain. In particular, the hub genes of the tan module included the heat shock protein (hsp) 10, hsp 60, hsp 70, and hsp 90 genes. Thirty-five antifreeze protein (afp) genes were identified according to the annotation results. Three afp genes were further identified among the hub genes. Six of these genes were selected for heterogeneous protein expression. One of them was expressed successfully. The thermal hysteresis activity (THA) analyses showed that the THA was 1.73°C. These results showed that the cytoskeleton, lipid metabolism, thermogenesis, HSPs and AFPs may play important roles in the cold resistance of E. pela.

Beauveria bassiana ERL836 and JEF-007 with similar virulence show different gene expression when interacting with cuticles of western flower thrips, Frankniella occidentalis

Background

Insect-killing fungal species, Beauveria bassiana, is as an environment-friendly pest management tool, and many isolates are on the track of industrialization. However, some of B. bassiana isolates show similar morphology and virulence against insect pests, and so it is hard to differentiate them. Herein we used two patented isolates, ERL836 and JEF-007, and investigated their virulence against western flower thrips, Frankliniella occidentalis, and further analyzed genome structures and transcriptional responses when interacting with cuticles of thrips to see possible differences on the initial step of fungal infection.

Results

The two isolates showed no significant differences in fungal growth, conidial production, and virulence against thrips, and they were structurally similar in genome. But, in transcription level, ERL836 appeared to infect thrips easily, while JEF-007 appeared to have more difficulty. In the GO analysis of ERL836 DEGs (differentially expressed genes), the number of up-regulated genes was much larger than that of down-regulated genes, when compared to JEF-007 DEGs (more genes down-regulated). Interestingly, in the enrichment analysis using shared DEGs between two infecting isolates, plasma membrane-mediated transporter activity and fatty acid degradation pathway including cytochrome P450 were more active in infecting ERL836.

Conclusion

The two B. bassiana isolates had similar morphology and virulence as well as genome structure, but in transcription level they differently interacted with the cuticle of western flower thrips. This comparative approach using shared DEG analysis could be easily applied to characterize the difference of the two B. bassiana isolates, JEF-007 and ERL836.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07253-y.

Potential Pathways and Genes Involved in Lac Synthesis and Secretion in Kerria chinensis (Hemiptera: Kerriidae) Based on Transcriptomic Analyses

Lac is a type of natural resin secreted by lac insects and is widely used in the military and other industries because of its excellent adhesion and insulation properties. The main ingredients of lac are lactones and lactides, which are formed from hydroxy fatty acids and sesquiterpene esters. In this study, we measured lac secretion rates by the insect Kerria chinensis at different developmental stages and identified lac secretion-minimum and lac secretion-active stages of the insect. We then analyzed transcriptomes of lac secretion-minimum and lac secretion-active stages of the insect. Based on expression profiles of genes in different stages of the insect, we identified pathways and genes that are potentially involved in lac synthesis and secretion in K. chinensis. Our study lays a foundation for future studies to reveal the molecular mechanisms and pathways of lac synthesis and secretion in this beneficial insect.

Edible insects as a food source: a review

Abstract

The practice of eating insects has been documented for a long time. With obvious economic and environmental advantages, edible insects are in the spotlight recently due to their enormous potential in multiple industries. A number of research has been conducted to maximize their value, but concerns have been raised for the food security. In this review, we compared the characteristics of edible insects with other traditional protein sources. The progress of modern entomophagy was introduced. The benefits and risks of eating insects were summarized. Additionally, the shortcoming of current production were discussed. Feasible ways and strategies were proposed to promote the consumption of edible insects.

Graphical abstract

Genome sequence of the Chinese white wax scale insect Ericerus pela: the first draft genome for the Coccidae family of scale insects

BACKGROUND

The Chinese white wax scale insect, Ericerus pela, is best known for producing wax, which has been widely used in candle production, casting, Chinese medicine, and wax printing products for thousands of years. The secretion of wax, and other unusual features of scale insects, is thought to be an adaptation to their change from an ancestral ground-dwelling lifestyle to a sedentary lifestyle on the higher parts of plants. As well as helping to improve its economic value, studies of E. pela might also help to explain the adaptation of scale insects. However, no genomic data are currently available for E. pela.

FINDINGS

To assemble the E. pela genome, 303.92 Gb of data were generated using Illumina and Pacific Biosciences sequencing, producing 277.22 Gb of clean data for assembly. The assembled genome size was 0.66 Gb, with 1,979 scaffolds and a scaffold N50 of 735 kb. The guanine + cytosine content was 33.80%. A total of 12,022 protein-coding genes were predicted, with a mean coding sequence length of 1,370 bp. Twenty-six fatty acyl-CoA reductase genes and 35 acyltransferase genes were identified. Evolutionary analysis revealed that E. pela and aphids formed a sister group and split ∼241.1 million years ago. There were 214 expanded gene families and 2,219 contracted gene families in E. pela.

CONCLUSION

We present the first genome sequence from the Coccidae family. These results will help to increase our understanding of the evolution of unique features in scale insects, and provide important genetic information for further research.

Chinese White Wax Solid Lipid Nanoparticles as a Novel Nanocarrier of Curcumin for Inhibiting the Formation of Staphylococcus aureus Biofilms

Chinese white wax solid lipid nanoparticles (cwSLNs) were prepared by high shear homogenization and ultrasound methods. Using an optimized formula, spherical cwSLNs with an average particle size of 401.9 ± 21.3 nm were obtained. The cwSLNs showed high entrapment efficiency, approximately 84.6%, for loading curcumin. The curcumin loaded cwSLNs (Cur-cwSLNs) exhibited sustained drug release properties. Notably, Cur-cwSLNs had a higher drug release rate at pH 4.5 than at pH 7.4, which suggested their applicability in an acidic environment. Cur-cwSLNs were able to inhibit the growth of Staphylococcus aureus and were more effective at reducing the biofilms produced by this bacterium compared to free curcumin. This study confirmed that cwSLNs may be novel carriers for increasing the bioavailability of curcumin with the potential to inhibit the formation of S. aureus biofilms.

Sexual Dimorphism in Wax Secretion Offers Ecological Adaptability During Ericerus pela (Hemiptera: Coccidae) Evolution

The scale insect, Ericerus pela Chavannes, shows a typical sexual dimorphism. Males and females are different not only in morphology, but also in their ability to secrete wax and ecological adaptability. Here we report the morphological and structural characteristics of wax glands on E. pela females and males. The differences in wax glands and wax secretion between females and males reflect their different needs for living habitats and different ecological strategies. Sciophilous male nymphs are with five types of wax glands, and the wax glands on the dorsum secrete a layer of wax filaments plausibly for protection against direct light irradiation. On the other hand, five types of wax glands were found on the abdomen of females. Heliophilous female nymphs hardly secrete any wax, but the wax glands located along the spiracle on the abdomen may help this insect to breathe. Female adults secrete wax filaments on eggs to protect them from predators and prevent themselves from sticking to each other. In summary, males appear to secreted wax for creating a shaded niche that fits their sciophilous life style, whereas females are likely to adopt an ecological strategy with thickened epidermis for heliophilous acclimatization and overwintering.

Policosanol fabrication from insect wax and optimization by response surface methodology

Background

Insect wax is a famous biological resource for the role in economic production in China. Insect wax is a good source of policosanol, which may is a candidate supplement in foodstuff and pharmaceuticals that has important physiological activities. Therefore, this work aims to investigate a high-yield and rapid method for policosanol fabrication from insect wax.

Results

The conditions for policosanol fabrication were optimized as follows: an oil bath temperature of 112.7°C and reductant dosage of 0.97 g (used for the reduction of 10.00 g of insect wax). The yield of policosanol reached 83.20%, which was 4 times greater than that of existing methods, such as saponification. The total content of policosanol obtained under the optimal conditions reached 87%. In other words, a high yield of policosanol was obtained from insect wax (723.84 mg/g), that was 55 times higher than that generated from beeswax-brown via saponification. The concentrations of metal residues in policosanol were within the limits of the European Union regulations and EFSA stipulation. The LD50 values for oral doses of insect wax and policosanol were both > 5 g/kg.

Conclusion

Policosanol was fabricated via solvent-free reduction from insect wax using LiAlH₄ at a high yield. The fabrication conditions were optimized. Policosanol and insect wax showed high security, which made them potential candidates as supplements in foods, pharmaceuticals and cosmetics. The rapid and high-yield method has great potential for commercial manufacturing of policosanol.

Toxicological characteristics of edible insects in China: A historical review

Edible insects are ideal food sources, which contain important nutrients and health-promoting compounds. With a rapid development of industrial insect farming, insect-derived food is a novel and emerging food industry. Edible insects have been traditionally consumed in various communities, while continuously gaining relevance in today's society; however, they currently remain underutilized. Although there are a large number of literature on edible insects, these literature primarily focus on the nutritional value edible insects. The toxicity assessment data of edible insects remain incomprehensive, especially for the new national standard that is currently in effect; and many data and conclusions are not accurately specified/reported. Therefore, we performed a literature review and summarized the data on the toxicological assessment of edible insects in China. The review first describes the research progress on safety toxicological assessment, and then offers references regarding the development of 34 edible insect species in China. These data can be a platform for the development of future toxicological assessment strategies, which can be carried out by a multidisciplinary team, possibly consisting of food engineers, agronomists, farmers, and so on, to improve the acceptability of edible insects.

Characterization and functional assay of a fatty acyl-CoA reductase gene in the scale insect, Ericerus pela Chavannes (Hemiptera: Coccoidae)

Ericerus pela Chavannes (Hemiptera: Coccoidae) is an economically important scale insect because the second instar males secrete a harvestable wax-like substance. In this study, we report the molecular cloning of a fatty acyl-CoA reductase gene (EpFAR) of E. pela. We predicted a 520-aa protein with the FAR family features from the deduced amino acid sequence. The EpFAR mRNA was expressed in five tested tissues, testis, alimentary canal, fat body, Malpighian tubules, and mostly in cuticle. The EpFAR protein was localized by immunofluorescence only in the wax glands and testis. EpFAR expression in High Five insect cells documented the recombinant EpFAR reduced 26-0:(S) CoA and to its corresponding alcohol. The data illuminate the molecular mechanism for fatty alcohol biosynthesis in a beneficial insect, E. pela.

A Lethal Fungus Infects the Chinese White Wax Scale Insect and Causes Dramatic Changes in the Host Microbiota

The Chinese white wax scale insect (Ericerus pela) is an economically valuable species with an important role in wax production. Recently, in a greenhouse in Kunming, we identified a genus of fungus that infects and kills E. pela females. This study sought to perform the molecular detection of entomopathogens and analyze the changes in the host microbiota after entomopathogen infection. We used library construction, high-throughput sequencing and real-time quantitative polymerase chain reaction (RT-qPCR) to identify the fungi infecting adult E. pela, to understand the changes in the host organism, and to determine the distribution of the entomopathogens. Cladosporium langeronii and C. sphaerospermum were the main pathogenic species that infected the E. pela females, and they were most prevalent in the dorsal cuticle. In vivo, after infection, the proportion of Cladosporium clearly increased. The infection had little influence on the fungal community but had a strong influence on the bacterial community. After infection, Arsenophonus was dominant, and numerous bacterial genera disappeared. However, Rickettsia, instead of Arsenophonus, became dominant in the Cladosporium-infected individuals that had also been infected with Rickettsia. We identified the species that infected E. pela females and determined the influence of infection on the host microorganisms.

Transcriptome analysis in different developmental stages of Batocera horsfieldi (Coleoptera: Cerambycidae) and comparison of candidate olfactory genes

The white-striped longhorn beetle Batocera horsfieldi (Coleoptera: Cerambycidae) is a polyphagous wood-boring pest that causes substantial damage to the lumber industry. Moreover olfactory proteins are crucial components to function in related processes, but the B. horsfieldi genome is not readily available for olfactory proteins analysis. In the present study, developmental transcriptomes of larvae from the first instar to the prepupal stage, pupae, and adults (females and males) from emergence to mating were built by RNA sequencing to establish a genetic background that may help understand olfactory genes. Approximately 199 million clean reads were obtained and assembled into 171,664 transcripts, which were classified into 23,380, 26,511, 22,393, 30,270, and 87, 732 unigenes for larvae, pupae, females, males, and combined datasets, respectively. The unigenes were annotated against NCBI’s non-redundant nucleotide and protein sequences, Swiss-Prot, Gene Ontology (GO), Pfam, Clusters of Eukaryotic Orthologous Groups (KOG), and KEGG Orthology (KO) databases. A total of 43,197 unigenes were annotated into 55 sub-categories under the three main GO categories; 25,237 unigenes were classified into 26 functional KOG categories, and 25,814 unigenes were classified into five functional KEGG Pathway categories. RSEM software identified 2,983, 3,097, 870, 2,437, 5,161, and 2,882 genes that were differentially expressed between larvae and males, larvae and pupae, larvae and females, males and females, males and pupae, and females and pupae, respectively. Among them, genes encoding seven candidate odorant binding proteins (OBPs) and three chemosensory proteins (CSPs) were identified. RT-PCR and RT-qPCR analyses showed that BhorOBP3, BhorCSP2, and BhorOBPC1/C3/C4 were highly expressed in the antenna of males, indicating these genes may may play key roles in foraging and host-orientation in B. horsfieldi. Our results provide valuable molecular information about the olfactory system in B. horsfieldi and will help guide future functional studies on olfactory genes.

In vivo evaluation of insect wax for hair growth potential

Insect wax is secreted by Ericerus pela Chavanness. It has been traditionally used to treat hair loss in China, but few reports have been published on the hair growth-promoting effect of insect wax. In this work, we examined the hair growth-promoting effects of insect wax on model animals. Different concentrations of insect wax were topically applied to the denuded backs of mice, and 5% minoxidil was applied topically as a positive control. We found that insect wax significantly promoted hair growth in a dose-dependent manner, 45% and 30% insect wax both induced hair to regrow, while less visible hair growth was observed in blank controls on the 16^th day. The experimental areas treated with 45% and 30% insect wax exhibited significant differences in hair scores compared to blank controls, and hair lengths in the 45% and 30% insect wax group was significantly longer than in blank controls on the 16^th and 20^th days. There were no new hair follicles forming in the treated areas, and the hair follicles were prematurely converted to the anagen phase from the telogen phase in experimental areas treated with 45% and 30% insect wax. Both 45% and 30% insect wax upregulated vascular endothelial growth factor expression. The results indicated that 45% and 30% insect wax showed hair growth-promoting potential approximately as potent as 5% minoxidil by inducing the premature conversion of telogen-to-anagen and by prolonging the mature anagen phase rather than increasing the number of hair follicles, which was likely related to the upregulation of VEGF expression. The dissociative policosanol in insect wax was considered the key ingredient most likely responsible for the hair growth promoting potential.

Transcriptome analysis in the beet webworm, Spoladea recurvalis (Lepidoptera: Crambidae)

The beet webworm, Spoladea recurvalis Fabricius, is a destructive pest on vegetable crops in tropics and subtropics; its main host plant is amaranth. It has become imperative to develop non-chemical methods to control S. recurvalis on amaranth. However, the lack of molecular information about this species has hindered the development of novel pest management strategies. In this study, high-throughput RNA sequencing covering de novo sequence assemblies, functional annotation of transcripts, gene function classification and enrichment was performed on S. recurvalis. Illumina sequencing generated a total of 120 435 transcript contigs ranging from 201 to 22 729 bases with a mean length of 688 bases. The assembled transcripts were subjected to Basic Local Alignment Search Tool-X (BLASTX) to obtain the annotations against non-redundant, Swiss-Prot, Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) protein databases. A subset of 58 225 transcript sequences returned hits from known proteins in the National Center for Biotechnology Information database, and the majority of the transcript sequences had the highest number of hits for Danaus plexippus (50.43%). A total of 1217 Gene Ontology-level 3 annotations were assigned to 51 805 transcripts, while 39 650 transcripts were predicted as functional protein-coding genes in the COG database and 20 037 transcripts were enriched to KEGG pathways. We identified 40 putative genes related to pheromone production and reception in S. recurvalis, with the expression of one gene between 0.29 and 1141.79 fragments per kilo base per million (FPKM) reads. The transcriptome sequence of S. recurvalis is a first step toward offering a comprehensive genomic resource which would enable better understanding of molecular mechanisms to enable development of effective pest management practices for this species.

Polysaccharide extracted from Chinese white wax scale ameliorates 2,4-dinitrochlorobenzene-induced atopic dermatitis-like symptoms in BALB/c mice

Atopic dermatitis (AD) is a common inflammatory skin disease with high rates of morbidity and is associated with erythema, pruritus, scaling of affected areas of skin. It is extremely important to introduce a therapeutic agent which has significant anti-inflammatory effect with less side-effect for treatment of AD. This study evaluated the effect of a natural compound from herbal extracts, the crude polysaccharide extracted from the white wax scale (CWPS), on AD-like mice. Repeated applications of 2,4-dinitrochlorobenzene (DNCB) were performed on ear and dorsal skin of BALB/c mice to induce AD-like symptoms and skin lesions. Oral administration of CWPS decreased serum IgE level and limited the infiltration of mast cells and eosinophils to the dermal tissues in the DNCB-induced AD mice. In addition, CWPS reduced Th1 and Th17 responses, leading to an attenuated cutaneous inflammatory response. Furthermore, in vitro study also demonstrated that CWPS limited T cell activation and cytokines (i.e. IFN-γ and IL-17) production induced by DNCB. We conclude that CWPS attenuates DNCB-induced AD-like skin lesion through modulating T cell-elicited immune responses and CD4+ T cell polarization, and could be exploited as a new therapeutic approach for AD.

Identification and evaluation of reference genes in the Chinese white wax scale insect Ericerus pela

BACKGROUND

The Chinese white wax scale insect, Ericerus pela, is a well-known resource insect. The females and males are dramatically distinct at each developmental stage. We sought to identify suitable reference genes to use as internal controls in molecular research on E. plea.

RESULTS

geNorm, RefFinder and Normfinder analyses showed that ßTub-2 was the best reference gene throughout different developmental stages; SdhA-1 was the most stable reference gene in different tissues, and ßTub-1 was the most reliable reference gene under treatment with different temperatures. The results also showed that the optimal number of reference genes for analyzing target gene expression levels in the three experimental conditions was two.

CONCLUSIONS

The identified reference genes are suitable reference genes for normalization in RT-qPCR of E. pela samples.

Molecular Characterization of Two Fatty Acyl-CoA Reductase Genes From Phenacoccus solenopsis (Hemiptera: Pseudococcidae)

Fatty acyl-CoA reductases (FARs) are key enzymes involved in fatty alcohol synthesis. Here, we cloned and characterized full-length cDNAs of two FAR genes from the cotton mealybug, Phenacoccus solenopsis. The results showed PsFAR I and PsFAR II cDNAs were 1,584 bp and 1,515 bp in length respectively. Both PsFAR I and PsFAR II were predicted to be located in the endoplasmic reticulum by Euk-mPLoc 2.0 approach. Both of them had a Rossmann folding region and a FAR_C region. Two conservative motifs were discovered in Rossmann folding region by sequence alignment including a NADPH combining motif, TGXXGG, and an active site motif, YXXXK. A phylogenetic tree made using MEGA 6.06 indicated that PsFAR I and PsFAR II were placed in two different branches. Gene expression analysis performed at different developmental stages showed that the expression of PsFar I is significantly higher than that of PsFar II in first and second instar nymphs and in male adults. Spirotetramat treatment at 125 mg/liter significantly increased the expression of PsFar I in third instar nymphs, but there was no effect in the expression of PsFar II Our results indicated these two FAR genes showed different expression patterns during insect development and after pesticide treatment, suggesting they play different roles in insect development and detoxification against pesticides.

Transcriptomic and proteomic analyses on the supercooling ability and mining of antifreeze proteins of the Chinese white wax scale insect

The Chinese white wax scale insect, Ericerus pela, can survive at extremely low temperatures, and some overwintering individuals exhibit supercooling at temperatures below -30°C. To investigate the deep supercooling ability of E. pela, transcriptomic and proteomic analyses were performed to delineate the major gene and protein families responsible for the deep supercooling ability of overwintering females. Gene Ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that genes involved in the mitogen-activated protein kinase, calcium, and PI3K-Akt signaling pathways and pathways associated with the biosynthesis of soluble sugars, sugar alcohols and free amino acids were dominant. Proteins responsible for low-temperature stress, such as cold acclimation proteins, glycerol biosynthesis-related enzymes and heat shock proteins (HSPs) were identified. However, no antifreeze proteins (AFPs) were identified through sequence similarity search methods. A random forest approach identified 388 putative AFPs in the proteome. The AFP gene ep-afp was expressed in Escherichia coli, and the expressed protein exhibited a thermal hysteresis activity of 0.97°C, suggesting its potential role in the deep supercooling ability of E. pela.

A de novo transcriptomic analysis to reveal functional genes in Apolygus lucorum

The widespread planting of genetically engineered cotton producing the Cry1Ac toxin has led to significantly reduced pesticide applications since 1997. However, consequently, the number of green mirid bugs (GMB), Apolygus lucorum (Meyer-Dür) has increased. So far the GMB, instead of the cotton bollworm Helicoverpa armigera (Hübner), has become the major pest in the transgenic Bt cotton field and has influenced cotton yield. Disproportionately, only a few studies on GMB at molecular level have been reported. Libraries from both third instar nymphs and adults were sequenced using Illumina technology, producing more than 106 million short reads and assembled into 63 029 unigenes of mean length 597 nt and N50 813 nt, ranging from 300 nt to 9771 nt. BLASTx analysis against Nr, Swissprot, GO and COG was performed to annotate these unigenes. As a result, 26 478 unigenes (42.01%) matched to known proteins and 107 immune-related, 320 digestive-related and 53 metamorphosis-related genes were detected in these annotated unigenes. Additionally, we profiled gene expression using mapping based differentially expressed genes (DEGs) strategy between the two developmental stages: nymph and adult. The results demonstrated that thousands of genes were significantly differentially expressed at different developmental stages. The transcriptome and gene expression data provided comprehensive and global gene resources of GMB. This transcriptome would improve our understanding of the molecular mechanisms of various underlying biological characteristics, including development, digestion and immunity in GMB. Therefore, these findings could help elucidate the intrinsic factors of the GMB resurgence, offering novel pest management targets for future transgenic cotton breeding.

Transcriptome analysis of sexually dimorphic Chinese white wax scale insects reveals key differences in developmental programs and transcription factor expression

The Chinese white wax scale insect, Ericerus pela, represents one of the most dramatic examples of sexual dimorphism in any insect species. In this study, we showed that although E. pela males display complete metamorphosis similar to holometabolous insects, the species forms the sister group to Acyrthosiphon pisum and cluster with hemimetabolous insects. The gene expression profile and Gene Ontology (GO) analyses revealed that the two sexes engaged in distinct developmental programs. In particular, female development appeared to prioritize the expression of genes related to cellular, metabolic, and developmental processes and to anatomical structure formation in nymphs. By contrast, male nymphal development is characterized by the significant down-regulation of genes involved in chitin, the respiratory system, and neurons. The wing and appendage morphogenesis, anatomical and tissue structure morphogenesis programs activated after male nymphal development. Transcription factors (that convey juvenile hormone or ecdysone signals, and Hox genes) and DNA methyltransferase were also differentially expressed between females and males. These results may indicate the roles that these differentially expressed genes play in regulating sexual dimorphism through orchestrating complex genetic programs. This differential expression was particularly prominent for processes linked to female development and wing development in males.

Protein profiles of Chinese white wax scale, Ericerus pela, at the male pupal stage by high-throughput proteomics

The Chinese white wax scale insect (Ericerus pela) is sexually dimorphic with holometabolous males and hemimetabolous females. Holometabolous insects were assumed to originate from hemimetabolous ancestors. Therefore, the male pupal stage is a major innovation compared with hemimetabolous female insects. Here, the protein profiles of the male pupae were obtained by high-throughput proteomics and analyzed using bioinformatics methods. A total of 1,437 peptides were identified and assigned to 677 protein groups. Most of the proteins had molecular weights below 40 kDa and isoelectric points from 4 to 7. Gene Ontology terms were assigned to 331 proteins, including metabolic process, developmental process, and cellular process. Kyoto Encyclopedia of Genes and Genomes annotations identified 142 pathways and most proteins were assigned to metabolism events. Pathways involved in cell growth and death, signal transduction, folding, and sorting and degradation were also identified. Six proteins that had undergone positive selection were classified into four groups, protein biosynthesis, protein degeneration, signal transduction, and detoxification. Many of the high-abundance proteins were enzymes involved in carbohydrate, lipid, and amino acid metabolism; signal transduction; degradation; and immunization, which indicated that metabolism, disruption, and development occurred intensely at the pupal stage. These processes are closely related to the physiological status of pupae. The results also suggested that these related proteins may be fundamental factors in the formation of pupae. This study describes pupal characterization at the molecular level and provides a basis for further physiological studies.

Cloning and expression analysis of four heat shock protein genes in Ericerus pela (Homoptera: Coccidae)

To explore the function of small heat shock protein genes (shsps) and hsp70 in Ericerus pela, we cloned the full-length cDNA sequences of hsp21.5, hsp21.7, hsp70, and hsc70 and the genomic sequence of hsc70. Open reading frames of the four hsps were 570, 564, 1,908, and 1,962 base pairs (bp), respectively, which encode proteins with calculated molecular mass of 21.5, 21.7, 69.8, and 71.6 kDa. Amino acid sequence analysis revealed the presence of the conserved Hsp motifs in all four proteins. The genomic DNA of hsc70 had four introns. ep-hsp21.5 was orthologous and ep-hsp21.7 was species specific. Expression of all four transcripts during heat or cold stress and development was examined by quantitative real-time polymerase chain reaction. All four hsps were upregulated during heat or cold stress in female adults, indicating a correlation between the four hsps and heat or cold-stress tolerance in female adults. ep-hsp21.7 and ep-hsp70 were upregulated during heat stress in male larvae, implying a correlation between the two hsps and heat-stress tolerance in male larvae. The four ep-hsps were also upregulated during the developmental process in males, and ep-hsp21.5, ep-hsp70, and ep-hsc70 were upregulated in females, which indicates their possible role in the developmental regulation of E. pela.

Transcriptome characteristics and six alternative expressed genes positively correlated with the phase transition of annual cambial activities in Chinese Fir (Cunninghamia lanceolata (Lamb.) Hook)

Background

The molecular mechanisms that govern cambial activity in angiosperms are well established, but little is known about these molecular mechanisms in gymnosperms. Chinese fir (Cunninghamia lanceolata (Lamb.) Hook), a diploid (2n  = 2x  = 22) gymnosperm, is one of the most important industrial and commercial timber species in China. Here, we performed transcriptome sequencing to identify the repertoire of genes expressed in cambium tissue of Chinese fir.

Methodology/Principal Findings

Based on previous studies, the four stage-specific cambial tissues of Chinese fir were defined using transmission electron microscopy (TEM). In total, 20 million sequencing reads (3.6 Gb) were obtained using Illumina sequencing from Chinese fir cambium tissue collected at active growth stage, with a mean length of 131 bp and a N50 of 90 bp. SOAPdenovo software was used to assemble 62,895 unigenes. These unigenes were further functionally annotated by comparing their sequences to public protein databases. Expression analysis revealed that the altered expression of six homologous genes (ClWOX1, ClWOX4, ClCLV1-like, ClCLV-like, ClCLE12, and ClPIN1-like) correlated positively with changes in cambial activities; moreover, these six genes might be directly involved in cambial function in Chinese fir. Further, the full-length cDNAs and DNAs for ClWOX1 and ClWOX4 were cloned and analyzed.

Conclusions

In this study, a large number of tissue/stage-specific unigene sequences were generated from the active growth stage of Chinese fir cambium. Transcriptome sequencing of Chinese fir not only provides extensive genetic resources for understanding the molecular mechanisms underlying cambial activities in Chinese fir, but also is expected to be an important foundation for future genetic studies of Chinese fir. This study indicates that ClWOX1 and ClWOX4 could be possible reverse genetic target genes for revealing the molecular mechanisms of cambial activities in Chinese fir.

Sequencing, de novo assembly, annotation and SSR and SNP detection of sabaigrass (Eulaliopsis binata) transcriptome

Eulaliopsis binata is one of the best fiber grass plants for its high fiber quality and production. Large scale trancriptome sequencing of E. binata was first performed using mixed leaf samples of 20 wild clusters. A total of 26,438,832 clean reads were generated and were assembled into 59,134 isogenes with an average length of 845bp. BLAST against the NCBI non-redundant protein, KEGG and GO databases has classified these isogenes into functional categories for understanding gene functions and regulation pathways. Only 15.0% of the assembled isogenes were similar to known proteins and 24.4% has no hits in the nr protein database. The total isogenes and 5306 highly expressed isogenes were performed by BLASTx with the MAIZEWALL, the cell wall navigator and the PlantTFDB databases. A total of 6681 simple sequence repeats (SSRs) and 147,177 single nucleotide polymorphisms (SNPs) were detected in the isogenes and 5723 pairs of SSR primers were designed.