Defense Regulatory Network Associated with circRNA in Rice in Response to Brown Planthopper Infestation.

The brown planthopper (BPH), Nilaparvata lugens (Stål), a rice-specific pest, has risen to the top of the list of significant pathogens and insects in recent years. Host plant-mediated resistance is an efficient strategy for BPH control. Nonetheless, BPH resistance in rice cultivars has succumbed to the emergence of distinct virulent BPH populations. Circular RNAs (circRNAs) play a pivotal role in regulating plant-environment interactions; however, the mechanisms underlying their insect-resistant functions remain largely unexplored. In this study, we conducted an extensive genome-wide analysis using high-throughput sequencing to explore the response of rice circRNAs to BPH infestations. We identified a total of 186 circRNAs in IR56 rice across two distinct virulence groups: IR-IR56-BPH (referring to IR rice infested by IR56-BPH) and IR-TN1-BPH, along with a control group (IR-CK) without BPH infestation. Among them, 39 circRNAs were upregulated, and 43 circRNAs were downregulated in the comparison between IR-IR56-BPH and IR-CK. Furthermore, in comparison with IR-CK, 42 circRNAs exhibited upregulation in IR-TN1-BPH, while 42 circRNAs showed downregulation. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that the targets of differentially expressed circRNAs were considerably enriched in a multitude of biological processes closely linked to the response to BPH infestations. Furthermore, we assessed a total of 20 randomly selected circRNAs along with their corresponding expression levels. Moreover, we validated the regulatory impact of circRNAs on miRNAs and mRNAs. These findings have led us to construct a conceptual model that circRNA is associated with the defense regulatory network in rice, which is likely facilitated by the mediation of their parental genes and competing endogenous RNA (ceRNA) networks. This model contributes to the understanding of several extensively studied processes in rice-BPH interactions.

Characterization of light-dependent rhythm of courtship vibrational signals in Nilaparvata lugens: essential involvement of cryptochrome genes.

BACKGROUND: Vibrational signal plays a crucial role in courtship communication in many insects. However, it remains unclear whether insect vibrational signals exhibit daily rhythmicity in response to changes in environmental cues. RESULTS: In this study, we observed daily rhythms of both female vibrational signals (FVS) and male vibrational signals (MVS) in the brown planthopper (BPH), Nilaparvata lugens (Stål), one of the most notorious rice pests across Asia. Notably, oscillations of FVS and MVS in paired BPHs were synchronized as part of male-female duetting interactions, displaying significant day-night rhythmicity. Furthermore, we observed light dependency of FVS emissions under different photoperiodic regimes (18 L:6 D and 6 L:18 D) and illumination intensity levels (>300 lx, 50 lx, and 25 lx). Subsequently, the potential role of circadian clock genes cryptochromes (Nlcry1 and Nlcry2) in regulating FVS daily oscillations was examined using gene knockdown via RNA interference. We observed sharp declines and disrupted rhythms in FVS frequencies when either of the Nlcrys was downregulated, with Nlcry2 knockdown showing a more prominent effect. Moreover, we recorded a novel FVS variant (with a dominant frequency of 361.76 ± 4.31 Hz) emitted by dsNlcry1-treated BPH females, which significantly diminished the impact of courtship stimuli on receptive males. CONCLUSION: We observed light-dependent daily rhythms of substrate-borne vibrational signals (SBVS) in BPH and demonstrated essential yet distinct roles of the two Nlcrys. These findings enhanced our understanding of insect SBVS and illustrated the potential of novel precision physical control strategies for disrupting mating behaviors in this rice pest. © 2023 Society of Chemical Industry.

Knockdown of the chromatin remodeling ATPase gene Brahma impairs the reproductive potential of the brown planthopper, Nilaparvata lugens.

