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

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.

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 Expression Analyses by qRT-PCR in Dendroctonus valens

Dendroctonus valens is the main pest of the genus Pinus. To facilitate gene expression analyses, suitable reference genes for adults and mature larvae of D. valens under different temperature conditions were determined. In particular, we obtained the sequences of candidate reference genes, ACT, TUB, SHDA, PRS18, 18S rRNA, and CYP4G55, from transcriptome data. Real-time quantitative PCR was used to analyze gene expression, and geNorm, NormFinder, and BestKeeper were used to evaluate expression stability. Under different temperature conditions, the expression levels of 18S rRNA, PRS18, and TUB were stable in adults, in which 18S rRNA > PRS18 > TUB. In mature larvae, the expression levels of TUB, 18S rRNA, and SDHA were stable, in which TUB > 18S rRNA > SDHA. The combination of 18S rRNA and PRS18 is recommended for studies of gene expression in adults and the combination of 18S rRNA and TUB is effective for studies of gene expression in mature larvae of D. valens under different temperature conditions.

Selection of Reference Genes for RT-qPCR Analysis Under Extrinsic Conditions in the Hawthorn Spider Mite, Amphitetranychus viennensis

Hawthorn spider mite, Amphitetranychus viennensis Zacher, is an economically important arthropod pest for fruit trees and woody ornamental plants. Extensive and repetitive use of synthetic acaricides has led to the development of resistance in A. viennensis. To understand the molecular basis of pesticide resistance, and to develop genetic-based control alternatives (e.g., RNAi-based biopesticides), a standardized protocol for real-time quantitative reverse transcription PCR (RT-qPCR) is needed. In the proceeding phase of this research, we screened for the internal references for RT-qPCR analysis from a pool of A. viennensis housekeeping genes under the intrinsic conditions, including developmental stage, sex, and diapause. Here, we continued our efforts to search for the reference genes under an array of extrinsic conditions, including temperature, humidity, photoperiod, host plant, and dietary RNAi. The stability of these candidate reference genes was investigated using geNorm, NormFinder, BestKeeper, and ΔCt method, respectively. Finally, RefFinder, a statistical platform integrating all four algorisms, provided a comprehensive list of genes for each extrinsic condition: (1) EF1A, α-tubulin and Actin3 were the best candidates for temperature, (2) GAPDH, 18S, and Actin3 were the most stable genes for humidity, (3) V-ATPase B, Actin3, and 18S were the top reference genes for photoperiod, (4) GAPDH, V-ATPase B, and α-tubulin were recommended for host plants, and (5) GAPDH, V-ATPase B, and RPS9 were the top choices for dietary RNAi. Overall, V-ATPase B, GAPDH, and Actin3 were the most commonly selected reference genes in A. viennensis regardless of the experimental conditions, including both intrinsic and extrinsic. Information present here lays the foundation for the genomic and functional genomic research in A. viennensis.

Identification and validation of potential reference gene for effective dsRNA knockdown analysis in Chilo partellus

Chilo partellus is an invasive polyphagous pest that has not been effectively managed with chemical pesticides. To select potential dsRNAs for use in an alternate control strategy, it is crucial to identify and evaluate stable reference genes for knockdown expression studies. This study evaluates the expression stability of seven candidate reference genes in C. partellus larvae fed on crude bacterially-expressed dsRNAs and purified dsRNAs at different time intervals, as well as the developmental stages and sexes. The expression stabilities of the reference genes were evaluated with different software programmes, such as BestKeeper, NormFinder, deltaCt, geNorm, and RefFinder. The overall results rank ELF as the most stably expressed reference gene when larvae were fed with crude bacteria-induced dsRNAs and purified dsRNA. However, Tubulin and HSP70 were more stable under different developmental stages and sexes. The expression levels of larvae that were fed crude bacteria-induced dsRNAs of Chitinase and Acetylcholinesterase were normalized with the four most stable reference genes (ELF, HSP70, V-ATPase and Tubulin) and the least stable reference gene (18S and HSP70) based on the geNorm algorithm. The least stable reference gene showed inconsistent knockdown expression, thereby confirming that the validation of a suitable reference gene is crucial to improve assay accuracy for dsRNA-targeted gene selection in C. partellus.

