Pathogen species specific sequence screening method, device and equipment and storage medium

Through multiple rounds of alignment and pan-genome analysis, the megablast and blastn programs were used to screen microbial species-specific sequences, solving the problems of database access restrictions and alignment interruptions in the prior art, and achieving rapid and accurate specific sequence screening.

CN119170103APending Publication Date: 2024-12-20THE FIRST AFFILIATED HOSPITAL OF XIAMEN UNIV

Patent Information

Application Number
CN202411195964.6
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2024-08-29
Publication Date
2024-12-20

AI Technical Summary

Technical Problem

The existing automated methods are limited by online database access restrictions in microbial species-specific sequence screening, and cannot flexibly respond to the needs of different species, and the comparison process is frequently and unexpectedly interrupted.

Method used

By reading the input file of the target species name, obtaining and preprocessing the genomic sequence, performing pan-genomic analysis, calling the megablast and blastn programs for multiple rounds of screening, excluding highly similar sequences and ignored core genomic sequences of non-target species, and obtaining specific sequence of sexual differences.

Benefits of technology

The problem of query sequence size limitation and interruption of alignment process when screening microbial species-specific sequences is solved, and fast and accurate specific sequence screening is achieved.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119170103A_ABST
    Figure CN119170103A_ABST
Patent Text Reader

Abstract

The invention provides a pathogen species specific sequence screening method, device and equipment and a storage medium, and the method comprises the following steps: reading an input file containing a target species name to obtain a genome sequence containing a target species, and preprocessing the genome sequence; then, extracting the number of the target species, and performing generic genome analysis on the genome sequence of the target species to obtain a core genome sequence; and finally, performing specific sequence screening on the target species according to the number of the target species and the core genome sequence, specifically, calling a BLAST + software megablast program to perform a first round of screening so as to exclude sequence fragments with high similarity with non-target species genomes, and performing a second round of screening so as to exclude sequence fragments with high similarity with non-target species genomes. And calling a BLAST + software blastn program to carry out a second round of screening so as to eliminate the core genome sequence ignored by the megablast program, and calling the blastn program to carry out a third round of screening so as to obtain a sequence having a specific difference with other non-target species in the genus. A target microbial species specific sequence screening process is constructed through local version BLAST +, and the problems that when the microbial species specific sequence is screened, the size of a query sequence is submitted online through NCBI, and the comparison process is frequently and accidentally interrupted are solved.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to the field of specific sequence screening, and in particular to a method, device, equipment and storage medium for screening pathogen species specific sequences. Background Art

[0002] Bacteria and other microorganisms are the main pathogens that cause many zoonotic infectious diseases. Rapid and accurate identification and detection of pathogens is crucial. Traditional microbial isolation, culture and immunological detection methods are time-consuming and complex. In recent years, molecular detection technologies based on specific genome sequences have received widespread attention.

[0003] However, manual screening of species-specific sequences is a time-consuming and laborious task. Existing automated methods are often limited by online database access restrictions or lack the flexibility to meet the needs of different species.

[0004] In view of this, this application is filed. Summary of the invention

[0005] The present invention discloses a pathogen species-specific sequence screening method, device, equipment and storage medium, aiming to solve the problems of the current lack of species-specific sequence screening methods and the limitation of the size of query sequences submitted online through NCBI when screening microbial species-specific sequences, as well as the frequent and unexpected interruptions of the comparison process.

[0006] The first embodiment of the present invention provides a method for screening pathogen species-specific sequences, comprising:

[0007] Reading an input file containing the name of a target species, obtaining a genome sequence containing the target species, and preprocessing the genome sequence;

[0008] Extracting the target species number, and performing pan-genome analysis on the genome sequence of the target species to obtain a core genome sequence;

[0009] The target species is subjected to specific sequence screening according to the target species number and the core genome sequence, wherein the specific sequence screening includes calling the megablast program for a first round of screening to exclude sequence fragments with high similarity to the genomes of non-target species, calling the blastn program for a second round of screening to exclude core genome sequences ignored by the megablast program, and calling the blastn program for a third round of screening to obtain sequences with specific differences from other non-target species in the genus.

[0010] Preferably, the step of obtaining the genome sequence of the target species is as follows:

