De Novo Transcriptome Analysis Reveals Abundant Gonad-specific Genes in the Ovary and Testis of Henosepilachna vigintioctopunctata

Comparative Transcriptome Analysis of Henosepilachna vigintioctomaculata Reveals Critical Pathways during Development

Henosepilachna vigintioctomaculata is distributed in several Asian countries. The larvae and adults often cause substantial economic losses to Solanaceae crops such as potato, tomato, eggplant, and Chinese boxthorn. Even though a chromosome-level genome has been documented, the expression profiles of genes involved in development are not determined. In this study, we constructed embryonic, larval, pupal, and adult transcriptomes, generated a comprehensive RNA-sequencing dataset including ~52 Gb of clean data, and identified 602,773,686 cleaned reads and 33,269 unigenes. A total of 18,192 unigenes were successfully annotated against NCBI nonredundant protein sequences, Swissprot, Eukaryotic Orthologous Groups, Gene Ontology (GO), or Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. There were 3580, 2040, 5160, 2496, 3008, and 3895 differentially expressed genes (DEGs) between adult/egg, egg/larval, larval/pupal, adult/pupal, egg/pupal, and adult/larval samples, respectively. GO and KEGG analyses of the DEGs highlighted several critical pathways associated with specific developing stages. This is the first comprehensive transcriptomic dataset encompassing all developmental stages in H. vigintioctomaculata. Our data may facilitate the exploitation of gene targets for pest control and can serve as a valuable gene resource for future molecular investigations.

Gonadal Transcriptomic Analysis Reveals Novel Sex-Related Genes in Bactrocera dorsalis

Bactrocera dorsalis (Hendel) (Diptera: Tephritidae) is one of the most devastating agricultural pests worldwide due to its high reproductive and invasive abilities. The elucidation of its gonadal developmental characteristics and the identification of sex-related genes will provide a useful genetic basis for reproductive-based pest control. Here, the gonadal transcriptome of B. dorsalis was sequenced, and novel gonad-specific expressed genes were analyzed. A total of 1338, 336, 35, and 479 differentially expressed genes (DEGs) were found in the testis (TE), ovary (OV), female accessory gland (FAG), and male accessory gland (MAG), respectively. Furthermore, 463 highly expressed gonad-specific genes were identified, with the TE having the highest number of specific highly expressed genes, at 402, followed by 51 in the OV, 9 in the MAG, and only 1 in the FAG. Strikingly, approximately half of highly expressed gonad-specific genes were uncharacterized. Then, it was found that 35, 17, 3, 2, and 1 of 202 uncharacterized highly expressed TE-specific genes encoded proteins that contained transmembrane domains, signal peptides, high-mobility group boxes, the zinc finger domain, and the BTB/POZ domain, respectively. Interestingly, approximately 40% of uncharacterized highly expressed gonad-specific genes encoding proteins were not predicted to possess functional motifs or domains. Finally, the spatiotemporal expression and sequence characterization of six novel highly expressed gonad-specific genes were analyzed. Altogether, our findings provide a valuable dataset for future functional analyses of sex-related genes and potential target sites for pest control.

SMRT Sequencing Technology Was Used to Construct the Batocera horsfieldi (Hope) Transcriptome and Reveal Its Features

Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi.

Silencing of dre4 Contributes to Mortality of Phyllotreta striolata

The striped flea beetle, Phyllotreta striolata, is one of the most destructive pests of Cruciferae crops worldwide. RNA interference (RNAi) is a promising alternative strategy for pest biological control, which overcomes the weakness of synthetic insecticides, such as pest resistance, food safety problems and toxicity to non-target insects. The homolog of Spt16/FACT, dre4 plays a critical role in the process of gene transcription, DNA repair, and DNA replication; however, the effects of dre4 silencing in P. striolata remain elusive. In this study, we cloned and characterized the full-length dre4 from P. striolata and silenced Psdre4 through microinjection and oral delivery; it was found that the silencing of dre4 contributed to the high mortality of P. striolata in both bioassays. Moreover, 1166 differentially regulated genes were identified after Psdre4 interference by RNA-seq analysis, which might have been responsible for the lethality. The GO analysis indicated that the differentially regulated genes were classified into three GO functional categories, including biological process, cellular component, and molecular function. The KEGG analysis revealed that these differentially regulated genes are related to apoptosis, autophagy, steroid hormone biosynthesis, cytochrome P450 and other signaling pathways. Our results suggest that Psdre4 is a fatal RNAi target and has significant potential for the development of RNA pesticides for P. striolata management.