Evaluation and Validation of Reference Genes for Quantitative Real-Time PCR in Helopeltis theivora Waterhouse (Hemiptera: Miridae)

Helopeltis theivora Waterhouse is a predominant sucking pest in many tropic economic crops, such as tea, cocoa and coffee. Quantitative real-time PCR (qRT-PCR) is one of the most powerful tools to analyze the gene expression level and investigate the mechanism of insect physiology at transcriptional level. Gene expression studies utilizing qRT-PCR have been applied to numerous insects so far. However, no universal reference genes could be used for H. theivora. To obtain accurate and reliable normalized data in H. theivora, twelve candidate reference genes were examined under different tissues, developmental stages and sexes by using geNorm, NormFinder, BestKeeper, Delta Ct and RefFinder algorithms, respectively. The results revealed that the ideal reference genes differed across the treatments, and the consensus rankings generated from stability values provided by these programs suggested a combination of two genes for normalization. To be specific, RPS3A and Actin were the best suitable reference genes for tissues, RPL13A and GAPDH were suitable for developmental stages, EF1α and RPL13A were suitable for sexes, and RPL13A and RPS3A were suitable for all samples. This study represents the first systematic analysis of reference genes for qRT-PCR experiments in H. theivora, and the results can provide a credible normalization for qRT-PCR data, facilitating transcript profiling studies of functional genes in this insect.

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.

Evaluation of Reference Genes for Normalization of RT-qPCR Gene Expression Data for Trichoplusia ni Cells During Antheraea pernyi (Lepidoptera: Saturniidae) Multicapsid Nucleopolyhedrovirus (AnpeNPV) Infection

Baculovirus infection impacts global gene expression in the host cell, including the expression of housekeeping genes. Evaluation of candidate reference genes during a viral infection will inform the selection of appropriate reference gene(s) for the normalization of expression data generated by Reverse Transcription Quantitative Real-timePolymerase Chain Reaction (RT-qPCR). Antheraea pernyi multicapsid nucleopolyhedrovirus (AnpeNPV) is able to infect the High Five cells (Tn-Hi5). In the present study, 10 candidate reference genes were evaluated in AnpeNPV-infected Tn-Hi5 cells. Gene expression data were analyzed using geNorm, NormFinder, BestKeeper, and RefFinder. The candidate genes were further validated as reliable reference genes for RT-qPCR by analyzing the expression of three target genes. The results of data analysis using four statistical methods showed that RPS18 was the least stable among the candidate reference genes tested. 18S rRNA and 28S rRNA were not suitable as reference genes for RT-qPCR analysis in AnpeNPV-infected Tn-Hi5 cells. Comprehensive ranking of the 10 candidate reference genes by RefFinder analysis indicated that Ann B, c45128_g1, and ACT were the top three genes. Normalization of the expression of three target genes using the candidate reference genes indicated the same expression pattern when Ann B and c45128_g1 were used as reference genes, with slight differences in the relative expression at each infection time point. Overall, Ann B and c45128_g1 were recommended to be more suitable than the most commonly used reference genes, such as ACT, GAPDH, and TUB, for RT-qPCR data normalization in AnpeNPV-infected Tn-Hi5 cells up to 48 hpi.

Selection of Reference Genes for the Normalization of RT-qPCR Data in Gene Expression Studies in Insects: A Systematic Review

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for quantifying expression levels of targeted genes during various biological processes in numerous areas of clinical and biological research. Selection of appropriate reference genes for RT-qPCR normalization is an elementary prerequisite for reliable measurements of gene expression levels. Here, by analyzing datasets published between 2008 and 2017, we summarized the current trends in reference gene selection for insect gene expression studies that employed the most widely used SYBR Green method for RT-qPCR normalization. We curated 90 representative papers, mainly published in 2013–2017, in which a total of 78 insect species were investigated in 100 experiments. Furthermore, top five journals, top 10 frequently used reference genes, and top 10 experimental factors have been determined. The relationships between the numbers of the reference genes, experimental factors, analysis tools on the one hand and publication date (year) on the other hand was investigated by linear regression. We found that the more recently the paper was published, the more experimental factors it tended to explore, and more analysis tools it used. However, linear regression analysis did not reveal a significant correlation between the number of reference genes and the study publication date. Taken together, this meta-analysis will be of great help to researchers that plan gene expression studies in insects, especially the non-model ones, as it provides a summary of appropriate reference genes for expression studies, considers the optimal number of reference genes, and reviews the average number of experimental factors and analysis tools per study.

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.

