Berberine Bridge Research Articles

Convergent evolution of berberine biosynthesis.

Berberine is an effective antimicrobial and antidiabetic alkaloid, primarily extracted from divergent botanical lineages, specifically Coptis (Ranunculales, early-diverging eudicot) and Phellodendron (Sapindales, core eudicot). In comparison with its known pathway in Coptis species, its biosynthesis in Phellodendron species remains elusive. Using chromosome-level genome assembly, coexpression matrix, and biochemical assays, we identified six key steps in berberine biosynthesis from Phellodendron amurense, including methylation, hydroxylation, and berberine bridge formation. Notably, we discovered a specific class of O-methyltransferases (NOMT) responsible for N-methylation. Structural analysis and mutagenesis of PaNOMT9 revealed its unique substrate-binding conformation. In addition, unlike the classical FAD-dependent berberine bridge formation in Ranunculales, Phellodendron uses a NAD(P)H-dependent monooxygenase (PaCYP71BG29) for berberine bridge formation, originating from the neofunctionalization of tryptamine 5-hydroxylase. Together, these findings reveal the convergence of berberine biosynthesis between Coptis and Phellodendron and signify the role of the convergent evolution in plant specialized metabolisms.

De novo production of protoberberine and benzophenanthridine alkaloids through metabolic engineering of yeast

Protoberberine alkaloids and benzophenanthridine alkaloids (BZDAs) are subgroups of benzylisoquinoline alkaloids (BIAs), which represent a diverse class of plant-specialized natural metabolites with many pharmacological properties. Microbial biosynthesis has been allowed for accessibility and scalable production of high-value BIAs. Here, we engineer Saccharomyces cerevisiae to de novo produce a series of protoberberines and BZDAs, including palmatine, berberine, chelerythrine, sanguinarine and chelirubine. An ER compartmentalization strategy is developed to improve vacuole protein berberine bridge enzyme (BBE) activity, resulting in >200% increase on the production of the key intermediate (S)-scoulerine. Another promiscuous vacuole protein dihydrobenzophenanthridine oxidase (DBOX) has been identified to catalyze two-electron oxidation on various tetrahydroprotoberberines at N7-C8 position and dihydrobenzophenanthridine alkaloids. Furthermore, cytosolically expressed DBOX can alleviate the limitation on BBE. This study highlights the potential of microbial cell factories for the biosynthesis of a diverse group of BIAs through engineering of heterologous plant enzymes.

Impact of nitrogen rates on biosynthesis pathways: A comparative study of diterpene synthases in clerodane diterpenoids and enzymes in benzylisoquinoline alkaloids

Tinospora sagittata is rich in secondary metabolites used in traditional medicine. However, environmental factors impact key enzymes in metabolite synthesis, highlighting the need for improved growth conditions. This study employs transcriptomics and metabolomics to assess nitrogen's impact on enzymes in secondary metabolites biosynthesis pathways. The gene expressions of berberine bridge enzymes (BBEs) like TsBBE2 had peak expression in low nitrogen treatments (A0 and A1) but were absent in higher nitrogen treatments (A2 and A3). Similar trends were observed for other enzymes such as (S)-scoulerine 9-O-methyltransferase (TsCMT3), Tetrahydroberberine oxidase (TsSTOX), and Columbamine O-methyltransferase (TsCoCOMT2-4) in response to nitrogen levels. In examining gene families related to diterpene synthases (diTPS), 1-deoxyxylulose 5-phosphate synthase (TsDXR1) expression increased with higher nitrogen fertilizer, while TsDXR2 peaked at maximal nitrogen levels. Geranylgeranyl diphosphate synthase (TsGGPP3 and TsGGPP5) decreased with nitrogen levels. (−)-kolavenyl diphosphate synthase (KPS) genes had higher expression in treatments, while ent-kaurene synthase (KSL) genes, especially TsKSL1 and TsKSL2, showed higher expression in control conditions with lower nitrogen fertilizer. Metabolite analysis confirmed more upregulated compounds in A3 compared to A0. These findings have practical implications for agriculture and pharmaceuticals, highlighting the link between nitrogen fertilization and specialized metabolism in medicinal plants.

MFS Transporter as the Molecular Switch Unlocking the Production of Cage-Like Acresorbicillinol C.

