SMRT sequencing of the full-length transcriptome of Odontotermes formosanus (Shiraki) under Serratia marcescens treatment

Unraveling metagenomics through long-read sequencing: a comprehensive review

The study of microbial communities has undergone significant advancements, starting from the initial use of 16S rRNA sequencing to the adoption of shotgun metagenomics. However, a new era has emerged with the advent of long-read sequencing (LRS), which offers substantial improvements over its predecessor, short-read sequencing (SRS). LRS produces reads that are several kilobases long, enabling researchers to obtain more complete and contiguous genomic information, characterize structural variations, and study epigenetic modifications. The current leaders in LRS technologies are Pacific Biotechnologies (PacBio) and Oxford Nanopore Technologies (ONT), each offering a distinct set of advantages. This review covers the workflow of long-read metagenomics sequencing, including sample preparation (sample collection, sample extraction, and library preparation), sequencing, processing (quality control, assembly, and binning), and analysis (taxonomic annotation and functional annotation). Each section provides a concise outline of the key concept of the methodology, presenting the original concept as well as how it is challenged or modified in the context of LRS. Additionally, the section introduces a range of tools that are compatible with LRS and can be utilized to execute the LRS process. This review aims to present the workflow of metagenomics, highlight the transformative impact of LRS, and provide researchers with a selection of tools suitable for this task.

A ferritin protein is involved in the development and reproduction of the whitefly, Bemisia tabaci

Ferritins are conserved iron-binding proteins that exist in most living organisms and play an essential role in the maintenance of cellular iron homeostasis. Although ferritin has been studied in many species, little is known about its role in the whitefly, Bemisia tabaci. In this study, we identified an iron-binding protein from B. tabaci and named it BtabFer1. The full-length cDNA of BtabFer1 is 1,043 bp and encodes a protein consisting of 224 amino acids with a deduced molecular weight of 25.26 kDa, and phylogenetic analysis shows that BtabFer1 is conserved among Hemiptera insects. The expression levels of BtabFer1 in different developmental stages and tissues were analyzed by real-time PCR, and results showed that BtabFer1 was ubiquitously expressed at all developmental stages and in all examined tissues. The RNAi-mediated knockdown of BtabFer1 caused a significant reduction in survival rate, egg production, and egg hatching rate of whiteflies. Knockdown of BtabFer1 also inhibited the transcription of genes in the juvenile hormone signal transduction pathway. Taken together, these results suggest that BtabFer1 plays a critical role in the development and reproduction of whiteflies. This study can broaden our understanding of ferritin in insect fecundity and development, as well as provide baseline data for future studies.

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.

Novel insights on genes and pathways involved in Pinus elliottii response to resinosis

Pinus elliottii, an important coniferous timber species, has recently become one of the most popular sources of resin in China. Resinosis is a common disease that may negatively affect pine tree growth and production. In this study, we used single-molecule real-time sequencing and Illumina RNA sequencing to generate an accurate transcriptome for P. elliottii. The transcriptome included 90,026 transcripts, 5160 long non-coding RNAs and 7710 transcription factors. We then analyzed RNA-sequencing, small RNA-sequencing and degradome data to identify genes, miRNAs and key miRNA-target pairs involved in response to resinosis in P. elliottii. We identified 1305 genes and 1151 miRNAs exhibiting significant differential expression in response to resinosis. According to the degradome sequencing analysis, 318 differentially expressed transcripts were targets of 14 differentially expressed miRNAs. Our study has provided resources for further functional characterization of genes and miRNAs involved in resinosis in P. elliottii, which should aid the future disease-resistance breeding of this species.

OfGNBP silencing enhances the toxicity of Serratia marcescens Bizio (SM1) to Odontotermes formosanus (Shiraki)

