Genome-wide gene expression profiling of the melon fly, Zeugodacus cucurbitae, during thirteen life stages

CRISPR/Cas9-mediated knockout of a male accessory glands-specific gene takeout1 decreases the fecundity of Zeugodacus cucurbitae female

BACKGROUND

The melon fly, Zeugodacus cucurbitae (Coquillett), is an invasive Tephritidae pest with robust fertility. The male accessory glands (MAGs) form a vital organ that ensures insect reproductive efficiency. Most of the secreted proteins by MAGs exhibit a male bias expression. Takeout, one of these proteins, is abundantly present in the MAGs of many insects.

RESULTS

In this study, we identified 32 takeout genes in Z. cucurbitae. The phylogenetic analysis and multiple sequence alignment results showed that Zctakeout1 is the most related homolog to the MAGs-specific takeout in Tephritidae. The real-time quantitative PCR results showed that Zctakeout1 was exclusively expressed in the male adult stage, and its expression level gradually increased with the increase in age and then remained stable at the sexually matured stage. The distribution among tissues demonstrated the specific expression of Zctakeout1 in the MAGs, and fluorescence immunohistochemical results confirmed the presence of Zctakeout1 in close proximity to binuclear cells of the mesoderm epidermal MAGs. In continuation, CRISPR/Cas9-mediated genome editing was employed, resulting in successfully generating a homozygous strain with an +8 bp insertion. The mating experiments with the Zctakeout1-/- males resulted in significant reductions in both the mating rate and egg production of females.

CONCLUSION

These findings prove that the MAGs-specific Zctakeout1 is essential in regulating fecundity in female Z. cucurbitae fruit flies. Our data suggests its utilization in future essential insect-specific gene-directed sterility insect technique (SIT) by the genetic manipulation to keep these important Tephritidae populations under control. © 2024 Society of Chemical Industry.

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.

Overexpression of miR-927-5p suppresses stalky expression and negatively reduces the spermatid production in Zeugodacus cucurbitae

BACKGROUND

The melon fly, Zeugodacus cucurbitae Coquillett, is one of the major pests attacking Cucurbitaceae crops. Identifying critical genes or proteins regulating fertility is essential for sustainable pest control and a research hotspot in insect physiology. MicroRNAs (miRNAs) are short RNAs that do not directly participate in protein translation, but instead function in post-transcriptional regulation of gene expression involved in male fertility.

RESULTS

We found that miR-927-5p is highly expressed in the testes and investigated its function in spermatogenesis in Z. cucurbitae. Fluorescence in situ hybridization (FISH) showed miR-927-5p in the transformation and maturation region of the testis, and overexpression of miR-927-5p reduced the number of sperms by 53%. In continuation, we predicted 12 target genes of miR-927-5p using bioinformatics combined with transcriptome sequencing data, and found that miR-927-5p targets the new gene Stalky in insects, which was validated by quantitative real-time PCR, RNA pull-down and dual luciferase reporter assays. FISH also confirmed the co-localization of miR-927-5p and the transcript Stalky_1 in the testis. Moreover, silencing of Stalky_1 by RNA interference reduced the number of sperms by 32% and reduced sperm viability by 39% in physiologically mature male adults. Meanwhile, the silencing of Stalky_1 also resulted in low hatchability.

CONCLUSION

Our work not only presents a new, so far unreported mechanism regulating spermatogenesis by miR-927-5p targeting a new unknown target, Stalky, which is providing new knowledge on the regulatory network of insect spermatogenesis, but also lays a foundation for the development of SIT against important tephritid fly pests. © 2024 Society of Chemical Industry.

Suppressing the expression of glutathione S-transferase gene GSTd10 increases the sensitivity of Zeugodacus cucurbitae against β-cypermethrin

