Complete genome sequence of a novel amalgavirus in sponge gourd, Luffa cylindrica.

A novel monopartite dsRNA virus, tentatively named "sponge gourd amalgavirus 1" (SGAV1), was discovered by high-throughput sequencing in sponge gourd (Luffa cylindrica) displaying mosaic symptoms in Jiashan County, Zhejiang Province, China. The genome of SGAV1 is 3,447 nucleotides in length and contains partially overlapping open reading frames (ORFs) encoding a putative replication factory matrix-like protein and a fusion protein, respectively. The fusion protein of SGAV1 shares 57.07% identity with the homologous protein of salvia miltiorrhiza amalgavirus 1 (accession no. DAZ91057.1). Phylogenetic analysis based on the RNA-dependent RNA polymerase (RdRp) protein suggests that SGAV1 belongs to the genus Amalgavirus of the family Amalgaviridae. Moreover, analysis of SGAV1-derived small interfering RNAs indicated that SGAV1 was actively replicating in the host plant. Semi-quantitative RT-PCR showed higher levels of SGAV1 expression in leaves than in flowers and fruits. This is the first report of a novel amalgavirus found in sponge gourd in China.

Complete genome sequence and genetic characterization of a novel segmented RNA virus infecting Nilaparvata lugens.

The brown planthopper (BPH), Nilaparvata lugens, is a significant agricultural pest capable of long-distance migration and transmission of viruses that cause severe disease in rice. In this study, we identified a novel segmented RNA virus in a BPH, and this virus exhibited a close relationship to members of a recently discovered virus lineage known as "quenyaviruses" within the viral kingdom Orthornavirae. This newly identified virus was named "Nilaparvata lugens quenyavirus 1" (NLQV1). NLQV1 consists of five positive-sense, single-stranded RNAs, with each segment containing a single open reading frame (ORF). The genomic characteristics and phylogenetic analysis support the classification of NLQV1 as a novel quenyavirus. Notably, all of the genome segments of NLRV contained the 5'-terminal sequence AUCUG. The characteristic virus-derived small interfering RNA (vsiRNA) profile of NLQV1 suggests that the antiviral RNAi pathway of the host BPH was activated in response to virus infection. These findings represent the first documented report of quenyaviruses in planthoppers, contributing to our understanding of quenyaviruses and expanding our knowledge of insect-specific viruses in planthoppers.

Identification and Characterization of Three Novel Solemo-like Viruses in the White-Backed Planthopper, Sogatella furcifera.

Agricultural insects play a crucial role in transmitting plant viruses and host a considerable number of insect-specific viruses (ISVs). Among these insects, the white-backed planthoppers (WBPH; Sogatella furcifera, Hemiptera: Delphacidae) are noteworthy rice pests and are responsible for disseminating the southern rice black-streaked dwarf virus (SRBSDV), a significant rice virus. In this study, we analyzed WBPH transcriptome data from public sources and identified three novel viruses. These newly discovered viruses belong to the plant-associated viral family Solemoviridae and were tentatively named Sogatella furcifera solemo-like virus 1-3 (SFSolV1-3). Among them, SFSolV1 exhibited a prevalent existence in different laboratory populations, and its complete genome sequence was obtained using rapid amplification of cDNA ends (RACE) approaches. To investigate the antiviral RNA interference (RNAi) response in WBPH, we conducted an analysis of virus-derived small interfering RNAs (vsiRNAs). The vsiRNAs of SFSolV1 and -2 exhibited typical patterns associated with the host's siRNA-mediated antiviral immunity, with a preference for 21- and 22-nt vsiRNAs derived equally from both the sense and antisense genomic strands. Furthermore, we examined SFSolV1 infection and distribution in WBPH, revealing a significantly higher viral load of SFSolV1 in nymphs' hemolymph compared to other tissues. Additionally, in adult insects, SFSolV1 exhibited higher abundance in male adults than in female adults.

Discovery of Two Novel Viruses of the Willow-Carrot Aphid, Cavariella aegopodii.

