Transcriptome and gene expression analysis of three developmental stages of the coffee berry borer, Hypothenemus hampei

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.

Combined PacBio Iso-Seq and Illumina RNA-Seq Analysis of the Tuta absoluta (Meyrick) Transcriptome and Cytochrome P450 Genes

Tuta absoluta (Meyrick) is a devastating invasive pest worldwide. The abamectin and chlorantraniliprole complex have become an alternative option for chemical control because they can enhance insecticidal activity and delay increased drug resistance. Notably, pests are inevitably resistant to various types of insecticides, and compound insecticides are no exception. To identify potential genes involved in the detoxification of abamectin and chlorantraniliprole complex in T. absoluta, PacBio SMRT-seq transcriptome sequencing and Illumina RNA-seq analysis of abamectin and chlorantraniliprole complex-treated T. absoluta were performed. We obtained 80,492 non-redundant transcripts, 62,762 (77.97%) transcripts that were successfully annotated, and 15,524 differentially expressed transcripts (DETs). GO annotation results showed that most of these DETs were involved in the biological processes of life-sustaining activities, such as cellular, metabolic, and single-organism processes. The KEGG pathway enrichment results showed that the pathways related to glutathione metabolism, fatty acid and amino acid synthesis, and metabolism were related to the response to abamectin and chlorantraniliprole complex in T. absoluta. Among these, 21 P450s were differentially expressed (11 upregulated and 10 downregulated). The qRT-PCR results for the eight upregulated P450 genes after abamectin and chlorantraniliprole complex treatment were consistent with the RNA-Seq data. Our findings provide new full-length transcriptional data and information for further studies on detoxification-related genes in T. absoluta.

Transcriptome analysis of aphids exposed to glandular trichomes in tomato reveals stress and starvation related responses

Understanding the responses of insect herbivores to plant chemical defences is pivotal for the management of crops and pests. However, the mechanisms of interaction are not entirely understood. In this study, we compared the whole transcriptome gene expression of the aphid Macrosiphum euphorbiae grown on two different varieties of tomato that differ in their inducible chemical defences. We used two isogenic lines of tomato with a shared genetic background that only differ in the presence of type IV glandular trichomes and their associated acylsucrose excretions. This works also reports a de novo transcriptome of the aphid M. euphorbiae. Subsequently, we identified a unique and distinct gene expression profile for the first time corresponding to aphid´s exposure to type IV glandular trichomes and acylsugars. The analysis of the aphid transcriptome shows that tomato glandular trichomes and their associated secretions are highly efficient in triggering stress-related responses in the aphid, and demonstrating that their role in plant defence goes beyond the physical impediment of herbivore activity. Some of the differentially expressed genes were associated with carbohydrate, lipid and xenobiotic metabolisms, immune system, oxidative stress response and hormone biosynthesis pathways. Also, the observed responses are compatible with a starvation syndrome. The transcriptome analysis puts forward a wide range of genes involved in the synthesis and regulation of detoxification enzymes that reveal important underlying mechanisms in the interaction of the aphid with its host plant and provides a valuable genomic resource for future study of biological processes at the molecular level using this aphid.

Phylogenetic Analysis of Small Hive Beetles From Native to Introduced Populations

The small hive beetle (SHB), a social parasite of beehives, is native to sub-Saharan Africa and has spread to America, Europe, and Australia. Recently, these beetles invaded China, causing widespread colony collapses in the honeybee, Apis cerana. In this study, single nucleotide polymorphisms (SNPs) were identified in the beetle genome from its native range (Africa), a region that was invaded by SHBs nearly 30 years ago (America), and more recent invasions (Asia). The beetles in the United States formed the earliest branch and show signs of two decades of gene flow and local adaptation to differentiate this population from the native ones. The beetles in China were deep branched and showed the highest fixation index when compared to the US populations. The number of SNPs in overexpressed genes was significantly higher than the transcriptome. Gene-expression profiles presented here distinguish the characters between adult and larvae SHBs.

Chromosome-level genome assembly, annotation, and phylogenomics of the gooseneck barnacle Pollicipes pollicipes

