Methods, equipment, media and products for risk assessment of avian influenza virus spillover

Abstract

The present application discloses a method, device, medium and product for assessing the risk of avian influenza virus spillover, and relates to the field of virus spillover risk assessment. The method for assessing the risk of avian influenza virus spillover includes obtaining the gene sequence of avian influenza virus; extracting the gene index information of multiple strains of typical viruses from the gene sequence; and based on the gene index information, assessing the risk of avian influenza virus spillover according to the region or the subtype / lineage to be assessed, and obtaining an assessment result. While improving the efficiency of spillover risk assessment, the present application can achieve early warning of spillover risk and high-throughput assessment.

Description

Technical Field

[0001] The present application relates to the field of virus spillover risk assessment, and in particular to an avian influenza virus spillover risk assessment method, equipment, medium and product. Background Art

[0002] Spillover refers to the transmission of pathogens from animal hosts to humans, and spillover risk assessment has important public health significance. There are many risk assessment tools in the existing technology, one of which is the second edition of the Influenza Pandemic Risk Assessment Tool (TIPRA) proposed in 2020, which expands the scope of assessment from viruses with human cases to animal influenza viruses with important public health significance. However, it does not clearly define the criteria for determining important public health significance, and most assessment reports are still for viruses with human cases. In addition, TIPRA still adopts the form of expert consultation and is carried out in the form of monthly risk assessment. Before the meeting, relevant information covering public health, animal health, virology and other aspects needs to be prepared. Although the reliability of the results is guaranteed, it is difficult to meet the timeliness requirements. Another is an influenza risk assessment tool based on expert consultation, among which the Influenza Risk Assessment Tool (IRAT) only evaluated 24 strains of viruses from 2011 to 2023. It can be seen that the current pandemic assessment tools are inefficient and difficult to achieve early warning and high-throughput assessment. Summary of the invention

[0003] The purpose of this application is to provide an avian influenza virus spillover risk assessment method, equipment, medium and product, which can improve the efficiency of spillover risk assessment and achieve early warning and high-throughput assessment of spillover risk.

[0004] To achieve the above objectives, this application provides the following solutions:

[0005] In a first aspect, the present application provides a method for assessing the risk of avian influenza virus spillover, comprising:

[0006] Obtain the genetic sequence of avian influenza virus;

[0007] Extracting gene index information of multiple typical virus strains from the gene sequence; the gene index information includes: nucleotide sequence, sampling year, fragment name or number, subtype or lineage name; the typical strains are selected based on set conditions;

[0008] Based on the gene indicator information, an avian influenza virus spillover risk assessment is performed according to the region or subtype / lineage to be assessed to obtain an assessment result.

[0009] Optionally, based on the gene indicator information, the risk assessment of avian influenza virus spillover is performed according to the region, and the process of obtaining the assessment result includes:

[0010] Based on the gene indicator information, homology screening is performed on avian influenza viruses in a set area to obtain a risk combination;

[0011] Count the strains to be calculated in the risk combination, and determine the annual number of risk strains found in the region and their host distribution;

[0012] The assessment results were obtained based on the annual number of findings and host distribution.

[0013] Optionally, based on the gene indicator information, homology screening is performed on avian influenza viruses in a set area to obtain a risk combination, specifically including:

[0014] Based on the gene indicator information, multiple sequence analysis is performed on the avian influenza virus to be evaluated in the set area according to MAFFT segmentation to obtain an aligned sequence;

[0015] Determine the genetic distance between each strain to be calculated and the typical strain in the aligned sequence by segmentation;

[0016] Screen the combinations whose genetic distance is less than the set value by segment, and obtain the first segment combination;

[0017] Determine the similarity of the first fragment combination, and screen the combinations whose similarity is greater than a set percentage by fragment to obtain a second fragment combination;

[0018] A combination of at least n fragments is screened from the second fragment combinations, and this combination is used as the risk combination.

[0019] Optionally, use the pairwiseAlignment() function of the Biostrings package in R language for multiple sequence analysis;

[0020] The pid() function is used to calculate the similarity of the first fragment combination.

[0021] Optionally, the set value is 0.1 and the set percentage is 90%.

