Genome sequence of the small brown planthopper, Laodelphax striatellus

Construction of cytomegalovirus promoter-driven gene expression system in Laodelphax striatellus

The small brown planthopper (SBPH, Laodelphax striatellus) is a significant rice pest, responsible for transmitting rice stripe virus (RSV) in a persistent and propagative manner. RSV is one of the most detrimental rice viruses, causing rice stripe disease, which results in considerable loss of rice grain yield. While RNA interference and gene knockout techniques have enabled gene downregulation in SBPH, no system currently exists for the overexpression of endogenous or exogenous genes. Consequently, the development of a protein expression system for SBPH is imperative to serve as a technical foundation for pest control and gene function investigations. This study aimed to construct an expression vector using the promoter of the constitutive-expressed tubulin gene of SBPH, and promoter of human cytomegalovirus (CMV). Fluorescence experiments demonstrated that both tubulin and CMV promoter could drive green fluorescent protein (GFP) expression in SBPH, and could also facilitate the expression of a nucleocapsid protein (NP) -GFP fusion protein containing viral NP with comparable efficiency. Through expression vector optimization, we have identified that the 3 tandem CMV promoters display a significantly higher promoter activity compared with both the 2 tandem CMV promoters and the single CMV promoter. In addition, the incorporation of Star polycation nanoparticles significantly enhanced the expression efficiency in SBPH. These results provide a promising technical platform for investigating gene functions in SBPH.

New genes helped acorn barnacles adapt to a sessile lifestyle

Barnacles are the only sessile lineages among crustaceans, and their sessile life begins with the settlement of swimming larvae (cyprids) and the formation of protective shells. These processes are crucial for adaptation to a sessile lifestyle, but the underlying molecular mechanisms remain poorly understood. While investigating these mechanisms in the acorn barnacle, Amphibalanus amphitrite, we discovered a new gene, bcs-6, which is involved in the energy metabolism of cyprid settlement and originated from a transposon by acquiring the promoter and cis-regulatory element. Unlike mollusks, the barnacle shell comprises alternate layers of chitin and calcite and requires another new gene, bsf, which generates silk-like fibers that efficiently bind chitin and aggregate calcite in the aquatic environment. Our findings highlight the importance of exploring new genes in unique adaptative scenarios, and the results will provide important insights into gene origin and material development.

Genomic and Transcriptomic Insights into the Genetic Basis of Foam Secretion in Rice Spittlebug Callitettix versicolor

Many animal species produce protective foams, the majority of which exhibit evolutionary adaptability. Although the function and composition of foams have been widely studied, the genetic basis of foam secretion remains unknown. Unlike most species that produce foam under specific situations, spittlebugs continuously secrete foams throughout all nymphal stages. Here, we capitalize on the rice spittlebug (Callitettix versicolor) to explore the genetic basis of foam secretion through genomic and transcriptomic approaches. Our comparative genomic analysis for C. versicolor and eight other insect species reveals 606 species-specific gene families and 66 expanded gene families, associated with carbohydrate and lipid metabolism. These functions are in accordance with the composition of foams secreted by spittlebugs. Transcriptomic analyses of malpighian tubules across developmental stages detected 3192 differentially expressed genes. Enrichment analysis of these genes highlights functions also revealed by our comparative genomic analysis and aligns with previous histochemical and morphological observations of foam secretion. This consistency suggests the important roles of these candidate genes in foam production. Our study not only provides novel insights into the genetic basis of foam secretion in rice spittlebugs but also contributes valuable knowledge for future evolutionary studies of spittlebugs and the development of pest control strategies for C. versicolor.

High-quality genome of the zoophytophagous stink bug, Nesidiocoris tenuis, informs their food habit adaptation

The zoophytophagous stink bug, Nesidiocoris tenuis, is a promising natural enemy of micro-pests such as whiteflies and thrips. This bug possesses both phytophagous and entomophagous food habits, enabling it to obtain nutrition from both plants and insects. This trait allows us to maintain its population density in agricultural fields by introducing insectary plants, even when the pest prey density is extremely low. However, if the bugs' population becomes too dense, they can sometimes damage crop plants. This dual character seems to arise from the food preferences and chemosensation of this predator. To understand the genomic landscape of N. tenuis, we examined the whole genome sequence of a commercially available Japanese strain. We used long-read sequencing and Hi-C analysis to assemble the genome at the chromosomal level. We then conducted a comparative analysis of the genome with previously reported genomes of phytophagous and hematophagous stink bugs to focus on the genetic factors contributing to this species' herbivorous and carnivorous tendencies. Our findings suggest that the gustatory gene set plays a pivotal role in adapting to food habits, making it a promising target for selective breeding. Furthermore, we identified the whole genomes of microorganisms symbiotic with this species through genomic analysis. We believe that our results shed light on the food habit adaptations of N. tenuis and will accelerate breeding efforts based on new breeding techniques for natural enemy insects, including genomics and genome editing.

Chromosome-level genome assembly of the Asian citrus psyllid, Diaphorina citri

Diaphorina citri is a global citrus pest. As a vector insect, it can transmit the causative agents of citrus huanglongbing, causing irreversible losses to the citrus industry. The acquisition of genomic information can provide a molecular genetic basis for effective control of D. citri. Here, the DNBSEQ™ , Oxford Nanopore Technologies, and Hi-C technologies are applied to generate a high-quality chromosome-level genome of D. citri. The genome size of D. citri was 523.78 Mb with a scaffold N50 of 47.05 Mb distributed on 13 chromosomes. A total of 250.64 Mb (47.85%) repeat sequences and 24 048 protein-coding genes were predicted. Genome resequencing of female and male individuals indicated that the sex chromosome system of D. citri is XO. Phylogenetic analysis demonstrated that D. citri and Pachypsylla venusta, which separated from their most recent common ancestor about 336.62 million years ago, were the most closely related. Additionally, we identified genes potentially involved in detoxification metabolism, pathogen transmission, and honeydew secretion for further investigation. The high-quality genome provides an important reference for developing effective management strategies of D. citri.

A plant virus manipulates the long-winged morph of insect vectors

Wing dimorphism of insect vectors is a determining factor for viral long-distance dispersal and large-area epidemics. Although plant viruses affect the wing plasticity of insect vectors, the potential underlying molecular mechanisms have seldom been investigated. Here, we found that a planthopper-vectored rice virus, rice stripe virus (RSV), specifically induces a long-winged morph in male insects. The analysis of field populations demonstrated that the long-winged ratios of male insects are closely associated with RSV infection regardless of viral titers. A planthopper-specific and testis-highly expressed gene, Encounter, was fortuitously found to play a key role in the RSV-induced long-winged morph. Encounter resembles malate dehydrogenase in the sequence, but it does not have corresponding enzymatic activity. Encounter is upregulated to affect male wing dimorphism at early larval stages. Encounter is closely connected with the insulin/insulin-like growth factor signaling pathway as a downstream factor of Akt, of which the transcriptional level is activated in response to RSV infection, resulting in the elevated expression of Encounter. In addition, an RSV-derived small interfering RNA directly targets Encounter to enhance its expression. Our study reveals an unreported mechanism underlying the direct regulation by a plant virus of wing dimorphism in its insect vectors, providing the potential way for interrupting viral dispersal.

NDUFA8 potentially rescues Wolbachia-induced cytoplasmic incompatibility in Laodelphax striatellus

The endosymbiont Wolbachia manipulates host reproduction by several strategies, one of the most important of which is cytoplasmic incompatibility (CI). CI can be rescued when Wolbachia-infected males mate with females infected with the same Wolbachia strain. However, the potential rescue mechanism of CI in the small brown planthopper Laodelphax striatellus is unclear. In this study, comparative transcriptome analysis was applied to explore the effect of Wolbachia on L. striatellus eggs. A total of 1387 differentially expressed genes were identified. RNA interference of 7 Wolbachia-upregulated key planthopper genes reduced egg reproduction, suggesting that Wolbachia might improve fecundity in L. striatellus by affecting these 7 genes. Suppressing the expression of another upregulated gene, NDUFA8 (encoding NADH dehydrogenase [ubiquinone] 1 α subcomplex subunit 8-like) by RNA interference significantly increased the mortality of early embryos without affecting the number of deposited eggs. Wolbachia infection upregulated the mRNA level of NDUFA8, and dsNDUFA8 treatment of Wolbachia-infected females recreated CI-like symptoms, suggesting that NDUFA8 is associated with the rescue phenotype. Because all L. striatellus populations worldwide are infected with Wolbachia, NDUFA8 is a potential pest control target.

Chromosome-level genome assembly and population genomic analysis provide novel insights into the immunity and evolution of Sogatella furcifera