Zeugodacus cucurbitae Coquillett (Diptera: Tephritidae) is an agriculturally and economically important pest worldwide that has developed resistance to β-cypermethrin. Glutathione S-transferases (GSTs) have been reported to be involved in the detoxification of insecticides in insects. We have found that both ZcGSTd6 and ZcGSTd10 were up-regulated by β-cypermethrin induction in our previous study, so we aimed to explore their potential relationship with β-cypermethrin tolerance in this study. The heterologous expression of ZcGSTd6 and ZcGSTd10 in Escherichia coli showed significantly high activities against 1-chloro-2,4-dinitrobenzene (CDNB). The kinetic parameters of ZcGSTd6 and ZcGSTd10 were determined by Lineweaver-Burk. The Vmax and Km of ZcGSTd6 were 0.50 μmol/min·mg and 0.3 mM, respectively. The Vmax and Km of ZcGSTd10 were 1.82 μmol/min·mg and 0.53 mM. The 3D modelling and molecular docking results revealed that β-cypermethrin exhibited a stronger bounding to the active site SER-9 of ZcGSTd10. The sensitivity to β-cypermethrin was significantly increased by 18.73% and 27.21%, respectively, after the knockdown of ZcGSTd6 and ZcGSTd10 by using RNA interference. In addition, the inhibition of CDNB at 50% (IC50) and the inhibition constants (Ki) of β-cypermethrin against ZcGSTd10 were determined as 0.41 and 0.33 mM, respectively. The Ki and IC50 of β-cypermethrin against ZcSGTd6 were not analysed. These results suggested that ZcGSTd10 could be an essential regulator involved in the tolerance of Z. cucurbitae to β-cypermethrin.

Two heat shock cognate 70 genes involved in spermatogenesis regulate the male fertility of Zeugodacus cucurbitae, as potential targets for pest control

The melon fly Zeugodacus cucurbitae Coquillett (Diptera: Tephritidae) is an agricultural quarantine pest threatening fruit and vegetable production. Heat shock cognate 70 (Hsc70), which is a homolog of the heat shock protein 70 (Hsp70), was first discovered in mice testes and plays an important role in spermatogenesis. In this study, we identified and cloned five Hsc70 genes from melon fly, namely ZcHsc70_1/2/3/4/5. Phylogenetic analysis showed that these proteins are closely related to Hsc70s from other Diptera insects. Spatiotemporal expression analysis showed that ZcHsc70_1 and ZcHsc70_2 are highly expressed in Z. cucurbitae testes. Fluorescence in situ hybridization further demonstrated that ZcHsc70_1 and ZcHsc70_2 are expressed in the transformation and maturation regions of testes, respectively. Moreover, RNA interference-based suppression of ZcHsc70_1 or ZcHsc70_2 resulted in a significant decrease of 74.61% and 63.28% in egg hatchability, respectively. Suppression of ZcHsc70_1 expression delayed the transformation of sperm cells to mature sperms. Meanwhile, suppression of ZcHsc70_2 expression decreased both sperm cells and mature sperms by inhibiting the meiosis of spermatocytes. Our findings show that ZcHsc70_1/2 regulates spermatogenesis and further affects the male fertility in the melon fly, showing potential as targets for pest control in sterile insect technique by genetic manipulation of males.

Comparative multi-tissue analyses identify testis-specific serine/threonine protein kinase (TSSK) genes involved in male fertility in the melon fly Zeugodacus cucurbitae

BACKGROUND

Zeugodacus cucurbitae is an agricultural pest species with robust reproductive capabilities capable of causing extensive damage. The advent of novel male fertility-related pest control strategies has been an area of active entomological research focused on the sterile insect technique (SIT) strategy.

RESULTS

RNA-sequencing analyses were conducted using 16 tissue samples from adult male Z. cucurbitae, leading to the identification of 5338 genes that were differentially expressed between the testes and three other analyzed tissue types. Of these genes, 808 exhibited high levels of testis expression. A quantitative polymerase chain reaction (qPCR) approach was used to validate the expression of ten of these genes selected at random, including ZcTSSK1 and ZcTSSK3, which are testis-specific serine/threonine protein kinase (TSSK) genes. Evaluation via a loss-of-function-based knockdown assay showed that the down-regulation of either of these two genes in males was associated with significantly decreased egg hatching rates. In situ hybridization analyses revealed the expression of both of these transcripts in the transformation zone, and significant decreases in Z. cucurbitae sperm numbers were observed following double-stranded RNA treatment. Together, these results suggested that inhibiting ZcTSSK1 and ZcTSSK3 expression was sufficient to alter male fertility in Z. cucurbitae.