[0011] The genome sequences of all strains of the target species were extracted from the local genome reference database. If the corresponding genome file did not exist locally, the corresponding genome sequence was downloaded from the Refseq database.

[0012] Preferably, preprocessing the genome sequence comprises:

[0013] The data with complete genome type was extracted, the compressed genome files were batch-decompressed, the sequence numbers of the RefSeq database were used as strain names, and the plasmid sequences contained in each genome were removed.

[0014] Preferably, the extraction of the target species number includes: classification numbers at the bacterial kingdom level and below, classification numbers at the genus level and below, classification numbers at the species level that are not classified below the genus level, and classification numbers at the species level and below.

[0015] Preferably, the pan-genome analysis of the genome sequence of the target species to obtain the core genome sequence is specifically:

[0016] Obtain analysis parameters and call the Panseq analysis tool to perform pan-genome analysis on the target species to generate a core genome sequence, wherein the analysis parameters include the length of the species genome sequence after cutting, the consistency range of the core genome sequence, and the proportion of the core genome sequence in all genomes of the target species.

[0017] Preferably, the target database for the first round of screening and comparison is a bacterial genome database, excluding the target species level;

[0018] The target database for the second round of screening was the bacterial genome database, excluding the target species level;

[0019] The target database for the third round of screening is all genomes at the genus level and below of the target species, but does not include genomes at the target species level and unclassified target species.

[0020] The second embodiment of the present invention provides a pathogen species-specific sequence screening device, comprising:

[0021] A preprocessing unit, used for reading an input file containing the name of a target species, obtaining a genome sequence containing the target species, and preprocessing the genome sequence;

[0022] A pan-genome analysis unit, used to extract the target species number and perform pan-genome analysis on the genome sequence of the target species to obtain a core genome sequence;

[0023] The specific sequence screening unit is used to perform specific sequence screening on the target species according to the target species number and the core genome sequence, wherein the specific sequence screening includes calling the megablast program to perform a first round of screening to exclude sequence fragments with high similarity to the genomes of non-target species, calling the blastn program to perform a second round of screening to exclude core genome sequences ignored by the megablast program, and calling the blastn program to perform a third round of screening to obtain sequences with specific differences from other non-target species in the genus.

[0024] The third embodiment of the present invention provides a pathogen species-specific sequence screening device, including a memory and a processor, wherein the memory stores a computer program, and the computer program can be executed by the processor to implement a pathogen species-specific sequence screening method as described in any one of the above.

[0025] The fourth embodiment of the present invention provides a computer-readable storage medium storing a computer program, wherein the computer program can be executed by a processor of a device where the computer-readable storage medium is located to implement a pathogen species-specific sequence screening method as described in any one of the above items.

[0026] Based on a pathogen species-specific sequence screening method, device, equipment and storage medium provided by the present invention, the genome sequence containing the target species is obtained by reading an input file containing the name of the target species, and the genome sequence is preprocessed; then, the target species number is extracted, and the genome sequence of the target species is subjected to pan-genome analysis to obtain the core genome sequence; finally, the target species is subjected to specific sequence screening according to the target species number and the core genome sequence, wherein the specific sequence screening includes calling the megablast program for the first round of screening to exclude highly similar sequence fragments with the genome of non-target species, calling the blastn program for the second round of screening to exclude the core genome sequence ignored by the megablast program, and calling the blastn program for the third round of screening to obtain sequences with specific differences from other non-target species in the genus. The target microbial species-specific sequence screening process is constructed by the local version of BLAST+, which solves the problem of the limitation of the query sequence size submitted online through NCBI when screening microbial species-specific sequences, and the problem of frequent unexpected interruptions in the alignment process. BRIEF DESCRIPTION OF THE DRAWINGS

[0027] Figure 1 It is a schematic flow chart of a pathogen species-specific sequence screening method provided in the first embodiment of the present invention;

[0028] Figure 2It is a schematic diagram of the specific process of pathogen species-specific sequence screening provided by the present invention;

[0029] Figure 3 It is a module schematic diagram of a pathogen species-specific sequence screening method and device provided in the second embodiment of the present invention. DETAILED DESCRIPTION

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

[0031] In order to better understand the technical solution of the present invention, the embodiments of the present invention are described in detail below with reference to the accompanying drawings.