[0022] Optionally, based on the gene indicator information, the avian influenza virus spillover risk assessment is performed according to the subtype / lineage to be assessed, and the process of obtaining the assessment result includes:

[0023] Based on the gene indicator information, homology screening is performed on the avian influenza virus of the subtype / lineage to be evaluated and a risk combination is obtained;

[0024] Based on the relationship between the indicators of the risk combination and the set threshold, an avian influenza virus spillover risk assessment is performed to obtain an assessment result; the indicators of the risk combination include the number of risk combinations, the number of risk strains and the number of typical strains.

[0025] In a second aspect, the present application provides a computer device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of any one of the above-mentioned methods for assessing the risk of avian influenza virus spillover.

[0026] In a third aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of any of the above-mentioned methods for assessing the risk of avian influenza virus spillover.

[0027] In a fourth aspect, the present application provides a computer program product, comprising a computer program, which, when executed by a processor, implements the steps of any one of the above-mentioned methods for assessing the risk of avian influenza virus spillover.

[0028] According to the specific embodiments provided in this application, this application discloses the following technical effects:

[0029] The present application provides a method, device, medium and product for assessing the risk of avian influenza virus spillover. By extracting the gene indicator information of multiple strains of typical viruses from the gene sequence, the indicators of avian influenza spillover risk assessment can be increased, and the width and breadth of the assessment system can be expanded. Based on the gene indicator information, the avian influenza virus spillover risk assessment is performed according to the region or the subtype / lineage to be assessed, and the assessment results are obtained, which can improve the efficiency of spillover risk assessment while achieving early warning of spillover risk and high-throughput assessment. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0031] Figure 1 This is an application environment diagram of an avian influenza virus spillover risk assessment method in one embodiment of the present application;

[0032] Figure 2 A schematic diagram of a flow chart of a method for assessing the risk of avian influenza virus spillover provided in one embodiment of the present application;

[0033] Figure 3A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0034] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0035] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0036] The avian influenza virus spillover risk assessment method provided in the embodiments of the present application can be applied to Figure 1 In the application environment shown. Among them, the terminal 102 communicates with the server 104 through the network. The data storage system can store the data that the server 104 needs to process. The data storage system can be set up separately, integrated on the server 104, or placed on the cloud or other servers. The terminal 102 can send the gene sequence of the avian influenza virus to the server 104. After the server 104 receives the gene sequence of the avian influenza virus, for the gene sequence of the avian influenza virus, the server 104 extracts the gene indicator information of multiple strains of typical viruses from the gene sequence. Based on the gene indicator information, the avian influenza virus spillover risk assessment is performed according to the region or subtype / lineage to be assessed to obtain an assessment result. The server 104 can feedback the obtained assessment results to the terminal 102. In addition, in some embodiments, the avian influenza virus spillover risk assessment method can also be implemented independently by the server 104 or the terminal 102. For example, the terminal 102 can directly perform spillover risk assessment on the gene sequence of the avian influenza virus, or the server 104 can obtain the gene sequence of the avian influenza virus from the data storage system and perform spillover risk assessment on the gene sequence of the avian influenza virus.

[0037] The terminal 102 may be, but is not limited to, various desktop computers, laptop computers, smart phones, tablet computers, IoT devices, and portable wearable devices. The IoT devices may be smart speakers, smart TVs, smart air conditioners, smart vehicle-mounted devices, etc. The portable wearable devices may be smart watches, smart bracelets, head-mounted devices, etc. The server 104 may be implemented as an independent server or a server cluster consisting of multiple servers, or may be a cloud server.

[0038] The gene sequences shared by GenBank and GISAID contain information such as subtype, time, space, and host, which provide data support for describing the historical epidemic situation of each subtype in detail and analyzing gene homology. H9N2 has been confirmed many times to participate in the recombination process of multiple human-infected avian influenza viruses such as H3N8, H7N9, and H10N8 as an internal gene donor, and can be used as a potential indicator for avian influenza spillover risk assessment. Based on this, in the description process of this application, H9N2 is used for explanation.

[0039] In an exemplary embodiment, Figure 2 As shown, a method for assessing the risk of avian influenza virus spillover is provided. The method is executed by a computer device, and can be executed by a computer device such as a terminal or a server alone, or by a terminal and a server together. In the embodiment of the present application, the method is applied to Figure 1 The server 104 in the example is used as an example to illustrate, including the following steps 200 to 202. Among them:

[0040] Step 200, obtain the gene sequence of avian influenza virus. For example, obtain the H9N2 gene sequence through public databases such as GenBank and GISAID, or obtain the following information through independent sampling and sequencing (conventional sampling, sequencing and statistical analysis of effective information):

[0041] (1) Internal gene sequences of each strain. Internal genes refer to the fragments other than HA gene and NA gene, including alkaline polymerase 1 (PB1), alkaline polymerase 2 (PB2), nucleoprotein (NP), NS gene, acidic polymerase (PA), and M protein.

[0042] (2) Sample collection time, specific to the year.

[0043] (3) Sample collection location, specific to the province.

[0044] (4) Sample collection host.

[0045] Step 201: Gene index information of multiple typical virus strains is extracted from the gene sequence. Gene index information includes: nucleotide sequence, sampling year, fragment name or number, subtype or lineage name. Typical strains are selected based on set conditions. For example, taking important H9N2 strains as typical strains, 12 important H9N2 strains are extracted as shown in Table 1.

[0046] Table 112 Important H9N2 strains

[0047]

[0048]

[0049] Step 202: Based on the gene indicator information, the avian influenza virus spillover risk assessment is performed according to the region or subtype / lineage to be assessed to obtain an assessment result. For example, in actual application, the avian influenza virus spillover risk assessment based on the internal gene homology of H9N2 can be divided into two forms: (1) specific region.

[0050] (2) Specific subtype or lineage.

[0051] Implementing the above steps 201 to 203 can improve the efficiency of spillover risk assessment while achieving early warning and high-throughput assessment of spillover risk.

[0052] In another exemplary embodiment of the present application, in order to provide a simple and quick way to assess the risk of avian influenza virus spillover in a specific area, step 202 provided above in the present application is replaced by the following steps (1) to (3):

[0053] (1) Based on the genetic indicator information, the avian influenza virus in the set area is homologously screened and the risk combination is obtained. Specifically:

[0054] (1.1) Based on the gene index information, the avian influenza virus to be evaluated in the set area is subjected to multiple sequence analysis according to MAFFT segmentation to obtain an aligned sequence. For example, sequences with a length greater than 90% are selected and the default parameters are used.

[0055] (1.2) Determine the genetic distance between each strain to be calculated and the typical strain in the aligned sequence by segment. For example, based on the aligned sequence in (1.1), the computer calculates the genetic distance between each strain to be calculated and the specific important strain of H9N2 by segment, and limits the sampling year of the strain to be calculated to be later than the specific important strain of H9N2. In actual application, the aligned sequence information can be read in through the read.dna() and phyDat() functions of the ape package of the R language, and the dist.ml() function is used to calculate the genetic distance between the two sequences.

[0056] (1.3) Segment by segment screening combinations with genetic distances less than a set value to obtain a first segment combination. For example, the computer segments screen combinations with genetic distances between the strain to be calculated and a specific important H9N2 strain less than 0.1.

[0057] (1.4) Determine the similarity of the first fragment combination, and screen the combinations with similarity greater than a set percentage by fragment, to obtain the second fragment combination. For example, the computer calculates the similarity of the combination after screening by (1.3) (i.e., the first fragment combination) by fragment. In actual application, the pairwiseAlignment() function of the Biostrings package of the R language can be used to perform global alignment analysis, and the pid() function can be used to calculate its similarity.

[0058] (1.5) A combination containing at least n (for example, 6) fragments is screened from the second fragment combination, and this combination is used as a risk combination. The risk combination obtained at this time has internal gene homology.

[0059] (2) Count the strains to be calculated in the risk combination and determine the annual number of risk strains found in the region and their host distribution. The purpose of counting the strains to be calculated in the risk combination is to achieve the deduplication of the names of the strains to be calculated involved in the risk combination.

[0060] For example, taking H9N2 as an example, through steps (1) and (2), H9N2 risk strains can be obtained, that is, all H9N2 strains (non-important H9N2 strains) are taken as strains to be calculated, and their internal gene homology with important H9N2 strains is calculated. H9N2 risk strains in specific areas are screened according to the sampling locations of H9N2 risk strains, and the annual number of findings and host distribution are counted according to the sampling time and sampling host.