Sogatella furcifera is a destructive agricultural pest causing large threats to rice production in China and Southeast Asian countries. Despite recent breakthroughs in long-read sequencing, high quality genomic data are very limited in S. furcifera. In present study, a chromosome-level assembly of the S. furcifera genome was completed (0.64 GB), comprising 15 chromosomes covered 95.04% of the estimated genome size, along with other 624 small scaffolds making up the remaining 4.96% of the genome of S. furcifera. A total of 24,669 protein-coding genes, 1211 long noncoding RNA and 7595 circular RNA transcripts were predicted in this study. Comparative genomic analysis revealed rapidly evolved genes were associated with multiple immune-related pathways in S. furcifera. Genome resequencing of 44 individuals from 12 geographic populations revealed frequent gene flow among populations. The systemic genomic analysis will provide more insights into the understanding of the immunity and evolutionary adaptation of S. furcifera.

Parallel Losses of Blue Opsin Correlate with Compensatory Neofunctionalization of UV-Opsin Gene Duplicates in Aphids and Planthoppers

Expanding on previous efforts to survey the visual opsin repertoires of the Hemiptera, this study confirms that homologs of the UV- and LW-opsin subfamilies are conserved in all Hemiptera, while the B-opsin subfamily is missing from the Heteroptera and subgroups of the Sternorrhyncha and Auchenorrhyncha, i.e., aphids (Aphidoidea) and planthoppers (Fulgoroidea), respectively. Unlike in the Heteroptera, which are characterized by multiple independent expansions of the LW-opsin subfamily, the lack of B-opsin correlates with the presence of tandem-duplicated UV-opsins in aphids and planthoppers. Available data on organismal wavelength sensitivities and retinal gene expression patterns lead to the conclusion that, in both groups, one UV-opsin paralog shifted from ancestral UV peak sensitivity to derived blue sensitivity, likely compensating for the lost B-opsin. Two parallel bona fide tuning site substitutions compare to 18 non-corresponding amino acid replacements in the blue-shifted UV-opsin paralogs of aphids and planthoppers. Most notably, while the aphid blue-shifted UV-opsin clade is characterized by a replacement substitution at one of the best-documented UV/blue tuning sites (Rhodopsin site 90), the planthopper blue-shifted UV-opsin paralogs retained the ancestral lysine at this position. Combined, the new findings identify aphid and planthopper UV-opsins as a new valuable data sample for studying adaptive opsin evolution.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

c-Jun regulates flotillin 2 transcription to benefit viral accumulation in insect vectors

The c-Jun N-terminal kinase (JNK) signaling pathway plays a critical role in viral infection in host cells. In addition to triggering immune reactions against pathogens, the JNK signaling pathway has also been found to benefit viral infection. Our previous work showed that JNK activation facilitated rice stripe virus (RSV) accumulation in the insect vector small brown planthopper, but the underlying mechanisms remain elusive. Here, we revealed a link between JNK activation and the transcriptional upregulation of the plasma membrane protein flotillin 2, which mediates RSV cell entry. c-Jun, a downstream substrate of JNKs, was identified as a transcription factor that targets the promoter of flotillin 2 at three binding sites. Phosphorylated c-Jun, especially at the serine 63 site, promoted the transcriptional activity of c-Jun on flotillin 2. JNK activation or inhibition affected c-Jun phosphorylation status and flotillin 2 expression. In the midguts of planthoppers, RSV infection significantly increased flotillin 2 expression and the phosphorylation level of JNKs and c-Jun. Manipulation of JNK status impacted viral acquisition in midgut cells. These findings reveal a new regulatory mechanism of the JNK signaling pathway and shed light on the virus-supportive effect of this pathway.

Chromosome-level Assembly, Dosage Compensation and Sex-biased Gene Expression in the Small Brown Planthopper, Laodelphax striatellus

In insects, sex chromosome differentiation often results in unequal gene dosages between sexes. Dosage compensation mechanisms evolve to balance gene expression, but the degree and mechanism of regulation often vary by insect species. In hemipteran species, the small brown planthopper (SBPH), Laodelphax striatellus, is an injurious crop pest, with a sex chromosome type XX in females and XO in males. This species offers the opportunity to study dosage compensation and sex-biased gene expression. In this study, we generated a chromosome-level genome of SBPH using Oxford Nanopore Technologies and high-throughput chromatin conformation capture (Hi-C) technology. We also sequenced RNA-seq data from 16 tissue samples to annotate the genome and analyze gene dosage compensation. We finally obtained a 510.2 megabases (Mb) genome with 99.12% of the scaffolds anchored on 15 chromosomes (14 autosomes and 1 X chromosome) and annotated 16,160 protein-coding genes based on full-length cDNA sequencing data. Furthermore, we found complete dosage compensation in all L. striatellus somatic tissues, but lack of dosage compensation in gonad tissue testis. We also found that female-biased genes were significantly enriched on the X chromosome in all tissues, whereas male-biased genes in gonad tissues were enriched on autosomes. This study not only provides a high-quality genome assembly but also lays a foundation for a better understanding of the sexual regulatory network in hemipteran insects.

Key role of exportin 6 in exosome-mediated viral transmission from insect vectors to plants

Exosomes play a key role in virus exocytosis and transmission. The exportin family is usually responsible for cargo nucleocytoplasmic trafficking, and they are frequently found in exosomes. However, the function of exportins sorted in exosomes remains unknown. Here, we successfully isolated "cup holder"-like exosomes from the saliva of ∼30,000 small brown planthoppers, which are vectors of rice stripe virus (RSV). RSV virions were packed in comparatively large exosomes. Four viral genomic RNAs at a certain ratio were identified in the saliva exosomes. The virions contained in the saliva exosomes were capable of replicating and causing disease in rice plants. Interference with each phase of the insect exosome system affected the transmission of RSV from the insect vectors to rice plants. Fragmented exportin 6 was coimmunoprecipitated with viral nucleocapsid protein in saliva and sorted to exosomes via interactions with the cargo sorting protein VPS37a. When the expression of exportin 6 was knocked down, the amounts of RSV secreted in saliva and rice plants were reduced by 60% and 74%, respectively. These results showed that exportin 6 acted as a vehicle for transporting RSV into exosomes to overcome the barrier of insect salivary glands for horizontal transmission. Exportin 6 would represent an ideal target that could be manipulated to control the outbreak of insect-borne viruses in the future.

Genome-Wide Identification and Expression Profiling of the Wnt Gene Family in Three Rice Planthoppers: Sogatella furcifera, Laodelphax striatellus, and Nilaparvata lugens

The Wnt gene family plays essential roles in regulating many developmental processes, including the maintenance of stem cells, cell division, and cell migration. The number of Wnt genes varies among species. Due to the diversity and importance of their functions, the Wnt gene family has gained extensive research interest in various animal species from invertebrates to vertebrates. However, knowledge of the Wnt gene family is limited in rice planthoppers. Three planthopper species, the white-backed planthopper (Sogatella furcifera Horvath), the small brown planthopper (Laodelphax striatellus Fallén) and the brown planthopper (Nilaparvata lugens Stål) (Hemiptera: Delphacidae), are devastating specialist pests of rice and cause serious damage to rice plants. To better study the evolution and function of the Wnt gene family in rice planthoppers, we identified 8 Wnt family genes in three rice planthoppers with both genomic and extensive transcriptomic resources available. We conducted a systematic analysis of the three kinds of rice planthoppers and analyzed the dynamic patterns of gene conservation, as well as Wnt gene loss and duplication. The expression profiles in different developmental stages of S. furcifera and different adult organs and tissues of L. striatellus provide preliminary functional implications for the Wnt genes in rice planthopper. This study presents the first genome-wide study of the Wnt gene family in rice planthoppers, and our findings provide insights into Wnt function and evolution in rice planthoppers.

Copper Chaperone for Superoxide Dismutase FoCCS1 in Frankliniella occidentalis May Be Associated with Feeding Adaptation after Host Shifting

Western flower thrips (Frankliniella occidentalis) pose a serious threat to the global vegetable and flower crop production. The regulatory mechanism for superoxide dismutase (SOD) in the feeding adaptation of F. occidentalis after host shifting remains unclear. In this study, the copper chaperone for SOD (CCS) and manganese SOD (MnSOD) genes in F. occidentalis were cloned, and their expression levels at different developmental stages was determined. The mRNA expression of FoCCS1 and FoMnSOD2 in F. occidentalis second-instar larvae and adult females of F₁, F₂, and F₃ generations was analyzed after shifting the thrips to kidney bean and broad bean plants, respectively. The F₂ and F₃ second-instar larvae and F₂ adult females showed significantly upregulated FoCCS1 mRNA expression after shifting to kidney bean plants. The F₁ second-instar larvae and F₂ adult females showed significantly upregulated FoCCS1 mRNA expression after shifting to broad bean plants. The RNA interference significantly downregulated the FoCCS1 mRNA expression levels and adult females showed significantly inhibited SOD activity after shifting to kidney bean and broad bean plants. F. occidentalis adult females subjected to RNA interference and released on kidney bean and broad bean leaves for rearing, respectively, significantly reduced the survival rate and fecundity. These findings suggest that FoCCS1 plays an active role in regulating the feeding adaptation ability of F. occidentalis after host shifting.