[0032] It should be clear that the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0033] The terms used in the embodiments of the present invention are only for the purpose of describing specific embodiments, and are not intended to limit the present invention. The singular forms "a", "said" and "the" used in the embodiments of the present invention and the appended claims are also intended to include plural forms, unless the context clearly indicates other meanings.

[0034] It should be understood that the term "and / or" used in this article is only a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in this article generally indicates that the associated objects before and after are in an "or" relationship.

[0035] The word "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determining" or "in response to detecting", depending on the context. Similarly, the phrases "if it is determined" or "if (stated condition or event) is detected" may be interpreted as "when it is determined" or "in response to determining" or "when detecting (stated condition or event)" or "in response to detecting (stated condition or event)", depending on the context.

[0036] The "first\second" mentioned in the embodiments is only to distinguish similar objects, and does not represent a specific order for the objects. It is understandable that the "first\second" can be interchanged with the specific order or sequence where permitted. It should be understood that the objects distinguished by "first\second" can be interchanged where appropriate, so that the embodiments described herein can be implemented in an order other than those illustrated or described herein.

[0037] The specific embodiments of the present invention are described in detail below with reference to the accompanying drawings.

[0038] The present invention discloses a pathogen species-specific sequence screening method, device, equipment and storage medium, aiming to solve the problem of limitation of query sequence size submitted online through NCB I when screening microbial species-specific sequences, and frequent unexpected interruption of the comparison process.

[0039] See also Figure 1 and Figure 2 The first embodiment of the present invention provides a pathogen species-specific sequence screening method, which can be performed by a pathogen species-specific sequence screening device (hereinafter referred to as the screening device), and in particular, is performed by one or more processors in the screening device to at least implement the following steps:

[0040] S101, reading an input file containing the name of a target species, obtaining a genome sequence containing the target species, and preprocessing the genome sequence;

[0041] In this embodiment, the screening device can be a terminal with data processing capabilities such as a laptop computer, a desktop computer, a server, etc., which can communicate with the database in real time. The screening device can be installed with a corresponding operating system and application software, and the functions required by this embodiment can be realized through the combination of the operating system and application software.

[0042] It should be noted that, in this embodiment, the screening process of pathogen species-specific sequences can be started by reading an input file containing the name of the target species. The input file can be a text file manually created by the user, which contains the specific name of the pathogen species to be screened, or relevant information can be entered through the user interface. The system obtains the genome sequences of all related strains of the species from the local genome reference database based on the input target species name. If the corresponding genome sequence does not exist in the local database, the system will automatically connect to the NCBI RefSeq database and download all relevant genome sequences of the target species from it. During the download process, the system can automatically select data with a complete genome type to ensure the accuracy and completeness of subsequent analysis.

[0043] After obtaining the genome sequence of the target species, these genome sequences are preprocessed. First, the genome sequence downloaded from the NCBI RefSeq database is decompressed. This can be done in an automated batch processing manner to efficiently process large-scale genome data. For the preprocessed genome sequence, the system will further remove the plasmid sequence. Plasmids are usually DNA molecules independent of bacterial chromosomes, which may carry drug resistance genes or virulence factors, but may interfere with the analysis of the core genome during the screening of species-specific sequences, so they are removed during the preprocessing process.

[0044] S102, extracting the target species number, and performing pan-genome analysis on the genome sequence of the target species to obtain a core genome sequence;

[0045] It should be noted that, in this embodiment, first, after reading the input file containing the name of the target species, the screening device automatically extracts the classification number of the species by querying the local or online database. They are the classification numbers at the bacterial kingdom level and below, the classification numbers at the genus level and below, the classification numbers at the species level that are not classified below the genus level, and the classification numbers at the species level and below. The classification numbers on the four evolutionary branches are used to set the final BLAST comparison reference database.

[0046] After extracting the target species number, the screening device performs pan-genome analysis on all genome sequences of the target species to obtain the core genome sequence. Pan-genome analysis is completed by calling the Panseq tool. It can automatically process and analyze the genome sequences of multiple strains and identify the core genome sequences shared by all strains. During the analysis process, the screening device will optimize the analysis process according to pre-set analysis parameters (such as the length of the genome sequence fragment, the consistency range of the core genome, and the proportion of the core genome in all genomes of the target species). These parameters can be adjusted according to the characteristics of the target species and specific research needs;