The immunity of insects plays a vital role in their survival. Our experiments found that lipopolysaccharide (LPS) and glucono-δ-lactone (GDL) could influence the virulence of Serratia marcescens Bizio (SM1) to Odontotermes formosanus (Shiraki) by affecting the immunity. Gram-negative binding proteins (GNBPs) are an important pattern recognition proteins that play a crucial role in the innate immune system. Therefore, two OfGNBPs were cloned in O. formosanus. The expression of OfGNBPs was significantly changed by LPS,SM1 and GDL, not prick. In addition, the immune-related gene expression, the phenoloxidase activity and antibacterial activity of donor termites and recipient termites were significantly induced by SM1. Furthermore, the knockdown of OfGNBP by RNA interference reduced not only individual immunity but also social immunity in O. formosanus, which increased the virulence of SM1 to O. formosanus. Importantly, dsOfGNBP alone also had good control effect on O. formosanus. In summary, we concluded that dsOfGNBPs are important termite immunosuppressants.

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.

Termicin silencing enhances the toxicity of Serratia marcescens Bizio (SM1) to Odontotermes formosanus (Shiraki)

Termites are often exposed to a variety of pathogens during their life cycle, which has led to the development of an innate immune system to resist these pathogens. Antimicrobial peptides (AMPs) play a crucial role in the innate immune system in insects. However, clear information on AMPs in termites has not been obtained. Therefore, exploring the function of AMPs in the subterranean termite Odontotermes formosanus (Shiraki) can lead to the development of novel termite control strategies that integrate RNA interference (RNAi) and pathogens. Here we first obtained two Oftermicins from O. formosanus and observed that the expression of these Oftermicin genes was significantly upregulated at the mRNA level after treatment with lipopolysaccharide (LPS) or Serratia marcescens Bizio (SM1). Interestingly, the expression of these Oftermicins increased not only in the donor termites but also in the recipient termites through transmission experiments. Bioassay experiments showed that the mortality of O. formosanus treated with SM1 after RNAi was significantly higher than that of other groups. In summary, dsOftermicins are important immunosuppressants for termite control and Oftermicins are optimal targets for termite control based on the combined use of RNAi and pathogens.

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.

Full-Length SMRT Transcriptome Sequencing and SSR Analysis of Bactrocera dorsalis (Hendel)

Simple Summary

In this study, a full-length transcriptome was analyzed with single-molecule real-time (SMRT) sequencing, which was first used to discover simple sequence repeat (SSR) genetic markers from B. dorsalis. Moreover, SSR markers from isoforms were screened for the identification of species diversity. These results could provide molecular biology methods for further population research.

Abstract

Bactrocera dorsalis (Hendel), as one of the most notorious and destructive invasive agricultural pests in the world, causes damage to over 250 different types of fruits and vegetables throughout tropical and subtropical areas. PacBio single-molecule real-time (SMRT) sequencing was used to generate the full-length transcriptome data of B. dorsalis. A total of 40,319,890 subreads (76.6 Gb, clean reads) were generated, including 535,241 circular consensus sequences (CCSs) and 386,916 full-length non-concatemer reads (FLNCs). Transcript cluster analysis of the FLNC reads revealed 22,780 high-quality reads (HQs). In total, 12,274 transcripts were functionally annotated based on four different databases. A total of 1978 SSR loci were distributed throughout 1714 HQ transcripts, of which 1926 were complete SSRs and 52 were complex SSRs. Among the total SSR loci, 2–3 nucleotide repeats were dominant, occupying 83.62%, of which di- and tri- nucleotide repeats were 39.38% and 44.24%, respectively. We detected 105 repeat motifs, of which AT/AT (50.19%), AC/GT (39.15%), CAA/TTG (32.46%), and ACA/TGT (10.86%) were the most common in di- and tri-nucleotide repeats. The repeat SSR motifs were 12–190 bp in length, and 1638 (88.02%) were shorter than 20 bp. According to the randomly selected microsatellite sequence, 80 pairs of primers were designed, and 174 individuals were randomly amplified by PCR using primers. The number of primers that had amplification products with clear bands and showed good polymorphism came to 41, indicating that this was a feasible way to explore SSR markers from the transcriptomic data of B. dorsalis. These results lay a foundation for developing highly polymorphic microsatellites for researching the functional genomics, population genetic structure, and genetic diversity of B. dorsalis.

