SMRT sequencing of the full-length transcriptome of the Rhynchophorus ferrugineus (Coleoptera: Curculionidae)

Full-length transcriptome annotation of a pyrosome, Pyrosoma atlanticum (Chordata, Thaliacea)

Pyrosomes represent a group of marine holozooplankton. They do not exist as individual entities but as intricate colonies comprising numerous zooids, well-known for their bioluminescent properties. Pyrosomes inhabit most oceanic waters spanning tropical to temperate regions. They serve as significant consumers of phytoplankton, playing a crucial role in transport of organic carbon and sedimentation of organic matter debris. Nonetheless, the knowledge about this group remains limited, particularly concerning genomic and gene aspects. To further investigate the pyrosome at the molecular level, we used Pyrosoma atlanticum, which is the best known of pyrosome, for investigation in this study. We performed PacBio Iso-Seq and Illumina RNA-seq to generate high-quality and full-length transcript data from P. atlanticum. The systematic gene functional annotation was performed by integrated data. The full-length transcriptome produced in this study represents the inaugural dataset of transcriptome within the class Thaliacea, serving as a reference for future investigations. Additionally, as a stem group for chordates, the information of pyrosomes can provide a valuable research foundation for the evolution of vertebrates.

Identification of testis development-related genes by combining Iso-Seq and RNA-Seq in Zeugodacus tau

Introduction: Zeugodacus tau (Walker) is an invasive pest. An effective method to control this pest is the sterile insect technique (SIT). To better apply this technique, it is necessary to understand testis development progression. Methods: Differentially expressed genes (DEGs) during testis development were analyzed by PacBio Iso-Seq and RNA-seq. Results: RNA-Seq library of Z. tau testes on day 1, 6, and 11 post eclosion were constructed. We identified 755 and 865 differentially expressed genes in the comparisons of T6 (testes on day 6) vs. T1 and T11 vs. T1, respectively. The KEGG pathway analysis showed that the DEGs were significantly enriched in retinol metabolism, vitamin B6 metabolism, and ascorbate and aldarate metabolism pathways. Knockdown of retinol dehydrogenase 12-like (rdh12-like), pyridoxal kinase (pdxk) and regucalcin (rgn), the representative gene in each of the above 3 pathways, reduced the hatching rate of Z. tau offspring. In addition, we identified 107 Drosophila spermatogenesis-related orthologous genes in Z. tau, of which innexin 2 (inx2) exhibited significantly up-regulated expression throughout testis development, and the knockdown of this gene reduced offspring hatching rate. Discussion: Our data indicated that rdh12-like, pdxk, rgn, and inx2 genes were related to testis development, and they were conserved in tephritid species. These results suggested that this gene might have the same function in tephritid. The findings provide an insight into testis development and spermatogenesis in tephritid species.

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.

Omics in the Red Palm Weevil Rhynchophorus ferrugineus (Olivier) (Coleoptera: Curculionidae): A Bridge to the Pest

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is the most devastating pest of palm trees worldwide. Mitigation of the economic and biodiversity impact it causes is an international priority that could be greatly aided by a better understanding of its biology and genetics. Despite its relevance, the biology of the RPW remains poorly understood, and research on management strategies often focuses on outdated empirical methods that produce sub-optimal results. With the development of omics approaches in genetic research, new avenues for pest control are becoming increasingly feasible. For example, genetic engineering approaches become available once a species's target genes are well characterized in terms of their sequence, but also population variability, epistatic interactions, and more. In the last few years alone, there have been major advances in omics studies of the RPW. Multiple draft genomes are currently available, along with short and long-read transcriptomes, and metagenomes, which have facilitated the identification of genes of interest to the RPW scientific community. This review describes omics approaches previously applied to RPW research, highlights findings that could be impactful for pest management, and emphasizes future opportunities and challenges in this area of research.

Tolerance to dietary linalool primarily involves co-expression of cytochrome P450s and cuticular proteins in Pagiophloeus tsushimanus (Coleoptera: Curculionidae) larvae using SMRT sequencing and RNA-seq

BACKGROUND

Pagiophloeus tsushimanus (Coleoptera: Curculionidae), an emerging forest pest exclusively infesting camphor trees, has recently caused severe ecological and economic damage in localized areas in China. Its population outbreak depends largely on the capacity to overcome the pressure of terpenoid-derived metabolites (e.g. linalool) from camphor trees. At present, the molecular basis of physiological adaptation of P. tsushimanus to dietary linalool is poorly understood, and there is no available reference genome or transcriptome.

RESULTS