[0047] The core genome sequences obtained through pan-genome analysis represent highly conserved gene fragments in all strains of the target species. These fragments are present in almost all strains within the species and are therefore highly species-specific. This high conservation makes the core genome sequence an important basis for distinguishing the target species from other species. At the same time, pan-genome analysis can also identify non-core genome sequences, that is, gene fragments that only exist in some strains.

[0048] S103, performing specific sequence screening on the target species according to the target species number and the core genome sequence, wherein the specific sequence screening includes calling the megablast program to perform a first round of screening to exclude sequence fragments with high similarity to the genomes of non-target species, calling the blastn program to perform a second round of screening to exclude core genome sequences ignored by the megablast program, and calling the blastn program to perform a third round of screening to obtain sequences with specific differences from other non-target species in the genus.

[0049] It should be noted that in this embodiment, in each round of screening results, the core genome fragments whose consistency with the genus-level alignment sequences of non-target species is higher than the set threshold and whose aligned genome length is greater than the set threshold are filtered out, and the core sequences are retained.

[0050] Specifically, in the present embodiment, the first round of screening: the core genome sequence of the target species is compared with the genome of the non-target species by the efficient comparison algorithm of the megablast program, and sequence fragments with high similarity in the non-target species are quickly identified and excluded. The main goal of this round of screening is to remove sequences that may cause false positive results, i.e., sequences that cannot be effectively distinguished between the target species and the non-target species. Through this round of screening, the sequence fragments retained are already target species specific to a great extent.

[0051] Further, because megablast may have omissions when processing short sequences or sequences with high variation, after the first round of screening, the system further calls the blastn program to carry out the second round of screening. The blastn program has higher sensitivity and can capture sequence fragments that megablast may ignore, particularly those short sequences with biological significance or sequences that are conservative in target species but have larger variations in non-target species. The second round of screening further optimizes the recognition process of specific sequences by in-depth comparison analysis. It should be noted that the screening conditions in the second round also include: expected value (for evaluating the statistical significance of the comparison result), word length (word length, the minimum length for initial exact matching), preferably, the parameters can be: "expected value" is set to 0.001, "word length" is set to 11, and "masking" and "filtering" options are turned off.

[0052] Next, the screening device conducts the third round of screening, still calling the blastn program, but this time the comparison database is limited to all genome sequences at the genus level of the target species. In this round, the screening device focuses on screening sequence fragments that are specific within the genus of the target species and have significant differences between species. The purpose of this process is to further ensure that the screened sequences are not only highly conserved within the target species, but also show sufficient differences among other species in the genus. The specific sequences finally obtained can effectively distinguish the target species from its closely related species in practical applications.

[0053] In this embodiment, the target database for the first round of screening and comparison is a bacterial genome database, excluding the target species level;

[0054] In the first round of screening, bacterial genome databases that do not include the target species at the genus level are selected for comparison, with the goal of excluding as many genera as possible that are widely present in higher taxonomic levels (such as families, orders, etc.). The goal of this round is to remove sequence fragments that are commonly present in all genera and higher levels other than the target species. These fragments may be shared among multiple species and cannot distinguish the target species from other species. Therefore, after excluding these sequences, the sequence fragments that remain have a relatively strong ability to distinguish within the genus of the target species.

[0055] The target database for the second round of screening was the bacterial genome database, excluding the target species level;

[0056] The second round of screening continued to use the bacterial genome database that did not include the target species at the genus level for comparison, with the goal of further filtering and optimizing the results after the first round of screening. In this round, the system used the blastn program, which has a higher comparison sensitivity and can identify sequences missed by the megablast program, especially those sequence fragments that are highly conserved within the species but still have a certain degree of specificity at a higher taxonomic level. Continuing to exclude sequences at the genus level of the target species is to further eliminate sequence fragments that may be shared with other species under the genus of the target species, ensuring that the sequences that are ultimately retained are more specific.

[0057] The target database for the third round of screening is all genomes at the genus level of the target species, and does not include genomes at the species level and unclassified species.