Sorbicillinoids are a class of fungal polyketides with diverse structures and distinguished bioactivities. Although remarkable progress has been achieved in their chemistry and biosynthesis, the efflux of sorbicillinoids is poorly understood. Here, we found MFS transporter AcsorT was responsible for the biosynthesis of sorbicillinoids in Acremonium chrysogenum. Combinatorial knockout and subcellular location demonstrated that the plasma membrane-associated AcsorT was responsible for the transportation of sorbicillinol and subsequent formation of oxosorbicillinol and acresorbicillinol C via the berberine bridge enzyme-like oxidase AcsorD in the periplasm. Homology modeling and site-directed mutation revealed that Tyr303 and Arg436 were the key residues of AcsorT, which was further explained by molecular dynamics simulation. Based on our study, it was suggested that AcsorT modulates sorbicillinoid production by coordinating its biosynthesis and export, and a transport model of sorbicillinoids was proposed in A. chrysogenum.

Biosynthetic Strategies of Berberine Bridge Enzyme-like Flavoprotein Oxidases toward Structural Diversification in Natural Product Biosynthesis.

Berberine bridge enzyme-like oxidases are often involved in natural product biosynthesis and are seen as essential enzymes for the generation of intricate pharmacophores. These oxidases have the ability to transfer a hydride atom to the FAD cofactor, which enables complex substrate modifications and rearrangements including (intramolecular) cyclizations, carbon-carbon bond formations, and nucleophilic additions. Despite the diverse range of activities, the mechanistic details of these reactions often remain incompletely understood. In this Review, we delve into the complexity that BBE-like oxidases from bacteria, fungal, and plant origins exhibit by providing an overview of the shared catalytic features and emphasizing the different reactivities. We propose four generalized modes of action by which BBE-like oxidases enable the synthesis of natural products, ranging from the classic alcohol oxidation reactions to less common amine and amide oxidation reactions. Exploring the mechanisms utilized by nature to produce its vast array of natural products is a subject of considerable interest and can lead to the discovery of unique biochemical activities.

A Secreted BBE-Like Enzyme Acting as a Drug-Binding Efflux Carrier Confers Microbial Self-Resistance to Mitomycin C.

The berberine bridge enzyme (BBE)-like flavoproteins have attracted continuous attention for their capability to catalyze various oxidative reactions. Here we demonstrate that MitR, a secreted BBE-like enzyme, functions as a special drug-binding efflux protein evolved from quinone reductase. Moreover, this protein provides self-resistance to its hosts toward the DNA-alkylating agent mitomycin C with a distinctive strategy, featured by independently performing drug binding and efflux.

A Berberine Bridge Enzyme-like Oxidase Mediates the Cage-like Acresorbicillinol C Biosynthesis.

Oxosorbicillinol and cage-like acresorbicillinol C are bioactive sorbicillinoids produced by Acremonium chrysogenum. We found that a berberine bridge enzyme-like oxidase AcsorD was responsible for their biosynthesis by gene deletion and heterologous expression. AcsorD catalyzed oxidation of sorbicillinol to form oxosorbicillinol in in vitro assays, which was successively condensed with sorbicillinol to form acresorbicillinol C spontaneously. Finally, site-directed mutation revealed that Tyr525 was the key residue in the catalysis of the oxidation reaction and unlocking cage-like acresorbicillinol C production.

Genome-Wide Characterization of Berberine Bridge Enzyme Gene Family in Wheat (Triticum aestivum L.) and the Positive Regulatory Role of TaBBE64 in Response to Wheat Stripe Rust.

Berberine bridge enzymes (BBEs), functioning as carbonate oxidases, enhance disease resistance in Arabidopsis and tobacco. However, the understanding of BBEs' role in monocots against pathogens remains limited. This study identified 81 TaBBEs with FAD_binding_4 and BBE domains. Phylogenetic analysis revealed a separation of the BBE gene family between monocots and dicots. Notably, RNA-seq showed TaBBE64's significant induction by both pathogen-associated molecular pattern treatment and Puccinia striiformis f. sp. tritici (Pst) infection at early stages. Subcellular localization revealed TaBBE64 at the cytoplasmic membrane. Knocking down TaBBE64 compromised wheat's Pst resistance, reducing reactive oxygen species and promoting fungal growth, confirming its positive role. Molecular docking and enzyme activity assays confirmed TaBBE64's glucose oxidation to produce H2O2. Since Pst relies on living wheat cells for carbohydrate absorption, TaBBE64's promotion of glucose oxidation limits fungal growth and resists pathogen infection. This study sheds light on BBEs' role in wheat resistance against biotrophic fungi.