The advancement of bioinformatics and sequencing technology has resulted in the identification of an increasing number of new RNA viruses. This study systematically identified the RNA virome of the willow-carrot aphid, Cavariella aegopodii (Hemiptera: Aphididae), using metagenomic sequencing and rapid amplification of cDNA ends (RACE) approaches. C. aegopodii is a sap-sucking insect widely distributed in Europe, Asia, North America, and Australia. The deleterious effects of C. aegopodii on crop growth primarily stem from its feeding activities and its role as a vector for transmitting plant viruses. The virome includes Cavariella aegopodii virga-like virus 1 (CAVLV1) and Cavariella aegopodii iflavirus 1 (CAIV1). Furthermore, the complete genome sequence of CAVLV1 was obtained. Phylogenetically, CAVLV1 is associated with an unclassified branch of the Virgaviridae family and is susceptible to host antiviral RNA interference (RNAi), resulting in the accumulation of a significant number of 22nt virus-derived small interfering RNAs (vsiRNAs). CAIV1, on the other hand, belongs to the Iflaviridae family, with vsiRNAs ranging from 18 to 22 nt. Our findings present a comprehensive analysis of the RNA virome of C. aegopodii for the first time, offering insights that could potentially aid in the future control of the willow-carrot aphid.

PTTH-Torso Signaling System Controls Developmental Timing, Body Size, and Reproduction through Regulating Ecdysone Homeostasis in the Brown Planthopper, Nilaparvata lugens.

In holometabolous insects, such as Drosophila and Bombyx, prothoracicotropic hormone (PTTH) is well established to be critical in controlling developmental transitions and metamorphosis by stimulating the biosynthesis of ecdysone in the prothoracic glands (PGs). However, the physiological role of PTTH and the receptor Torso in hemimetabolous insects remains largely unexplored. In this study, homozygous PTTH- and Torso-null mutants of the brown planthopper (BPH), Nilaparvata lugens, were successfully generated by employing clustered regularly interspaced short palindromic repeats/CRISPR-associated 9 (CRISPR-Cas9). Further characterization showed that both NlPTTH-/- and NlTorso-/- mutants exhibited prolonged nymphal duration and increased final adult size. Enzyme-linked immunosorbent assay (ELISA) revealed that NlPTTH-/- and NlTorso-/- mutants exhibited a significant reduction in 20-hydroxyecdysone (20E) in fifth-instar nymphs at 48 h post-ecdysis compared to Wt controls. Furthermore, our results indicated that both NlPTTH-/- and NlTorso-/- mutants had shortened lifespan, reduced female fecundity, and reduced egg hatching rates in adults. These findings suggest a conserved role for the PTTH-Torso signaling system in the regulation of developmental transitions by stimulating ecdysone biosynthesis in hemimetabolous insects.

Musashi orchestrates melanism in Laodelphax striatellus.

In insects, melanism, a fundamental pigmentation process, is of significant importance in evolutionary biology due to its complex genetic foundation. We investigated the role of the RNA-binding gene Musashi (msi) in melanism in Laodelphax striatellus, a Hemiptera species. We identified a single L. striatellus msi homolog, Lsmsi, encoding a 357 amino acid protein with 2 RNA recognition motifs. RNA interference-mediated knockdown of LsMsi resulted in complete body melanism and increased cuticular permeability. Additionally, we found the involvement of G protein-coupled receptor A42 and tyrosine hydroxylase (Th) in L. striatellus melanism. Knockdown of LsTh lightened the epidermis, showing dehydration signs, while LsA42 knockdown enhanced LsTh expression, leading to melanism. Surprisingly, Lsmsi knockdown decreased both LsA42 and LsTh expression, which was expected to cause whitening but resulted in melanism. Further, we found that Lsmsi influenced downstream genes like phenoloxidase homolog LsPo and dopa decarboxylase (Ddc) homolog LsDdc in the tyrosine-mediated melanism pathway. Extending to Nilaparvata lugens and Sogatella furcifera, we demonstrated the conserved role of msi in melanism among Delphacidae. Given MSI proteins' roles in cancer and tumors in vertebrates, our study is the first to link msi in insects to Delphacidae body color melanization via the tyrosine-mediated pathway, offering fresh perspectives on the genetic basis of insect melanism and msi functions.

Complete genome analysis of a novel narnavirus in sweet viburnum (Viburnum odoratissimum).

Trees and shrubs provide important ecological services. However, few studies have surveyed the virome in trees and shrubs. In this study, we discovered a new positive-sense RNA virus originating from Viburnum odoratissimum, which we named "Vo narna-like virus". The complete genome of Vo narna-like virus is 3,451 nt in length with an open reading frame (ORF) encoding the RNA-dependent RNA polymerase (RdRP) protein. Phylogenetic analysis placed this virus within the betanarnavirus clade, sharing 53.63% amino acid sequence identity with its closest relative, Qingdao RNA virus 2. The complete sequence of the virus was confirmed by rapid amplification of cDNA ends (RACE) and Sanger sequencing. Small interfering RNA (siRNA) analysis indicated that this virus interacts with the RNA interference (RNAi) pathway of V. odoratissimum. This is the first report of a narnavirus in V. odoratissimum.