The brown planthopper (BPH), Nilaparvata lugens (Stål), is one of the most destructive pests in rice-growing regions of Asia. Extensive studies have suggested that SWI/SNF chromatin remodeling ATPase Brahma (BRM) plays multiple roles in the insect model Drosophila. Yet much less is known about the physiological properties for NlBRM. In the present study, the cloned full-length cDNA of NlBRM was 5637 bp and contained an ORF of 5292 bp encoding a 194.53 kD protein. The spatiotemporal dynamics of NlBRM was investigated by qPCR, which showed that it was abundantly expressed in the egg and ovary. Then significant downregulation of NlBRM by dsRNA injection had a relatively greater impact on female survival than male. Moreover, the number of oviposition marks of the NlBRM-RNAi females were declined by 61.11% - 73.33% compared with the controls during the subsequent 5 days after dsRNA injection. Meanwhile, the number of newly hatched BPH nymphs also decreased correspondingly by 93.56% - 100%. Phenotypic analysis revealed that none of normally banana-shaped eggs were discernable in the ovaries of NlBRM-deficient females, where mRNA expression of N. lugens vitellogenin gene was also reduced. Our results demonstrated that NlBRM played a crucial role in ovarian development and fecundity of BPH, likely by regulating the vitellogenin gene in vivo, which could be as a promising target for parental RNAi-based control of this serious rice pest.

Phenotypic and transcriptomic responses of two Nilaparvata lugens populations to the Mudgo rice containing Bph1.

The Bph1 gene was the first reported brown planthopper (BPH, Nilaparvata lugens) resistance gene in Mudgo rice and was widely used as a commercial cultivar for controlling BPH infestations. However, rapid adaptations of BPH on the Mudgo rice resulted in its resistance breakdown and the emergence of virulent BPH populations. Thus, specific BPH populations and rice varieties can serve as good model systems for studying the roles of different bio-compounds and proteins in the insect-plant interactions. Although our understandings have been improved on the complexity of BPH and rice interactions, the underlying molecular mechanisms remain largely unknown. Here we analyzed the feeding performances and the transcriptomic responses of two BPH populations (Mugdo-BPH and TN1-BPH) during compatible (Mudog-BPH feeding on Mudgo rice) and incompatible (TN1-BPH feeding on Mudgo rice) interactions. The electrical penetration graph (EPG) results indicated that the BPH feeding and performances during the incompatible interaction are significantly affected in terms of decreased honeydew, loss of weight, decreased phloem sap ingestion (N4 waveform), but increased non-penetration (NP waveform) phase. Abundance of glucose and trehalose was reduced in BPH during the incompatible interaction. Transcriptomic surveys of insects in both interactions revealed that genes involved in cuticle formation, detoxification, metabolite transport, digestion, RNA processing, lipid or fatty acid metabolism, and proteolysis were significantly down-regulated during the incompatible interaction, whereas genes involved in insulin signaling were significantly upregulated. Knockdown of four genes, including the sugar transporter NlST45, the serine and arginine-rich protein NlSRp54, the cytochrome P450 gene NlCYP6AY1, and the cuticle protein NlCPR70 through RNA-interference revealed thess genes are important for BPH survival. Overall, the results of this study will be helpful for the future researches on BPH virulence shifts.

Molecular Characterization of Ca2+/Calmodulin-Dependent Protein Kinase II Isoforms in Three Rice Planthoppers-Nilaparvata lugens, Laodelphax striatellus, and Sogatella furcifera.

This study reports the identification of splice variants for the calcium/calmodulin-dependent protein kinase II (CaMKII) gene from Nilaparvata lugens, Laodelphax striatellus, and Sogatella furcifera. CaMKII is a multifunctional serine/threonine protein kinase that transduces Ca2+ signals in cells to control a range of cellular processes in the nervous system and muscular tissue. Sequence analysis showed that CaMKII was 99.0% identical at the amino acid level among three rice planthoppers, with the exception of a variable region located in the association domain. Four kinds of 20-81 amino acid "inserts" were found in the variable region. The phylogenetic tree of the deduced amino acid sequences showed that the NlCaMKII isoforms were more closely related to the LsCaMKII isoforms and were slightly distinct from SfCaMKII. CaMKII-E was the dominant type among the five main isoforms. CaMKII genes were constitutively expressed in various nymphal and adult stages and in tested tissues with the predominant transcription occurring in the head. There was no major tissue specificity of isoform expression, but the expression pattern and relative abundance of isoforms varied when compared with the RT-PCR between tissues. In addition, RNAi in N. lugens with dsRNA at a concentration of 200 ng nymph-1 induced a mortality of 77.7% on the 10th day and a reduction in the mRNA expression level of 67.2%. Unlike the holometabolous insect Helicoverpa armigera, the knockdown of NlCaMKII did not suppress the expression of 20E response genes, such as ECR, USP1, and HR3, in N. lugens. These results indicate that the role of CaMKII in hemimetabolous insects may be different from that in holometabolous insects.