[0058] In the third round of screening, the target database for screening is limited to all genomes at the genus level of the target species, but does not include genomes at the target species level and unclassified species. The main purpose of this round of screening is to further refine the screening within the genus of the target species to identify sequences that do not exist or have significant differences in other species within the genus. This means that the third round of screening focuses on distinguishing the differences in gene sequences between the target species and its closely related species within the genus. By excluding genomic data at the species level and unclassified species, the specificity of the screening sequence can be further improved, ensuring that the final sequence obtained is highly conserved within the target species and has sufficient differences in its closely related species to avoid cross-reactions.

[0059] See also Figure 3 The second embodiment of the present invention provides a pathogen species-specific sequence screening device, comprising:

[0060] A preprocessing unit 201 is used to read an input file containing the name of a target species, obtain a genome sequence containing the target species, and preprocess the genome sequence;

[0061] A pan-genome analysis unit 202 is used to extract the target species number and perform a pan-genome analysis on the genome sequence of the target species to obtain a core genome sequence;

[0062] The specific sequence screening unit 203 is used to perform specific sequence screening on the target species according to the target species number and the core genome sequence, wherein the specific sequence screening includes calling the megablast program to perform a first round of screening to exclude sequence fragments with high similarity to the genomes of non-target species, calling the blastn program to perform a second round of screening to exclude core genome sequences ignored by the megablast program, and calling the blastn program to perform a third round of screening to obtain sequences with specific differences from other non-target species in the genus.

[0063] The third embodiment of the present invention provides a pathogen species-specific sequence screening device, including a memory and a processor, wherein the memory stores a computer program, and the computer program can be executed by the processor to implement a pathogen species-specific sequence screening method as described in any one of the above.

[0064] The fourth embodiment of the present invention provides a computer-readable storage medium storing a computer program, wherein the computer program can be executed by a processor of a device where the computer-readable storage medium is located to implement a pathogen species-specific sequence screening method as described in any one of the above items.

[0065] Based on a pathogen species-specific sequence screening method, device, equipment and storage medium provided by the present invention, the genome sequence containing the target species is obtained by reading an input file containing the name of the target species, and the genome sequence is preprocessed; then, the target species number is extracted, and the genome sequence of the target species is subjected to pan-genome analysis to obtain a core genome sequence; finally, the target species is subjected to specific sequence screening according to the target species number and the core genome sequence, wherein the specific sequence screening includes calling the megablast program for a first round of screening to exclude highly similar sequence fragments with the genome of non-target species, calling the blastn program for a second round of screening to exclude core genome sequences ignored by the megablast program, and calling the blastn program for a third round of screening to obtain sequences with specific differences from other non-target species in the genus. The problem of the limitation of the size of the query sequence submitted online through NCBI when screening microbial species-specific sequences and the frequent and unexpected interruption of the alignment process is solved.

[0066] Exemplarily, the computer program described in the third and fourth embodiments of the present invention may be divided into one or more modules, and the one or more modules are stored in the memory and executed by the processor to complete the present invention. The one or more modules may be a series of computer program instruction segments capable of completing specific functions, and the instruction segments are used to describe the execution process of the computer program in the device for implementing a battery self-modeling device based on cyclic charging and discharging. For example, the device described in the second embodiment of the present invention.

[0067] The processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor, etc. The processor is the control center of the method for screening a pathogen species-specific sequence, and uses various interfaces and lines to connect the various parts of the method for screening a pathogen species-specific sequence.

[0068] The memory can be used to store the computer program and / or module, and the processor realizes various functions of a pathogen species-specific sequence screening method by running or executing the computer program and / or module stored in the memory, and calling the data stored in the memory. The memory can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, an application required for at least one function (such as a sound playback function, a text conversion function, etc.), etc.; the data storage area can store data created according to the use of the mobile phone (such as audio data, text message data, etc.), etc. In addition, the memory can include a high-speed random access memory, and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device, or other volatile solid-state storage devices.

[0069] Wherein, if the implemented module is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, the present invention implements all or part of the processes in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Wherein, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium. It should be noted that the content contained in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable media do not include electric carrier signals and telecommunication signals.

[0070] It should be noted that the device embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. In addition, in the accompanying drawings of the device embodiments provided by the present invention, the connection relationship between the modules indicates that there is a communication connection between them, which may be specifically implemented as one or more communication buses or signal lines. A person of ordinary skill in the art may understand and implement it without paying any creative effort.