BACKGROUND

The barnacles are a group of >2,000 species that have fascinated biologists, including Darwin, for centuries. Their lifestyles are extremely diverse, from free-swimming larvae to sessile adults, and even root-like endoparasites. Barnacles also cause hundreds of millions of dollars of losses annually due to biofouling. However, genomic resources for crustaceans, and barnacles in particular, are lacking.

RESULTS

Using 62× Pacific Biosciences coverage, 189× Illumina whole-genome sequencing coverage, 203× HiC coverage, and 69× CHi-C coverage, we produced a chromosome-level genome assembly of the gooseneck barnacle Pollicipes pollicipes. The P. pollicipes genome is 770 Mb long and its assembly is one of the most contiguous and complete crustacean genomes available, with a scaffold N50 of 47 Mb and 90.5% of the BUSCO Arthropoda gene set. Using the genome annotation produced here along with transcriptomes of 13 other barnacle species, we completed phylogenomic analyses on a nearly 2 million amino acid alignment. Contrary to previous studies, our phylogenies suggest that the Pollicipedomorpha is monophyletic and sister to the Balanomorpha, which alters our understanding of barnacle larval evolution and suggests homoplasy in a number of naupliar characters. We also compared transcriptomes of P. pollicipes nauplius larvae and adults and found that nearly one-half of the genes in the genome are differentially expressed, highlighting the vastly different transcriptomes of larvae and adult gooseneck barnacles. Annotation of the genes with KEGG and GO terms reveals that these stages exhibit many differences including cuticle binding, chitin binding, microtubule motor activity, and membrane adhesion.

CONCLUSION

This study provides high-quality genomic resources for a key group of crustaceans. This is especially valuable given the roles P. pollicipes plays in European fisheries, as a sentinel species for coastal ecosystems, and as a model for studying barnacle adhesion as well as its key position in the barnacle tree of life. A combination of genomic, phylogenetic, and transcriptomic analyses here provides valuable insights into the evolution and development of barnacles.

Gut transcriptome of two bark beetle species stimulated with the same kairomones reveals molecular differences in detoxification pathways

Dendroctonus bark beetles are the most destructive agents in coniferous forests. These beetles come into contact with the toxic compounds of their host's chemical defenses throughout their life cycle, some of which are also used by the insects as kairomones to select their host trees during the colonization process. However, little is known about the molecular mechanisms by which the insects counteract the toxicity of these compounds. Here, two sibling species of bark beetles, D. valens and D. rhizophagus, were stimulated with vapors of a blend of their main kairomones (α-pinene, β-pinene and 3-carene), in order to compare the transcriptional response of their gut. A total of 48 180 unigenes were identified in D. valens and 43 704 in D. rhizophagus, in response to kairomones blend. The analysis of differential gene expression showed a transcriptional response in D. valens (739 unigenes, 0.58–10.36 Log2FC) related to digestive process and in D. rhizophagus (322 unigenes 0.87–13.08 Log2FC) related to xenobiotics metabolism. The expression profiles of detoxification genes mainly evidenced the up-regulation of COEs and GSTs in D. valens, and the up-regulation of P450s in D. rhizophagus. Results suggest that terpenes metabolism comes accompanied by an integral hormetic response, result of compensatory mechanisms, including the activation of other metabolic pathways, to ensure the supply of energy and the survival of organisms which is specific for each species, according to its life history and ecological strategy.

CARACTERIZACIÓN Y DISMINUCIÓN EN LA EXPRESIÓN DE UNA QUITINA SINTASA MEDIADA POR ARNi EN Hypothenemus hampei (CURCULIONIDAE)