Differential Responses of OsMPKs in IR56 Rice to Two BPH Populations of Different Virulence Levels.

The conserved mitogen-activated protein kinase (MAPK) cascades play vital roles in plant defense responses against pathogens and insects. In the current study, the expression profiles of 17 OsMPKs were determined in the TN1 and IR56 rice varieties under the infestation of brown planthopper (BPH), one of the most destructive hemimetabolous rice pests. The virulent IR56 BPH population (IR56-BPH) and the avirulent TN1 BPH population (TN-BPH) were used to reveal the roles of OsMPKs in the compatible (IR56-BPH infested on the TN1 and IR56 rice varieties, and TN1-BPH infested on the TN1 rice variety) and the incompatible (TN1-BPH infested on the IR56 rice variety) interaction. The statistical analysis revealed that rice variety, BPH population type, and infestation period have significant effects on the transcription of OsMPKs. Out of these genes, five OsMPKs (OsMPK1, OsMPK3, OsMPK7, OsMPK14, and OsMPK16) were found to exhibit upregulated expression only during incompatible interaction. Six OsMPKs (OsMPK4, OsMPK5, OsMPK8, OsMPK9, OsMPK12, and OsMPK13) were associated with both incompatible and compatible interactions. The transcription analysis of salicylic acid, jasmonic acid, and ethylene phytohormone signaling genes revealed their roles during the rice⁻BPH interactions. The upregulated expression of OsC4H, OsCHS, and OsCHI in the incompatible interaction implied the potential defense regulatory roles of phenylpropanoids. In both varieties, the elevated transcript accumulations of OsGST and OsSOD, and the increased enzyme activities of POD, SOD, and GST at 1 day post-infestation (dpi), but not at 3 dpi, indicated that reactive oxygen species (ROS) signaling might be an early event in rice⁻BPH interactions. Furthermore, upregulated transcription of OsLecRK3 and OsLecRK4 was found only during an incompatible interaction, suggesting their involvement in the BPH resistance response in the IR56 rice variety. Lastly, based on the findings of this study, we have proposed a model of interactions of IR56 rice with TN1-BPH and IR56-BPH that depicts the resistance and susceptibility reactions, respectively.

The Roles of E93 and Kr-h1 in Metamorphosis of Nilaparvata lugens.

Metamorphosis is a crucial process in insect development. Ecdysone-induced protein 93 (E93) is a determinant that promotes adult metamorphosis in both hemimetabolous and holometabolous insects. Krüppel-homolog 1 (Kr-h1), an early juvenile hormone (JH)-inducible gene, participates in JH signaling pathway controlling insect metamorphosis. In the current study, an E93 cDNA (NlE93) and two Kr-h1 cDNA variants (NlKr-h1-a and NlKr-h1-b) were cloned from Nilaparvata lugens (Stål), one of the most destructive hemimetabolous insect pests on rice. Multiple sequence alignment showed that both NlE93 and NlKr-h1 share high identity with their orthologs from other insects. The expression patterns revealed that decreasing NlKr-h1 mRNA levels were correlated with increasing NlE93 mRNA levels and vice versa. Moreover, RNA interference (RNAi) assays showed that the knockdown of one of the two genes resulted in significantly upregulated expression of the other. Correspondingly, phenotypical observation of the RNAi insects revealed that depletion of NlE93 prevented nymph-adult transition (causing a supernumerary nymphal instar), while depletion of NlKr-h1 triggered precocious formation of incomplete adult features. The results suggest that Nlkr-h1 and NlE93 are mutual repressors, fitting into the MEKRE93 pathway. The balance between these two genes plays a critical role in the metamorphosis of N. lugens determining the proper timing for activating metamorphosis during the nymphal stage.

Knockdown of a putative argininosuccinate lyase gene reduces arginine content and impairs nymphal development in Nilaparvata lugens.