Laodelphax striatellus saliva mucin enables the formation of stylet sheathes to facilitate its feeding and rice stripe virus transmission

BACKGROUND

Laodelphax striatellus transmits rice stripe virus (RSV) during sap feeding on the rice plant. The insect saliva proteins have direct and indirect roles in mediating RSV transmission; however, the function of most saliva proteins remains unclear.

RESULTS

In this study, we sequenced L. striatellus saliva proteins using shotgun liquid chromatography-electrospray ionization-tandem mass spectrometry. We identified 41 secreted saliva proteins, among which a saliva mucin-like protein, designated LssaMP, was the most abundant. In silico analysis revealed the sequence conservation among planthoppers. We revealed that the LssaMP gene is specifically expressed in the salivary glands and the protein is secreted as a component of gel saliva. Using LssaMP-specific double-stranded RNA (dsRNA) to silence gene expression, we revealed that LssaMP is required for formation of the salivary sheath, an important structure for sap feeding. Disrupting LssaMP expression resulted in inefficient formation of the feeding structure, thereby stopping insects from secreting watery saliva and acquiring sufficient nutrients from the phloem sap. We confirmed that RSV is mainly released via the watery saliva, which passes through the salivary sheathes into the plant phloem. An insufficient feeding structure results in decreased release of watery saliva, as well as the arboviruses.

CONCLUSION

This study clarified the function of an insect saliva protein in mediating insect feeding, as well as arbovirus transmission. © 2022 Society of Chemical Industry.

The Adipokinetic Peptides of Hemiptera: Structure, Function, and Evolutionary Trends

The Hemiptera comprise the most species-rich order of the hemimetabolous insects. Members of a number of superfamilies, most notably especially the more basal ones such as white flies, psyllids and aphids, belong to the most destructive agricultural insects known worldwide. At the other end of the phylogenetic tree are hemipterans that are notorious medical pests (e.g. kissing bugs). Most of the hemipteran species are good flyers, and lipid oxidation plays a pivotal role to power the contraction of flight muscles and, in aquatic water bugs, also deliver the ATP for the extensive swimming action of the leg muscles. Mobilization of stored lipids (mostly triacylglycerols in the fat body) to circulating diacylglycerols in the hemolymph is regulated by a set of small neuropeptides, the adipokinetic hormones (AKHs). We searched the literature and publicly available databases of transcriptomes and genomes to present here AKH sequences from 191 hemipteran species. Only few of these peptides were sequenced via Edman degradation or mass spectrometry, and even fewer were characterized with molecular biology methods; thus, the majority of the AKHs we have identified by bioinformatics are merely predicted sequences at this stage. Nonetheless, a total of 42 AKH primary sequences are assigned to Hemiptera. About 50% of these structures occur also in other insect orders, while the remaining 50% are currently unique for Hemiptera. We find 9 novel AKHs not shown to be synthesized before in any insect. Most of the hemipteran AKHs are octapeptides (28) but there is an impressive number of decapeptides (12) compared to other speciose orders such as Diptera and Lepidoptera. We attempt to construct a hypothetical molecular peptide evolution of hemipteran AKHs and find quite a bit of overlapping with current phylogenetic ideas of the Hemiptera. Lastly, we discuss the possibility to use the sequence of the aphid AKH as lead peptide for the research into a peptide mimetic fulfilling criteria of a green insecticide.

ATP-binding cassette transporters ABCF2 and ABCG9 regulate rice black-streaked dwarf virus infection in its insect vector, Laodelphax striatellus (Fallén)

The majority of plant viral disease is transmitted and spread by insect vectors in the field. The small brown planthopper, Laodelphax striatellus (Fallén), is the only efficient vector for rice black-streaked dwarf virus (RBSDV), a devastating plant virus that infects multiple grain crops, including rice, maize, and wheat. Adenosine triphosphate (ATP)-binding cassette (ABC) transporters participate in various biological processes. However, little is known about whether ABC transporters affect virus infection in insects. In this study, RBSDV accumulation was significantly reduced in L. striatellus after treatment with verapamil, an effective inhibitor of ABC transporters. Thirty-four ABC transporter genes were identified in L. striatellus and expression analysis showed that LsABCF2 and LsABCG9 were significantly upregulated and downregulated, respectively, after RBSDV infection. LsABCF2 and LsABCG9 were expressed during all developmental stages, and LsABCG9 was highly expressed in the midgut of L. striatellus. Knockdown of LsABCF2 promoted RBSDV accumulation, while knockdown of LsABCG9 suppressed RBSDV accumulation in L. striatellus. Our data showed that L. striatellus might upregulate the expression of LsABCF2 and downregulate LsABCG9 expression to suppress RBSDV infection. These results will contribute to understanding the effects of ABC transporters on virus transmission and provide theoretical basis for virus management in the field.

GrpE is involved in mitochondrial function and is an effective target for RNAi-mediated pest and arbovirus control

Laodelphax striatellus is a sap-feeding pest and the main insect vector of rice stripe virus (RSV). There is an urgent need to identify molecular targets to control this insect pest and plant arboviruses. In this study, we identified a L. striatellus gene (named LsGrpE) encoding a GroP-E-like protein. We found that the LsGrpE protein localized to mitochondria. Using gene-specific dsRNA to interfere with the expression of LsGrpE led to a significant increase in insect mortality, and most of the surviving insects could not develop into adults. Further analyses revealed that LsGrpE deficiency caused mitochondrial dysfunction and inhibited the insulin pathway, resulting in diabetes-like symptoms such as elevated blood sugar, inactive behaviour, developmental delay, and death. In addition, LsGrpE deficiency significantly reduced the RSV titre in surviving L. striatellus, and indirectly prevented viral vertical transmission by reducing the number of adults. We generated transgenic rice plants expressing LsGrpE-specific dsRNA, and the dsRNA was acquired by L. striatellus during feeding, resulting in increased insect mortality and the prevention of arboviral transmission. This study clarifies the function of LsGrpE and demonstrates that LsGrpE can be used as a molecular target of plant-generated dsRNA to resist this sap-feeding pest, a17nd therefore, its transmitted arboviruses.

Genome-wide SNPs reveal the fine-scale population structure of Laodelphax striatellus in China using double-digest restriction site-associated DNA sequencing

The small brown planthopper (SBPH), Laodelphax striatellus (Fallén) is one of the most destructive rice pests and has caused serious economic losses in China. To clarify the genetic differentiation and population genetic structure of this insect pest, we investigated the genomic polymorphisms, genetic differentiation, and phylogeography of 31 SBPH populations from 28 sampling sites from three climatic zones of China using double-digest restriction site-associated DNA sequencing (ddRADseq). In total, 2,813,221,369 high-quality paired-end reads from 306 individuals and 1925 single nucleotide polymorphisms (SNPs) were obtained. Low levels of genetic diversity and significant genetic differentiation were observed among the SBPH populations, and three genetic clusters were detected in China. Neutrality tests and bottleneck analysis provided strong evidence for recent rapid expansion with a severe bottleneck in most populations. Our work provides new insights into the genetics of the SBPH and will contribute to the development of effective management strategies for this pest.

Metagenome-assembled genomes of phytoplankton microbiomes from the Arctic and Atlantic Oceans

BACKGROUND

Phytoplankton communities significantly contribute to global biogeochemical cycles of elements and underpin marine food webs. Although their uncultured genomic diversity has been estimated by planetary-scale metagenome sequencing and subsequent reconstruction of metagenome-assembled genomes (MAGs), this approach has yet to be applied for complex phytoplankton microbiomes from polar and non-polar oceans consisting of microbial eukaryotes and their associated prokaryotes.

RESULTS

Here, we have assembled MAGs from chlorophyll a maximum layers in the surface of the Arctic and Atlantic Oceans enriched for species associations (microbiomes) with a focus on pico- and nanophytoplankton and their associated heterotrophic prokaryotes. From 679 Gbp and estimated 50 million genes in total, we recovered 143 MAGs of medium to high quality. Although there was a strict demarcation between Arctic and Atlantic MAGs, adjacent sampling stations in each ocean had 51-88% MAGs in common with most species associations between Prasinophytes and Proteobacteria. Phylogenetic placement revealed eukaryotic MAGs to be more diverse in the Arctic whereas prokaryotic MAGs were more diverse in the Atlantic Ocean. Approximately 70% of protein families were shared between Arctic and Atlantic MAGs for both prokaryotes and eukaryotes. However, eukaryotic MAGs had more protein families unique to the Arctic whereas prokaryotic MAGs had more families unique to the Atlantic.