“Stepwise or concerted or a hybrid mechanism for berberine bridge enzyme?”: A computational mechanistic study by QM-MM and QM methods

Berberine bridge enzyme (BBE) is a plant-based amine oxidase that catalyzes conversion of (S)-reticuline into (S)-scoulerine using flavin as cofactor in a stereospecific way in an alkaloid biosynthesis pathway. Based on the active site enzyme variants, a concerted mechanism was proposed involving hydride transfer, proton transfer, and substrate cyclization processes in a single step. In this mechanism, Glu417 residue acts as the catalytic base which deprotonates the phenolic proton of the substrate while a hydride ion transfers from the substrate to flavin and the substrate cyclizes. However, based on solvent and substrate deuterium kinetic isotope effect studies, it was proposed that the oxidation process occurs in a stepwise fashion in which a hydride ion transfer from substrate to flavin first, then cyclization of the substrate occurs together with the proton transfer process. In this study, we formulated computational models to elucidate the oxidation mechanism of (S)-reticuline into (S)-scoulerine using the crystal structure of enzyme complexed with (S)-reticuline. Both QM and QM-MM calculations revealed that a hybrid of concerted and stepwise mechanisms might be operative during the catalysis. It was found that a concerted hydride-proton transfer processes occurs forming a reactive intermediate which subsequently cyclize without an energy barrier as a decoupled step.

Mapping of prehaustorial resistance against wheat leaf rust in einkorn (Triticum monococcum), a progenitor of wheat

Wheat leaf rust (Puccinia triticina) is one of the most significant fungal diseases of wheat, causing substantial yield losses worldwide. Infestation is currently being reduced by fungicide treatments and mostly vertical resistance. However, these measures often break down when the fungal virulence pattern changes, resulting in a breakdown of vertical resistances. In contrast, the prehaustorial resistance (phr) that occurs in the einkorn–wheat leaf rust interaction is race-independent, characterized by an early defense response of plants during the prehaustorial phase of infestation. Einkorn (Triticum monococcum) is closely related to Triticum urartu as a progenitor of wheat and generally shows a high level of resistance against leaf rust of wheat. Hence, einkorn can serve as a valuable source to improve the level of resistance to the pathogen in future wheat lines. In particular, einkorn accession PI272560 is known to exhibit a hypersensitive prehaustorial effector triggered immune reaction, preventing the infection of P. triticina. Remarkably, this effector-triggered immune reaction turned out to be atypical as it is non-race-specific (horizontal). To genetically dissect the prehaustorial resistance (phr) in PI272560, a biparental F2 population of 182 plants was established after crossing PI272560 with the susceptible T. boeoticum accession 36554. Three genetic maps comprising 2,465 DArT-seq markers were constructed, and a major QTL was detected on chromosome 5A. To locate underlying candidate genes, marker sequences flanking the respective QTL were aligned to the T. urartu reference genome and transcriptome data available from the parental accessions were used. Within the QTL interval of approximately 16.13 million base pairs, the expression of genes under inoculated and non-inoculated conditions was analyzed via a massive analysis of cDNA (MACE). Remarkably, a single gene located 3.4 Mbp from the peak marker within the major QTL was upregulated (20- to 95-fold) after the inoculation in the resistant accession in comparison to the susceptible T. boeoticum accession. This gene belongs to a berberine bridge enzyme-like protein that is suspected to interact on the plant surface with glycoside hydrolases (GH) secreted by the fungus and to induce a hypersensitive defense reaction in the plant after fungal infections.

Berberine bridge enzyme-like oxidases of cellodextrins and mixed-linked β-glucans control seed coat formation.