CONCLUSION

These transcriptional sequencing results provide a foundation for further efforts to clarify the regulators of Z. cucurbitae male fertility. These preliminary analyses of the functions of ZcTSSK family genes as regulators of spermatogenesis underscore their importance in the processes integral to male fecundity and their potential as targets for pest control efforts centered on the genetic manipulation of males. © 2023 Society of Chemical Industry.

Genome-wide screening and expression of glutathione S-transferase genes reveal that GSTe4 contributes to sensitivity against β-cypermethrin in Zeugodacus cucurbitae

Glutathione S-transferases (GSTs) are an essential multifunctional protein family with common detoxifying enzymes. In this study, 34 GST genes were identified from the melon fly, Zeugodacus cucurbitae, one of the most destructive pests worldwide. These GSTs include 32 cytosolic genes and two microsomal genes. Furthermore, these cytosolic GSTs were classified into six classes: 11 delta, 13 epsilon, three theta, one sigma, two zeta, and two omega. Most of these showed dynamic expression during the developmental stage, some of which showed stage-specific expression. The expression in various adult tissues showed that most of them were expressed in anti-stress-related tissues. The transcriptional response of the delta and epsilon families was determined when Z. cucurbitae was exposed to three insecticides, abamectin, dinotefuran, and β-cypermethrin. Seven genes were significantly up-regulated by abamectin exposure. Moreover, five and four genes were significantly up-regulated with dinotefuran and β-cypermethrin exposure, respectively, demonstrating their involvement in the detoxification of these such toxic substances in Z. cucurbitae. One example of these genes, ZcGSTe4, was randomly selected to explore its function in response to β-cypermethrin exposure. Over-expressed ZcGSTe4 in E. coli showed significant tolerance to β-cypermethrin, and RNAi-mediated suppression of ZcGSTe4 also increased the sensitivity of melon fly to this agent. This study provides a foundation for further studies on the mechanism of detoxification metabolism in the melon fly.

A Unified Protocol for CRISPR/Cas9-Mediated Gene Knockout in Tephritid Fruit Flies Led to the Recreation of White Eye and White Puparium Phenotypes in the Melon Fly

Tephritid fruit flies are among the most invasive and destructive agricultural pests worldwide. Over recent years, many studies have implemented the CRISPR/Cas9 genome-editing technology to dissect gene functions in tephritids and create new strains to facilitate their genetics, management, and control. This growing literature allows us to compare diverse strategies for delivering CRISPR/Cas9 components into tephritid embryos, optimize procedures, and advance the technology to systems outside the most thoroughly studied species within the family. Here, we revisit five years of CRISPR research in Tephritidae and propose a unified protocol for candidate gene knockout in fruit flies using CRISPR/Cas9. We demonstrated the efficiency of our protocol by disrupting the eye pigmentation gene white eye (we) in the melon fly, Zeugodacus cucurbitae (Coquillett) (Diptera: Tephritidae). High rates of somatic and germline mutagenesis were induced by microinjecting pre-assembled Cas9-sgRNA complexes through the chorion of embryos at early embryogenesis, leading to the rapid development of new mutant lines. We achieved comparable results when targeting the we orthologue in the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), illustrating the reliability of our methods when transferred to other related species. Finally, we functionally validated the recently discovered white pupae (wp) loci in the melon fly, successfully recreating the white puparium phenotype used in suppression programs of this and other major economically important tephritids. This is the first demonstration of CRISPR-based genome-editing in the genus Zeugodacus, and we anticipate that the procedures described here will contribute to advancing genome-editing in other non-model tephritid fruit flies.

Recapitulation of the embryonic transcriptional program in holometabolous insect pupae