CONCLUSION

Our study provides a genomic context to complex phytoplankton microbiomes to reveal that their community structure was likely driven by significant differences in environmental conditions between the polar Arctic and warm surface waters of the tropical and subtropical Atlantic Ocean. Video Abstract.

Comprehensive evaluation of structural variant genotyping methods based on long-read sequencing data

Background

Structural variants (SVs) play a crucial role in gene regulation, trait association, and disease in humans. SV genotyping has been extensively applied in genomics research and clinical diagnosis. Although a growing number of SV genotyping methods for long reads have been developed, a comprehensive performance assessment of these methods has yet to be done.

Results

Based on one simulated and three real SV datasets, we performed an in-depth evaluation of five SV genotyping methods, including cuteSV, LRcaller, Sniffles, SVJedi, and VaPoR. The results show that for insertions and deletions, cuteSV and LRcaller have similar F1 scores (cuteSV, insertions: 0.69–0.90, deletions: 0.77–0.90 and LRcaller, insertions: 0.67–0.87, deletions: 0.74–0.91) and are superior to other methods. For duplications, inversions, and translocations, LRcaller yields the most accurate genotyping results (0.84, 0.68, and 0.47, respectively). When genotyping SVs located in tandem repeat region or with imprecise breakpoints, cuteSV (insertions and deletions) and LRcaller (duplications, inversions, and translocations) are better than other methods. In addition, we observed a decrease in F1 scores when the SV size increased. Finally, our analyses suggest that the F1 scores of these methods reach the point of diminishing returns at 20× depth of coverage.

Conclusions

We present an in-depth benchmark study of long-read SV genotyping methods. Our results highlight the advantages and disadvantages of each genotyping method, which provide practical guidance for optimal application selection and prospective directions for tool improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08548-y.

Whole genome sequencing of spotted stem borer, Chilo partellus, reveals multiple genes encoding enzymes for detoxification of insecticides

Spotted stem borer, Chilo partellus, is the most important constraint for increasing the production and productivity of maize and sorghum, the two major coarse cereals in Asia and Africa. The levels of resistance to this pest in the cultivated germplasm are low to moderate, and hence, farmers have to use insecticides for effective control of this pest. However, there is no information on the detoxification mechanisms in C. partellus, which is one of the constraints for deployment of appropriate insecticides to control this pest. The ability to detoxify insecticides varies across insect populations, and hence, we sequenced different populations of C. partellus to identify and understand detoxification mechanisms to devise appropriate strategies for deployment of different insecticides for controlling this pest. Larval samples were sequenced from three different cohorts of C. partellus using the Illumina HiSeq 2500 platform. The data were subjected to identify putative genes that are involved in detoxification on insecticides in our cohort insect species. These studies resulted in identification of 64 cytochrome P450 genes (CYP450s), and 36 glutathione S-transferases genes (GSTs) encoding metabolic detoxification enzymes, primarily responsible for xenobiotic metabolism in insects. A total of 183 circadian genes with > 80% homolog and 11 olfactory receptor genes that mediate chemical cues were found in the C. partellus genome. Also, target receptors related to insecticide action, 4 acetylcholinesterase (AChE), 14 γ-aminobutyric acid (GABA), and 15 nicotinic acetylcholine (nAChR) receptors were detected. This is the first report of whole genome sequencing of C. partellus useful for understanding mode of action of different insecticides, and mechanisms of detoxification and designing target-specific insecticides to develop appropriate strategies to control C. partellus for sustainable crop production.

Integrating whole-genome sequencing with multi-omic data reveals the impact of structural variants on gene regulation in the human brain

Structural variants (SVs), genomic rearrangements of >50 bp, are an important source of genetic diversity and have been linked to many diseases. However, it remains unclear how they modulate human brain function and disease risk. Here, we report 170,996 SVs discovered using 1,760 short-read whole genomes from aged adults and Alzheimer’s disease individuals. By applying quantitative trait locus (SV-xQTL) analyses, we quantified the impact of cis-acting SVs on histone modifications, gene expression, splicing, and protein abundance in post-mortem brain tissues. More than 3,200 SVs were associated with at least one molecular phenotype. We found reproducibility of 65–99% SV-eQTLs across cohorts and brain regions. SV associations with mRNA and proteins shared the same direction of effect in more than 87% of SV-gene pairs. Mediation analysis showed ~8% of SV-eQTLs mediated by histone acetylation, and ~11% by splicing. Additionally, associations of SVs with progressive supranuclear palsy identified previously known and novel SVs.

Flotillin 2 Facilitates the Infection of a Plant Virus in the Gut of Insect Vector

Most plant viruses require insect vectors for transmission. One of the key steps for the transmission of persistent-circulative plant viruses is overcoming the gut barrier to enter epithelial cells. To date, little has been known about viral cofactors in gut epithelial cells of insect vectors. Here, we identified flotillin 2 as a plasma membrane protein that facilitates the infection of rice stripe virus (RSV) in its vector, the small brown planthopper. Flotillin 2 displayed a prominent plasma membrane location in midgut epithelial cells. The nucleocapsid protein of RSV and flotillin 2 colocalized on gut microvilli, and a nanomolar affinity existed between the two proteins. Knockout of flotillin 2 impeded the entry of virions into epithelial cells, resulting in a 57% reduction of RSV levels in planthoppers. The knockout of flotillin 2 decreased disease incidence in rice plants fed by viruliferous planthoppers from 40% to 11.7%. Furthermore, flotillin 2 mediated the infection of southern rice black-streaked dwarf virus in its vector, the white-backed planthopper. This work implies the potential of flotillin 2 as a target for controlling the transmission of rice stripe disease.

IMPORTANCE Plant viral diseases are a major threat to world agriculture. The transmission of 80% of plant viruses requires vector insects, and 54% of vector-borne plant viruses are persistent-circulative viruses, which must overcome the barriers of gut cells with the help of proteins on the cell surface. Here, we identified flotillin 2 as a membrane protein that mediates the cell entry of rice stripe virus in its vector insect, small brown planthopper. Flotillin 2 displays a prominent cellular membrane location in midgut cells and can specifically bind to virions. The loss of flotillin 2 impedes the entry of virions into the midgut cells of vector insects and substantially suppresses viral transmission to rice. Therefore, flotillin 2 may be a promising target gene for manipulation in vector insects to control the transmission of rice stripe disease and perhaps that of other rice virus diseases in the future.

A Review on Transcriptional Responses of Interactions between Insect Vectors and Plant Viruses

This review provides a synopsis of transcriptional responses pertaining to interactions between plant viruses and the insect vectors that transmit them in diverse modes. In the process, it attempts to catalog differential gene expression pertinent to virus–vector interactions in vectors such as virus reception, virus cell entry, virus tissue tropism, virus multiplication, and vector immune responses. Whiteflies, leafhoppers, planthoppers, and thrips are the main insect groups reviewed, along with aphids and leaf beetles. Much of the focus on gene expression pertinent to vector–virus interactions has centered around whole-body RNA extraction, whereas data on virus-induced tissue-specific gene expression in vectors is limited. This review compares transcriptional responses in different insect groups following the acquisition of non-persistent, semi-persistent, and persistent (non-propagative and propagative) plant viruses and identifies parallels and divergences in gene expression patterns. Understanding virus-induced changes in vectors at a transcriptional level can aid in the identification of candidate genes for targeting with RNAi and/or CRISPR editing in insect vectors for management approaches.

Population genomics provides insights into lineage divergence and local adaptation within the cotton bollworm

The cotton bollworm Helicoverpa armigera is a cosmopolitan pest and its diverse habitats plausibly contribute to the formation of diverse lineages. Despite the significant threat it poses to economic crops worldwide, its evolutionary history and genetic basis of local adaptation are poorly understood. In this study, we de novo assembled a high-quality chromosome-level reference genome of H. a. armigera (contig N50 = 7.34 Mb), with 99.13% of the HaSCD2 assembly assigned into 31 chromosomes (Z-chromosome + 30 autosomes). We constructed an ultra-dense variation map across 14 cotton bollworm populations and identified a novel lineage in northwestern China. Historical inference showed that effective population size changes coincided with global temperature fluctuation. We identified nine differentiated genes in the three H. armigera lineages (H. a. armigera, H. a. conferta, and the new northwestern Chinese lineage), of which per and clk genes are involved in circadian rhythm. Selective sweep analyses identified a series of GO categories related to climate adaptation, feeding behavior and insecticide tolerance. Our findings reveal fundamental knowledge of the local adaptation of different cotton bollworm lineages and will guide the formulation of cotton bollworm management measures at different scales.

