This application claims the benefit of and priority to Chinese patent application 202310065167.5 filed Feb. 6, 2023, which is incorporated herein by reference in its entirety.
Antibiotics were once a “secret weapon” for humans to fight against many diseases. In the late 19th and early 20th centuries, the discovery of a series of antibiotics greatly increased human lifespan. In recent years, with the continuous application of antibiotics, drug abuse has gradually emerged, which leads to an increase in clinical antibiotic resistance and adverse reactions, and brings a heavy burden to the global economy. Effectively controlling the abuse of antibiotics in healthcare is an important link in addressing the global issue of antibiotic resistance.
Pathogenic microorganisms refer to microorganisms that can invade the human body, and cause infections or even infectious diseases, which are also known as pathogens. Pathogenic microorganisms mainly include bacteria, viruses, fungi, parasites, mycoplasma, chlamydia, rickettsia, spirochetes, etc. There are various types of microbial samples. Intestinal samples include feces, mucous membranes, etc. Liquid samples include urine, blood, cerebrospinal fluid, saliva, sputum, alveolar lavage fluid, amniotic fluid, etc. Swab samples include samples from oral cavity, reproductive tract, skin, etc. Other samples include tissues, liver, eyes, placenta, etc.
The Latin scientific name of Genus Klebsiella is Klebsiella Trevisan, and its systematic classification level is Genus. It is a straight bacterium with a diameter of 0.3-1.0 μm. 0.6-6.0 in length μM, and it's single, paired, or short chain arrangement. Currently reported strains of this genus include Klebsiella pneumoniae, Klebsiella aerogenes, Klebsiella oxytoca, Klebsiella quasipneumoniae, Klebsiella variicola, Klebsiella michiganensis, etc.
Among them, Klebsiella pneumoniae, as a model strain of the genus Genus Klebsiella, is widely present in the environment, and is a common opportunistic pathogen that easily colonizes the respiratory and intestinal tracts of patients, and causes infections in multiple parts of the digestive tract, respiratory tract, blood, etc. It is one of the pathogenic bacteria that cause human pneumonia and is also one of the common drug-resistant bacteria in hospitals. According to a study by the Second Military Medical University, the resistance rate against meropenem of carbapenem resistant Klebsiella pneumoniae isolated from 2014 to 2017 to was 62.5% (252/403).
Meropenem is a widely antibacterial and injectable antibiotic, which is used to treat various infections including meningitis and pneumonia. It is a type of β Lactam antibiotics and belongs to the category of carbapenems. Meropenem has a wide antibacterial spectrum, strong antibacterial activity, and is stable to βlactamase.
Bacterial drug sensitivity testing is currently the most commonly used method for detecting bacterial resistance in clinic and laboratory both domestically and internationally. There are methods such as paper disc method, agar dilution method, broth dilution method, and concentration gradient method. Except for the paper disc method, all other methods can obtain relatively accurate minimum inhibitory concentration (MIC) of drugs. It's required firstly to obtain pure culture in the bacterial drug sensitivity test, which is not suitable for difficult to cultivate and non cultured bacteria, and takes a long time. Sometimes, it is difficult to meet the current clinical needs for rapid diagnosis and targeted treatment of severe and emergency infections. The traditional detection and identification methods of pathogenic microorganisms fail to meet the comprehensive needs of wide coverage, speed, and accuracy. The diagnosis and treatment of infectious diseases are mainly based on empirical and directional methods. Clinical doctors and patients urgently need innovative detection methods to identify infectious pathogens more comprehensively, accurately, and quickly, which can assist in diagnosis and reasonable standardized medication treatment, shorten treatment courses, reduce mortality rates, and reduce medical costs.
With the promotion of emerging technologies such as PCR technology, whole genome sequencing technology, microfluidic technology, VITEK-2 compact fully automated bacterial identification/drug sensitivity system, the exploration of new technologies for bacterial resistance detection is gradually deepening, and various new technologies and methods for bacterial resistance detection are becoming increasingly mature. Although the VITEK-2 compact fully automated bacterial identification/drug sensitivity system is simple and fast, its accuracy in bacterial identification/drug sensitivity evaluation is influenced by the sample status and bacterial culture, and its usage cost is relatively high.
Therefore, there is an urgent need to develop a system and method in this field that can quickly, accurately, and cost-effectively predict the sensitivity of Klebsiella strains against meropenem.
In response to the aforementioned shortcomings and requirements of prior art in this field, the present invention aims to provide a system and method of predicting sensitivity of Klebsiella against MeropeneM.
Technical solution of the present invention is as follows.
A system of predicting sensitivity of Klebsiella against MeropeneM, comprises a computing unit. The computing unit comprises a computer readable storage medium on which is stored a computer program. An Exp (−k) power value calculation method is implemented when the computer program is executed by a processor. The Exp(−k) power value calculation method comprises following computing steps:
In formula I,
In above formula I, the first pair of parentheses includes C1 minus 0.223 divided by 0.757; the second pair of parentheses includes C2 minus 0.951 divided by 0.090; the third pair of parentheses includes C3 minus 0.952 divided by 0.103; the fourth pair of parentheses includes C4 minus 0.469 divided by 1.090; and the fifth pair of parentheses includes C5 minus 0.196 divided by 0.670.
In above formula I, k=0−1.44+0.708×(C1 minus 0.223 divided by 0.757)−0.660×(C2 minus 0.951 divided by 0.090)+0.088×(C3 minus 0.952 divided by 0.103)+3.048×(C4 minus 0.469 divided by 1.090)−0.460×(C5 minus 0.196 divided by 0.670).
The system of predicting sensitivity of Klebsiella against MeropeneM, also comprises a result output unit. The computing unit transmits the calculated Exp(−k) power value to the result output unit. The result output unit recognizes Exp(−k) power value and outputs result;
The system of predicting sensitivity of Klebsiella against MeropeneM, also comprises an experiment unit and a data input unit. The experiment unit and data input unit are communicated through a data-path. The experiment unit outputs experiment results which are transmitted to the data input unit and transformed to independent variables, The data input unit and computing unit are communicated through the data-path, Independent variables are transmitted to the computing unit.
The independent variables include: values of C1, C2, C3, C4, C5;
A method of predicting sensitivity of Klebsiella against MeropeneM, comprises:
In formula I,
In above formula I, the first pair of parentheses includes C1 minus 0.223 divided by 0.757; the second pair of parentheses includes C2 minus 0.951 divided by 0.090; the third pair of parentheses includes C3 minus 0.952 divided by 0.103; the fourth pair of parentheses includes C4 minus 0.469 divided by 1.090; and the fifth pair of parentheses includes C5 minus 0.196 divided by 0.670.
In above formula I, k=0−1.44+0.708×(C1 minus 0.223 divided by 0.757)−0.660×(C2 minus 0.951 divided by 0.090)+0.088×(C3 minus 0.952 divided by 0.103)+3.048×(C4 minus 0.469 divided by 1.090)−0.460×(C5 minus 0.196 divided by 0.670).
A predicting result corresponding to Exp(−k) power value<1 is the candidate Klebsiella strain sensitive to MeropeneM, and a predicting result corresponding to Exp(−k) power value≥1 is the candidate Klebsiella strain resistant against MeropeneM.
One aspect of the present invention proposes a method of predicting sensitivity of Klebsiella against MeropeneM.
In this invention, after routine processing of the obtained microbial samples, necessary steps such as DNA extraction and sequencing can be carried out. Through bioinformatics process analysis, the state of the relevant features of the Klebsiella prediction system in the samples can be obtained. The feature state information can be imported into the system to predict the drug sensitivity of the samples. Compared to traditional methods, it has advantages such as simple operation, short detection time, and accurate species identification.
In order to effectively evaluate the performance of a prediction system, it is necessary to establish a dataset that is not involved in the establishment of the prediction system, and evaluate the accuracy of the prediction system on this dataset. This independent dataset is called the test set. The evaluation methods for system prediction effectiveness include F1 score, Precision, Recall, and confusion matrix.
The method of the present invention also has the following advantages. The present invention utilizes a test set to evaluate the accuracy of the system. The average accuracy of the method is 0.994, F1 score is 0.984, and recall score is 0.968. On the one hand, the present invention is less affected by subjective factors from such as operators, and has good detection stability; on the other hand, it achieves rapid and accurate identification of infectious pathogens and prediction of drug sensitivity of the test samples, meanwhile, it assists in diagnosis and rational & standardized medication treatment, has high throughput, and reduce medical costs.
In order to facilitate the understanding of the present invention, a more comprehensive description of the present invention will be provided in the embodiments below.
Unless otherwise defined, all technical and scientific terms used in this article have the same meanings as those commonly understood by a person skilled in the art of the present invention. The terms used in the specification of the present invention in this article are only for the purpose of describing specific embodiments and are not intended to limit the present invention.
The reagents used in the following embodiments are commercially available unless otherwise specified.
The 160 samples used in the experimental example of the present invention are pure cultures of Klebsiella strains isolated from clinical blood culture, from Beijing Union Medical College Hospital of the Chinese Academy of Medical Sciences.
All tested strains (strains) were identified as Klebsiella (scientific name: Genus Klebsiella, Latin name: Klebsiella Trevisan, systematic classification level: Genus) by mass spectrometry MALDI-TOF MS.
On the Illumina Novaseq NGS sequencing platform, these strains include 122 Klebsiella pneumoniae strains, 18 Klebsiella aerogenes strains, 7 Klebsiella oxytoca strains, 6 Klebsiella quasipneumoniae strains, 4 Klebsiella variicola strains, and 3 Klebsiella michiganensis strains, all of them are reported strains of the Klebsiella genus.
The above strains or strains can be obtained from common cases of Klebsiella pneumoniae pneumonia or from the applicant's laboratory. The applicant promises to distribute strains to the public within 20 years from the application date of the present invention for verifying the technical effects of the present invention.
This group of examples provides a system of predicting sensitivity of Klebsiella against MeropeneM. All examples of this group possesses the following common features: as shown in
In formula I,
In above formula I, the first pair of parentheses includes C1 minus 0.223 divided by 0.757; the second pair of parentheses includes C2 minus 0.951 divided by 0.090; the third pair of parentheses includes C3 minus 0.952 divided by 0.103; the fourth pair of parentheses includes C4 minus 0.469 divided by 1.090; and the fifth pair of parentheses includes C5 minus 0.196 divided by 0.670.
In above formula I, k=0−1.44+0.708×(C1 minus 0.223 divided by 0.757)−0.660×(C2 minus 0.951 divided by 0.090)+0.088×(C3 minus 0.952 divided by 0.103)+3.048×(C4 minus 0.469 divided by 1.090)−0.460×(C5 minus 0.196 divided by 0.670).
In some examples of this invention, said natural constant e=2.718281828459045.
In more specific examples, above genes are all known genes reported in the art, specifically as follows:
Klebsiella pneumoniae KpnE gene is Klebsiella pneumoniae KpnE gene recorded in “An unequivocal superbug: PDR Klebsiella pneumoniae with an arsenal of resistance and virulence factor genes”.
KPC-1 gene is KPC-1 gene recorded in “Novel Carbapenem-Hydrolyzing b-Lactamase, KPC-1, from a Carbapenem-Resistant Strain of Klebsiella pneumoniae”.
floR gene is floR gene recorded in “Spread of the florfenicol resistance floR gene among clinical Klebsiella pneumoniae isolates in China”.
In further examples, as shown in
Preferably, Exp(−k) power value calculated by the computing unit is transmitted to the result output unit.
In more further example, as shown in
In more specific example, said data-path is data transmission carrier well-known to a person skilled in the arts of computer and electronics. The data path can be selected from wired or wireless form, for example, it can be a wired path, line, wireless path, WiFi connection, wireless channel, etc.
Preferably, said independent variables include: values of C1, C2, C3, C4, C5;
Preferably, said experiment results comprise: the number of mphA, marA, Klebsiella pneumoniae KpnE, KPC-1, floR gene copies in the candidate Klebsiella strain.
The copy number of known genes in known strains can be routinely obtained by a person skilled in the arts of molecular biology and bioinformatics through conventional techniques such as sequencing and bioinformatics analysis. The mphA, marA, Klebsiella pneumoniae KpnE, KPC-1, and floR genes involved in the experiment results output by the experiment unit of the prediction system of the present invention are all reported genes in the art, and their gene information and primary structural sequences can be queried through the NCBI website or other known bioinformatics databases. By conducting whole genome sequencing of the candidate Klebsiella strain, the numbers of each of the aforementioned genes copies in the strain can be obtained.
In other specific examples, the number of mphA, marA, Klebsiella pneumoniae KpnE, KPC-1, floR gene copies in the candidate Klebsiella strain are obtained through a second-generation high-throughput sequencing method.
In more specific examples, the number of gene copies in the candidate Klebsiella strain=depth of gene contigs/depth of genome contigs;
The second-generation high-throughput sequencing method has a conventional technical meaning that is well-known to a person skilled in the art, while obtaining gene copy numbers using the second-generation high-throughput sequencing method is a conventional technical means that is well-known to a person skilled in the art. In some specific embodiments, the specific method for calculating the gene copy number is as follows:
Sequencing of strains is conducted by using second-generation high-throughput sequencing methods. The average sequencing depth is about 150×, and the approximate sequencing amount for the Klebsiella genus is about 1G. Using the depth of contigs obtained and calculated by SPAdes (v3.13.0) assembly software during the assembly process as the standard, the longest contigs segment is defined as the genome segment. Prokka software (1.14.6) is used to conduct gene predicting on contigs, and all gene cds and protein sequences on contigs are obtained. blat (v. 36) software and Diamond (v2.0.4.142) software are respectively used to compare the cds and protein sequences in the CARD database, sequences with a similarity greater than 90% are positive sequences, and annotation results for all resistance genes are obtained. The number of all gene copies on contigs are calculated using formula II as follows:
Formula II: the number of all gene copies on said contigs=depth of said contigs/depth of genome contigs
If a gene has two or more genome copies on different or the same contigs, the final number of gene copies is equal to the sum of all calculated number of the gene copies. Examples of calculation methods are as follows:
Assuming that the KPC-1 gene has 2 copies on one contigs and no copies on other contigs, the number of the KPC-1 gene copies is:
Assuming that the KPC-1 gene has one copy on one contig1 and contig2, but no copies on other contigs, the number of the KPC-1 gene copies is:
In more specific embodiments, the result output unit, experiment unit, and data input unit are all set as computer readable storage media on which computer programs are stored.
In some embodiments, when the computer program on the computer readable storage medium of the result output unit is executed by the processor, a method of comparing the Exp (−k) power value with 1 is implemented and the result is output;
That method of comparing the Exp (−k) power value with 1 is implemented and the result is output refers to:
In other embodiments, a method for calculating the number of gene copies is implemented when the computer program on the computer readable storage medium of the experiment unit is executed by the processor;
The method for calculating the number of gene copies is a conventional technical means well known to a person skilled in the art, and the specific steps are as follows:
In some embodiments, the computer program on the computer readable storage medium of the data input unit is executed by the processor to achieve dimensionless processing of the number of gene copies.
The dimensionless processing involves removing data dimensions or data units from the number of gene copies to obtain dimensionless values. Generally speaking, the data dimension or unit of the number of gene copies is: copy, number, or copies. In other embodiments, as shown in
This group examples provides a method of predicting sensitivity of Klebsiella against meropenem. All examples of this group possesses the following common features:
In formula I,
In above formula I, the first pair of parentheses includes C1 minus 0.223 divided by 0.757; the second pair of parentheses includes C2 minus 0.951 divided by 0.090; the third pair of parentheses includes C3 minus 0.952 divided by 0.103; the fourth pair of parentheses includes C4 minus 0.469 divided by 1.090; and the fifth pair of parentheses includes C5 minus 0.196 divided by 0.670.
In above formula I, k=0−1.44+0.708×(C1 minus 0.223 divided by 0.757)−0.660×(C2 minus 0.951 divided by 0.090)+0.088×(C3 minus 0.952 divided by 0.103)+3.048×(C4 minus 0.469 divided by 1.090)−0.460×(C5 minus 0.196 divided by 0.670).
A predicting result corresponding to Exp(−k) power value<1 is the candidate Klebsiella strain sensitive to MeropeneM, and a predicting result corresponding to Exp(−k) power value≥1 is the candidate Klebsiella strain resistant against MeropeneM.
In above formula I, e, as a mathematical constant, is the base of a natural logarithmic function, also known as a natural constant, natural base, or Euler number. It is an infinite non recurring decimal, which has the conventional technical meaning commonly understood by a ordinary technical person skilled in the art of mathematics. Its value is approximately: e=2.71828182845904523536.
In some examples of this invention, the value of said natural constant e is 2.718281828459045;
In some specific examples, the number of mphA, marA, Klebsiella pneumoniae KpnE, KPC-1, floR gene copies in the candidate Klebsiella strain are obtained through a second-generation high-throughput sequencing method.
In more specific examples, the number of gene copies in the candidate Klebsiella strain-depth of gene contigs/depth of genome contigs;
The second-generation high-throughput sequencing method has a conventional technical meaning that is well-known to a person skilled in the art, while obtaining the number of gene copies by using the second-generation high-throughput sequencing method is a conventional technical means that is well-known to a person skilled in the art.
In some specific embodiments, the specific method for calculating the number of gene copies is as follows:
Sequencing of strains is conducted by using second-generation high-throughput sequencing methods. The average sequencing depth is about 150×, and the approximate sequencing amount for the Klebsiella genus is about 1G. Using the depth of contigs obtained and calculated by SPAdes (v3.13.0) assembly software during the assembly process as the standard, the longest contigs segment is defined as the genome segment. Prokka software (1.14.6) is used to conduct gene predicting on contigs, and all gene cds and protein sequences on contigs are obtained. blat (v. 36) software and Diamond (v2.0.4.142) software are respectively used to compare the cds and protein sequences in the CARD database, sequences with a similarity greater than 90% are positive sequences, and annotation results for all resistance genes are obtained. The number of all gene copies on contigs are calculated using formula II as follows:
Formula II: the number of all gene copies on said contigs=depth of said contigs/depth of genome contigs
If a gene has two or more genome copies on different or the same contigs, the final number of gene copies is equal to the sum of all calculated number of the gene copies. Examples of calculation methods are as follows:
Assuming that the KPC-1 gene has 2 copies on one contigs and no copies on other contigs, the number of the KPC-1 gene copies is:
Assuming that the KPC-1 gene has one copy on one contig1 and contig2, but no copies on other contigs, the number of the KPC-1 gene copies is:
The prediction system of the present invention was evaluated using 160 clinical samples, and the comparison between the broth micro dilution classification results and the system prediction results of 160 clinical samples is shown in Table 1. In the table below, S represents sensitivity and R represents resistance.
The confusion matrix generated by the test result data is shown in Table 2:
Assuming TP (True Positive) represents the number of true positive cases, FP (False Positive) represents the number of false positive cases, FN (False Negative) represents the number of false negative cases, and TN (True Negative) represents the number of true negative cases. Precision refers to the proportion of positive samples in the positive case determined by the classifier. The recall rate refers to the proportion of predicted positive cases to the total positive cases. Accuracy refers to the proportion of correct judgments made by the classifier on the entire sample. F1 score is the harmonic mean of accuracy and recall, with a maximum of 1 and a minimum of 0. The calculation results of each indicator are as follows:
The above examples only express the embodiments of the present invention, and their description is more specific and detailed, but they cannot be understood as a limitation on the scope of the invention patent. It should be pointed out that for an ordinary technical person skilled in the art, several deformations and improvements can be made without departing from the concept of the present invention, all of which fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention patent should be based on the claims.
Number | Date | Country | Kind |
---|---|---|---|
202310065167.5 | Feb 2023 | CN | national |