Nilaparvata lugens is a typical phloem feeder. Rice phloem is high in simple sugars and very low in essential amino acids. Nilaparvata lugens harbors an ascomycete Entomomyces delphacidicola that hypothetically biosynthesizes several amino acids to meet the nutrition requirement of the planthopper. Among these amino acids, here, we focused on arginine biosynthesis. A complete cDNA of an E. delphacidicola gene, arginine-succinate lyase, EdArg4, the last step in arginine biosynthesis, was obtained. RNAi-mediated suppression of EdArg4 reduced arginine content in the hemolymph, and decreased the expression of several arginine biosynthesis genes. Silencing of EdArg4 delayed nymphal development and led to nymphal lethality. About 20% of the EdArg4 RNAi surviving adults were deformed. The most obvious defect was wider and larger abdomen. The EdArg4 RNAi-treated planthoppers had thickened wings and enlarged antennae, legs, and anal tubes and a few adults did not normally emerge. Arginine deficiency in the EdArg4 RNAi planthoppers repressed nitric oxide signaling, determined at the transcriptional level. We infer that E. delphacidicola biosynthesizes essential arginine to compensate for nutrition deficiency in N. lugens.

Reference genes for quantitative real-time PCR analysis in symbiont Entomomyces delphacidicola of Nilaparvata lugens (Stål).

Nilaparvata lugens (Stål) (Hemiptera: Delphacidae) is a major rice pest that harbors an endosymbiont ascomycete fungus, Entomomyces delphacidicola str. NLU (also known as yeast-like symbiont, YLS). Driving by demand of novel population management tactics (e.g. RNAi), the importance of YLS has been studied and revealed, which greatly boosts the interest of molecular level studies related to YLS. The current study focuses on reference genes for RT-qPCR studies related to YLS. Eight previously unreported YLS genes were cloned, and their expressions were evaluated for N. lugens samples of different developmental stages and sexes, and under different nutritional conditions and temperatures. Expression stabilities were analyzed by BestKeeper, geNorm, NormFinder, ΔCt method and RefFinder. Furthermore, the selected reference genes for RT-qPCR of YLS genes were validated using targeted YLS genes that respond to different nutritional conditions (amino acid deprivation) and RNAi. The results suggest that ylsRPS15p/ylsACT are the most suitable reference genes for temporal gene expression profiling, while ylsTUB/ylsACT and ylsRPS15e/ylsGADPH are the most suitable reference gene choices for evaluating nutrition and temperature effects. Validation studies demonstrated the advantage of using endogenous YLS reference genes for YLS studies.

A Genome-Wide Identification and Analysis of the Basic Helix-Loop-Helix Transcription Factors in Brown Planthopper, Nilaparvata lugens.

The basic helix-loop-helix (bHLH) transcription factors in insects play essential roles in multiple developmental processes including neurogenesis, sterol metabolism, circadian rhythms, organogenesis and formation of olfactory sensory neurons. The identification and function analysis of bHLH family members of the most destructive insect pest of rice, Nilaparvata lugens, may provide novel tools for pest management. Here, a genome-wide survey for bHLH sequences identified 60 bHLH sequences (NlbHLHs) encoded in the draft genome of N. lugens. Phylogenetic analysis of the bHLH domains successfully classified these genes into 40 bHLH families in group A (25), B (14), C (10), D (1), E (8) and F (2). The number of NlbHLHs with introns is higher than many other insect species, and the average intron length is shorter than those of Acyrthosiphon pisum. High number of ortholog families of NlbHLHs was found suggesting functional conversation for these proteins. Compared to other insect species studied, N. lugens has the highest number of bHLH members. Furthermore, gene duplication events of SREBP, Kn(col), Tap, Delilah, Sim, Ato and Crp were found in N. lugens. In addition, a putative full set of NlbHLH genes is defined and compared with another insect species. Thus, our classification of these NlbHLH members provides a platform for further investigations of bHLH protein functions in the regulation of N. lugens, and of insects in general.

ATP phosphoribosyltransferase from symbiont Entomomyces delphacidicola invovled in histidine biosynthesis of Nilaparvata lugens (Stål).