A cytorhabdovirus-based expression vector in Nilaparvata lugens, Laodelphax striatellus, and Sogatella furcifera

The brown planthopper (BPH, Nilaparvata lugens), the small brown planthopper (SBPH, Laodelphax striatellus), and the white-backed planthopper (WBPH, Sogatella furcifera) are problematic insect pests and cause severe yield losses through phloem sap-sucking and virus transmission. Barley yellow striate mosaic virus (BYSMV), a plant cytorhabdovirus, has been developed as versatile expression platforms in SBPHs and cereal plants. However, bio-safe overexpression vectors based on recombinant BYSMV (rBYSMV) remain to be developed and applied to the three kinds of planthoppers. Here, we found that rBYSMV was able to infect SBPHs, BPHs and WBPHs through microinjection with crude extracts from rBYSMV-infected barley leaves. To ensure bio-safety of the rBYSMV vectors, we generated an rBYSMV mutant by deleting the accessory protein P3, a putative viral movement protein. As expected, the resulting mutant abolished viral systemic infection in barley plants but had no effects on BYSMV infectivity in insect vectors. Subsequently, we used the modified rBYSMV vector to overexpress iron transport peptide (ITP) in the three kinds of planthoppers and revealed the potential functions of ITP. Overall, our results provide bio-safe overexpression platforms to facilitate functional genomics studies of planthoppers.

Alternative splicing landscape of small brown planthopper and different response of JNK2 isoforms to rice stripe virus infection

Alternative splicing (AS) is a frequent posttranscriptional regulatory event occurring in response to various endogenous and exogenous stimuli in most eukaryotic organisms. However, little is known about the effects of insect-transmitted viruses on AS events in insect vectors. The present study used third-generation sequencing technology and RNA sequencing (RNA-Seq) to evaluate the AS response in the small brown planthopper Laodelphax striatellus to rice stripe virus (RSV). The full-length transcriptome of L. striatellus was obtained using single-molecule real-time sequencing technology (SMRT). Posttranscriptional regulatory events, including AS, alternative polyadenylation, and fusion transcripts, were analyzed. A total of 28,175 nonredundant transcript isoforms included 24,950 transcripts assigned to 8,500 annotated genes of L. striatellus, and 5,000 of these genes (58.8%) had AS events. RNA-Seq of the gut samples of insects infected by RSV for 8 d identified 3,458 differentially expressed transcripts (DETs); 2,185 of these DETs were transcribed from 1,568 genes that had AS events, indicating that 31.4% of alternatively spliced genes responded to RSV infection of the gut. One of the c-Jun N-terminal kinase (JNK) genes, JNK2, experienced exon skipping, resulting in three transcript isoforms. These three isoforms differentially responded to RSV infection during development and in various organs. Injection of double-stranded RNAs targeting all or two isoforms indicated that three or at least two JNK2 isoforms facilitated RSV accumulation in planthoppers. These results implied that AS events could participate in the regulation of complex relationships between viruses and insect vectors. Importance Alternative splicing (AS) is a regulatory mechanism that occurs after gene transcription. AS events can enrich protein diversity to promote the reactions of the organisms to various endogenous and exogenous stimulations. It is not known how insect vectors exploit AS events to cope with transmitted viruses. The present study used third-generation sequencing technology to obtain the profile of AS events in the small brown planthopper Laodelphax striatellus, which is an efficient vector for rice stripe virus (RSV). The results indicated that 31.4% of alternatively spliced genes responded to RSV infection in the gut of planthoppers. One of the c-Jun N-terminal kinase (JNK) genes, JNK2, produced three transcript isoforms by AS. These three isoforms showed different responses to RSV infection, and at least two isoforms facilitated viral accumulation in planthoppers. These results implied that AS events could participate in the regulation of complex relationships between viruses and insect vectors.

Vitellogenin from planthopper oral secretion acts as a novel effector to impair plant defenses

Vitellogenin (Vg) is a well-known nutritious protein involved in reproduction in nearly all oviparous animals, including insects. Recently, Vg has been detected in saliva proteomes of several piercing-sucking herbivorous arthropods, including the small brown planthopper (Laodelphax striatellus, SBPH). Its function, however, remains unexplored. We investigated the molecular mechanism underlying SBPH orally secreted Vg-mediated manipulation of plant-insect interaction by RNA interference, phytohormone and H₂ O₂ profiling, protein-protein interaction studies and herbivore bioassays. A C-terminal polypeptide of Vg (VgC) in SBPH, when secreted into rice plants, acted as a novel effector to attenuate host rice defenses, which in turn improved insect feeding performance. Silencing Vg reduced insect feeding and survival on rice. Vg-silenced SBPH nymphs consistently elicited higher H₂ O₂ production, a well-established defense mechanism in rice, whereas expression of VgC in planta significantly hindered hydrogen peroxide (H₂ O₂ ) accumulation and promoted insect performance. VgC interacted directly with the rice transcription factor OsWRKY71, a protein which is involved in induction of H₂ O₂ accumulation and plant resistance to SBPH. These findings indicate a novel effector function of Vg: when secreted into host rice plants, this protein effectively weakened H₂ O₂ -mediated plant defense through its association with a plant immunity regulator.

Chromosome-level genome assembly of the mirid predator Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), an important natural enemy in the rice ecosystem

Though the genomes of many rice herbivorous pests have recently been well characterized, little is known about the genome of their natural enemies. Here, by using the Illumina and PacBio platforms, we sequenced and assembled the whole genome of the mirid species Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), which is an economically and ecologically important natural enemy in the rice ecosystem acting as a dominant predator for planthoppers and leafhoppers in the field. Through Hi-C scaffolding, 1615 scaffolds with a total size of 338.08 Mb were successfully anchored onto 13 chromosomes. The assembled genome size was 345.75 Mb with a final scaffold N50 of 27.58 Mb. Approximately 107.51 Mb of sequences accounting for 31.10% of the genome were identified as repeat elements, and 14,644 protein-coding genes were annotated. Phylogenetic analysis showed that C. lividipennis clustered with other Hemipteran species and diverged from Apolygus lucorum about 66.7 million years ago. Gene families related to detoxification, environmental adaptation and digestion were analysed comparatively with other Hemipteran species, but no significant expansion or contraction was found in C. lividipennis. We also observed male meiosis in C. lividipennis, which showed a typical post-reduction of sex chromosomes and a karyotype of 2n = 22 + XY. As the first natural-enemy genome in the rice ecosystem, the genomic resource of C. lividipennis not only expands our understanding of the multitrophic interactions (host plant-prey-predator), but also provides a genomic basis for better understanding this dominant predator and therefore promotes sustainable rice pest management and food grain production.

Improved draft reference genome for the Glassy-winged Sharpshooter (Homalodisca vitripennis), a vector for Pierce's disease

Homalodisca vitripennis (Hemiptera: Cicadellidae), known as the glassy-winged sharpshooter, is a xylem feeding leafhopper and an important agricultural pest as a vector of Xylella fastidiosa, which causes Pierce’s disease in grapes and a variety of other scorch diseases. The current H. vitripennis reference genome from the Baylor College of Medicine's i5k pilot project is a 1.4-Gb assembly with 110,000 scaffolds, which still has significant gaps making identification of genes difficult. To improve on this effort, we used a combination of Oxford Nanopore long-read sequencing technology combined with Illumina sequencing reads to generate a better assembly and first-pass annotation of the whole genome sequence of a wild-caught Californian (Tulare County) individual of H. vitripennis. The improved reference genome assembly for H. vitripennis is 1.93-Gb in length (21,254 scaffolds, N50 = 650 Mb, BUSCO completeness = 94.3%), with 33.06% of the genome masked as repetitive. In total, 108,762 gene models were predicted including 98,296 protein-coding genes and 10,466 tRNA genes. As an additional community resource, we identified 27 orthologous candidate genes of interest for future experimental work including phenotypic marker genes like white. Furthermore, as part of the assembly process, we generated four endosymbiont metagenome-assembled genomes, including a high-quality near complete 1.7-Mb Wolbachia sp. genome (1 scaffold, CheckM completeness = 99.4%). The improved genome assembly and annotation for H. vitripennis, curated set of candidate genes, and endosymbiont MAGs will be invaluable resources for future research of H. vitripennis.

Comparative Transcriptome Analysis of Chemoreception Organs of Laodelphax striatellus in Response to Rice Stripe Virus Infection