Evaluation of potential reference genes for real-time qPCR analysis in a biparental beetle, Lethrus apterus (Coleoptera: Geotrupidae)

Hormones play an important role in the regulation of physiological, developmental and behavioural processes. Many of these mechanisms in insects, however, are still not well understood. One way to investigate hormonal regulation is to analyse gene expression patterns of hormones and their receptors by real-time quantitative polymerase chain reaction (RT-qPCR). This method, however, requires stably expressed reference genes for normalisation. In the present study, we evaluated 11 candidate housekeeping genes as reference genes in samples of Lethrus apterus, an earth-boring beetle with biparental care, collected from a natural population. For identifying the most stable genes we used the following computational methods: geNorm, NormFinder, BestKeeper, comparative delta Ct method and RefFinder. Based on our results, the two body regions sampled (head and thorax) differ in which genes are most stably expressed. We identified two candidate reference genes for each region investigated: ribosomal protein L7A and RP18 in samples extracted from the head, and ribosomal protein L7A and RP4 extracted from the muscles of the thorax. Additionally, L7A and RP18 appear to be the best reference genes for normalisation in all samples irrespective of body region. These reference genes can be used to study the hormonal regulation of reproduction and parental care in Lethrus apterus in the future.

Selection and validation of reference genes for qRT-PCR analysis during biological invasions: The thermal adaptability of Bemisia tabaci MED

The Bemisia tabaci Mediterranean (MED) cryptic species has been rapidly invading to most parts of the world owing to its strong ecological adaptability, which is considered as a model insect for stress tolerance studies under rapidly changing environments. Selection of a suitable reference gene for quantitative stress-responsive gene expression analysis based on qRT-PCR is critical for elaborating the molecular mechanisms of thermotolerance. To obtain accurate and reliable normalization data in MED, eight candidate reference genes (β-act, GAPDH, β-tub, EF1-α, GST, 18S, RPL13A and α-tub) were examined under various thermal stresses for varied time periods by using geNorm, NormFinder and BestKeeper algorithms, respectively. Our results revealed that β-tub and EF1-α were the best reference genes across all sample sets. On the other hand, 18S and GADPH showed the least stability for all the samples studied. β-act was proved to be highly stable only in case of short-term thermal stresses. To our knowledge this was the first comprehensive report on validation of reference genes under varying temperature stresses in MED. The study could expedite particular discovery of thermotolerance genes in MED. Further, the present results can form the basis of further research on suitable reference genes in this invasive insect and will facilitate transcript profiling in other invasive insects.

Validation of reference genes for normalization of gene expression by qRT-PCR in a resveratrol-producing entophytic fungus (Alternaria sp. MG1)

Alternaria sp. MG1, an endophytic fungus isolated from Vitis vinifera, can independently produce resveratrol, indicating that this species contains the key genes for resveratrol biosynthesis. Identification of these key genes is essential to understand the resveratrol biosynthesis pathway in this strain, which is currently unknown in microorganisms. qRT-PCR is an efficient and widely used method to identify the key genes related to unknown pathways at the level of gene expression. Verification of stable reference genes in this strain is essential for qRT-PCR data normalization, although results have been reported for other Alternaria sp. strains. In this study, nine candidate reference genes including TUBA, EF1, EF2, UBC, UFD, RPS5, RPS24, ACTB and 18S were evaluated for expression stability in a diverse set of six samples representing different growth periods. We compared cell culture conditions and an optimized condition for resveratrol production. The comparison of the results was performed using four statistical softwares. A combination of TUBA and EF1 was found to be suitable for normalization of Alternaria sp. MG1 in different developmental stages, and 18S was found to be the least stable. The reference genes verified in this study will facilitate further research to explore gene expression and molecular mechanisms as well as the improvement of secondary metabolite yields in Alternaria sp. MG1. To our knowledge, this is the first validation of reference genes in Alternaria with the capability to produce resveratrol. Additionally, these results provide useful guidelines for the selection of reference genes in other Alternaria species.