Exploration of Candidate Genes Involved in the Biosynthesis, Regulation and Recognition of the Male-Produced Aggregation Pheromone of Halyomorpha halys

Simple Summary The brown marmorated stink bug, Halyomorpha halys, is an invasive insect pest native to Asia that was accidently introduced into North America and Europe. The male-produced aggregation pheromone of H. halys has great application potential in population monitoring and control of this pest. In this study, we carried out a comprehensive bioinformatics analysis to identify the candidate genes involved in the biosynthesis of the male-specific sesquiterpene aggregation pheromone of H. halys. We identified a novel gene, designated as HhTPS1, which is related to the synthesis of the precursor molecule of the aggregation pheromone. Moreover, potential P450 reductase and transcription factors related to the aggregation pheromone biosynthesis pathway were detected. In addition, two olfactory-related genes, chemosensory protein 5 (HhCSP5) and olfactory receptor 85b (HhOr85b), were also identified, which are assumed to be involved in the olfactory recognition of this pheromone. Our results lay a solid foundation for further elucidating the biosynthesis of this aggregation pheromone and its behavioral regulation of H. halys. Abstract The aggregation pheromone of the brown marmorated stink bug, Halyomorpha halys (Stål), is produced by adult males, and plays an important role in the behavioral regulation of H. halys. However, information on the molecular mechanisms underlying this pheromone’s biosynthesis is limited. In this study, HhTPS1, a key candidate synthase gene in the aggregation pheromone biosynthesis pathway of H. halys, was identified. Then, through weighted gene co-expression network analysis, the candidate P450 enzyme genes in the biosynthetic downstream of this pheromone and the related candidate transcription factor in this pathway were also identified. In addition, two olfactory-related genes, HhCSP5 and HhOr85b, involved in the recognition of the aggregation pheromone of H. halys, were detected. We further identified the key amino acid sites of HhTPS1 and HhCSP5 that interact with substrates by using molecular docking analysis. This study provides basic information for further investigations into the biosynthesis pathways and recognition mechanisms of aggregation pheromones in H. halys. It also provides key candidate genes for bioengineering bioactive aggregation pheromones necessary for the development of technologies for the monitoring and control of H. halys.


Introduction
Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive, polyphagous insect pest. Adults and nymphs of H. halys damage many important crops, genes involved in the biosynthesis of the aggregation pheromone of H. halys have not yet been fully explored, and transcription factors regulating the expression of these biosynthesis genes have not been reported. The process of olfactory recognition in insects begins with the adsorption of mostly hydrophobic odorant molecules on the surface of the antenna. From here, the molecules have to pass through the pores and pore tubules of the olfactory hair cuticle to enter the sensillum. The odorant molecule has to then transverse the hydrophilic fluid filling the cuticular hair (the sensillum lymph), before reaching the dendritic membrane of the olfactory receptor neurons [24]. Detailed analyses of sensillum lymph from different insect species have led to the discovery of high concentrations of soluble proteins that bind hydrophobic compounds, including pheromones, as well as general odorant molecules. These small globular proteins have been already classified as odorant/pheromone-binding proteins (O/PBPs), chemosensory proteins (CSPs), and Niemann-Pick-type C2 proteins (NPC2), among others. When the solubilized odorant binds to odorant receptors in the dendritic membrane of the olfactory receptor cell, the cell is depolarized and generates action potentials that transmit the olfactory signal to the antennal lobe, eventually stimulating the corresponding physiological or behavioral responses [24][25][26]. Although several olfactory genes have previously been described from H. halys, such as CSP, OBP and OR [27,28], other olfactory genes may regulate the recognition of aggregation pheromones in H. halys that have not yet been reported.
In this study, diverse bioinformatics approaches were used to identify key enzymecoding genes involved in the biosynthesis pathway of aggregation pheromone in H. halys, and the transcription factors that regulate the expression of these enzymes were analyzed. In addition, the olfactory genes that may be responsible for the recognition of aggregation pheromones in H. halys were also identified. This study identifies candidate genes for the production of aggregation pheromones of H. halys, aiding the development of new management and control technologies for H. halys.