Plants have evolved various resistance mechanisms to cope with biotic stresses that threaten their survival. The BBE23 member (At5g44360/BBE23) of the Arabidopsis berberine bridge enzyme-like (BBE-l) protein family (Arabidopsis thaliana) has been characterized in this paper in parallel with the closely related and previously described CELLOX (At4g20860/BBE22). In addition to cellodextrins, both enzymes, renamed here as CELLODEXTRIN OXIDASE 2 and 1 (CELLOX2 and CELLOX1), respectively, oxidize the mixed-linked β-1→3/β-1→4-glucans (MLGs), recently described as capable of activating plant immunity, reinforcing the view that the BBE-l family includes members that are devoted to the control of the homeostasis of potential cell wall-derived damage-associated molecular patterns (DAMPs). The 2 putatively paralogous genes display different expression profiles. Unlike CELLOX1, CELLOX2 is not expressed in seedlings or adult plants and is not involved in immunity against Botrytis cinerea. Both are instead expressed in a concerted manner in the seed coat during development. Whereas CELLOX2 is expressed mainly during the heart stage, CELLOX1 is expressed at the immediately later stage, when the expression of CELLOX2 decreases. Analysis of seeds of cellox1 and cellox2 knockout mutants shows alterations in the coat structure: the columella area is smaller in cellox1, radial cell walls are thicker in both cellox1 and cellox2, and the mucilage halo is reduced in cellox2. However, the coat monosaccharide composition is not significantly altered, suggesting an alteration of the organization of the cell wall, thus reinforcing the notion that the architecture of the cell wall in specific organs is determined not only by the dynamics of the synthesis/degradation of the main polysaccharides but also by its enzymatic oxidation.

Cannabinoid Biosynthesis Using Noncanonical Cannabinoid Synthases.

We report enzymes from the berberine bridge enzyme (BBE) superfamily that catalyze the oxidative cyclization of the monoterpene moiety in cannabigerolic acid (CBGA) to form cannabielsoin (CBE). The enzymes are from a variety of organisms and are previously uncharacterized. Out of 232 homologues chosen from the enzyme superfamily, four orthologues were shown to accept CBGA as a substrate and catalyze the biosynthesis of CBE. The four enzymes discovered in this study were recombinantly expressed and purified in Pichia pastoris. These enzymes are the first report of heterologous expression of BBEs that did not originate from the Cannabis plant that catalyze the production of cannabinoids using CBGA as substrate. This study details a new avenue for discovering and producing natural and unnatural cannabinoids.

Flavofun: Exploration of fungal flavoproteomes

Fungi produce a plethora of natural products exhibiting a fascinating diversity of chemical structures with an enormous potential for medical applications. Despite the importance of understanding the scope of natural products and their biosynthetic pathways, a systematic analysis of the involved enzymes has not been undertaken. In our previous studies, we examined the flavoprotein encoding gene pool in archaea, eubacteria, the yeast Saccharomyces cerevisiae, Arabidopsis thaliana, and Homo sapiens. In the present survey, we have selected the model fungus Neurospora crassa as a starting point to investigate the flavoproteomes in the fungal kingdom. Our analysis showed that N. crassa harbors 201 flavoprotein-encoding genes amounting to 2% of the total protein-encoding genome. The majority of these flavoproteins (133) could be assigned to primary metabolism, termed the “core flavoproteome”, with the remainder of flavoproteins (68) serving in, as yet unidentified, reactions. The latter group of “accessory flavoproteins” is dominated by monooxygenases, berberine bridge enzyme-like enzymes, and glucose-methanol-choline-oxidoreductases. Although the exact biochemical role of most of these enzymes remains undetermined, we propose that they are involved in activities closely associated with fungi, such as the degradation of lignocellulose, the biosynthesis of natural products, and the detoxification of harmful compounds in the environment. Based on this assumption, we have analyzed the accessory flavoproteomes in the fungal kingdom using the MycoCosm database. This revealed large differences among fungal divisions, with Ascomycota, Basidiomycota, and Mucoromycota featuring the highest average number of genes encoding accessory flavoproteins. Moreover, a more detailed analysis showed a massive accumulation of accessory flavoproteins in Sordariomycetes, Agaricomycetes, and Glomeromycotina. In our view, this indicates that these fungal classes are proliferative producers of natural products and also interesting sources for flavoproteins with potentially useful catalytic properties in biocatalytic applications.

Molecular dynamics simulations and kinetic measurements provide insights into the structural requirements of substrate size-dependent specificity of oligogalacturonide oxidase 1 (OGOX1)