Effect of acute ammonia toxicity on inflammation, oxidative stress and apoptosis in head kidney macrophage of Pelteobagrus fulvidraco and the alleviation of curcumin

Ammonia is one of the most major pollutant and stress factors of aquaculture systems, and has seriously endangered fish health. However, few studies have been performed on mechanisms of the detrimental impact of ammonia stress and mitigation in fish. A study was carried out to investigate the response of genes involved in inflammation, antioxidation, polarization and apoptosis in head kidney macrophages to acute ammonia toxicity, and the alleviation effect of curcumin. The cells were divided into six groups, as follows: The control group composed of untreated macrophages (CON), the experimental groups, consisting of macrophages treated with 0.23 mg L^-1 ammonia (AM), 45 μmol L^-1 curcumin (CUR), 0.23 mg L^-1 ammonia and 5 μmol L^-1 curcumin (5A), 0.23 mg L^-1 ammonia and 25 μmol L^-1 curcumin (25A), 0.23 mg L^-1 ammonia and 45 μmol L^-1 curcumin (45A). The cells were pretreated with different concentrations of curcumin for 1 h and then incubated with ammonia for 24 h. The results showed that ammonia poisoning could increase ROS levels, up-regulate the expression of antioxidant enzymes (SOD and GPx), inflammatory cytokines (IL-1, IL-6 and TNF-α) and inflammatory mediators (NF-κB p65 and COX-2), decrease cell viability, down-regulate the expression of M2 marker (Arg-1) and anti-apoptosis (Bcl-2), but curcumin could alleviate the adverse effect of ammonia toxicity. Overall, these results have important implications for understanding of the mechanism of ammonia toxicity and the mitigating effect of curcumin in fish.

Full-length SMRT transcriptome sequencing and microsatellite characterization in Paulownia catalpifolia

Paulownia catalpifolia is an important, fast-growing timber species known for its high density, color and texture. However, few transcriptomic and genetic studies have been conducted in P. catalpifolia. In this study, single-molecule real-time sequencing technology was applied to obtain the full-length transcriptome of P. catalpifolia leaves treated with varying degrees of drought stress. The sequencing data were then used to search for microsatellites, or simple sequence repeats (SSRs). A total of 28.83 Gb data were generated, 25,969 high-quality (HQ) transcripts with an average length of 1624 bp were acquired after removing the redundant reads, and 25,602 HQ transcripts (98.59%) were annotated using public databases. Among the HQ transcripts, 16,722 intact coding sequences, 149 long non-coding RNAs and 179 alternative splicing events were predicted, respectively. A total of 7367 SSR loci were distributed throughout 6293 HQ transcripts, of which 763 complex SSRs and 6604 complete SSRs. The SSR appearance frequency was 28.37%, and the average distribution distance was 5.59 kb. Among the 6604 complete SSR loci, 1-3 nucleotide repeats were dominant, occupying 97.85% of the total SSR loci, of which mono-, di- and tri-nucleotide repeats were 44.68%, 33.86% and 19.31%, respectively. We detected 112 repeat motifs, of which A/T (42.64%), AG/CT (12.22%), GA/TC (9.63%), GAA/TTC (1.57%) and CCA/TGG (1.54%) were most common in mono-, di- and tri-nucleotide repeats, respectively. The length of the repeat SSR motifs was 10-88 bp, and 4997 (75.67%) were ≤ 20 bp. This study provides a novel full-length transcriptome reference for P. catalpifolia and will facilitate the identification of germplasm resources and breeding of new drought-resistant P. catalpifolia varieties.