Many vector-borne viruses possess the ability to manipulate vector behaviors to facilitate their transmission. There is evidence that the mechanism of this phenomenon has been described in part as direct manipulation through regulating vector chemosensation. Rice stripe virus (RSV) is transmitted by the small brown planthopper, Laodelphax striatellus (Fallen), in a persistent, circulative-propagative manner. The effect of RSV infection on the olfactory system of L. striatellus has not been fully elucidated. Here, we employed transcriptomic sequencing to analyze gene expression profiles in antennae, legs and heads (without antennae) from L. striatellus females and males with/without RSV infection. Comparisons of the differentially expressed genes (DEGs) among antennae, legs and heads indicated that tissue-specific changes in the gene expression profile were greater than sex-specific changes. A total of 17 olfactory related genes were differentially expressed in viruliferous antennae as compared to nonviruliferous antennae, including LstrOBP4/9, LstrCSP1/2/5, LstrGR28a/43a/43a-1, LstrIR1/2/NMDA1, LstrOR67/85e/56a/94 and LstrSNMP2/2-2. There are 23 olfactory related DEGs between viruliferous and nonviruliferous legs, including LstrOBP2/3/4/12/13, LstrCSP13/5/10, LstrIR1/2/Delta2/Delta2-1/kainate2/NMDA2, LstrOR12/21/31/68 and LstrORco. A low number of olfactory related DEGs were found between viruliferous and nonviruliferous heads, including LstrCSP1, LstrOBP2, LstrOR67 and LstrSNMP2-2. Among these DEGs, the expression patterns of LstrOBP2, LstrOBP3 and LstrOBP9 in three tissues was validated by quantitative real-time PCR. The demonstration of overall changes in the genes in L. striatellus' chemoreception organs in response to RSV infection would not only improve our understanding of the effect of RSV on the olfactory related genes of insect vectors but also provide insights into developing approaches to control the plant virus transmission and spread as well as pest management in the future.

RNA N6 -methyladenosine modification suppresses replication of rice black streaked dwarf virus and is associated with virus persistence in its insect vector

N⁶ methylation of adenosine (m⁶ A) was recently discovered to play a role in regulating the life cycle of various viruses by modifying viral and host RNAs. However, different studies on m⁶ A effects on the same or different viruses have revealed contradictory roles for m⁶ A in the viral life cycle. In this study, we sought to define the role of m⁶ A on infection by rice black streaked dwarf virus (RBSDV), a double-stranded RNA virus, of its vector small brown planthopper (SBPH). Infection by RBSDV decreased the level of m⁶ A in midgut cells of SBPHs. We then cloned two genes (LsMETTL3 and LsMETTL14) that encode m⁶ A RNA methyltransferase in SBPHs. After interference with expression of the two genes, the titre of RBSDV in the midgut cells of SBPHs increased significantly, suggesting that m⁶ A levels were negatively correlated with virus replication. More importantly, our results revealed that m⁶ A modification might be the epigenetic mechanism that regulates RBSDV replication in its insect vector and maintains a certain virus threshold required for persistent transmission.

Regulation of RNA Interference Pathways in the Insect Vector Laodelphax striatellus by Viral Proteins of Rice Stripe Virus

RNA interference (RNAi), especially the small interfering RNA (siRNA) and microRNA (miRNA) pathways, plays an important role in defending against viruses in plants and insects. However, how insect-transmitted phytoviruses regulate the RNAi-mediated antiviral response in vector insects has barely been uncovered. In this study, we explored the interaction between rice stripe virus (RSV) and the miRNA and siRNA pathways of the small brown planthopper, which is a vector insect. The transcript and protein levels of key genes in the two RNAi pathways did not change during the RSV infection process. When the expression of insect Ago1, Ago2, or Translin was silenced by the injection of double-stranded RNAs targeting these genes, viral replication was promoted with Ago2 silencing but inhibited with Translin silencing. Protein-protein binding assays showed that viral NS2 and RNA-dependent RNA polymerase interacted with insect Ago2 and Translin, respectively. When NS2 was knocked down, the transcript level of Ago2 increased and viral replication was inhibited. Therefore, viral NS2 behaved like an siRNA suppressor in vector insects. This protein-binding regulation of insect RNAi systems reflects a complicated and diverse coevolution of viruses with their vector insects.

Cytochrome P450 Monooxygenases CYP6AY3 and CYP6CW1 Regulate Rice Black-Streaked Dwarf Virus Replication in Laodelphax striatellus (Fallén)

The small brown planthopper, Laodelphax striatellus (Fallén), is an important agricultural pest that causes significant losses by sucking and transmitting multiple plant viruses, such as rice black-streaked dwarf virus (RBSDV). Insecticides are commonly used to control planthoppers and cause the induction or overexpression of cytochrome P450 monooxygenases (P450s) from the CYP3 and CYP4 clades after insecticide application. However, little is known about the roles of insecticides and P450s in the regulation of viral replication in insects. In this study, RBSDV-infected L. striatellus were injected with imidacloprid, deltamethrin, pymetrozine, and buprofezin, respectively. The insecticide treatments caused a significant decrease in RBSDV abundance in L. striatellus. Treatment of piperonyl butoxide (PBO), an effective inhibitor of P450s, significantly increased the RBSDV abundance in L. striatellus. Fourteen P450 candidate genes in the CYP3 clade and 21 in the CYP4 clade were systematically identified in L. striatellus, and their expression patterns were analyzed under RBSDV infection, in different tissues, and at different developmental stages. Among the thirty-five P450 genes, the expression level of CYP6CW1 was the highest, while CYP6AY3 was the lowest after RBSDV infection. Knockdown of CYP6CW1 and CYP6AY3 significantly increased the virus abundance and promoted virus replication in L. striatellus. Overall, our data reveal that CYP6CW1 and CYP6AY3 play a critical role in the regulation of virus replication in L.striatellus.

Characterization of protein-protein interactions between rice viruses and vector insects

Planthoppers are the most notorious rice pests, because they transmit various rice viruses in a persistent-propagative manner. Protein-protein interactions (PPIs) between virus and vector are crucial for virus transmission by vector insects. However, the number of known PPIs for pairs of rice viruses and planthoppers is restricted by low throughput research methods. In this study, we applied DeNovo, a virus-host sequence-based PPI predictor, to predict potential PPIs at a genome-wide scale between three planthoppers and five rice viruses. PPIs were identified at two different confidence thresholds, referred to as low and high modes. The number of PPIs for the five planthopper-virus pairs ranged from 506 to 1985 in the low mode and from 1254 to 4286 in the high mode. After eliminating the "one-too-many" redundant interacting information, the PPIs with unique planthopper proteins were reduced to 343-724 in the low mode and 758-1671 in the high mode. Homologous analysis showed that 11 sets and 31 sets of homologous planthopper proteins were shared by all planthopper-virus interactions in the two modes, indicating that they are potential conserved vector factors essential for transmission of rice viruses. Ten PPIs between small brown planthopper and rice stripe virus (RSV) were verified using glutathione-S-transferase (GST)/His-pull down or co-immunoprecipitation assay. Five of the ten PPIs were proven positive, and three of the five SBPH proteins were confirmed to interact with RSV. The predicted PPIs provide new clues for further studies of the complicated relationship between rice viruses and their vector insects.

iTRAQ-based quantitative proteomics suggests mitophagy involvement after Rice black-streaked dwarf virus acquisition in insect vector small brown planthopper Laodelphax striatellus Fallén