Identification and Functional Analysis of the fruitless Gene in a Hemimetabolous Insect, Nilaparvata lugens.

The fruitless (fru) gene functions as a crucial "tuner" in male insect courtship behavior through distinct expression patterns. In Nilaparvata lugens, our previous research showed doublesex (dsx) influencing male courtship songs, causing mating failures with virgin females. However, the impact of fru on N. lugens mating remains unexplored. In this study, the fru homolog (Nlfru) in N. lugens yielded four spliceosomes: Nlfru-374-a/b, Nlfru-377, and Nlfru-433, encoding proteins of 374aa, 377aa, and 433aa, respectively. Notably, only Nlfru-374b exhibited male bias, while the others were non-sex-specific. All NlFRU proteins featured the BTB conserved domain, with NlFRU-374 and NlFRU-377 possessing the ZnF domain with different sequences. RNAi-mediated Nlfru or its isoforms' knockdown in nymph stages blocked wing-flapping behavior in mating males, while embryonic knockdown via maternal RNAi resulted in over 80% of males losing wing-flapping ability, and female receptivity was reduced. Nlfru expression was Nldsx-regulated, and yet courtship signals and mating success were unaffected. Remarkably, RNAi-mediated Nlfru knockdown up-regulated the expression of flightin in macropterous males, which regulated muscle stiffness and delayed force response, suggesting Nlfru's involvement in muscle development regulation. Collectively, our results indicate that Nlfru functions in N. lugens exhibit a combination of conservation and species specificity, contributing insights into fru evolution, particularly in Hemiptera species.

The genomic history and global migration of a windborne pest.

Many insect pests, including the brown planthopper (BPH), undergo windborne migration that is challenging to observe and track. It remains controversial about their migration patterns and largely unknown regarding the underlying genetic basis. By analyzing 360 whole genomes from around the globe, we clarify the genetic sources of worldwide BPHs and illuminate a landscape of BPH migration showing that East Asian populations perform closed-circuit journeys between Indochina and the Far East, while populations of Malay Archipelago and South Asia undergo one-way migration to Indochina. We further find round-trip migration accelerates population differentiation, with highly diverged regions enriching in a gene desert chromosome that is simultaneously the speciation hotspot between BPH and related species. This study not only shows the power of applying genomic approaches to demystify the migration in windborne migrants but also enhances our understanding of how seasonal movements affect speciation and evolution in insects.

Molecular and biological characterization of a bunyavirus infecting the brown planthopper (Nilaparvata lugens).

A negative-strand symbiotic RNA virus, tentatively named Nilaparvata lugens Bunyavirus (NLBV), was identified in the brown planthopper (BPH, Nilaparvata lugens). Phylogenetic analysis indicated that NLBV is a member of the genus Mobuvirus (family Phenuiviridae, order Bunyavirales). Analysis of virus-derived small interfering RNA suggested that antiviral immunity of BPH was successfully activated by NLBV infection. Tissue-specific investigation showed that NLBV was mainly accumulated in the fat-body of BPH adults. Moreover, NLBV was detected in eggs of viruliferous female BPHs, suggesting the possibility of vertical transmission of NLBV in BPH. Additionally, no significant differences were observed for the biological properties between NLBV-infected and NLBV-free BPHs. Finally, analysis of geographic distribution indicated that NLBV may be prevalent in Southeast Asia. This study provided a comprehensive characterization on the molecular and biological properties of a symbiotic virus in BPH, which will contribute to our understanding of the increasingly discovered RNA viruses in insects.

Cross-kingdom RNA interference mediated by insect salivary microRNAs may suppress plant immunity.