Herein, we constructed the transcriptome profiling of P. tsushimanus larvae reared on linalool-infused diets using RNA sequencing and single-molecule real-time sequencing. A total of 20,325 high-quality full-length transcripts were identified as a reference transcriptome, of which 14,492 protein-coding transcripts including 130 transcription factors (TFs), and 5561 long non-coding RNAs (lncRNAs) were detected. Also, 30 alternative splicing events and 8049 simple sequence repeats were captured. Gene ontology enrichment of differential expressed transcripts revealed that overall up-regulation of both cytochrome P450s (CYP450s) and cuticular proteins (CPs), was the primary response characteristic against dietary linalool. Other physiological effects possibly caused by linalool exposure, such as increase in Reactive Oxygen Species (ROS) and hormetic stimulation, were compensated by a handful of induced genes encoding antioxidases, heat shock proteins (HSPs), juvenile hormone (JH) epoxide hydrolases, and digestive enzymes. Additionally, based on co-expression networks analysis, a diverse array of hub lncRNAs and TFs co-expressed with CYP450s and CPs were screened as the potential gene regulators. Temporal expression of candidate transcripts determined by quantitative real-time PCR also indicated a cooperative relationship between the inductions of CYP450s and CPs upon exposure to linalool.

CONCLUSIONS

Our present study provides an important transcriptome resource of P. tsushimanus, and lays a valuable foundation for understanding how this specialist pest copes with chemical challenges in its specific host environments.

The transcriptomic response of Hyphantria cunea (Drury) to the infection of Serratia marcescens Bizio based on full-length SMRT transcriptome sequencing

Hyphantria cunea (Drury) is a globally important forest pest. We found that the Serratia marcescens Bizio strain SM1 had insecticidal activity against H. cunea, but the transcriptomic response of H. cunea to SM1 were not clear. Therefore, we performed full-length sequencing of the transcriptomes of H. cunea larvae infected with SM1 and the control group. A total of 1,183 differentially expressed genes (DEGs) were identified by comparing the group infected with SM1 and the control group, including 554 downregulated genes and 629 upregulated genes. We found many downregulated genes in metabolic pathways. Furthermore, some of these downregulated genes were involved in cellular immunity, melanization, and detoxification enzymes, which showed that SM1 weakened H. cunea immunity. In addition, genes in the juvenile hormone synthesis pathway were upregulated, which was detrimental to the survival of H. cunea. This research analyzed the transcriptomic response of H. cunea to SM1 by high-throughput full-length transcriptome sequencing. The results provide useful information to explore the relationship between S. marcescens and H. cunea, and theoretical support for the application of S. marcescens and the control of H. cunea in the future.

The combination of SMRT sequencing and Illumina sequencing highlights organ-specific and age-specific expression patterns of miRNAs in Sika Deer

Due to the lack of high-quality Sika Deer (Cervus nippon) transcriptome and sRNAome across multiple organs or development stages, it is impossible to comprehensively analyze the mRNA and miRNA regulatory networks related to growth, development and immunity response. In this study, we used single molecule-real time sequencing (SMRT-seq) and Illumina sequencing methods to generate transcriptome and sRNAome from ten tissues and four age groups of Sika Deer to help us understand molecular characteristics and global miRNA expression profiles. The results showed that a total of 240,846 consensus transcripts were generated with an average length of 2,784 bp. 4,329 Transcription factors (TFs), 109,000 Simple Sequence Repeats (SSRs) and 18,987 Long non-coding RNAs (LncRNAs) were identified. Meanwhile, 306 known miRNAs and 143 novel miRNAs were obtained. A large number of miRNAs showed organ-specific and age-specific differential expression patterns. In particular, we found that the organ-specific miRNAs were enriched in the brain, some of which shared only between the brain and adrenal. These miRNAs were involved in maintaining specific functions within the brain and adrenal. By constructing miRNA96mRNA interaction networks associated with Sika Deer immunity, we found that miRNAs (miR-148a, miR-26a, miR-214, let-7b, etc.) and mRNAs (CD6, TRIM38, C3, CD163, etc.) might play an important role in the immune response of Sika Deer spleen. Together, our study generated an improved transcript annotation for Sika Deer by SMRT-seq and revealed the role of miRNA in regulating the growth, development and immunity response of Sika Deer.

A full-length transcriptome and gene expression analysis of three detoxification gene families in a predatory stink bug, Picromerus lewisi