La broca del café, Hypothenemus hampei, es un insecto plaga que causa daños significativos al grano de café y grandes pérdidas económicas a los productores en todo el mundo. Al igual que otros insectos, la broca de café requiere de quitina sintasas (CHS) para la biosíntesis de la quitina, componente principal del exoesqueleto del insecto, y de vital importancia para su crecimiento y desarrollo. En este estudio, el gen CHS1 de la broca del café (HhCHS1) fue identificado, caracterizado y posteriormente silenciado mediante el uso de ARNi, mecanismo que permite degradar el ARNm e interrumpir la expresión de proteínas de interés en un organismo. Los perfiles de expresión del gen HhCHS1, medidos por RT-qPCR, mostraron niveles de expresión diferencial en las diferentes etapas del desarrollo del insecto. Los niveles más altos de expresión se encontraron en larvas de segundo estadio (L2) y machos adultos. El ARNcd administrado por vía oral, a concentraciones de 2 µg/100 µL, generó un silenciamiento efectivo del gen HhCHS1 (84 %) después de 7 días de tratamiento. Estos resultados sugieren que el gen HhCHS1 desempeña un papel importante en el desarrollo del insecto, y que, por ende, podría usarse como objetivo para desarrollar nuevas estrategias de manejo de este insecto plaga, mediante el uso de ARNi.  Palabras Clave: ARNcd, café, insecto plaga, quitina-sintasa, RT-qPCR.

Trends and Applications of Omics Technologies to Functional Characterisation of Enzymes and Protein Metabolites Produced by Fungi

Identifying and adopting industrial applications for proteins and enzymes derived from fungi strains have been at the focal point of several studies in recent times. To facilitate such studies, it is necessary that advancements and innovation in mycological and molecular characterisation are concomitant. This review aims to provide a detailed overview of the necessary steps employed in both qualitative and quantitative research using the omics technologies that are pertinent to fungi characterisation. This stems from the understanding that data provided from the functional characterisation of fungi and their metabolites is important towards the techno-economic feasibility of large-scale production of biological products. The review further describes how the functional gaps left by genomics, internal transcribe spacer (ITS) regions are addressed by transcriptomics and the various techniques and platforms utilised, including quantitive reverse transcription polymerase chain reaction (RT-qPCR), hybridisation techniques, and RNA-seq, and the insights such data provide on the effect of environmental changes on fungal enzyme production from an expressional standpoint. The review also offers information on the many available bioinformatics tools of analysis necessary for the analysis of the overwhelming data synonymous with the omics approach to fungal characterisation.

A coffee berry borer (Hypothenemus hampei) genome assembly reveals a reduced chemosensory receptor gene repertoire and male-specific genome sequences

Coffee berry borer-CBB (Hypothenemus hampei) is a globally important economic pest of coffee (Coffea spp.). Despite current insect control methods for managing CBB, development of future control strategies requires a better understanding of its biology and interaction with its host plant. Towards this objective, we performed de novo CBB genome and transcriptome sequencing, improved CBB genome assembly and predicted 18,765 protein-encoding genes. Using genome and transcriptome data, we annotated the genes associated with chemosensation and found a reduced gene repertoire composed by 67 odorant receptors (ORs), 62 gustatory receptors (GRs), 33 ionotropic receptors (IRs) and 29 odorant-binding proteins (OBPs). In silico transcript abundance analysis of these chemosensory genes revealed expression enrichment in CBB adults compared with larva. Detection of differentially expressed chemosensory genes between males and females is likely associated with differences in host-finding behavior between sexes. Additionally, we discovered male-specific genome content and identified candidate male-specific expressed genes on these scaffolds, suggesting that a Y-like chromosome may be involved in the CBB's functional haplodiploid mechanism of sex determination.

The Developmental Transcriptome of Bagworm, Metisa plana (Lepidoptera: Psychidae) and Insights into Chitin Biosynthesis Genes