Plant viruses trigger numerous responses in their insect vectors. Using iTRAQ-based quantitative proteomics analysis, early responses of the insect vector, the small brown planthopper (Laodelphax striatellus Fallén, SBPH), after acquiring Rice black-streaked dwarf virus (RBSDV) at 3 days and 5 days post first access to diseased plants (padp) were revealed. A total of 582 differentially abundant proteins (DAPs) in SBPH with a fold change >1.500 or <0.667 (p-value < 0.05) were identified. The proteomic analysis in SBPH at 3 days padp revealed 106 highly abundant proteins and 193 of low abundance, while 5 days padp revealed 214 highly abundant proteins and 182 of low abundance. Among them, 51 highly abundant proteins and 42 of low abundance were shown consistently at both 3 days and 5 days padp. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis mapping and Gene Ontology (GO) term classification suggested impairment of mitochondria in SBPH after RBSDV acquisition, and the 77 out of 582 differentially abundant SBPH proteins analyzed by the STRING program revealed the interaction network of the mitochondrial DAPs, showing an overall down-regulation of mitochondrial proteins including the electron transport chain proteins and mitochondrial ribosome proteins. The high abundance of Parkin at 5 days padp suggests that activation of mitophagy induced degradation of mitochondria occurred. Further verification of autophagy/mitophagy-related genes by reverse-transcription quantitative RT-PCR (RT-qPCR) in SBPH after RBSDV acquisition showed up-regulation of the autophagy receptors Optineurin (OPTN), Sequestosome-1 (SQSTM1, also known as p62) and Tax1-binding protein 1 (TAX1BP1) which targets ubiquitinated damaged mitochondria during mitophagy. The phosphorylation of the three autophagy receptors may be up-regulated through an increase of transcription level TRAF-associated NFκB activator (TANK)-binding kinase 1 (TBK1). As a result, an overall reduction in the abundance of mitochondrial proteins was observed and the selective autophagic degradation was up-regulated through increased transcription level of OPTN, p62/SQSTM1, TAX1BP1 and TBK1. Therefore, acquisition of RBSDV associated with up-regulated autophagy and selective mitochondrial degradation in SBPH suggest prevention of mitochondrial-mediated apoptosis and extension of the vector life span. BIOLOGICAL SIGNIFICANCE: RBSDV causes severe yield loss in rice plants. RBSDV is transmitted efficiently only through SBPH. It is important to understand how RBSDV infects SBPH in a persistent, circulative and propagative manner. However, there has been no study on the interaction between RBSDV and SBPH at the early acquisition stage using a proteomics approach. In this study, we combined iTRAQ technique and LC-MS/MS to analyze the vector proteomics at both the initial and latent infection stages after RBSDV acquisition and verified the results by RT-qPCR. Our results revealed that significantly low DAPs were involved in various pathways, including biosynthesis of secondary metabolites, ribosomes, carbon metabolism, biosynthesis of amino acids and TCA cycle. Further clustering of the DAPs revealed significant changes in SBPH mitochondria, including decreased proteins in mitochondrial ribosomes and electron transport chain complex I, II and V. On the other hand, there was a high abundance of Parkin, suggesting the occurrence of mitochondria damage and subsequent Parkin-mediated mitophagy for clearance of impaired mitochondria. Moreover, the decreased level of PMPCB in terms of gene expression and protein abundance suggested decreased PINK1 turnover, promoting Parkin/PINK1-mediated mitophagy. Further analysis on autophagy/mitophagy-related gene transcription level indicated up-regulation of OPTN, p62/SQSTM1, TAX1BP1 and TBK1, promoting selective autophagy in SBPH after RBSDV acquisition. These findings provided new insights into the effects of RBSDV on SBPH after early acquisition by selective degradation of mitochondria, especially on reprogramming of energy metabolism and decreased mitochondria biogenesis, to prevent apoptosis and prolong the life span of SBPH post virus acquisition.

The Draft Genome of Yellow Stem Borer, an Agriculturally Important Pest, Provides Molecular Insights into Its Biology, Development and Specificity Towards Rice for Infestation

Simple Summary

Yellow stem borer (YSB), is the most destructive and widely occurring pest that attacks rice throughout the growing season. Rice (Oryza sativa L.) is a major staple cereal worldwide, providing essential caloric requirements for more than half of the world’s population. Annual losses to rice borers are approximately 5–10%, but losses in individual fields may reach up to 50–60%. The use of traditional pest management strategies in controlling YSB is somewhat challenging due to its unique internal feeding habit. Genome sequence information of economically important crop pests is important for designing or developing pest-resistant rice varieties. In an approach to achieve this, we present our first-ever study on the draft genome sequence of YSB. The information provided from our current study might be useful in developing genome-based approaches for the management of pest species.

Abstract

Yellow stem borer (YSB), Scirpophaga incertulas (Walker) (Lepidoptera: Crambidae), a major monophagous insect pest of rice, causes significant yield losses. The rice–YSB interaction is very dynamic, making it difficult for management. The development of resistant lines has been unsuccessful as there are no effective resistant sources in the germplasm. Genome information is necessary for a better understanding of interaction with rice in terms of its recognition, response, and infestation mechanism. The draft genome of YSB is predicted to have 46,057 genes with an estimated size of 308 Mb, being correlated with the flow cytometry analysis. The existence of complex metabolic mechanisms and genes related to specific behavior was identified, being conditioned by a higher level of regulation. We deciphered the possible visual, olfactory, and gustatory mechanisms responsible for its evolution as a monophagous pest. Comparative genomic analysis revealed that YSB is unique in the way it has evolved. The obvious presence of high-immunity-related genes, well-developed RNAi machinery, and diverse effectors provides a means for developing genomic tools for its management. The identified 21,696 SSR markers can be utilized for diversity analysis of populations across the rice-growing regions. We present the first draft genome of YSB. The information emanated paves a way for biologists to design novel pest management strategies as well as for the industry to design new classes of safer and specific insecticide molecules.

Novel alternative splicing of GABA receptor RDL exon 9 from Laodelphax striatellus modulates agonist potency

The resistance to dieldrin gene (RDL) encodes the primary subunit of the insect ionotropic γ-aminobutyric acid (GABA) receptor (GABAR), which is the target of phenylpyrazole and isoxazoline insecticides. The splice variants in exons 3 and 6 of RDL, which have been widely explored in many insects, modulate the agonist potency of the homomeric RDL GABAR and potentially play an important role in the development of insects. In the present study, four splice variants of exon 9 were identified in RDL of the small brown planthopper, Laodelphax striatellus (LsRDL), resulting in LsRDL-9a, LsRDL-9a', LsRDL-9b, and LsRDL-9c. LsRDL-9a has one more amino acid (E, glutamic acid) compared with LsRDL-9a', and LsRDL-9b lacked two amino acids and had seven different amino acids compared with LsRDL-9c. Two-electrode voltage-clamp recording on LsRDLs expressed in Xenopus oocytes showed that alternative splicing of exon 9 has significant impact on LsRDL sensitivity to the agonists GABA and β-alanine, whereas no significant difference was observed in the potencies of the non-competitive antagonists (NCAs) ethiprole and fluralaner on the splice variants. Our results suggest that alternative splicing of RDL exon 9 broadens functional capabilities of the GABAR in L. striatellus by influencing the action of GABA.

A lateral oviduct secreted protein plays a vital role for egg movement through the female reproductive tract in the brown planthopper

The oviduct serves as a delivery tube for mature eggs ovulated from ovaries to egg-laying sites. Oviduct secreted components play important roles in ovulation and fertilization in mammals, however, no oviduct secreted protein has been characterized in an insect to date. Here, we identified a gene highly expressed in the lateral oviduct of the adult females in the brown planthopper (BPH), Nilaparvata lugens, the most destructive rice insect pest. Western blotting and immunofluorescence analyses revealed that the gene encodes a protein that is specifically expressed in the lateral oviduct as a component of the gel-like material secreted by the oviduct epithelial cells into the lumen of the swollen part of the lateral oviducts. The protein was tentatively named N. lugens oviduct secreted protein (Nlodsp). RNA interference (RNAi) against NlOdsp transcripts caused a failure of the lateral oviducts to deliver oocytes to the common oviduct that was, by consequence, plugged by 1-2 oocytes. Moreover, although oocytes in the Nlodsp-deficient ovariole were not released to the oviduct, they continued to develop, finally resulting in the presence of several matured oocytes in an ovariole. These defects evidently declined female fecundity. Together, our results demonstrate that NlOdsp plays an essential role in egg transport through the oviduct during ovulation. This work deepens our understanding of insect reproductive system and provides a potential target gene for RNAi-based insect pest control.

Long Non-Coding RNAs in Insects

Only a small subset of all the transcribed RNAs are used as a template for protein translation, whereas RNA molecules that are not translated play a very important role as regulatory non-coding RNAs (ncRNAs). Besides traditionally known RNAs (ribosomal and transfer RNAs), ncRNAs also include small non-coding RNAs (sncRNAs) and long non-coding RNAs (lncRNAs). The lncRNAs, which were initially thought to be junk, have gained a great deal attention because of their regulatory roles in diverse biological processes in animals and plants. Insects are the most abundant and diverse group of animals on this planet. Recent studies have demonstrated the role of lncRNAs in almost all aspects of insect development, reproduction, and genetic plasticity. In this review, we describe the function and molecular mechanisms of the mode of action of different insect lncRNAs discovered up to date.

Laodelphax striatellus Atg8 facilitates Rice stripe virus infection in an autophagy-independent manner

Rice stripe virus (RSV) is the causative agent of rice stripe disease and is completely dependent on insect vectors for its plant-to-plant transmission. Laodelphax striatellus is the major insect vector for RSV. In this study, we explored the interactions between RSV infection and L. striatellus autophagy, a potential intrinsic antiviral mechanism in insects. We found that L. striatellus autophagic activity did not affect RSV infection; however, the autophagy-related-8 (Atg8) gene significantly enhanced virus infection. During RSV initial infection within the L. striatellus midgut, silencing of Atg8 expression significantly decreased the phosphorylation of c-Jun N-terminal kinase (p-JNK); however, when RSV infection is absent, silencing of Atg8 did not alter p-JNK levels. These results indicated that Atg8 might activate the JNK machinery by allowing more virus infection into cells. We further revealed that Atg8-deficiency significantly decreased RSV accumulation on the surface of the insect midgut epithelial cells, suggesting a receptor trafficking function of the γ-aminobutyric acid receptor-associated protein family. Using the RSV ovary entry as a model, in which vitellogenin receptor (VgR) mediates RSV cell entry, we clarified that Atg8-deficiency decreased the abundance of VgR localizing on the cytomembrane and disturbed the attachment of RSV in the germarium zones. Collectively, these results revealed an autophagy-independent function of L. striatellus Atg8 that enhances RSV initial infection by increasing virus attachment on the infection sites.