Communication between insects and plants relies on the exchange of bioactive molecules that traverse the species interface. Although proteinic effectors have been extensively studied, our knowledge of other molecules involved in this process remains limited. In this study, we investigate the role of salivary microRNAs (miRNAs) from the rice planthopper Nilaparvata lugens in suppressing plant immunity. A total of three miRNAs were confirmed to be secreted into host plants during insect feeding. Notably, the sequence-conserved miR-7-5P is specifically expressed in the salivary glands of N. lugens and is secreted into saliva, distinguishing it significantly from homologues found in other insects. Silencing miR-7-5P negatively affects N. lugens feeding on rice plants, but not on artificial diets. The impaired feeding performance of miR-7-5P-silenced insects can be rescued by transgenic plants overexpressing miR-7-5P. Through target prediction and experimental testing, we demonstrate that miR-7-5P targets multiple plant genes, including the immune-associated bZIP transcription factor 43 (OsbZIP43). Infestation of rice plants by miR-7-5P-silenced insects leads to the increased expression of OsbZIP43, while the presence of miR-7-5P counteracts this upregulation effect. Furthermore, overexpressing OsbZIP43 confers plant resistance against insects which can be subverted by miR-7-5P. Our findings suggest a mechanism by which herbivorous insects have evolved salivary miRNAs to suppress plant immunity, expanding our understanding of cross-kingdom RNA interference between interacting organisms.

Silencing NlFAR7 destroyed the pore canals and related structures of the brown planthopper.

Fatty acyl-CoA reductase (FAR) is one of the key enzymes, which catalyses the conversion of fatty acyl-CoA to the corresponding alcohols. Among the FAR family members in the brown planthopper (Nilaparvata lugens), NlFAR7 plays a pivotal role in both the synthesis of cuticular hydrocarbons and the waterproofing of the cuticle. However, the precise mechanism by which NlFAR7 influences the formation of the cuticle structure in N. lugens remains unclear. Therefore, this paper aims to investigate the impact of NlFAR7 through RNA interference, transmission electron microscope, focused ion beam scanning electron microscopy (FIB-SEM) and lipidomics analysis. FIB-SEM is employed to reconstruct the three-dimensional (3D) architecture of the pore canals and related cuticle structures in N. lugens subjected to dsNlFAR7 and dsGFP treatments, enabling a comprehensive assessment of changes in the cuticle structures. The results reveal a reduction in the thickness of the cuticle and disruptions in the spiral structure of pore canals, accompanied by widened base and middle diameters. Furthermore, the lipidomics comparison analysis between dsNlFAR7- and dsGFP-treated N. lugens demonstrated that there were 25 metabolites involved in cuticular lipid layer synthesis, including 7 triacylglycerols (TGs), 5 phosphatidylcholines (PCs), 3 phosphatidylethanolamines (PEs) and 2 diacylglycerols (DGs) decreased, and 4 triacylglycerols (TGs) and 4 PEs increased. In conclusion, silencing NlFAR7 disrupts the synthesis of overall lipids and destroys the cuticular pore canals and related structures, thereby disrupting the secretion of cuticular lipids, thus affecting the cuticular waterproofing of N. lugens. These findings give significant attention with reference to further biochemical researches on the substrate specificity of FAR protein, and the molecular regulation mechanisms during N. lugens life cycle.

Salivary proteins potentially derived from horizontal gene transfer are critical for salivary sheath formation and other feeding processes.

Herbivorous insects employ an array of salivary proteins to aid feeding. However, the mechanisms behind the recruitment and evolution of these genes to mediate plant-insect interactions remain poorly understood. Here, we report a potential horizontal gene transfer (HGT) event from bacteria to an ancestral bug of Eutrichophora. The acquired genes subsequently underwent duplications and evolved through co-option. We annotated them as horizontal-transferred, Eutrichophora-specific salivary protein (HESPs) according to their origin and function. In Riptortus pedestris (Coreoidea), all nine HESPs are secreted into plants during feeding. The RpHESP4 to RpHESP8 are recently duplicated and found to be indispensable for salivary sheath formation. Silencing of RpHESP4-8 increases the difficulty of R. pedestris in probing the soybean, and the treated insects display a decreased survivability. Although silencing the other RpHESPs does not affect the salivary sheath formation, negative effects are also observed. In Pyrrhocoris apterus (Pyrrhocoroidea), five out of six PaHESPs are secretory salivary proteins, with PaHESP3 being critical for insect survival. The PaHESP5, while important for insects, no longer functions as a salivary protein. Our results provide insight into the potential origin of insect saliva and shed light on the evolution of salivary proteins.

Inhibition of Rice Stripe Virus Accumulation by Polyubiquitin-C in Laodelphax striatellus.