[0061] (3) The assessment results are obtained based on the annual number of discoveries and the distribution of hosts. For example, the computer assesses the spillover risk by comparing multiple locations or setting thresholds. Specifically:

[0062] a) Comparison of multiple regions: For example, a total of 15 H9N2 risk strains were found in region A from 2001 to 2021, mainly distributed from 2001 to 2005. In terms of host distribution, they were mainly found in poultry, and also in pigs and wild birds. A total of 2 H9N2 risk strains were found in region B from 2001 to 2021, which were found in poultry in 2002 and 2009 respectively. In comparison, the H9N2 epidemic intensity in region A is higher, and there is a higher risk of avian influenza virus spillover.

[0063] b) Setting thresholds: You can also set corresponding thresholds based on the above indicators to directly conduct spillover risk assessments. For example, you can set ① more than 5 H9N2 risk strains were found between 2001 and 2021. ② at least 2 of the poultry, wild birds, pigs, and humans have records of H9N2 risk strains. ① and ② must be met at the same time. According to the above standards, area A can be judged as high risk and area B as low risk.

[0064] The above text uses H5N1 (specific subtype or lineage), specific region (spatial range, spatial resolution), time range (2001-2021), and time resolution (years) as examples for detailed description.

[0065] In another exemplary embodiment of the present application, in order to provide a simple and rapid way to assess the risk of spillover of a specific avian influenza subtype or lineage, the above step 202 performs an assessment of the risk of spillover of avian influenza virus according to the subtype / lineage to be assessed based on the gene indicator information, and the process of obtaining the assessment result may include the following steps (1) to (2). Among them:

[0066] (1) Based on the gene index information, the avian influenza virus of the subtype / lineage to be evaluated is homologously screened and a risk combination is obtained. The method for obtaining the risk combination is as described in steps (1.1) to (1.5) in the above embodiment.

[0067] (2) Based on the relationship between the risk combination indicators and the set thresholds, an avian influenza virus spillover risk assessment is performed to obtain an assessment result. The risk combination indicators include the number of risk combinations, the number of risk strains, and the number of typical strains. For example, a computer is used to compare the number of risk combinations, the number of risk strains, the number of H9N2 strains, and other indicators of different subtypes or lineages with H9N2, or to set thresholds, to determine the spillover risk of a specific subtype or lineage. Among them:

[0068] (2.1) The number of risk combinations of the risk combinations in step (1) is obtained by counting.

[0069] (2.2) Extract the names of the strains to be calculated involved in the risk combination, remove duplicates and count them to obtain the number of risk strains.

[0070] (2.3) Extract the names of important H9N2 strains involved in the risk combination, remove duplicates and count them to obtain the number of H9N2 strains.

[0071] (2.4) Comparison of different subtypes or lineages: For example, the number of risk combinations for H5N1 is 4069, the number of risk strains is 1282, and the number of H9N2 strains is 12; the number of risk combinations for H5N2 is 928, the number of risk strains is 164, and the number of H9N2 strains is 12. The number of risk combinations and the number of risk strains for H5N1 are both greater than those for H5N2, so the spillover risk of H5N1 is greater than that of H5N2.

[0072] (2.5) Setting thresholds. For example, if the threshold for the number of risk combinations is set to 50, the threshold for the number of risk strains is set to 10, and the threshold for the number of H9N2 strains is set to 12, then all three indicators of H5N1 will reach this threshold, and the spillover risk is high.

[0073] Based on the above description, the present application can quickly determine the internal gene homology between two strains through the similarity and genetic distance of 6 internal genes, and determine the spillover risk of a specific subtype or lineage by taking the number of strains with high homology to H9N2 important strains in a specific subtype or lineage, the number of H9N2 important strains involved and the number of combinations formed as examples. The risk of avian influenza virus spillover in a certain area can be assessed by the prevalence of H9N2 strains with high homology to H9N2 important strains in that area.

[0074] In summary, the two avian influenza virus spillover risk assessment methods provided in this application are in parallel relationship, and both can independently implement spillover risk assessment. In the implementation process, different implementation forms of avian influenza spillover risk can be evaluated based on H9N2 internal gene homology.