Identification of Farnesyl Diphosphate (FDP) and Terpene Synthase (TPS) Genes in H. halys
The genomic coding sequences and protein sequences files of H. halys (assembly Hhal_2.0) were downloaded from NCBI. The FDP and trans-IDS-type TPS genes of H. halys were identified using the blastp method and the FDP and TPS sequences from seven other insect species, including M. histrionica, Phyllotreta striolata and Ips pini, with E values < 10 −10 . Names of IDS followed the reporting study [29]. The orthologous gene clusters between H. halys and harlequin bug (M. histrionica) were identified using the OrthoVenn2 (https://orthovenn2.bioinfotoolkits.net/) web tool [30]. Parameter settings were: E-value = 1× 10 -10 and Inflation value = 1.5. For M. histrionica, the open reading frames were obtained The process of olfactory recognition in insects begins with the adsorption of mostly hydrophobic odorant molecules on the surface of the antenna. From here, the molecules have to pass through the pores and pore tubules of the olfactory hair cuticle to enter the sensillum. The odorant molecule has to then transverse the hydrophilic fluid filling the cuticular hair (the sensillum lymph), before reaching the dendritic membrane of the olfactory receptor neurons [24]. Detailed analyses of sensillum lymph from different insect species have led to the discovery of high concentrations of soluble proteins that bind hydrophobic compounds, including pheromones, as well as general odorant molecules. These small globular proteins have been already classified as odorant/pheromone-binding proteins (O/PBPs), chemosensory proteins (CSPs), and Niemann-Pick-type C2 proteins (NPC2), among others. When the solubilized odorant binds to odorant receptors in the dendritic membrane of the olfactory receptor cell, the cell is depolarized and generates action potentials that transmit the olfactory signal to the antennal lobe, eventually stimulating the corresponding physiological or behavioral responses [24][25][26]. Although several olfactory genes have previously been described from H. halys, such as CSP, OBP and OR [27,28], other olfactory genes may regulate the recognition of aggregation pheromones in H. halys that have not yet been reported.
In this study, diverse bioinformatics approaches were used to identify key enzymecoding genes involved in the biosynthesis pathway of aggregation pheromone in H. halys, and the transcription factors that regulate the expression of these enzymes were analyzed. In addition, the olfactory genes that may be responsible for the recognition of aggregation pheromones in H. halys were also identified. This study identifies candidate genes for the production of aggregation pheromones of H. halys, aiding the development of new management and control technologies for H. halys.

Identification of Farnesyl Diphosphate (FDP) and Terpene Synthase (TPS) Genes in H. halys
The genomic coding sequences and protein sequences files of H. halys (assembly Hhal_2.0) were downloaded from NCBI. The FDP and trans-IDS-type TPS genes of H. halys were identified using the blastp method and the FDP and TPS sequences from seven other insect species, including M. histrionica, Phyllotreta striolata and Ips pini, with E values < 10 −10 . Names of IDS followed the reporting study [29]. The orthologous gene clusters between H. halys and harlequin bug (M. histrionica) were identified using the OrthoVenn2 (https://orthovenn2.bioinfotoolkits.net/, accessed on 14 February 2022) web tool [30]. Parameter settings were: E-value = 1× 10 -10 and Inflation value = 1.5. For M. histrionica, the open reading frames were obtained by TransDecoder v5.5.0 [31] from the unigenes identified in transcriptome assembly (Genbank ID:GECQ00000000.1) [21]. MEGA6.0 [32] was used to compare the IDS protein sequences of H. halys with the IDS protein sequences of seven other insects, and to construct a maximum likelihood phylogenetic tree, with the bootstrap value set at 1000 and the cut-off value set at 50%. Transcriptomes of 15 samples from different developmental life stages of H. halys (2nd instar, 4th instar), sexes (male and female), and tissues (antenna, salivary gland, principal salivary gland, and accessory salivary gland) were downloaded from the NCBI SRA database (Table S1). The raw data was processed with fastp [33], including adapter trimming and quality filtering, to obtain high-quality clear data, which were then aligned to the assembly Hhal_2.0 with HISAT2 [34]. We then calculated the count value of each gene expressed in the transcriptome with featureCount [35]. Then, we converted the count values for each gene into TPM expression values using Tbtools [36].

Weighted Gene Co-Expression Network Analysis
To identify co-expressed gene modules of the HhTPS gene, we used the same 15 samples mentioned above. TPM values were used as standardized gene expression levels. The genes with zero variance, and genes with more than 10% missing sample material were filtered using the R (V.3.6.1) software (https://cran.r-project.org/bin/windows/base/old/, accessed on 8 February 2022) WGCNA (V.1.70-3) package [37], which builds a weighted gene co-expression network. To do this, first, the weight value was calculated using the function PicksoftThreshold, selecting the appropriate weight value and allowing the network to match the unmetered network distribution, and then selecting the optimal soft valve according to the analysis results. An automatic block-wise multiple-set network construction module [37] was used to build a co-expression network. We set the parameters at minmodulesize = 30, maxBlockSize = 10,000, and other parameters according to the default settings. For the co-expressed genes identified in this way, gene function annotation and enrichment analyses were performed using the KOBAS online tool [38], and the corrected p-value of < 0.01 was set as the threshold.

Molecular Docking of HhTPS1 and HhCSP 5 in H. halys
The amino acid sequences of HhTPS1 and HhCSP5 were used for blastp searches in the NCBI PDB and the RCSB PDB databases. The PDB ID: 2rah [39] and PDB ID: 2gvs [40] were selected as templates for HhTPS1 and HhCSP5, respectively. MODELLER 9v15 was used for protein homologous modeling [41]. The modeling results were evaluated with the Structural Analysis and Verification Server (https://saves.mbi.ucla.edu/, accessed on 14 February 2022). The PROCHECK server (http://servicesn.mbi.ucla.edu/PROCHECK/, accessed on 14 February 2022) was used for evaluating the quality of the modelled HhTPS1 and HhCSP5 proteins' structures concerning energy and stereo-chemical geometry [42]. Using the ERRAT program from SAVES (Structure Analysis and Verification Server), nonbonded interactions were analyzed, and the quality factors for the modelled HhTPS1 and HhCSP5 proteins were recorded in percentages [43]. Compatibility of the atomic model (3D) with its amino acid sequence (1D) was determined with the help of the VERIFY-3D program from SAVES and the results were obtained (3D-1D Profile) in percentages [44].

Identification of Key FDP and TPS Genes in the Aggregation Pheromone Biosynthesis Pathway
Through homology comparison, we identified seven IDS genes from the genome of H. halys, and a phylogenetic analysis showed that the HhIDS2 and FDP genes of the harlequin bug (Genbank ID:AVZ23978.1) were in the same clan ( Figure 2). The protein sequence homology of these two genes was 90%. The HhIDS1 and HhIDS3-7 genes were in the same clan as the TPS gene of the harlequin bug (Genbank ID:AVZ23977.1), and the HhIDS7 gene was in the same subclan as the TPS gene of the harlequin bug ( Figure 2). The protein sequence homology of these two genes was 80%.

Identification of Key FDP and TPS Genes in the Aggregation Pheromone Biosynthesis Pathway
Through homology comparison, we identified seven IDS genes from the genome of H. halys, and a phylogenetic analysis showed that the HhIDS2 and FDP genes of the harlequin bug (Genbank ID:AVZ23978.1) were in the same clan ( Figure 2). The protein sequence homology of these two genes was 90%. The HhIDS1 and HhIDS3-7 genes were in the same clan as the TPS gene of the harlequin bug (Genbank ID:AVZ23977.1), and the HhIDS7 gene was in the same subclan as the TPS gene of the harlequin bug ( Figure 2). The protein sequence homology of these two genes was 80%. The OrthoVenn2 was used to identify orthologous genes among H. halys and M. histrionica, and a total of 19,589 orthologous clusters were identified (Table S2). Among these clusters, only two IDS genes were detected that were shared between H. halys and M. histrionica. One orthologous cluster was XP_014276183.1 (HhIDS7) in H. halys and GECQ01420512.1 (TPS, Genbank ID:AVZ23977.1) in M. histrionica. Another orthologous cluster was XP_024219390.1(HhIDS2) in H. halys and GECQ01414919.1 (FDP, Genbank ID:AVZ23978.1) M. histrionica (Table S2). According to the sequence homology (>99%) of these IDS genes with that of recently reported genes (HhTPS1 Genbank ID: MG870387; HhTPS2 Genbank ID: MG917093 and HhFPPS Genbank ID: MG870389) [29], the HhIDS7, HhIDS1 and HhIDS2 genes were named as HhTPS1, HhTPS2 and HhFPPS, respectively. The results of expression profile analysis showed that HhTPS2 and HhIDS3-6 were almost not expressed in the second instar, fourth instar, adult male, and adult female of H. halys, while HhFPPS was expressed in the second instar, fourth instar, adult male, and adult female with a TPS value between 13 and 18, with the highest expression level in the fourth instar. HhTPS1 was also expressed in the second instar, fourth instar, adult male, and adult female with a TPS value between 6 and 22, and with the highest level of expression in the adult male and adult female ( Figure 2). The OrthoVenn2 was used to identify orthologous genes among H. halys and M. histrionica, and a total of 19,589 orthologous clusters were identified (Table S2). Among these clusters, only two IDS genes were detected that were shared between H. halys and M. histrionica. One orthologous cluster was XP_014276183.1 (HhIDS7) in H. halys and GECQ01420512.1 (TPS, Genbank ID:AVZ23977.1) in M. histrionica. Another orthologous cluster was XP_024219390.1(HhIDS2) in H. halys and GECQ01414919.1 (FDP, Genbank ID:AVZ23978.1) M. histrionica (Table S2). According to the sequence homology (>99%) of these IDS genes with that of recently reported genes (HhTPS1 Genbank ID: MG870387; HhTPS2 Genbank ID: MG917093 and HhFPPS Genbank ID: MG870389) [29], the HhIDS7, HhIDS1 and HhIDS2 genes were named as HhTPS1, HhTPS2 and HhFPPS, respectively. The results of expression profile analysis showed that HhTPS2 and HhIDS3-6 were almost not expressed in the second instar, fourth instar, adult male, and adult female of H. halys, while HhFPPS was expressed in the second instar, fourth instar, adult male, and adult female with a TPS value between 13 and 18, with the highest expression level in the fourth instar. HhTPS1 was also expressed in the second instar, fourth instar, adult male, and adult female with a TPS value between 6 and 22, and with the highest level of expression in the adult male and adult female (Figure 2). Based on the above evolutionary relationship analysis and expression analysis, we concluded that HhTPS1 is a key TPS gene for the synthesis of sesquipiperitol as a precursor for the aggregation pheromone in H. halys. HhTPS1 may be responsible for catalyzing (E,E)-FDP generation (1S,6S,7R) sesquipiperitol. HhFPPS is the FDP gene responsible

Co-Expression Network Analysis to Identify the HhTPS1 Co-Expressed Genes
RNA sequencing was used to analyze the gene expression partners of all H. halys genes across the 15 samples types studied. Samples were clustered based on their gene expression level. As shown in Figure S1, the clustering results based on average distance and a hierarchical clustering algorithm showed that the materials used for sequencing had high repeatability in each process and the gene expression patterns in the same process were similar and clustered together. After screening the soft valve values (β), we selected the β value of nine. A weighted gene co-expression network was successfully constructed by the WGCNA method, and the network was divided into 19 modules. Among these modules, the TPS gene was detected in the red module, which suggests that the module may participate in the aggregation pheromone biosynthesis and action process of this pest. The number of genes in the red module was 1185 (Table S3). KEGG functional enrichment analysis of the genes of this module found that these genes were enriched in 25 pathways, including metabolic pathways, oxidative phosphorylation, carbon metabolism, peroxisome, and metabolism of xenobiotics by cytochrome P450 pathways, among others (Table S4). GO functional enrichment analysis showed that these genes were enriched in 76 GO terms, including the mitochondrion, mitochondrial respiratory chain complex I, mitochondrial electron transport, NADH to ubiquinone, plasma membrane, NADH dehydrogenase activity, and transferase activity, and that for transferring hexosyl groups, among others. (Table S5).
Further analysis of the red module (Table S3) identified candidate genes that may be involved in the regulation and recognition of aggregation pheromones in H. halys (Table 1). We identified seven P450 genes, including three CYP6a family members, three CYP6j family members, and one CYP4C family member. It is worth noting that these seven genes may be responsible for aggregation pheromone downstream biosynthesis. Seventeen transcription factors were also identified, including ten zinc finger proteins, two coiled-coil-helix-coiledcoil-helix domain-containing proteins, one putative homeodomain transcription factor, one RING finger protein, and three other types of transcription factors. These transcription factors may be responsible for the regulation of the expression of FDP, TPS, and P450 genes in the aggregation pheromone biosynthesis pathway. In addition, two olfactory-related genes were identified: one chemosensory protein 5 (HhCSP5) and one odorant receptor 85b (HhOr85b) ( Table 1).

Discussion
The aggregation pheromone of H. halys consists of sesquiterpenes, which play an important role in the ecology of H. halys [5,6]. However, these compounds have a complex structure and are an enantiomer mixture of two sesquiterpene alcohols, with a low chemical synthesis yield and potential for serious environmental pollution. Therefore, although pheromone-baited traps are effective in attracting H. halys, the cost of their production is too high for practical use. The average annual cost of using H. halys aggregation pheromone traps in the United States is about $900 US dollars per acre over the season [5]. The cost of chemical synthesis of aggregation pheromones in H. halys limits the large-scale application of this pheromone trapping technology.
Compared with these costs of chemical synthesis, biosynthesis technology can potential to greatly reduce the cost of this aggregation pheromone, making the large-scale application of this pheromone feasible. Currently, the biosynthesis of a variety of sesquiterpenes through yeast cells, transgenic plants, or other technologies has been successfully implemented [46,47].
Biosynthesis technology is an effective tool to solve the problem of the high cost of chemical synthesis of the aggregation pheromone of H. halys. In this study, we identified the genes for the key enzymes (HhFPPS and HhTPS1) in the aggregation pheromone biosynthesis pathway, which have an orthologous relationship with the FDP and trans-IDStype TPS genes of the aggregation pheromone biosynthesis pathway in the harlequin bug, with more than 80% protein sequence homology. Combined with expression profiling and molecular docking results, HhTPS1 was further confirmed as the key enzyme gene for the biosynthesis of aggregation pheromones in H. halys (Figures 2 and S2). The in vitro synthase activity assay of the HhFPPS and HhTPS1 genes has recently been confirmed [29]. The identification of these genes is useful for future potential large-scale biosynthesis of the aggregation pheromone of H. halys using metabolic engineering techniques, which would The highest binding energy and highest interaction energy between HhCSP5 and (3R,6S,7R,10S)-10,11-epoxy-1-bisabolen-3-ol was 5.2411 and −22.0535 kcal/mol, respectively. There was Alkyl force between LYS320 of HhCSP5 and (3R,6S,7R,10S)-10,11-epoxy-1-bisabolen-3-ol, and there was a conventional hydrogen bond between GLU317 and the (3R,6S,7R,10S)-10,11-epoxy-1-bisabolen-3-ol ( Figure 3B). These amino acids were the key amino acid residues of the interaction between HhCSP5 and aggregation pheromone.

Discussion
The aggregation pheromone of H. halys consists of sesquiterpenes, which play an important role in the ecology of H. halys [5,6]. However, these compounds have a complex structure and are an enantiomer mixture of two sesquiterpene alcohols, with a low chemical synthesis yield and potential for serious environmental pollution. Therefore, although pheromone-baited traps are effective in attracting H. halys, the cost of their production is too high for practical use. The average annual cost of using H. halys aggregation pheromone traps in the United States is about $900 US dollars per acre over the season [5]. The cost of chemical synthesis of aggregation pheromones in H. halys limits the large-scale application of this pheromone trapping technology.
Compared with these costs of chemical synthesis, biosynthesis technology can potential to greatly reduce the cost of this aggregation pheromone, making the large-scale application of this pheromone feasible. Currently, the biosynthesis of a variety of sesquiterpenes through yeast cells, transgenic plants, or other technologies has been successfully implemented [46,47].
Biosynthesis technology is an effective tool to solve the problem of the high cost of chemical synthesis of the aggregation pheromone of H. halys. In this study, we identified the genes for the key enzymes (HhFPPS and HhTPS1) in the aggregation pheromone biosynthesis pathway, which have an orthologous relationship with the FDP and trans-IDS-type TPS genes of the aggregation pheromone biosynthesis pathway in the harlequin bug, with more than 80% protein sequence homology. Combined with expression profiling and molecular docking results, HhTPS1 was further confirmed as the key enzyme gene for the biosynthesis of aggregation pheromones in H. halys (Figures 2 and S2). The in vitro synthase activity assay of the HhFPPS and HhTPS1 genes has recently been confirmed [29]. The identification of these genes is useful for future potential large-scale biosynthesis of the aggregation pheromone of H. halys using metabolic engineering techniques, which would make widespread application of this pheromone inexpensive. In the future, the in vivo function of the HhTPS1 gene should be further confirmed using RNAi, and CRISPR/Cas9.
Through molecular docking methods, the key amino acid sites of HhTPS1 that catalyze (E,E)-FDP were identified ( Figure S2). The potentially important roles of ARG132, ASP258, and LYS272 in the interaction between HhTPS1 and (E,E)-FDP have been supported by a recent study [29]. These amino acid sites are conducive to our understanding of the molecular interaction between HhTPS1 and (E,E)-FDP. In the future, the directional evolution of HhTPS1 enzymes could be realized through mutation of these amino acid sites to obtain the synthases with more catalytic activity and efficiency.
Through co-expression analysis, seven relevant P450 genes were identified (Table 1). These P450 genes are candidate genes responsible for the production of aggregation pheromones by P450 epoxidase catalyzation of sesquipiperitol or zingiberenol in the aggregation pheromone biosynthesis pathway [12,20]. Whether these reactions are catalyzed by a single or two separate P450 epoxidases remains to be determined [12]. So far, the P450 genes in the aggregation pheromone biosynthesis pathway of H. halys are still unknown and have not been reported in other stink bug species, such as harlequin bug M. histrionica, although some candidate P450 genes have been predicted in M. histrionica [21]. The P450 genes identified in this study are therefore candidate genes for understanding the downstream biosynthesis of the aggregation pheromone of H. halys and its relatives. These candidate genes can be validated by strategies such as RNAi, CRISPR/Cas9, or yeast or insect-cell expression. In the future, HhTPS1 and P450 co-expression by Saccharomyces cerevisiae might allow useful strains to be engineered for large scale production. Such an engineered fungal strain could be used in the industrial production of the aggregation pheromone.
The transcription factors identified in this study may be involved in the regulation of the HhFPPS, HhTPS1 or HhP450 genes in the aggregation pheromone biosynthesis pathway of H. halys. The function of these transcription factors needs to be verified by RNAi, gene editing, yeast-one-hybrid assays, or electrophoretic mobility shift assays. Based on the importance of aggregation pheromones in the behavior regulation of H. halys, these transcription factors are important candidate target genes for the future management and control of H. halys by gene editing or RNAi technology.
In the current study, we also identified a chemosensory protein and an olfactory receptor in H. halys that may recognize the aggregation pheromones. For these two candidate genes, in vivo functional studies can be carried out through RNAi or gene editing, followed by assays on the physiology or behavior responses of H. halys using the electroantennography technique and Y-tube olfactometers. For in vitro functional studies, the fluorescence-based binding assay could be conducted for HhCSP5, and for HhOr85b, we recommend use of the two-electrode voltage-clamp recording (in Xenopus oocyte) expression system.
In this study, we also detected the key amino acids of HhCSP5 that bind to the two main components of the aggregation pheromones. The interaction between this chemosensory protein (HhCSP5) and the aggregation pheromone is useful to design and develop precise environmentally friendly behavior regulation substances through a reverse chemical ecology approach to control H. halys populations. This method has been successfully applied in the screening of a variety of insect behavior-modifying substances [48][49][50]. Using this method, we anticipate being able to obtain H. halys behavior-regulating substances that might have similar functions to the natural aggregation pheromone while being cheaper to synthesize and more stable during use. The docking analysis of the interaction between HhOr85b and the aggregation pheromone was not carried out, because we could not detect an appropriate template with reasonable homology. In insects, the recognition of pheromones often involves multiple odorant-binding proteins and chemoreceptors [49,[51][52][53]. There-fore, in addition to the HhCSP5 identified in this study, there may be other OBPs or CSPs responsible for the recognition and transport of aggregation pheromones in H. halys.
In conclusion, we identified HhTPS1, a candidate enzyme for aggregation pheromone biosynthesis in H. halys, and we also found some P450s, transcriptional factors, and olfactoryrelated genes that are co-expressed with HhTPS1. Altogether, these genes provide a solid foundation for understanding the biosynthesis pathways and recognition mechanisms of aggregation pheromones in H. halys. These genes also provide key candidate targets for the management and control of H. halys through CRISPR/Cas9 or RNAi technology.