Coordination between terminal variation of the viral genome and insect microRNAs regulates rice stripe virus replication in insect vectors

Maintenance of a balance between the levels of viral replication and selective pressure from the immune systems of insect vectors is one of the prerequisites for efficient transmission of insect-borne propagative phytoviruses. The mechanism regulating the adaptation of RNA viruses to insect vectors by genomic variation remains unknown. Our previous study demonstrated an extension of the 3’-untranslated terminal region (UTR) of two genomic segments of rice stripe virus (RSV). In the present study, a reverse genetic system for RSV in human cells and an insect vector, the small brown planthopper Laodelphax striatellus, was used to demonstrate that the 3’-terminal extensions suppressed viral replication in vector insects by inhibiting promoter activity due to structural interference with the panhandle structure formed by viral 3’- and 5’-UTRs. The extension sequence in the viral RNA1 segment was targeted by an endogenous insect microRNA, miR-263a, which decreased the inhibitory effect of the extension sequence on viral promoter activity. Surprisingly, the expression of miR-263a was negatively regulated by RSV infection. This elaborate coordination between terminal variation of the viral genome and endogenous insect microRNAs controls RSV replication in planthopper, thus reflecting a distinct strategy of adaptation of phytoviruses to insect vectors.

Mutations frequently happen when insect-transmitted RNA viruses circulate between insect vectors and plant or mammalian hosts. However, the significance of these mutations for viral fitness in the two distinct organisms is poorly understood. We discovered that a high proportion of rice stripe virus (RSV) had terminally extended genomes when the virus infected insect vectors. In the present study, we found that the extension sequence suppressed viral replication in insect vectors by impairing a special structure formed by the two ends of the viral genomes. An endogenous insect small RNA was able to bind the extension sequence to relieve the inhibitory effect. However, the expression of this small RNA was reduced in the presence of RSV to ultimately maintain the inhibitory effect of the extension sequence. This elaborate coordination between virus and vector enables a limited level of RSV replication that does not produce serious damage to vectors, thus reflecting a distinct strategy of adaptation of insect-transmitted plant viruses.

The nucleocapsid protein of rice stripe virus in cell nuclei of vector insect regulates viral replication

Rice stripe virus (RSV) transmitted by the small brown planthopper causes severe rice yield losses in Asian countries. Although viral nuclear entry promotes viral replication in host cells, whether this phenomenon occurs in vector cells remains unknown. Therefore, in this study, we systematically evaluated the presence and roles of RSV in the nuclei of vector insect cells. We observed that the nucleocapsid protein (NP) and viral genomic RNAs were partially transported into vector cell nuclei by utilizing the importin α nuclear transport system. When blocking NP nuclear localization, cytoplasmic RSV accumulation significantly increased. In the vector cell nuclei, NP bound the transcription factor YY1 and affected its positive regulation to FAIM. Subsequently, decreased FAIM expression triggered an antiviral caspase-dependent apoptotic reaction. Our results reveal that viral nuclear entry induces completely different immune effects in vector and host cells, providing new insights into the balance between viral load and the immunity pressure in vector insects.

Filling in the gaps: A reevaluation of the Lygus hesperus peptidome using an expanded de novo assembled transcriptome and molecular cloning

Peptides are the largest and most diverse class of molecules modulating physiology and behavior. Previously, we predicted a peptidome for the western tarnished plant bug, Lygus hesperus, using transcriptomic data produced from whole individuals. A potential limitation of that analysis was the masking of underrepresented genes, in particular tissue-specific transcripts. Here, we reassessed the L. hesperus peptidome using a more comprehensive dataset comprised of the previous transcriptomic data as well as tissue-specific reads produced from heads and accessory glands. This augmented assembly significantly improves coverage depth providing confirmatory transcripts for essentially all of the previously identified families and new transcripts encoding a number of new peptide precursors corresponding to 14 peptide families. Several families not targeted in our initial study were identified in the expanded assembly, including agatoxin-like peptide, CNMamide, neuropeptide-like precursor 1, and periviscerokinin. To increase confidence in the in silico data, open reading frames of a subset of the newly identified transcripts were amplified using RT-PCR and sequence validated. Further PCR-based profiling of the putative L. hesperus agatoxin-like peptide precursor revealed evidence of alternative splicing with near ubiquitous expression across L. hesperus development, suggesting the peptide serves functional roles beyond that of a toxin. The peptides predicted here, in combination with those identified in our earlier study, expand the L. hesperus peptidome to 42 family members and provide an improved platform for initiating molecular and physiological investigations into peptidergic functionality in this non-model agricultural pest.

Bamboo Transposon Research: Current Status and Perspectives

Bamboo, a fast-growing non-timber forest plant with many uses, is a valuable species for green development. However, bamboo flowering is very infrequent, extending, in general, for up to 120 years. Ecologically, bamboo species are generally better adapted to various environments than other grasses. Therefore, the species deserves a special status in what could be called Ecological Bioeconomy. An understanding of the genetic processes of bamboo can help us sustainably develop and manage bamboo forests. Transposable elements (TEs), jumping genes or transposons, are major genetic elements in plant genomes. The rapid development of the bamboo reference genome, at the chromosome level, reveals that TEs occupy over 63.24% of the genome. This is higher than found in rice, Brachypodium, and sorghum. The bamboo genome contains diverse families of TEs, which play a significant role in bamboo's biological processes including growth and development. TEs provide important clues for understanding the evolution of the bamboo genome. In this chapter, we briefly describe the current status of research on TEs in the bamboo genome, their regulation, and transposition mechanisms. Perspectives for future research are also provided.

Characterization and analysis of the transcriptome in Opisina arenosella from different developmental stages using single-molecule real-time transcript sequencing and RNA-seq

Opisina arenosella is one of the main pests harming coconut trees. To date, there have been few studies on the molecular genetics, biochemistry and physiology of O. arenosella at the transcriptional level, and there are no available reference genomes. Here, Illumina RNA sequencing combined with PacBio single-molecule real-time analysis was applied to study the transcriptome of this pest at different developmental stages, providing reference data for transcript expression analysis. Twelve samples of O. arenosella from different stages of development were sequenced using Illumina RNA sequencing, and the pooled RNA samples were sequenced with PacBio technology (Iso-Seq). A full-length transcriptome with 41,938 transcripts was captured, and the N50 and N90 lengths were 3543 bp and 1646 bp, respectively. A total of 36,925 transcripts were annotated in public databases, 6493 of which were long noncoding RNAs, while 2510 represented alternative splicing events. There were significant differences in the gene expression profiles at different developmental stages, with high levels of differential gene expression associated with growth, development, carbohydrate metabolism and immunity. This work provides resources and information for the study of the transcriptome and gene function of O. arenosella and provides a valuable foundation for understanding the changes in gene expression during development.

Regulatory Functions of Nilaparvata lugens GSK-3 in Energy and Chitin Metabolism

Glucose metabolism is a biologically important metabolic process. Glycogen synthase kinase (GSK-3) is a key enzyme located in the middle of the sugar metabolism pathway that can regulate the energy metabolism process in the body through insulin signaling. This paper mainly explores the regulatory effect of glycogen synthase kinase on the metabolism of glycogen and trehalose in the brown planthopper (Nilaparvata lugens) by RNA interference. In this paper, microinjection of the target double-stranded GSK-3 (dsGSK-3) effectively inhibited the expression of target genes in N. lugens. GSK-3 gene silencing can effectively inhibit the expression of target genes (glycogen phosphorylase gene, glycogen synthase gene, trehalose-6-phosphate synthase 1 gene, and trehalose-6-phosphate synthase 2 gene) in N. lugens and trehalase activity, thereby reducing glycogen and glucose content, increasing trehalose content, and regulating insect trehalose balance. GSK-3 can regulate the genes chitin synthase gene and glucose-6-phosphate isomerase gene involved in the chitin biosynthetic pathway of N. lugens. GSK-3 gene silencing can inhibit the synthesis of chitin N. lugens, resulting in abnormal phenotypes and increased mortality. These results indicated that a low expression of GSK-3 in N. lugens can regulate the metabolism of glycogen and trehalose through the insulin signal pathway and energy metabolism pathway, and can regulate the biosynthesis of chitin, which affects molting and wing formation. The relevant research results will help us to more comprehensively explore the molecular mechanism of the regulation of energy and chitin metabolism of insect glycogen synthase kinases in species such as N. lugens.