Silencing uridine diphosphate N-acetylglucosamine pyrophosphorylase gene impairs larval development in Henosepilachna vigintioctopunctata

BACKGROUND

Uridine diphosphate-N-acetylglucosamine (UDP-GlcNAc) diphosphorylase (UAP) catalyzes the formation of UDP-GlcNAc, the precursor for the production of chitin in ectodermally derived epidermal cells and midgut, for GlcNAcylation of proteins and for generation of glycosyl-phosphatidyl-inositol anchors in all tissues in Drosophila melanogaster.

RESULTS

Here, we identified a putative HvUAP gene in Henosepilachna vigintioctopunctata. Knockdown of HvUAP at the second-, third- and fourth-instar stages impaired larval development. Most resultant HvUAP hypomorphs showed arrested development at the third-, fourth-instar larval or prepupal stages, and became paralyzed, depending on the age when treated. Some HvUAP-silenced larvae had weak and soft scoli. A portion of HvUAP-depleted beetles formed misshapen pupae. No HvUAP RNA interference pupae successfully emerged as adults. Dissection and microscopic observation revealed that knockdown of HvUAP affected gut growth and food ingestion, reduced cuticle thickness, and negatively affected the formation of newly generated cuticle layers during ecdysis. Furthermore, HvUAP deficiency inhibited development of the tracheal respiratory system and thinned tracheal taenidia.

CONCLUSION

The phenotypical defects in HvUAP hypomorphs suggest that HvUAP is involved in the production of chitin. Moreover, our findings will enable the development of a double-stranded RNA-based pesticide to control H. vigintioctopunctata. © 2021 Society of Chemical Industry.

RNAi suppression of the nuclear receptor FTZ-F1 impaired ecdysis, pupation, and reproduction in the 28-spotted potato ladybeetle, Henosepilachna vigintioctopunctata

Fushi-tarazu factor 1 (FTZF1) is an ecdysone-inducible transcription factor that plays a vital role during the metamorphosis in insects. In this study, we functionally characterized HvFTZ-F1 in H. vigintioctopunctata, a dreadful solanaceous crop pest, by using a dietary RNA interference technique. The HvFTZ-F1 expression levels were elevated in the 1st and 2nd-instars before molting and declined immediately after ecdysis. The HvFTZ-F1 silencing led to high mortality in the 1st instars, while the expression of the osmosis-regulative gene, HvAQPAn.G, was significantly increased in the 1st instars. HvFTZ-F1 silencing downregulated the Halloween and 20E-related genes, decreased the ecdysteroids titer, suppressed the expression of pigmentation-related genes, and reduced the catecholamines titer. In the 4th instars, HvFTZ-F1 silencing caused 100% mortality by arresting the development at the prepupal stage and preventing new abdominal cuticle formation. In the female adults, HvFTZ-F1 silencing caused an evident decrease in fecundity, prolonged the pre-oviposition period, reduced the number of eggs and hatching rate, severely atrophied the ovaries. Moreover, the 20E-related genes and the dopamine synthesis genes were suppressed in the dsHvFTZ-F1-treated females. Overall, our results revealed that HvFTZ-F1 regulates ecdysis, pupation, and reproduction in H. vigintioctopunctata, thereby could be a promising molecular target for the development of RNAi-based biopesticides to control H. vigintioctopunctata.

SMRT sequencing of the full-length transcriptome of Gekko gecko