Histidine is an essential amino acid assumed to be synthesized by an obligatory yeast-like symbiont (Entomomyces delphacidicola str. NLU) in Nilaparvata lugens, an important rice pest. The adenosine-triphosphate phosphoribosyltransferase (ATP-PRTase) facilities the committed first step of the histidine biosynthesis pathway. In the current study, a putative ATP-PRTase was cloned and verified to be of E. delphacidicola origin (EdePRTase). The expression of the gene was spatial and temporal universal with a profile that matched the distribution of the fungal symbiont. RNA interference aided the knockdown of the EdePRTase-suppressed EdePRTase expression by 32-48 %. Hemolymph histidine level was also reduced followed by significant reduction of adult body weight. However, other performance characters including nymph development, survival, and adult sex ratio were not adversely affected by the knockdown. Furthermore, forced histidine exposure (through injection or feeding) significantly inhibited the EdePRTase mRNA levels at higher concentrations, but significantly increased EdePRTase expression levels at lower concentrations (feeding only). The significance of these findings support that the EdePRTase is from symbiont E. delphacidicola, and its involvement in histidine biosynthesis of N. lugens was discussed. The results provide a better understanding of EdePRTase and the encoded functional ATP-PRTase enzyme regulation in N. lugens and insects in general.

Ran Involved in the Development and Reproduction Is a Potential Target for RNA-Interference-Based Pest Management in Nilaparvata lugens.

Ran (RanGTPase) in insects participates in the 20-hydroxyecdysone signal transduction pathway in which downstream genes, FTZ-F1, Krüppel-homolog 1 (Kr-h1) and vitellogenin, are involved. A putative Ran gene (NlRan) was cloned from Nilaparvata lugens, a destructive phloem-feeding pest of rice. NlRan has the typical Ran primary structure features that are conserved in insects. NlRan showed higher mRNA abundance immediately after molting and peaked in newly emerged female adults. Among the examined tissues ovary had the highest transcript level, followed by fat body, midgut and integument, and legs. Three days after dsNlRan injection the NlRan mRNA abundance in the third-, fourth-, and fifth-instar nymphs was decreased by 94.3%, 98.4% and 97.0%, respectively. NlFTZ-F1 expression levels in treated third- and fourth-instar nymphs were reduced by 89.3% and 23.8%, respectively. In contrast, NlKr-h1 mRNA levels were up-regulated by 67.5 and 1.5 folds, respectively. NlRan knockdown significantly decreased the body weights, delayed development, and killed >85% of the nymphs at day seven. Two apparent phenotypic defects were observed: (1) Extended body form, and failed to molt; (2) The cuticle at the notum was split open but cannot completely shed off. The newly emerged female adults from dsNlRan injected fifth-instar nymphs showed lower levels of NlRan and vitellogenin, lower weight gain and honeydew excretion comparing with the blank control, and no offspring. Those results suggest that NlRan encodes a functional protein that was involved in development and reproduction. The study established proof of concept that NlRan could serve as a target for dsRNA-based pesticides for N. lugens control.

Identification and Function Analysis of enolase Gene NlEno1 from Nilaparvata lugens (Stål) (Hemiptera:Delphacidae).

The enolase [EC 4.2.1.11] is an essential enzyme in the glycolytic pathway catalyzing the conversion of 2-phosphoglycerate (2-PGE) to phosphoenolpyruvate (PEP). In this study, a full-length cDNA encoding α-enolase was cloned from rice brown planthopper (Nilaparvata lugens) and is provisionally designated as NlEno1. The cDNA sequence of NlEno1 was 1,851 bp with an open reading frame (ORF) of 1,305 bp and encoding 434 amino acids. The deduced protein shares high identity of 80-87% with ENO1-like protein from Hemiptera, Diptera, and Lepidoptera speices. The NlEno1 showed the highest mRNA expression level in hemolymph, followed by fat body, salivary gland, ovaries and egg, and showed trace mRNA levels in testis. The mRNA of NlEno1 showed up-regulated level in virulent N. lugens population Mudgo, IR56 and IR42 when compared with TN1 population. Injection of double-stranded RNA (dsRNA) of NlEno1 into the adults significantly down-regulated the NlEno1 mRNA level along with decreased eggs and offspring. Moreover, injection of NlEno1-dsRNA decreased mRNA level of Vitellogenin (Vg) gene. These results showed that the NlEno1, as a key glycolytic enzyme, may play roles in regulation of fecundity and adaptation of N. lugens to resistant rice varieties.