Many hosts utilize the ubiquitin system to defend against viral infection. As a key subunit of the ubiquitin system, the role of polyubiquitin in the viral infection of insects is unclear. Here, we identified the full-length cDNA of the polyubiquitin-C (UBC) gene in Laodelphax striatellus, the small brown planthopper (SBPH). LsUBC was expressed in various tissues and was highly expressed in salivary glands, midgut, and reproductive systems. Furthermore, the LsUBC expression profiles in the developmental stages showed that LsUBC was ubiquitously expressed in seven developmental stages and was highest expressed in female adults with SBPH. qRT-PCR analyses indicated that rice stripe virus (RSV) infection promoted the LsUBC expression. Knockdown of LsUBC mRNA via RNA interference increased RSV accumulation. These findings suggest that LsUBC inhibits RSV accumulation in L. striatellus.

Maintenance of persistent transmission of a plant arbovirus in its insect vector mediated by the Toll-Dorsal immune pathway.

Throughout evolution, arboviruses have developed various strategies to counteract the host's innate immune defenses to maintain persistent transmission. Recent studies have shown that, in addition to bacteria and fungi, the innate Toll-Dorsal immune system also plays an essential role in preventing viral infections in invertebrates. However, whether the classical Toll immune pathway is involved in maintaining the homeostatic process to ensure the persistent and propagative transmission of arboviruses in insect vectors remain unclear. In this study, we revealed that the transcription factor Dorsal is actively involved in the antiviral defense of an insect vector (Laodelphax striatellus) by regulating the target gene, zinc finger protein 708 (LsZN708), which mediates downstream immune-related effectors against infection with the plant virus (Rice stripe virus, RSV). In contrast, an antidefense strategy involving the use of the nonstructural-protein (NS4) to antagonize host antiviral defense through competitive binding to Dorsal from the MSK2 kinase was employed by RSV; this competitive binding inhibited Dorsal phosphorylation and reduced the antiviral response of the host insect. Our study revealed the molecular mechanism through which Toll-Dorsal-ZN708 mediates the maintenance of an arbovirus homeostasis in insect vectors. Specifically, ZN708 is a newly documented zinc finger protein targeted by Dorsal that mediates the downstream antiviral response. This study will contribute to our understanding of the successful transmission and spread of arboviruses in plant or invertebrate hosts.

Molecular Characterization of an Isolate of Bean Common Mosaic Virus First Identified in Gardenia Using Metatranscriptome and Small RNA Sequencing.