Holometabolous insects are predominantly motionless during metamorphosis, when no active feeding is observed and the body is enclosed in a hardened cuticle. These physiological properties as well as undergoing processes resemble embryogenesis, since at the pupal stage organs and systems of the imago are formed. Therefore, recapitulation of the embryonic expression program during metamorphosis could be hypothesized. To assess this hypothesis at the transcriptome level, we have performed a comprehensive analysis of the developmental datasets available in the public domain. Indeed, for most datasets, the pupal gene expression resembles the embryonic rather than the larval pattern, interrupting gradual changes in the transcriptome. Moreover, changes in the transcriptome profile during the pupa-to-imago transition are positively correlated with those at the embryo-to-larvae transition, suggesting that similar expression programs are activated. Gene sets that change their expression level during the larval stage and revert it to the embryonic-like state during the metamorphosis are enriched with genes associated with metabolism and development.

Vitelline Membrane Protein 26 Mutagenesis, Using CRISPR/Cas9, Results in Egg Collapse in Plutella xylostella

Vitelline membrane proteins (VMPs) are the main proteins that form the inner shell (vitelline membrane layer) of insect eggs and are an integral part of egg formation and embryo development. Here, we characterized the molecular structure and expression patterns of the VMP26 gene and analyzed its reproductive functions in diamondback moth, Plutella xylostella (L.), a worldwide migratory pest of cruciferous plants. The PxVMP26 gene was shown to be a single exon gene that contained an open reading frame of 852 base pairs (bp) encoding 283 amino acids. Both qPCR and western blot analyses showed that PxVMP26 was specifically expressed in female adults and was significantly highly expressed in the ovary. Further anatomical analysis indicated that the expression level of PxVMP26 in the ovarian tube with an incomplete yolk was significantly higher than that in the ovarian tube with a complete yolk. CRISPR/Cas9-induced PxVMP26 knockout successfully created two homozygous strains with 8- and 46-bp frameshift mutations. The expression deficiency of the PxVMP26 protein was detected in the mutant strains using immunofluorescence and western blot. No significant difference was found in the number of eggs laid within three days between wild and mutant individuals, but there was a lower egg hatchability. The loss of the PxVMP26 gene changed the mean egg size, damaged the structure of the vitelline membrane, and increased the proportion of abnormal eggs due to water loss, resulting in egg collapse. This first analysis of the roles of the VMP gene in the oocyte formation and embryonic development of P. xylostella, using CRISPR/Cas9 technology, provides a basis for screening new genetic control targets of P. xylostella.

Expression and Role of Vitellogenin Genes in Ovarian Development of Zeugodacus cucurbitae

Vitellogenin (Vg) genes encode the major egg yolk protein precursor in arthropods. In this study, four Vgs were identified in Zeugodacus cucurbitae (Coquillett). Sequence analysis showed that four ZcVgs had the conserved Vg domain. Phylogenetic analysis indicated that four ZcVgs were homologous to the Vgs of Tephritidae insects. The temporal and spatial expression patterns of ZcVgs were analyzed by quantitative real-time polymerase chain reaction (RT-qPCR), and the four ZcVgs showed high expression levels in female adults, especially in the fat body. The expression of ZcVg1 and ZcVg3 was down-regulated by a low dosage (0.5 μg) of 20-hydroxyecdysone (20E), and ZcVg2, ZcVg3, and ZcVg4 were up-regulated by a high dosage (1.0 and 2.0 μg) of 20E. The expression of ZcVg1 and ZcVg2 was up-regulated by 5 μg of juvenile hormone (JH), while all of the ZcVgs were down-regulated by a low and high dosage of JH. Expression of ZcVgs was down-regulated after 24 h of starvation and recovered to normal after nutritional supplementation. After micro-injection of the gene-specific double-stranded RNA, the ZcVgs’ expression was significantly suppressed, and ovarian development was delayed in Z. cucurbitae females. The results indicate that RNA interference of reproduction-related genes is a potential pest control method that works by manipulating female fertility.

lnc94638 is a testis-specific long non-coding RNA involved in spermatozoa formation in Zeugodacus cucurbitae (Coquillett)