The predatory stink bug P. Lewisi shows potential for Integrated Pest Management programs for controlling Lepidoptera pest insects in crops and forests. The importance of this insect for biological control has stimulated several studies into its biology and ecology. However, P. lewisi has little genetic information available. In the present study, PacBio single-molecule real-time (SMRT) sequencing and Illumina RNA-seq sequencing technologies were used to reveal the full-length transcriptome profiling and tissue-specific expression patterns of P. lewisi. A total of 12,997 high-quality transcripts with an average length of 2,292 bp were obtained from different stages of P. lewisi using SMRT sequencing. Among these, 12,101 were successfully annotated in seven public databases. A total of 67 genes of cytochrome P450 monooxygenases, 43 carboxylesterase genes, and 18 glutathione S-transferase genes were identified, most of which were obtained with full-length ORFs. Then, tissue-specific expression patterns of 5th instar nymphs were analyzed using Illumina sequencing. Several candidate genes related to detoxification of insecticides and other xenobiotics as well as the degradation of odors, were identified in the guts and antennae of P. lewisi. The current study offered in-depth knowledge to understand the biology and ecology of this beneficial predator and related species.

SMRT Sequencing Reveals Candidate Genes and Pathways With Medicinal Value in Cipangopaludina chinensis

Cipangopaludina chinensis is an economically important aquatic snail with high medicinal value. However, molecular biology research on C. chinensis is limited by the lack of a reference genome, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. chinensis. Herein, we conducted the first full-length transcriptome analysis of C. chinensis using PacBio single-molecule real-time (SMRT) sequencing technology. We identified a total of 26,312 unigenes with an average length of 2,572 bp, of which the largest number of zf-c2h2 transcription factor families (120,18.24%) were found, and also observed that the majority of the 8,058 SSRs contained 4-7 repeat units, which provided data for subsequent work on snail genetics Subsequently, 91.86% (24,169) of the genes were successfully annotated to the four major databases, while the highest homology was observed with Pomacea canaliculata. Functional annotation revealed that the majority of transcripts were enriched in metabolism, signal transduction and Immune-related pathways, and several candidate genes involved in drug metabolism and immune response were identified (e.g., CYP1A1, CYP2J, CYP2U1, GST, ,PIK3, PDE3A, PRKAG). This study lays a foundation for future molecular biology research and provides a reference for studying genes associated with the medicinal value of C. chinensis.

Full-Length Transcriptome Analysis of the Halophyte Nitraria sibirica Pall

Background: Nitraria sibirica Pall. is one of the pioneer tree species in saline–alkali areas due to its extreme salt tolerance. However, the lack of information on its genome limits the further exploration of the molecular mechanisms in N. sibirica under salt stress. Methods: In this study, we used single-molecule real-time (SMRT) technology based on the PacBio Iso-Seq platform to obtain transcriptome data from N. sibirica under salt treatment for the first time, which is helpful for our in-depth analysis of the salt tolerance and molecular characteristics of N. sibirica. Results: Our results suggested that a total of 234,508 circular consensus sequences (CCSs) with a mean read length of 2121 bp were obtained from the 19.26 Gb raw data. Furthermore, based on transcript cluster analysis, 93,713 consensus isoforms were obtained, including 92,116 high-quality isoforms. After removing redundant sequences, 49,240 non-redundant transcripts were obtained from high-quality isoforms. A total of 37,261 SSRs, 1816 LncRNAs and 47,314 CDSs, of which 40,160 carried complete ORFs, were obtained. Based on our transcriptome data, we also analyzed the coding genes of H+-PPase, and the results of both bioinformatics and functional analyses indicated that the gene prediction via full-length transcripts obtained by SMRT technology is reliable and effective. In summary, our research data obtained by SMRT technology provides more reliable and accurate information for the further analysis of the regulatory network and molecular mechanism of N. sibirica under salt stress.

The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.

Haplotype-resolved genome assembly enables gene discovery in the red palm weevil Rhynchophorus ferrugineus

The red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is an economically-important invasive species that attacks multiple species of palm trees around the world. A better understanding of gene content and function in R. ferrugineus has the potential to inform pest control strategies and thereby mitigate economic and biodiversity losses caused by this species. Using 10x Genomics linked-read sequencing, we produced a haplotype-resolved diploid genome assembly for R. ferrugineus from a single heterozygous individual with modest sequencing coverage (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}∼ 62x). Benchmarking against conserved single-copy Arthropod orthologs suggests both pseudo-haplotypes in our R. ferrugineus genome assembly are highly complete with respect to gene content, and do not suffer from haplotype-induced duplication artifacts present in a recently published hybrid assembly for this species. Annotation of the larger pseudo-haplotype in our assembly provides evidence for 23,413 protein-coding loci in R. ferrugineus, including over 13,000 predicted proteins annotated with Gene Ontology terms and over 6000 loci independently supported by high-quality Iso-Seq transcriptomic data. Our assembly also includes 95% of R. ferrugineus chemosensory, detoxification and neuropeptide-related transcripts identified previously using RNA-seq transcriptomic data, and provides a platform for the molecular analysis of these and other functionally-relevant genes that can help guide management of this widespread insect pest.