Tokay Gecko (Gekko gecko) is a rare and endangered medicinal animal in China. Its dry body has been used as an anti-asthmatic agent for two thousand years. To date, the genome and transcriptome of this species remain poorly understood. Here, we adopted single molecule real-time (SMRT) sequencing to obtain full-length transcriptome data and characterized the transcriptome structure. We identified 882,273 circular consensus (CCS) reads, including 746,317 full-length nonchimeric (FLNC) reads. The transcript cluster analysis revealed 212,964 consensus sequences, including 203,994 high-quality isoforms. In total, 111,372 of 117,888 transcripts were successfully annotated against eight databases (Nr, eggNOG, Swiss-Prot, GO, COG, KOG, Pfam and KEGG). Furthermore, 23,877 alternative splicing events, 169,128 simple sequence repeats (SSRs), 10,437 lncRNAs and 7,932 transcription factors were predicted across all transcripts. To our knowledge, this report is the first to document the G. gecko transcriptome using SMRT sequencing. The full-length transcript data might accelerate transcriptome research and lay the foundation for further research on G. gecko.

Single-molecule long-read sequencing reveals the potential impact of posttranscriptional regulation on gene dosage effects on the avian Z chromosome

Background

Mammalian sex chromosomes provide dosage compensation, but avian lack a global mechanism of dose compensation. Herein, we employed nanopore sequencing to investigate the genetic basis of gene expression and gene dosage effects in avian Z chromosomes at the posttranscriptional level.

Results

In this study, the gonad and head skin of female and male duck samples (n = 4) were collected at 16 weeks of age for Oxford nanopore sequencing. Our results revealed a dosage effect and local regulation of duck Z chromosome gene expression. Additionally, AS and APA achieve tissue-specific gene expression, and male-biased lncRNA regulates its Z-linked target genes, with a positive regulatory role for gene dosage effects on the duck Z chromosome. In addition, GO enrichment and KEGG pathway analysis showed that the dosage effects of Z-linked genes were mainly associated with the cellular response to hormone stimulus, melanin biosynthetic, metabolic pathways, and melanogenesis, resulting in sex differences.

Conclusions

Our data suggested that post transcriptional regulation (AS, APA and lncRNA) has a potential impact on the gene expression effects of avian Z chromosomes. Our study provides a new view of gene regulation underlying the dose effects in avian Z chromosomes at the RNA post transcriptional level.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08360-8.

RNAi of vacuolar-type H+-ATPase genes causes growth delay and molting defect in Henosepilachna vigintioctopunctata

Henosepilachna vigintioctopunctata is one of the most serious insect pests to a large number of nightshades and cucurbits. RNA interference (RNAi) triggered by double-stranded RNA (dsRNA) offers a reduced risk approach to control the beetle. Identification of amenable target genes and determination of appropriate life stage for dsRNA treatment are two critical steps in order to improve RNAi efficiency. In the present paper, we identified three vATPase genes, namely HvvATPaseC, HvvATPaseE and HvvATPaseH. We found that the three transcripts were widely expressed in the eggs, first- to fourth-instar larvae, prepupae, pupae and adults. They were abundantly transcribed in the hindgut and Malpighian tubules, in contrast to the epidermis and fat body. Three days' ingestion of dsvATPaseC, dsvATPaseE and dsvATPaseH by the fourth-instar larvae significantly decreased corresponding transcript level by 90.1, 88.9 and 97.2%, greatly reduced larval fresh weight by 28.0, 29.9 and 28.0%, and caused 66.7, 100 and 78.7% larval lethality respectively. Comparably, 3 days' exposure of the third-instar larvae to dsvATPaseC significantly reduced HvvATPaseC mRNA level by 89.5%, decreased approximately 80% of the larval fresh weight, and killed 100% of the treated larvae. Therefore, the three vATPase genes, especially HvvATPaseE, are potential amenable target genes and young larvae are more susceptible to dsRNA. Our findings will enable the development of the dsRNA-based pesticide to control H. vigintioctopunctata.