Bagworm, Metisa plana (Lepidoptera: Psychidae) is a ubiquitous insect pest in the oil palm plantations. M. plana infestation could reduce the oil palm productivity by 40% if it remains untreated over two consecutive years. Despite the urgency to tackle this issue, the genome and transcriptome of M. plana have not yet been fully elucidated. Here, we report a comprehensive transcriptome dataset from four different developmental stages of M. plana, comprising of egg, third instar larva, pupa and female adult. The de novo transcriptome assembly of the raw data had produced a total of 193,686 transcripts, which were then annotated against UniProt, NCBI non-redundant (NR) database, Gene Ontology, Cluster of Orthologous Group, and Kyoto Encyclopedia of Genes and Genomes databases. From this, 46,534 transcripts were annotated and mapped to 146 known metabolic or signalling KEGG pathways. The paper further identified 41 differentially expressed transcripts encoding seven genes in the chitin biosynthesis pathways, and their expressions across each developmental stage were further analysed. The genetic diversity of M. plana was profiled whereby there were 21,516 microsatellite sequences and 379,895 SNPs loci found in the transcriptome of M. plana. These datasets add valuable transcriptomic resources for further study of developmental gene expression, transcriptional regulations and functional gene activities involved in the development of M. plana. Identification of regulatory genes in the chitin biosynthesis pathway may also help in developing an RNAi-mediated pest control management by targeting certain pathways, and functional studies of the genes in M. plana.

Transcriptome and gene expression analysis of Rhynchophorus ferrugineus (Coleoptera: Curculionidae) during developmental stages

Background

Red palm weevil, Rhynchophorus ferrugineus Olivier, is one of the most destructive pests harming palm trees. However, genomic resources for R. ferrugineus are still lacking, limiting the ability to discover molecular and genetic means of pest control.

Methods

In this study, PacBio Iso-Seq and Illumina RNA-seq were used to generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th-instar larva, adult) to increase the understanding of the life cycle and molecular characteristics of the pest.

Results

Sequencing generated 625,983,256 clean reads, from which 63,801 full-length transcripts were assembled with N50 of 3,547 bp. Expression analyses revealed 8,583 differentially expressed genes (DEGs). Moreover, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these DEGs were mainly related to the peroxisome pathway which associated with metabolic pathways, material transportation and organ tissue formation. In summary, this work provides a valuable basis for further research on the growth and development, gene expression and gene prediction, and pest control of R. ferrugineus.

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

Background

Red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is one of the most destructive insects for palm trees in the world. However, its genome resources are still in the blank stage, which limits the study of molecular and growth development analysis.

Methods

In this study, we used PacBio Iso-Seq and Illumina RNA-seq to first generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th larva, female and male) to increase our understanding of the life cycle and molecular characteristics of R. ferrugineus.

Results

A total of 63,801 nonredundant full-length transcripts were generated with an average length of 2,964 bp from three developmental stages, including the 7th instar larva, pupa, female adult and male adult. These transcripts showed a high annotation rate in seven public databases, with 54,999 (86.20%) successfully annotated. Meanwhile, 2,184 alternative splicing (AS) events, 2,084 transcription factors (TFs), 66,230 simple sequence repeats (SSR) and 9,618 Long noncoding RNAs (lncRNAs) were identified. In summary, our results provide a new source of full-length transcriptional data and information for the further study of gene expression and genetics in R. ferrugineus.

Characterization of Dendrolimus houi Lajonquiere (Lepidoptera: Lasiocampidae) Transcriptome across All Life Stages

Dendrolimus houi Lajonquiere is a phytophagous caterpillar infesting many economically important coniferous tree species in China, causing serious economic and ecological environment losses. Based on previous research, it has one generation per year in South China and East China in contrast to two generations per year in Yunnan province in southwestern China. The species is potentially resilient to climatic extremes in these regions with the eggs and 1st instar larvae surviving in the winter (5 °C), older instar larvae and pupae surviving high temperatures in the summer (35 °C), suggesting some temperature stress tolerance during different developmental stages. However, little is known in this species at the genetic and genomic level. In this study, we used high throughput sequencing to obtain transcriptome data from different developmental stages (eggs, 1st-3rd instar larvae, 4th-5th instar larvae, 6th-7th instar larvae, pupae, male and female adults), which were collected from Fujian province. In total, we obtained approximately 90 Gb of data, from which 33,720 unigenes were assembled and 17,797 unigenes were annotated. We furtherly analyzed the differentially expressed genes (DGEs) across all stages, the largest number between the eggs and 1st instar larvae stage and gene expression varied significantly in different developmental stages. Furthermore, 4138 SSR genes and 114,977 SNP loci were screened from transcriptome data. This paper will be a foundation for further study towards improved integrated pest management strategies for this species.