[0071] The above is only a preferred specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by a person skilled in the art within the technical scope disclosed by the present invention should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims

1. A method for screening pathogen species-specific sequences, characterized in that: include: Reading an input file containing the name of a target species, obtaining a genome sequence containing the target species, and preprocessing the genome sequence; Extracting the target species number, and performing pan-genome analysis on the genome sequence of the target species to obtain a core genome sequence; The target species is subjected to specific sequence screening according to the target species number and the core genome sequence, wherein the specific sequence screening includes calling the megablast program for a first round of screening to exclude sequence fragments with high similarity to the genomes of non-target species, calling the blastn program for a second round of screening to exclude core genome sequences ignored by the megablast program, and calling the blastn program for a third round of screening to obtain sequences with specific differences from other non-target species in the genus.

2. A pathogen species-specific sequence screening method according to claim 1, characterized in that: The step of obtaining the genome sequence of the target species is as follows: The genome sequences of all strains of the target species were extracted from the local genome reference database. If the corresponding genome file did not exist locally, the corresponding genome sequence was downloaded from the Refseq database.

3. A pathogen species-specific sequence screening method according to claim 1, characterized in that: Preprocessing the genome sequence includes: The data with complete genome type were extracted, the compressed genome files were batch-decompressed, the sequence numbers of the RefSeq database were used as strain names, and the plasmid sequences contained in each genome were removed.

4. The method for screening pathogen species-specific sequences according to claim 1, characterized in that: The extracted target species numbers include: classification numbers at the bacterial kingdom level and below, classification numbers at the genus level and below, classification numbers at the species level that are not classified below the genus level, and classification numbers at the species level and below.

5. The method for screening pathogen species-specific sequences according to claim 1, characterized in that: The pan-genome analysis of the genome sequence of the target species is performed to obtain the core genome sequence, specifically: Obtain analysis parameters and call the Panseq analysis tool to perform pan-genome analysis on the target species to generate a core genome sequence, wherein the analysis parameters include the length of the species genome sequence after cutting, the consistency range of the core genome sequence, and the proportion of the core genome sequence in all genomes of the target species.

6. The method for screening pathogen species-specific sequences according to claim 1, characterized in that: The target database for the first round of screening and comparison is a bacterial genome database, excluding the target species level; The target database for the second round of screening was the bacterial genome database, excluding the target species level; The target database for the third round of screening is all genomes at the genus level and below of the target species, but does not include genomes at the target species level and unclassified target species.

7. A pathogen species-specific sequence screening device, characterized in that: include: A preprocessing unit, used for reading an input file containing the name of a target species, obtaining a genome sequence containing the target species, and preprocessing the genome sequence; A pan-genome analysis unit, used to extract the target species number and perform pan-genome analysis on the genome sequence of the target species to obtain a core genome sequence; A specific sequence screening unit is used to perform specific sequence screening on the target species according to the target species number and the core genome sequence, wherein the specific sequence screening includes calling the megablast program for a first round of screening to exclude sequence fragments with high similarity to the genomes of non-target species, calling the blastn program for a second round of screening to exclude core genome sequences ignored by the megablast program, and calling the blastn program for a third round of screening to obtain sequences with specific differences from other non-target species in the genus.

8. A pathogen species-specific sequence screening device, characterized in that: It comprises a memory and a processor, wherein the memory stores a computer program, and the computer program can be executed by the processor to implement a pathogen species-specific sequence screening method as described in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that: A computer program is stored, and the computer program can be executed by a processor of the device where the computer-readable storage medium is located to implement a pathogen species-specific sequence screening method as described in any one of claims 1 to 6.

Citation Information

Patent Citations

  • Method and system for screening specific sequences of pathogenic species

    CN115719616A

  • Method for qualitative and quantitative detection of microorganism in human body

    US20190048393A1

Cited By

  • Low-abundance pathogen detection method and system based on metagenome sequencing data

    CN120032716A

  • Specific detection primer of bifidobacterium longum subsp. Longum BBMN68 and application thereof

    CN122060887A

  • Specific detection primer of bifidobacterium longum subsp. longum bbmn68 and application thereof

    CN122060887B