RNA interference targeting ecdysone receptor blocks the larval-pupal transition in Henosepilachna vigintioctopunctata

Henosepilachna vigintioctopunctata is a serious insect pest which attacks a large number of nightshades and cucurbits in Asian countries, Brazil and Australia. Prolonged application of traditional pesticides has caused environmental pollution and exerted deleterious effects on human health. Finding new approaches with high target specificity and low environmental contamination has become an urgent task. RNA interference (RNAi) induced by double-stranded RNA (dsRNA) is expected to be applicable to managing this pest. Here we evaluated the effects of Escherichia coli-expressed dsRNAs targeting ecdysone receptor (EcR) gene via dietary delivery in laboratory and foliar spraying in a greenhouse. The target transcript was successfully knocked down when the 4th-instar larvae had fed on potato foliage dipped with dsEcR in a laboratory bioassay. Around 85% of the HvEcR RNAi larvae remained as prepupae or became abnormal pupae, and failed to emerge into adults. Ingestion of dsEcR-immersed foliage by the 3rd-instar larvae effectuated a comparable RNAi response and brought about more severe defects: all the resultant larvae arrested development, remained as prepupae and finally died. For assay in the greenhouse, a dsEcR-contained E. coli suspension was directly sprayed to the foliage of greenhouse-growing potato plants and the 3rd- and 4th-instar larvae were transferred to the leaves. High RNAi efficacy was obtained and identical RNAi phenotypes were observed in treated larvae. In addition, spraying dsEcR reduced leaf damage. Our results indicate a possibility of practical application of dsEcR as an environmentally friendly RNA pesticide to control H. vigintioctopunctata larvae.

Research Progress on Oviposition-Related Genes in Insects

Oviposition-related genes have remained a consistent focus of insect molecular biology. Previous research has gradually clarified our mechanistic understanding of oviposition-related genes, including those related to oviposition-gland-related genes, oogenesis-related genes, oviposition-site-selection-related genes, and genes related to ovulation and hatching. Moreover, some of this research has revealed how the expression of single oviposition-related genes affects the expression of related genes, and more importantly, how individual node genes function to link the expression of upstream and downstream genes. However, the research to date is not sufficient to completely explain the overall interactions among the genes of the insect oviposition system. Through a literature review of a large number of studies, this review provides references for future research on oviposition-related genes in insects and the use of RNAi or CRISPR/Cas9 technology to verify the functions of oviposition-related genes and to prevent and control harmful insects.

Functional divergence of white genes in Henosepilachna vigintioctopunctata revealed by RNA interference

Henosepilachna vigintioctopunctata is a serious pest of Solanaceae and Cucurbitaceae in many Asian countries. RNA interference (RNAi) can effectively reduce transcript abundance in this beetle, offering opportunities to explore the biological function of specific genes. The white gene encodes a half-type ATP-binding cassette transporter that plays an essential role in tryptophan, guanine and uric acid transport across membranes. Mutations that disrupt the function of white are known to cause eye pigmentation phenotypes in many insect species. Here, we found evidence for five white gene paralogues present in H. vigintioctopunctata transcriptome datasets sequenced from a range of developmental stages. We individually knocked down each of the five white genes through the injection of corresponding double-stranded RNAs (dsRNAs) to the fourth-instar larvae to determine whether functional divergence has occurred. We found that injecting 1 μg dswhite3 caused compound eye colour of pupae and adults to develop as red/brown and brown, respectively, compared with black eyes in control beetles. Injection of 2 μg dswhite3 increased RNAi efficacy and produced a clearer eye colour phenotype. At both doses, the ocular diaphragm (a ring of black pigment surrounding each eye) did not change in the white3 RNAi hypomorphs. Moreover, our data revealed that injection of dswhite2 at the fourth-instar larval stage impaired the climbing ability of both male and female adults. Our results confirmed, for the first time, functional divergence of duplicated white genes in an insect species.