[0075] The present application also provides an application scenario, which applies the above-mentioned avian influenza virus overflow risk assessment method. Specifically: the avian influenza virus overflow risk assessment method provided in this embodiment can be applied in the content distribution scenario. The internal distribution scenario includes the content production link, the content processing link and the content distribution link; the gene sequence of the avian influenza virus enters the content processing link from the content production link, obtains the corresponding content features through human-computer collaboration, and enters the downstream content distribution link. The avian influenza virus overflow risk assessment method provided in this embodiment belongs to the overflow risk assessment link in the content processing link. Specifically, in the content processing link process for the gene sequence of the avian influenza virus, the overflow risk assessment can be performed on the gene sequence of the avian influenza virus based on the collaborative method of machine processing and manual processing.

[0076] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Figure 3 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. Among them, the processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store overflow risk assessment data. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a method for assessing the overflow risk of avian influenza virus is implemented.

[0077] Those skilled in the art will understand that Figure 3 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components. In an exemplary embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above-mentioned method embodiments when executing the computer program.

[0078] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0079] In an exemplary embodiment, a computer program product is provided, including a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0080] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.

[0081] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM).

[0082] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but is not limited thereto.

[0083] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0084] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A method for assessing the risk of avian influenza virus spillover, characterized in that: The avian influenza virus spillover risk assessment method comprises: Obtain the genetic sequence of avian influenza virus; Extracting gene index information of multiple typical virus strains from the gene sequence; the gene index information includes: nucleotide sequence, sampling year, fragment name or number, subtype or lineage name; typical strains are selected based on set conditions; Based on the gene indicator information, an avian influenza virus spillover risk assessment is performed according to the region or subtype / lineage to be assessed to obtain an assessment result; Based on the genetic indicator information, the risk assessment of avian influenza virus spillover is conducted according to the region, and the process of obtaining the assessment results includes: (1) Based on the gene indicator information, homology screening is performed on avian influenza viruses in a set area and risk combinations are obtained, specifically: (1.1) Based on the gene indicator information, multiple sequence analysis is performed on the avian influenza virus to be evaluated in the set area according to MAFFT segmentation to obtain an aligned sequence; (1.2) determining the genetic distance between each strain to be calculated in the aligned sequence and the typical strain by segmentation; (1.3) Screening the combinations whose genetic distance is less than the set value by segment to obtain the first segment combination; (1.4) determining the similarity of the first fragment combination, and screening the combinations whose similarity is greater than a set percentage by fragment, to obtain a second fragment combination; (1.5) screening the second fragment combinations to obtain a combination having at least n fragments, and using the combination as the risk combination; (2) Count the strains to be calculated in the risk combination and determine the annual number of risk strains found in the region and their host distribution; (3) The assessment results are based on the annual number of discoveries and host distribution.

2. The method for assessing the risk of avian influenza virus spillover according to claim 1, characterized in that: Multiple sequence analysis was performed using the pairwiseAlignment() function of the Biostrings package in R; The pid() function is used to calculate the similarity of the first fragment combination.

3. The method for assessing the risk of avian influenza virus spillover according to claim 1, characterized in that: The set value is 0.1 and the set percentage is 90%.

4. The method for assessing the risk of avian influenza virus spillover according to claim 1, characterized in that: Based on the gene indicator information, the risk assessment of avian influenza virus spillover is performed according to the subtype / lineage to be assessed, and the process of obtaining the assessment result includes: Based on the gene indicator information, homology screening is performed on the avian influenza virus of the subtype / lineage to be evaluated and a risk combination is obtained; Based on the relationship between the indicators of the risk combination and the set threshold, an avian influenza virus spillover risk assessment is performed to obtain an assessment result; the indicators of the risk combination include the number of risk combinations, the number of risk strains and the number of typical strains.

5. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the avian influenza virus spillover risk assessment method according to any one of claims 1 to 4.

6. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for assessing the risk of avian influenza virus spillover according to any one of claims 1 to 4 is implemented.

7. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the method for assessing the risk of avian influenza virus spillover according to any one of claims 1 to 4 is implemented.

Images