Oligogalacturonides (OGs) are pectin fragments released from the breakdown of the homogalacturonan during pathogenesis that act as Damage-Associated Molecular Patterns. OG-oxidase 1 (OGOX1) is an Arabidopsis berberine bridge enzyme-like (BBE-l) oligosaccharide oxidase that oxidizes OGs, impairing their elicitor activity and concomitantly releasing H2O2. The OG-oxidizing activity of OGOX1 is markedly pH-dependent, with optimum pH around 10, and is higher towards OGs with a degree of polymerization higher than two. Here, the molecular determinants of OGOX1 responsible for the binding of OGs with different lengths have been investigated through molecular dynamics simulations and enzyme kinetics studies. OGOX1 was simulated in complex with OGs with different degree of polymerization such as di-, tri-, tetra- and penta-galacturonide, in water solution at alkaline pH. Our simulations revealed that, among the four OGOX1/OG combinations, the penta-galacturonide (OG5) showed the best conformation in the active site to be efficiently oxidized by OGOX1. The optimal conformation can be stabilized by salt-bridges formed between the carboxyl groups of OG5 and five positively charged amino acids of OGOX1, highly conserved in all OGOX paralogs. Our results suggest that these interactions limit the mobility of OG5 as well as longer OGs, contributing to maintain the terminal monomer of OGs in the optimal orientation in order to be oxidized by the enzyme. In accordance with these results, the enzyme efficiency (Kcat/KM) of OGOX1 on OG5 (40.04) was found to be significantly higher than that on OG4 (13.05) and OG3 (0.6).

Berberine Bridge Enzyme-like Oligosaccharide Oxidases Act as Enzymatic Transducers Between Microbial Glycoside Hydrolases and Plant Peroxidases.

Oligogalacturonide (OG)-oxidase 1 (OGOX1) and cellodextrin (CD)-oxidase (CELLOX) are plant berberine bridge enzyme-like oligosaccharide oxidases that oxidize OGs and CDs, cell-wall fragments with the nature of damage-associated molecular patterns. The oxidation of OGs and CDs attenuates their elicitor activity and concomitantly releases H2O2. By using a multiple enzyme-based assay, we demonstrate that the H2O2 generated downstream of the combined action between a fungal polygalacturonase and OGOX1 or an endoglucanase and CELLOX can be directed by plant peroxidases (PODs) either towards a reaction possibly involved in plant defense, such as the oxidation of monolignol or a reaction possibly involved in a developmental event, such as the oxidation of auxin (indole-3-acetic acid), pointing to OGOX1 and CELLOX as enzymatic transducers between microbial glycoside hydrolases and plant PODs. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A Berberine Bridge Enzyme-Like Protein, GmBBE-like43, Confers Soybean's Coordinated Adaptation to Aluminum Toxicity and Phosphorus Deficiency.

Phosphorus (P) deficiency and aluminum (Al) toxicity often coexist and are two major limiting factors for crop production in acid soils. The purpose of this study was to characterize the function of GmBBE-like43, a berberine bridge enzyme-like protein-encoding gene, in soybean (Glycine max) adaptation to Al and low P stresses. Present quantitative real-time PCR (qRT-PCR) assays confirmed the phosphate (Pi)-starvation enhanced and Al-stress up-regulated expression pattern of GmBBE-like43 in soybean roots. Meanwhile, the expression of a GmBBE-like43-GFP chimera in both common bean hairy roots and tobacco leaves demonstrated its cell wall localization. Moreover, both transgenic Arabidopsis and soybean hairy roots revealed the function of GmBBE-like43 in promoting root growth under both Al and low P stresses. GmBBE-like43-overexpression also resulted in more H2O2 production on transgenic soybean hairy root surface with oligogalacturonides (OGs) application and antagonized the effects of Al on the expression of two SAUR-like genes. Taken together, our results suggest that GmBBE-like43 might be involved in the soybean's coordinated adaptation to Al toxicity and Pi starvation through modulation of OGs-oxidation in the cell wall.

Biosynthesis of Annullatin D in Penicillium roqueforti Implies Oxidative Lactonization between Two Hydroxyl Groups Catalyzed by a BBE-like Enzyme