Long non-coding RNAs (lncRNAs) generally display tissue-specific distributions, and testis-specific lncRNAs form the highest proportion of lncRNAs in many species. Here, we presented a detailed analysis of testis-specific lncRNAs in the melon fly, Zeugodacus cucurbitae, a highly destructive insect pest of cucurbitaceous and other related crops. Most testis-specific lncRNAs were found to be long intergenic non-coding RNAs (lincRNA). The size distribution of these lncRNAs ranged between 600 and 1000 nucleotides. Testis-specific lncRNAs that harboured one isoform number and two exons were the most abundant. Compared to other male tissues, the testis had more highly expressed lncRNAs. The quantitative real-time polymerase chain reaction results of 10 randomly selected testis-specific lncRNAs showed expression patterns consistent with RNA-seq data. Further analysis of the most highly expressed testis-specific lncRNA, lnc94638, was undertaken. Fluorescent in situ hybridization assays localized lnc94638 to the apical region of the testis that contains mature spermatozoa. RNA interference-mediated knockdown of lnc94638 expression reduced spermatozoa numbers and impaired the fertility of Z. cucurbitae male. This study provides a catalogue of testis-specific lncRNAs, shows that the testis-specific lnc94638 is involved in spermatogenesis and has the potential to be used for treating male sterility.

Time-Course Transcriptome Analysis Reveals Global Gene Expression Profiling and Dynamic Developmental Signatures across Complete Life Cycle of Bombyx mori

Background: The silkworm (Bombyx mori) is an important lepidopteran model insect worldwide which undergoes a complete metamorphosis developmental process. Although genome sequencing has been long performed, no transcriptome data covering the complete life cycle are available. Methods: Herein, a total of 10 samples were collected consecutively at four different developmental stages, including eggs of 24 h after oviposition (Ed) and eggs of 24 h after artificial egg-hatching (E); larvae from fist to fifth instar (L1–L5); early and late pupa (P4 and P8); and adult moth (M), were subjected to Illumina RNA-Seq and time-course analysis. Results: The summations of the gene expression of the silkworm ten developmental stages show: at Ed stage, eggs develop towards diapause status, the total gene expression level is relatively low; at E stage, after artificial egg-hatching, the expression level improves rapidly; during larval stages from L1–L5, the expression level rises gradually and reaches a peak at L5 stage; during pupae and moth stages, the total gene expression decline stage by stage. The results revealed a dynamical gene expression profile exhibiting significant differential expressions throughout the silkworm life cycle. Moreover, stage-specific key genes were identified at different developmental stages, suggesting their functions mainly characterized in maintaining insect development and immunity homeostasis or driving metamorphosis. GO annotation and KEGG enrichment analysis further revealed the most significantly enriched and fundamentally biological processes during silkworm growth. Conclusion: Collectively, our omics data depicted the first comprehensive landscape of dynamic transcriptome throughout complete developmental processes of B. mori. Our findings also provide valuable references and novel insights into understanding the molecular developmental remodeling events for other Lepidoptera species.

Identification and Expression Characterization of ATP-Binding Cassette (ABC) Transporter Genes in Melon Fly

Simple Summary

The melon fly, Zeugodacus cucurbitae, is an important agricultural pest. At present, chemical pesticide treatment is the main method for field control, but this promotes pesticide resistance by Z. cucurbitae, because of its frequent use. ABC transporters are involved in detoxification metabolism, but few studies have yet considered their expression in melon fly. In this study, we identified the ABC transporters genes at a genome-wide level in melon fly, and analysed their spatiotemporal expression patterns, as well as changes in expression after insecticides treatments. A total of 49 ABC transporters were identified, and their expression levels varied at different developmental stages and between tissues. After three insecticides treatment, ZcABCB7 and ZcABCC2 were up-regulated. After β-cypermethrin induction, tissues were dissected at 12, 24 and 48 h, and the expression levels of a number of ABC genes were highly expressed within the fat body. From these results, we conclude that ZcABCB7 and ZcABCC2 may be involved in detoxification metabolism, and that the fat body is the main tissue that plays this role.

Abstract