Gardenia (Gardenia jasminoides) is a popular and economically vital plant known for its ornamental and medicinal properties. Despite its widespread cultivation, there has been no documentation of plant viruses on gardenia yet. In the present study, gardenia leaves exhibiting symptoms of plant viral diseases were sampled and sequenced by both metatranscriptome and small RNA sequencing. As a consequence, bean common mosaic virus (BCMV) was identified in gardenia for the first time and named BCMV-gardenia. The full genome sequence of BCMV-gardenia is 10,054 nucleotides (nt) in length (excluding the poly (A) at the 3' termini), encoding a large polyprotein of 3,222 amino acids. Sequence analysis showed that the N-termini of the polyprotein encoded by BCMV-gardenia is less conserved when compared to other BCMV isolates, whereas the C-termini is the most conserved. Maximum likelihood phylogenetic analysis showed that BCMV-gardenia was clustered closely with other BCMV isolates identified outside the leguminous plants. Our results indicated that the majority of BCMV-gardenia virus-derived small interfering RNAs (vsiRNAs) were 21 nt and 22 nt, with 21 nt being more abundant. The first nucleotide at the 5' termini of vsiRNAs derived from BCMV-gardenia preferred U and A. The ratio of vsiRNAs derived from sense (51.1%) and antisense (48.9%) strands is approaching, and the distribution of vsiRNAs along the viral genome is generally even, with some hot spots forming in local regions. Our findings could provide new insights into the diversity, evolution, and host expansion of BCMV and contribute to the prevention and treatment of this virus.

Comparative transcriptomic analysis of salivary glands between the zoophytophagous Cyrtorhinus lividipennis and the phytozoophagous Apolygus lucorum.

BACKGROUND: Saliva plays a crucial role in shaping the feeding behavior of insects, involving processes such as food digestion and the regulation of interactions between insects and their hosts. Cyrtorhinus lividipennis serves as a predominant natural enemy of rice pests, while Apolygus lucorum, exhibiting phytozoophagous feeding behavior, is a destructive agricultural pest. In this study, a comparative transcriptome analysis, incorporating the published genomes of C.lividipennis and A.lucorum, was conducted to reveal the role of salivary secretion in host adaptation. RESULTS: In contrast to A.lucorum, C.lividipennis is a zoophytophagous insect. A de novo genome analysis of C.lividipennis yielded 19,706 unigenes, including 16,217 annotated ones. On the other hand, A.lucorum had altogether 20,111 annotated genes, as obtained from the published official gene set (20,353 unigenes). Functional analysis of the top 1,000 salivary gland (SG)-abundant genes in both insects revealed that the SG was a dynamically active tissue engaged in protein synthesis and secretion. Predictions of other tissues and signal peptides were compared. As a result, 94 and 157 salivary proteins were identified in C.lividipennis and A.lucorum, respectively, and were categorized into 68 and 81 orthogroups. Among them, 26 orthogroups were shared, potentially playing common roles in digestion and detoxification, including several venom serine proteases. Furthermore, 42 and 55 orthogroups were exclusive in C.lividipennis and A.lucorum, respectively, which were exemplified by a hyaluronidase in C.lividipennis that was associated with predation, while polygalacturonases in A.lucorum were involved in mesophyll-feeding patterns. CONCLUSIONS: Findings in this study provide a comprehensive insight into saliva secretions in C.lividipennis and A.lucorum via a transcriptome approach, reflecting the intricate connections between saliva secretions and feeding behaviors. It is found that conserved salivary secretions are involved in shaping the overlapping feeding patterns, while a plethora of unique salivary secretions may drive the evolution of specific feeding behaviors crucial for their survival. These results enhance our understanding of the feeding mechanisms in different insects from the perspective of saliva and contribute to future environmentally friendly pest control by utilizing predatory insects.

Complete genome sequence of a novel robigovirus infecting Mentha arvensis.

The complete genomic sequence of a novel robigovirus, provisionally named "Mentha arvensis robigovirus 1" (MARV1), was determined by combining next-generation sequencing (NGS), reverse transcription polymerase chain reaction (RT-PCR), and rapid amplification of cDNA ends (RACE) PCR. The complete genomic sequence of this new virus is 7617 nucleotides in length, excluding the 3' poly(A) tail. The MARV1 genome encodes a putative replicase, "triple gene block" proteins, and a coat protein. Phylogenetic analysis demonstrated that MARV1 is a member of the genus Robigovirus, with closest relationships to African oil palm ringspot virus (AOPRV). Furthermore, MARV1-derived small interfering RNAs (siRNAs) showed typical patterns of plant-virus-derived siRNAs produced by the host antiviral RNA interference pathway. This is the first report of a plant virus of the genus Robigovirus in M. arvensis.

Biosynthesis of eriodictyol in citrus waster by endowing P450BM3 activity of naringenin hydroxylation.

The flavonoid naringenin is abundantly present in pomelo peels, and the unprocessed naringenin in wastes is not friendly for the environment once discarded directly. Fortunately, the hydroxylated product of eriodictyol from naringenin exhibits remarkable antioxidant and anticancer properties. The P450s was suggested promising for the bioconversion of the flavonoids, but less naturally existed P450s show hydroxylation activity to C3' of the naringenin. By well analyzing the catalytic mechanism and the conformations of the naringenin in P450, we proposed that the intermediate Cmpd I ((porphyrin)Fe = O) is more reasonable as key conformation for the hydrolyzation, and the distance between C3'/C5' of naringenin to the O atom of CmpdI determines the hydroxylating activity for the naringenin. Thus, the "flying kite model" that gradually drags the C-H bond of the substrate to the O atom of CmpdI was put forward for rational design. With ab initio design, we successfully endowed the self-sufficient P450-BM3 hydroxylic activity to naringenin and obtained mutant M5-5, with kcat, Km, and kcat/Km values of 230.45 min-1, 310.48 µM, and 0.742 min-1 µM-1, respectively. Furthermore, the mutant M4186 was screened with kcat/Km of 4.28-fold highly improved than the reported M13. The M4186 also exhibited 62.57% yield of eriodictyol, more suitable for the industrial application. This study provided a theoretical guide for the rational design of P450s to the nonnative compounds. KEY POINTS: •The compound I is proposed as the starting point for the rational design of the P450BM3 •"Flying kite model" is proposed based on the distance between O of Cmpd I and C3'/C5' of naringenin •Mutant M15-5 with 1.6-fold of activity than M13 was obtained by ab initio modification.