Annullatins from Cordyceps annullata are alkylated aromatic polyketides including annullatin D with a fused dihydrobenzofuran lactone ring system. Here, we report the identification of a silent biosynthetic gene cluster for annullatins from Penicillium roqueforti by heterologous expression in Aspergillus nidulans, gene deletion, and feeding experiments as well as by biochemical characterization. The polyketide core structure is consecutively modified by hydroxylation and prenylation. A berberine bridge enzyme-like protein catalyzes the final step, an oxidative lactonization between two hydroxyl groups, to form (2S, 9S)-annullatin D.

Endophytic bacterium Bacillus aryabhattai induces novel transcriptomic changes to stimulate plant growth.

In nature, plants interact with a wide range of microorganisms, and most of these microorganisms could induce growth through the activation of important molecular pathways. The current study evaluated whether the endophytic bacterium Bacillus aryabhattai encourages plant growth and the transcriptional changes that might be implicated in this effect. The endophytic bacterium promotes the growth of Arabidopsis and tobacco plants. The transcriptional changes in Arabidopsis plants treated with the bacterium were also identified, and the results showed that various genes, such as cinnamyl alcohol dehydrogenase, apyrase, thioredoxin H8, benzaldehyde dehydrogenase, indoleacetaldoxime dehydratase, berberine bridge enzyme-like and gibberellin-regulated protein, were highly expressed. Also, endophytic bacterial genes, such as arginine decarboxylase, D-hydantoinase, ATP synthase gamma chain and 2-hydroxyhexa-2,4-dienoate hydratase, were activated during the interaction. These findings demonstrate that the expression of novel plant growth-related genes is induced by interaction with the endophytic bacterium B. aryabhattai and that these changes may promote plant growth in sustainable agriculture.

Prunus persica apoplastic proteome analysis reveals candidate proteins involved in the resistance and defense against Taphrina deformans

Taphrina deformans is the fungus responsible for the peach leaf curl disease. To gain insight into the molecular mechanisms involved in plant resistance and response to the fungus, apoplastic differentially abundant proteins (DAPs) in a resistant (DR) and/or in a susceptible genotype (FL) were identified after 12 and 96 h post inoculation (hpi) and compared to those at 0 hpi. The Prunus persica apoplastic proteome was assessed by LC-MS/MS analysis. Altogether 332 proteins were identified, and their molecular and biological functions were classified. In both genotypes, major changes occurred at 96 hpi when the fungus had achieved the filamentous form. However, at 96 hpi, DR exhibited a greater number of increased proteins than FL. DAPs were enriched in biotic stress response, with most of the proteins belonging to the pathogenesis related (PR)-type. PRs exhibited the greatest fold changes of induction in DR. While PRs acting on pathogen cell wall (PR2, PR3 and PR4) were increased in both susceptible and resistant genotypes, others were exclusively induced in DR, such as some isoforms of PR5, defensin and PR17. Proteins exclusively induced in DR upon T.deformans inoculation such as four berberine bridge enzymes, two snakins and a GDS-lipase were identified. Moreover, upon inoculation cuticle was thickened to a greater extent in DR than in FL. This work reveals the active role of the apoplast against T. deformans and not only contributes to the elucidation of responses involved in resistance to leaf curl disease but also improves the knowledge on peach defenses against pathogens.

The Agronomic Traits, Alkaloids Analysis, FT-IR and 2DCOS-IR Spectroscopy Identification of the Low-Nicotine-Content Nontransgenic Tobacco Edited by CRISPR-Cas9.

In this study, the agricultural traits, alkaloids content and Fourier transform infrared spectroscopy (FT-IR) and two-dimensional correlation infrared spectroscopy (2DCOS-IR) analysis of the tobacco after Berberine Bridge Enzyme-Like Proteins (BBLs) knockout were investigated. The knockout of BBLs has limited effect on tobacco agricultural traits. After the BBLs knockout, nicotine and most alkaloids are significantly reduced, but the content of myosmine and its derivatives increases dramatically. In order to identify the gene editing of tobacco, principal component analysis (PCA) was performed on the FT-IR and 2DCOS-IR spectroscopy data. The results showed that FT-IR can distinguish between tobacco roots and leaves but cannot classify the gene mutation tobacco from the wild one. 2DCOS-IR can enhance the characteristics of the samples due to the increased apparent resolution of the spectra. Using the autopeaks in the synchronous map for PCA analysis, we successfully identified the mutants with an accuracy of over 90%.