The ATP-binding cassette (ABC) transporter is a protein superfamily that transports specific substrate molecules across lipid membranes in all living species. In insects, ABC transporter is one of the major transmembrane protein families involved in the development of xenobiotic resistance. Here, we report 49 ABC transporter genes divided into eight subfamilies (ABCA-ABCH), including seven ABCAs, seven ABCBs, 10 ABCCs, two ABCDs, one ABCE, three ABCFs, 16 ABCGs, and three ABCHs according to phylogenetic analysis in Zeugodacus cucurbitae, a highly destructive insect pest of cucurbitaceous and other related crops. The expressions level of 49 ABC transporters throughout various developmental stages and within different tissues were evaluated by quantitative transcriptomic analysis, and their expressions in response to three different insecticides were evaluated by quantitative real-time polymerase chain reaction (qRT-PCR). These ABC transporter genes were widely expressed at developmental stages but most highly expressed in tissues of the midgut, fat body and Malpighian tube. When challenged by exposure to three insecticides, abamectin, β-cypermethrin, and dinotefuran, the expressions of ZcABCB7 and ZcABCC2 were significantly up-regulated. ZcABCB1, ZcABCB6, ZcABCB7, ZcABCC2, ZcABCC3, ZcABCC4, ZcABCC5, and ZcABCC7 were significantly up-regulated in the fat body at 24 h after β-cypermethrin exposure. These data suggest that ZcABCB7 and ZcABCC2 might play key roles in xenobiotic metabolism in Z. cucurbitae. Collectively, these data provide a foundation for further analysis of ABCs in Z. cucurbitae.

The PacBio Full-Length Transcriptome of the Tea Aphid as a Reference Resource

The tea aphid, Aphis aurantii, has become one of the destructive pests in tea plantations in the tropics and subtropics. Very few functional studies have so far focused on the developmental and reproductive biology at a molecular level, because of the lack of comprehensive genetic information. Full-length transcriptomes represent a very highly efficient approach to obtain reference gene sequences in non-model insects. In the present study, the transcriptome of A. aurantii was comprehensively sequenced using PacBio Iso-Seq technology. A total of 46.8 Gb nucleotides and 15,938 non-redundant full-length transcripts were obtained, 13,498 (84.69%) of which were annotated into seven databases. Of these transcripts, 2,029 alternative splicing events and 15,223 simple sequence repeats were detected. Among these transcripts, 4,571 (28.68%) and 11,367 (71.32%) were long non-coding RNAs (lncRNAs) and protein-coding genes, respectively. Five hundred and ninety transcription factors were detected. The first full-length transcriptome represents a significant increase in the known genetic information of A. aurantii. It will assist the future functional study of genes involved in its development and reproduction.

Genome-wide analysis of long non-coding RNAs in adult tissues of the melon fly, Zeugodacus cucurbitae (Coquillett)

Background

Long non-coding RNAs (lncRNAs) are involved in many fundamental biological processes, such as transcription regulation, protein degradation, and cell differentiation. Information on lncRNA in the melon fly, Zeugodacus cucurbitae (Coquillett) is currently limited.

Results

We constructed 24 RNA-seq libraries from eight tissues (midgut, Malpighian tubules, fat body, ovary, and testis) of Z. cucurbitae adults. A total of 3124 lncRNA transcripts were identified. Among those, 1464 were lincRNAs, 1037 were intronic lncRNAs, 301 were anti-sense lncRNAs, and 322 were sense lncRNAs. The majority of lncRNAs contained two exons and one isoform. Differentially expressed lncRNAs were analyzed between tissues, and Malpighian tubules versus testis had the largest number. Some lncRNAs exhibited strong tissue specificity. Specifically expressed lncRNAs were identified and filtered in tissues of female and male Z. cucurbitae based on their expression levels. Four midgut-specific lncRNAs were validated by quantitative real-time polymerase chain reaction (RT-qPCR), and the data were consistent with RNA-seq data. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of targets of midgut-specific lncRNAs indicated an enrichment of the metabolic process.

Conclusions

This was the first systematic identification of lncRNA in the melon fly. Expressions of lncRNAs in multiple adult tissues were evaluated by quantitative transcriptomic analysis. These qualitative and quantitative analyses of lncRNAs, especially the tissue-specific lncRNAs in Z. cucurbitae, provide useful data for further functional studies.