Patent Application
20180148762
The present invention relates to a method of determining an infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment, a method of selecting a treatment of a patient suffering from an infection with a potentially resistant Shigella strain, and a method of determining an antimicrobial drug, e.g. antibiotic, resistance profile for bacterial microorganisms of Shigella species, as well as computer program products used in these methods.
Antibiotic resistance is a form of drug resistance whereby a sub-population of a microorganism, e.g. a strain of a bacterial species, can survive and multiply despite exposure to an antibiotic drug. It is a serious and health concern for the individual patient as well as a major public health issue. Timely treatment of a bacterial infection requires the analysis of clinical isolates obtained from patients with regard to antibiotic resistance, in order to select an efficacious therapy. Generally, for this purpose an association of the identified resistance with a certain microorganism (i.e. ID) is necessary.
Antibacterial drug resistance (ADR) represents a major health burden. According to the World Health Organization's antimicrobial resistance global report on surveillance, ADR leads to 25,000 deaths per year in Europe and 23,000 deaths per year in the US. In Europe, 2.5 million extra hospital days lead to societal cost of 1.5 billion euro. In the US, the direct cost of 2 million illnesses leads to 20 billion dollar direct cost. The overall cost is estimated to be substantially higher, reducing the gross domestic product (GDP) by up to 1.6%.
Shigella species are Gram-negative, facultative anaerobic rods belonging to the family of Enterobacteriaceae. The genus is divided into four serogroups with multiple serotypes: A (S dysenteriae); B (S flexneri); C (S boydii); and D (S sonnei). Shigella species are a major cause of diarrhoea and dysentery (diarrhoea with blood and mucus in the stools) throughout the world. The susceptibility of Shigella species to various antibiotics (e.g. ampicillin, trimetropin-sulfamethoxazole) decreased significantly from the 1970-80 to the 1990-2000.
In general the mechanisms for resistance of bacteria against antimicrobial treatments rely to a very substantial part on the organism's genetics. The respective genes or molecular mechanisms are either encoded in the genome of the bacteria or on plasmids that can be interchanged between different bacteria. The most common resistance mechanisms include:
In addition, some pathogens show natural resistance against drugs. For example, an organism can lack a transport system for an antibiotic or the target of the antibiotic molecule is not present in the organism.
Pathogens that are in principle susceptible to drugs can become resistant by modification of existing genetic material (e.g. spontaneous mutations for antibiotic resistance, happening in a frequency of one in about 100 mio bacteria in an infection) or the acquisition of new genetic material from another source. One example is horizontal gene transfer, a process where genetic material contained in small packets of DNA can be transferred between individual bacteria of the same species or even between different species. Horizontal gene transfer may happen by transduction, transformation or conjugation.
Generally, testing for susceptibility/resistance to antimicrobial agents is performed by culturing organisms in different concentration of these agents.
In brief, agar plates are inoculated with patient sample (e.g. urine, sputum, blood, stool) overnight. On the next day individual colonies are used for identification of organisms, either by culturing or using mass spectroscopy. Based on the identity of organisms new plates containing increasing concentration of drugs used for the treatment of these organisms are inoculated and grown for additional 12-24 hours. The lowest drug concentration which inhibits growth (minimal inhibitory concentration—MIC) is used to determine susceptibility/resistance for tested drugs. The process takes at least 2 to 3 working days during which the patient is treated empirically. A significant reduction of time-to-result is needed especially in patients with life-threatening disease and to overcome the widespread misuse of antibiotics.
Recent developments include PCR based test kits for fast bacterial identification (e.g. Biomerieux Biofire Tests, Curetis Unyvero Tests). With these test the detection of selected resistance loci is possible for a very limited number of drugs, but no correlation to culture based AST is given. Mass spectroscopy is increasingly used for identification of pathogens in clinical samples (e.g. Bruker Biotyper), and research is ongoing to establish methods for the detection of susceptibility/resistance against antibiotics.
For some drugs such it is known that at least two targets are addressed, e.g. in case of Ciprofloxacin (drug bank ID 00537; http://www.drugbank.ca/drugs/DB00537) targets include DNA Topoisomerase IV, DNA Topoisomerase II and DNA Gyrase. It can be expected that this is also the case for other drugs although the respective secondary targets have not been identified yet. In case of a common regulation, both relevant genetic sites would naturally show a co-correlation or redundancy.
It is known that drug resistance can be associated with genetic polymorphisms. This holds for viruses, where resistance testing is established clinical practice (e.g. HIV genotyping). More recently, it has been shown that resistance has also genetic causes in bacteria and even higher organisms, such as humans where tumors resistance against certain cytostatic agents can be linked to genomic mutations.
Wozniak et al. (BMC Genomics 2012, 13(Suppl 7):S23) disclose genetic determinants of drug resistance in Staphylococcus aureus based on genotype and phenotype data. Stoesser et al. disclose prediction of antimicrobial susceptibilities for Escherichia coli and Klebsiella pneumoniae isolates using whole genomic sequence data (J Antimicrob Chemother 2013; 68: 2234-2244).
Chewapreecha et al (Chewapreecha et al (2014) Comprehensive Identification of single nucleotid polymorphisms associated with beta-lactam resistance within pneumococcal mosaic genes. PLoS Genet 10(8): e1004547) used a comparable approach to identify mutations in gram-positive Streptococcus Pneumonia.
The fast and accurate detection of infections with Shigella species and the prediction of response to anti-microbial therapy represent a high unmet clinical need.
This need is addressed by the present invention.
The present inventors addressed this need by carrying out whole genome sequencing of a large cohort of Shigella clinical isolates and comparing the genetic mutation profile to classical culture based antimicrobial susceptibility testing with the goal to develop a test which can be used to detect bacterial susceptibility/resistance against antimicrobial drugs using molecular testing.
The inventors performed extensive studies on the genome of bacteria of Shigella species either susceptible or resistant to antimicrobial, e.g. antibiotic, drugs. Based on this information, it is now possible to provide a detailed analysis on the resistance pattern of Shigella strains based on individual genes or mutations on a nucleotide level. This analysis involves the identification of a resistance against individual antimicrobial, e.g. antibiotic, drugs as well as clusters of them. This allows not only for the determination of a resistance to a single antimicrobial, e.g. antibiotic, drug, but also to groups of antimicrobial drugs, e.g. antibiotics such as lactam or quinolone antibiotics, or even to all relevant antibiotic drugs.
Therefore, the present invention will considerably facilitate the selection of an appropriate antimicrobial, e.g. antibiotic, drug for the treatment of a Shigella infection in a patient and thus will largely improve the quality of diagnosis and treatment.
According to a first aspect, the present invention discloses a diagnostic method of determining an infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment, which can be also described as a method of determining an antimicrobial drug, e.g. antibiotic, resistant Shigella infection of a patient, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least two genes from the group of genes listed in Table 1 or Table 2 below, wherein the presence of said at least two mutations is indicative of an infection with an antimicrobial drug resistant, e.g. antibiotic resistant, Shigella strain in said patient.
An infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment herein means an infection of a patient with Shigella species wherein it is unclear if the Shigella species is susceptible to treatment with a specific antimicrobial drug or if it is resistant to the antimicrobial drug.
In step b) above, as well as corresponding steps, at least one mutation in at least two genes is determined, so that in total at least two mutations are determined, wherein the two mutations are in different genes.
According to a second aspect, the present invention relates to a method of selecting a treatment of a patient suffering from an infection with a potentially resistant Shigella stain, e.g. from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least two genes from the group of genes listed in Table 1 or Table 2 above, wherein the presence of said at least two mutations is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs; and
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection.
A third aspect of the present invention relates to a method of determining an antimicrobial drug, e.g. antibiotic, resistance profile for bacterial microorganisms of Shigella species, comprising:
obtaining or providing a first data set of gene sequences of a plurality of clinical isolates of Shigella species;
providing a second data set of antimicrobial drug, e.g. antibiotic, resistance of the plurality of clinical isolates of Shigella species;
aligning the gene sequences of the first data set to at least one, preferably one, reference genome of Shigella, and/or assembling the gene sequence of the first data set, at least in part;
analyzing the gene sequences of the first data set for genetic variants to obtain a third data set of genetic variants;
correlating the third data set with the second data set and statistically analyzing the correlation; and
determining the genetic sites in the genome of Shigella associated with antimicrobial drug, e.g. antibiotic, resistance.
In addition, the present invention relates in a fourth aspect to a method of determining an antimicrobial drug, e.g. antibiotic, resistance profile for a bacterial microorganism belonging to the species Shigella comprising the steps of
a) obtaining or providing a sample containing or suspected of containing the bacterial microorganism;
b) determining the presence of a mutation in at least one gene of the bacterial microorganism as determined by the method according to the third aspect of the present invention;
wherein the presence of a mutation is indicative of a resistance to an antimicrobial, e.g. antibiotic, drug.
Furthermore, the present invention discloses in a fifth aspect a diagnostic method of determining an infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment, which can, like in the first aspect, also be described as method of determining an antimicrobial drug, e.g. antibiotic, resistant Shigella infection of a patient, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing a bacterial microorganism belonging to the species Shigella from the patient;
b) determining the presence of at least one mutation in at least one gene of the bacterial microorganism belonging to the species Shigella as determined by the method according to the third aspect of the present invention, wherein the presence of said at least one mutation is indicative of an antimicrobial drug, e.g. antibiotic, resistant Shigella infection in said patient.
Also disclosed is in a sixth aspect a method of selecting a treatment of a patient suffering from an infection with a potentially resistant Shigella strain, e.g. from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing a bacterial microorganism belonging to the species Shigella from the patient;
b) determining the presence of at least one mutation in at least one gene of the bacterial microorganism belonging to the species Shigella as determined by the method according to the third aspect of the present invention, wherein the presence of said at least one mutation is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs; and
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection.
A seventh aspect of the present invention relates to a method of acquiring, respectively determining, an antimicrobial drug, e.g. antibiotic, resistance profile for a bacterial microorganisms of Shigella species, comprising:
obtaining or providing a first data set of gene sequences of a clinical isolate of Shigella species;
providing a second data set of antimicrobial drug, e.g. antibiotic, resistance of a plurality of clinical isolates of Shigella species;
aligning the gene sequences of the first data set to at least one, preferably one, reference genome of Shigella, and/or assembling the gene sequence of the first data set, at least in part;
analyzing the gene sequences of the first data set for genetic variants to obtain a third data set of genetic variants of the first data set;
correlating the third data set with the second data set and statistically analyzing the correlation; and
determining the genetic sites in the genome of Shigella of the first data set associated with antimicrobial drug, e.g. antibiotic, resistance.
According to an eighth aspect, the present invention discloses a computer program product comprising executable instructions which, when executed, perform a method according to the third, fourth, fifth, sixth or seventh aspect of the present invention.
Further aspects and embodiments of the invention are disclosed in the dependent claims and can be taken from the following description, FIGURES and examples, without being limited thereto.
The enclosed drawings should illustrate embodiments of the present invention and convey a further understanding thereof. In connection with the description they serve as explanation of concepts and principles of the invention. Other embodiments and many of the stated advantages can be derived in relation to the drawings. The elements of the drawings are not necessarily to scale towards each other. Identical, functionally equivalent and acting equal features and components are denoted in the FIGURES of the drawings with the same reference numbers, unless noted otherwise.
Unless defined otherwise, technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
An “antimicrobial drug” in the present invention refers to a group of drugs that includes antibiotics, antifungals, antiprotozoals, and antivirals. According to certain embodiments, the antimicrobial drug is an antibiotic.
The term “nucleic acid molecule” refers to a polynucleotide molecule having a defined sequence. It comprises DNA molecules, RNA molecules, nucleotide analog molecules and combinations and derivatives thereof, such as DNA molecules or RNA molecules with incorporated nucleotide analogs or cDNA.
The term “nucleic acid sequence information” relates to information which can be derived from the sequence of a nucleic acid molecule, such as the sequence itself or a variation in the sequence as compared to a reference sequence.
The term “mutation” relates to a variation in the sequence as compared to a reference sequence. Such a reference sequence can be a sequence determined in a predominant wild type organism or a reference organism, e.g. a defined and known bacterial strain or substrain. A mutation is for example a deletion of one or multiple nucleotides, an insertion of one or multiple nucleotides, or substitution of one or multiple nucleotides, duplication of one or a sequence of multiple nucleotides, translocation of one or a sequence of multiple nucleotides, and, in particular, a single nucleotide polymorphism (SNP).
In the context of the present invention a “sample” is a sample which comprises at least one nucleic acid molecule from a bacterial microorganism. Examples for samples are: cells, tissue, body fluids, biopsy specimens, blood, urine, saliva, sputum, plasma, serum, cell culture supernatant, swab sample and others. According to certain embodiments, the sample is a patient sample (clinical isolate).
New and highly efficient methods of sequencing nucleic acids referred to as next generation sequencing have opened the possibility of large scale genomic analysis. The term “next generation sequencing” or “high throughput sequencing” refers to high-throughput sequencing technologies that parallelize the sequencing process, producing thousands or millions of sequences at once. Examples include Massively Parallel Signature Sequencing (MPSS), Polony sequencing, 454 pyrosequencing, Illumina (Solexa) sequencing, SOLiD sequencing, Ion semiconductor sequencing, DNA nanoball sequencing, Helioscope™ single molecule sequencing, Single Molecule SMRT™ sequencing, Single Molecule real time (RNAP) sequencing, Nanopore DNA sequencing, Sequencing By Hybridization, Amplicon Sequencing, GnuBio.
Within the present description the term “microorganism” comprises the term microbe. The type of microorganism is not particularly restricted, unless noted otherwise or obvious, and, for example, comprises bacteria, viruses, fungi, microscopic algae and protozoa, as well as combinations thereof. According to certain aspects, it refers to one or more Shigella species, particularly Shigella boydii, Shigella dysenteriae, Shigella flexneri, Shigella sonnei and/or other Shigella species.
A reference to a microorganism or microorganisms in the present description comprises a reference to one microorganism as well a plurality of microorganisms, e.g. two, three, four, five, six or more microorganisms.
A vertebrate within the present invention refers to animals having a vertebrae, which includes mammals—including humans, birds, reptiles, amphibians and fishes. The present invention thus is not only suitable for human medicine, but also for veterinary medicine.
According to certain embodiments, the patient in the present methods is a vertebrate, more preferably a mammal and most preferred a human patient.
Before the invention is described in exemplary detail, it is to be understood that this invention is not limited to the particular component parts of the process steps of the methods described herein as such methods may vary. It is also to be understood that the terminology used herein is for purposes of describing particular embodiments only, and is not intended to be limiting. It must be noted that, as used in the specification and the appended claims, the singular forms “a,” “an” and “the” include singular and/or plural referents unless the context clearly dictates otherwise. For example, the term “a” as used herein can be understood as one single entity or in the meaning of “one or more” entities. It is also to be understood that plural forms include singular and/or plural referents unless the context clearly dictates otherwise. It is moreover to be understood that, in case parameter ranges are given which are delimited by numeric values, the ranges are deemed to include these limitation values.
Regarding the dosage of the antimicrobial, e.g. antibiotic, drugs, it is referred to the established principles of pharmacology in human and veterinary medicine. For example, Forth, Henschler, Rummel “Allgemeine und spezielle Pharmakologie und Toxikologie”, 9th edition, 2005, pp. 781-919, might be used as a guideline. Regarding the formulation of a ready-to-use medicament, reference is made to “Remington, The Science and Practice of Pharmacy”, 22nd edition, 2013, pp. 777-1070.
Assembling of a gene sequence can be carried out by any known method and is not particularly limited.
According to certain embodiments, mutations that were found using alignments can also be compared or matched with alignment-free methods, e.g. for detecting single base exchanges, for example based on contigs that were found by assemblies. For example, reads obtained from sequencing can be assembled to contigs and the contigs can be compared to each other.
According to a first aspect, the present invention relates to a diagnostic method of determining an infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment, which can also be described as method of determining an antimicrobial drug, e.g. antibiotic, resistant Shigella infection of a patient, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least two genes from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, and pdxA, or from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and SSON53_06085, wherein the presence of said at least two mutations is indicative of an infection with an antimicrobial, e.g. antibiotic, resistant Shigella strain in said patient.
In this method, as well as the other methods of the invention, the sample can be provided or obtained in any way, preferably non-invasive, and can be e.g. provided as an in vitro sample or prepared as in vitro sample.
According to certain aspects, mutations in at least two, three, four, five, six, seven, eight, nine or ten genes are determined in any of the methods of the present invention, e.g. in at least two genes or in at least three genes. Instead of testing only single genes or mutants, a combination of several variant positions can improve the prediction accuracy and further reduce false positive findings that are influenced by other factors. Therefore, it is in particular preferred to determine the presence of a mutation in 2, 3, 4, 5, 6, 7, 8 or 9 (or more) genes selected from Table 1 or 2.
For the above genes, i.e. the genes also denoted in Tables 1 and 2, the highest probability of a resistance to at least one antimicrobial drug, e.g. antibiotic, could be observed, with p-values smaller than 10−30, particularly smaller than 10−40, indicating the high significance of the values (n=470; α=0.05). Details regarding Tables 1 and 2 can be taken from Tables 3 and 4 (4a, 4b, 4c) disclosed in the Examples. Having at least two genes with mutations determined, a high probability of an antimicrobial drug, e.g. antibiotic, resistance could be determined. The genes in Table 1 thereby represent the 50 best genes for which a mutation was observed in the genomes of Shigella species, whereas the genes in Table 2 represent the 50 best genes for which a cross-correlation could be observed for the antimicrobial drug, e.g. antibiotic, susceptibility testing for Shigella species as described below.
According to certain embodiments, the obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient in this method—as well as the other methods of the invention—can comprise the following:
A sample of a vertebrate, e.g. a human, e.g. is provided or obtained and nucleic acid sequences, e.g. DNA or RNA sequences, are recorded by a known method for recording nucleic acid, which is not particularly limited. For example, nucleic acid can be recorded by a sequencing method, wherein any sequencing method is appropriate, particularly sequencing methods wherein a multitude of sample components, as e.g. in a blood sample, can be analyzed for nucleic acids and/or nucleic acid fragments and/or parts thereof contained therein in a short period of time, including the nucleic acids and/or nucleic acid fragments and/or parts thereof of at least one microorganism of interest, particularly of at least one Shigella species. For example, sequencing can be carried out using polymerase chain reaction (PCR), particularly multiplex PCR, or high throughput sequencing or next generation sequencing, preferably using high-throughput sequencing. For sequencing, preferably an in vitro sample is used.
The data obtained by the sequencing can be in any format, and can then be used to identify the nucleic acids, and thus genes, of the microorganism, e.g. of Shigella species, to be identified, by known methods, e.g. fingerprinting methods, comparing genomes and/or aligning to at least one, or more, genomes of one or more species of the microorganism of interest, i.e. a reference genome, etc., forming a third data set of aligned genes for a Shigella species—discarding additional data from other sources, e.g. the vertebrate. Reference genomes are not particularly limited and can be taken from several databases. Depending on the microorganism, different reference genomes or more than one reference genomes can be used for aligning. Using the reference genome—as well as also the data from the genomes of the other species, e.g. Shigella species—mutations in the genes for each species and for the whole multitude of samples of different species, e.g. Shigella species, can be obtained.
For example, it is useful in genome-wide association studies to reference the points of interest, e.g. mutations, to one constant reference for enhanced standardization. In case of the human with a high consistency of the genome and 99% identical sequences among individuals this is easy and represents the standard, as corresponding reference genomes are available in databases. In case of organisms that trigger infectious diseases (e.g. bacteria and viruses) this is much more difficult, though. One possibility is to fall back on a virtual pan genome which contains all sequences of a certain genus. A further possibility is the analysis of all available references, which is much more complex. Therein all n references from a database (e.g. RefSeq) are extracted and compared with the newly sequenced bacterial genomes k. After this, matrices (% of mapped reads, % of covered genome) are applied to estimate which reference is best suited to all new bacteria. However, n×k complete alignments are carried out. Having a big number of references, though, stable results can be obtained, as is the case for Shigella.
According to certain embodiments, the genomes of Shigella species are referenced to one reference genome. However, it is not excluded that for other microorganisms more than one reference genome is used. In the present methods, the reference genome of Shigella is NC 016822 as annotated at the NCBI according to certain embodiments.
The reference sequence was obtained from Shigella strain NC_016822 (http://www.genome.jp/dbget-bin/www_bget?refseq+NC_016822), and can be found as SEQ ID NO 1 in the sequence listing.
LOCUS NC_016822 4988504 bp DNA circular CON 7 Feb. 2015
DEFINITION Shigella sonnei 53G main chromosome, complete genome.
ACCESSION NC_016822
VERSION NC_016822.1 GI:377520096
DBLINK BioProject: PRJNA224116
KEYWORDS RefSeq; complete genome; complete replicon.
SOURCE Shigella sonnei 53G
REFERENCE 1
REFERENCE 2 (bases 1 to 4988504)
Institute, Pathogen Sequencing Unit, Wellcome Trust Genome Campus,
Hinxton, Cambridge, Cambridgeshire CB10 1SA, UNITED KINGDOM
Alternatively or in addition, the gene sequence of the first data set can be assembled, at least in part, with known methods, e.g. by de-novo assembly or mapping assembly. The sequence assembly is not particularly limited, and any known genome assembler can be used, e.g. based on Sanger, 454, Solexa, Illumina, SOLid technologies, etc., as well as hybrids/mixtures thereof.
According to certain embodiments, the data of nucleic acids of different origin than the microorganism of interest, e.g. Shigella species, can be removed after the nucleic acids of interest are identified, e.g. by filtering the data out. Such data can e.g. include nucleic acids of the patient, e.g. the vertebrate, e.g. human, and/or other microorganisms, etc. This can be done by e.g. computational subtraction, as developed by Meyerson et al. 2002. For this, also aligning to the genome of the vertebrate, etc., is possible. For aligning, several alignment-tools are available. This way the original data amount from the sample can be drastically reduced.
Also after such removal of “excess” data, fingerprinting and/or aligning, and/or assembly, etc. can be carried out, as described above, forming a third data set of aligned and/or assembled genes for a Shigella species.
Using these techniques, genes with mutations of the microorganism of interest, e.g. Shigella species, can be obtained for various species.
When testing these same species for antimicrobial drug, e.g. antibiotic, susceptibility of a number of antimicrobial drugs, e.g. antibiotics, e.g. using standard culturing methods on dishes with antimicrobial drug, e.g. antibiotic, intake, as e.g. described below, the results of these antimicrobial drug, e.g. antibiotic, susceptibility tests can then be cross-referenced/correlated with the mutations in the genome of the respective microorganism, e.g. Shigella. Using several, e.g. 50 or more than 50, 100 or more than 100, 200 or more than 200, 300 or more than 300, 400 or more than 400, or 450 or more than 450 different species of a microorganism, e.g. different Shigella species, statistical analysis can be carried out on the obtained cross-referenced data between mutations and antimicrobial drug, e.g. antibiotic, susceptibility for these number of species, using known methods.
Regarding culturing methods, samples can be e.g. cultured overnight. On the next day individual colonies can be used for identification of organisms, either by culturing or using mass spectroscopy. Based on the identity of organisms new plates containing increasing concentration of antibiotics used for the treatment of these organisms are inoculated and grown for additional 12-24 hours. The lowest drug concentration which inhibits growth (minimal inhibitory concentration—MIC) can be used to determine susceptibility/resistance for tested antibiotics.
Correlation of the nucleic acid/gene mutations with antimicrobial drug, e.g. antibiotic, resistance can be carried out in a usual way and is not particularly limited. For example, resistances can be correlated to certain genes or certain mutations, e.g. SNPs, in genes. After correlation, statistical analysis can be carried out.
In addition, statistical analysis of the correlation of the gene mutations with antimicrobial drug, e.g. antibiotic, resistance is not particularly limited and can be carried out, depending on e.g. the amount of data, in different ways, for example using analysis of variance (ANOVA) or Student's t-test, for example with a sample size n of 50 or more, 100 or more, 200 or more, 300 or more, 400 or more or 450 or more, and a level of significance (α-error-level) of e.g. 0.05 or smaller, e.g. 0.05, preferably 0.01 or smaller. A statistical value can be obtained for each gene and/or each position in the genome as well as for all antibiotics tested, a group of antibiotics or a single antibiotic. The obtained p-values can also be adapted for statistical errors, if needed.
For statistically sound results a multitude of individuals should be sampled, with n=50 or more, 100 or more, 200 or more, 300 or more, 400 or more or 450 or more, and a level of significance (α-error-level) of e.g. 0.05 or smaller, e.g. 0.05, preferably 0.01 or smaller. According to certain embodiments, particularly significant results can be obtained for n=200 or more, 300 or more, 400 or more or 450 or more.
For statistically sound results a multitude of individuals should be sampled, with n=50, 100, 200, 300, 400 or 450, and a level of significance (α-error-level) of e.g. 0.05 or smaller, e.g. 0.05, preferably 0.01 or smaller. According to certain embodiments, particularly significant results can be obtained for n=200, 300, 400 or 450. According to certain embodiments, p-values can be adjusted for multiple testing using approaches as bonferroni adjustment or by controlling the false discovery rate ((FDR)
After the above procedure has been carried out for more than 450, e.g. 470, individual species of Shigella, the data disclosed in Tables 1 and 2 were obtained for the statistically best correlations between gene mutations and antimicrobial drug, e.g. antibiotic, resistances. Thus, mutations in these genes were proven as valid markers for antimicrobial drug, e.g. antibiotic, resistance.
According to a further aspect, the present invention relates in a second aspect to a method of selecting a treatment of a patient suffering from an infection with a potentially resistant Shigella stain, e.g. from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least two genes from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, and pdxA, or from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and SSON53_06085, wherein the presence of said at least two mutations is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs; and
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection.
In this method, the steps a) of obtaining or providing a sample and b) of determining the presence of at least one mutation are as in the method of the first aspect.
The identification of the at least one or more antimicrobial, e.g. antibiotic, drug in step c) is then based on the results obtained in step b) and corresponds to the antimicrobial, e.g. antibiotic, drug(s) that correlate(s) with the mutations. Once these antimicrobial drugs, e.g. antibiotics, are ruled out, the remaining antimicrobial drugs, e.g. antibiotic drugs/antibiotics, can be selected in step d) as being suitable for treatment.
In the description, references to the first and second aspect also apply to the 19th, 20th, 21st and 22nd aspect, referring to the same genes, unless clear from the context that they don't apply.
According to certain embodiments, the antimicrobial drug, e.g. antibiotic, in the method of the first or second aspect, as well as in the other methods of the invention, is at least one selected from the group of β-lactams, β-lactam inhibitors, quinolines and derivatives thereof, aminoglycosides, polyketides, respectively tetracyclines, and folate synthesis inhibitors.
In the methods of the invention the resistance of Shigella to one or more antimicrobial, e.g. antibiotic, drugs can be determined according to certain embodiments.
Results with particular significance and low p-values can be obtained in the methods of the first and second aspect when the antimicrobial, e.g. antibiotic, drug is selected from lactam antibiotics, and particularly when the presence of a mutation in the following genes is determined: SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, and/or pdxA, or SSON53_12105, SSON53_12455, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and/or SSON53_06085.
For the lactam antibiotics, the p-values are that low for these genes that a statistically significant determination of antibiotic susceptibility is possible in particular.
According to certain embodiments of the first and/or second aspect of the invention the antimicrobial, e.g. antibiotic, drug is selected from polyketide antibiotics, preferably tetracycline antibiotics, and the presence of a mutation in the following genes is determined: metH, SSON53_24790, SSON53_08770, SSON53_24780, SSON53_15280, SSON53_10655, SSON53_13555, and/or SSON53_06415, or metH, SSON53_24790, SSON53_08770, SSON53_24780, SSON53_15280, SSON53_10655, SSON53_13555, and/or SSON53_06415.
According to certain embodiments, the antimicrobial drug is an antibiotic/antibiotic drug.
According to certain embodiments of the first and/or second aspect of the invention, determining the nucleic acid sequence information or the presence of a mutation comprises determining the presence of a single nucleotide at a single position in a gene. Thus the invention comprises methods wherein the presence of a single nucleotide polymorphism or mutation at a single nucleotide position is detected.
According to certain embodiments, the antibiotic drug in the methods of the present invention is selected from the group consisting of Amoxicillin/K Clavulanate (AUG), Ampicillin (AM), Aztreonam (AZT), Cefazolin (CFZ), Cefepime (CPE), Cefotaxime (CFT), Ceftazidime (CAZ), Ceftriaxone (CAX), Cefuroxime (CRM), Cephalotin (CF), Ciprofloxacin (CP), Ertapenem (ETP), Gentamicin (GM), Imipenem (IMP), Levofloxacin (LVX), Meropenem (MER), Piperacillin/Tazobactam (P/T), Ampicillin/Sulbactam (A/S), Tetracycline (TE), Tobramycin (TO), and Trimethoprim/Sulfamethoxazole (T/S).
The inventors have surprisingly found that mutations in certain genes are indicative not only for a resistance to one single antimicrobial, e.g. antibiotic, drug, but to groups containing several drugs, sometimes even for drugs from different drug classes.
According to certain embodiments of the first and/or second aspect of the invention, the gene is from Table 1, the antibiotic drug is selected from lactam antibiotics and a mutation in at least one of the following genes is detected with regard to reference genome NC_016822: SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA.
According to certain embodiments of the first and/or second aspect of the invention, the gene is from Table 1, the antibiotic drug is selected from polyketide, preferably tetracycline antibiotics and a mutation in at least one of the following genes is detected with regard to reference genome NC_016822: metH, SSON53_24790, SSON53_08770, SSON53_24780, SSON53_15280, SSON53_10655, SSON53_13555, SSON53_06415.
According to certain embodiments of the first and/or second aspect of the invention, the gene is from Table 2, the antibiotic drug is selected from lactam antibiotics and a mutation in at least one of the following genes is detected with regard to reference genome NC_016822: SSON53_12105, SSON53_12455, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, SSON53_06085.
According to certain embodiments of the first and/or second aspect of the invention, the gene is from Table 2, the antibiotic drug is selected from polyketide, preferably tetracycline antibiotics and a mutation in at least one of the following genes is detected with regard to reference genome NC_016822: metH, SSON53_24790, SSON53_08770, SSON53_24780, SSON53_15280, SSON53_10655, SSON53_13555, SSON53_06415.
For specific antimicrobial drugs, e.g. antibiotics, specific positions in the above genes can be determined where a high statistical significance is observed. The inventors found that, apart from the above genes indicative of a resistance against antibiotics, also single nucleotide polymorphisms (=SNP's) may have a high significance for the presence of a resistance against defined antibiotic drugs. The analysis of these polymorphisms on a nucleotide level may further improve and accelerate the determination of a drug resistance to antimicrobial drugs, e.g. antibiotics, in Shigella.
According to certain embodiments of the first and/or second aspect of the invention, the gene is from Table 1, the antibiotic drug is selected from lactam antibiotics and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 2241215, 2309519, 2233949, 2314492, 4595667, 87512, 4502174, 2543074, 1159191, 2325584, 1303652, 4707429, 3134050, 4648709, 1655745, 949314, 1427818, 2636621, 1429397, 1690889, 4856482, 1723070, 4857482, 4705326, 862961, 2876708, 887353, 980321, 4429947, 4431976, 3376607, 4565385, 750613, 1781189, 4821412, 2533864, 1511130, 1513422, 1901825, 1990704, 2506615, 2539010, 2803779, 2877526, 3239464, 545532, 1257335, 906989, 3237959, 60014.
According to certain embodiments of the first and/or second aspect of the invention, the gene is from Table 1, the antibiotic drug is selected from polyketide, preferably tetracycline antibiotics and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 4595667, 4707429, 1655745, 4705326, 2876708, 1990704, 2539010, 1257335.
According to certain embodiments of the first and/or second aspect of the invention, the gene is from Table 2, the antibiotic drug is selected from lactam antibiotics and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 2241215, 2309519, 2314492, 4595667, 87512, 4502174, 2543074, 1159191, 2325584, 4707429, 3134050, 1655745, 949314, 1427818, 2636621, 1429397, 1690889, 4856482, 1723070, 4857482, 4705326, 862961, 2876708, 980321, 4429947, 4431976, 3376607, 4565385, 750613, 1781189, 4821412, 2533864, 1511130, 1513422, 1901825, 1990704, 2506615, 2539010, 2803779, 2877526, 3239464, 545532, 1257335, 906989, 3237959, 60014, 921872, 983367, 1182426, 1200860.
According to certain embodiments of the first and/or second aspect of the invention, the gene is from Table 2, the antibiotic drug is selected from polyketide, preferably tetracycline antibiotics and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 4595667, 4707429, 1655745, 4705326, 2876708, 1990704, 2539010, 1257335.
According to certain embodiments of the first and/or second aspect of the invention, the antibiotic drug is CF and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 2241215, 2309519, 2314492, 4595667, 87512, 4502174, 2543074, 1159191, 2325584, 4707429, 3134050, 1655745, 949314, 1427818, 2636621, 1429397, 1690889, 4856482, 1723070, 4857482, 4705326, 862961, 2876708, 980321, 4429947, 4431976, 3376607, 4565385, 750613, 1781189, 4821412, 2533864, 1511130, 1513422, 1901825, 1990704, 2506615, 2539010, 2803779, 2877526, 3239464, 545532, 1257335, 906989, 3237959, 60014, 921872, 983367, 1182426, 1200860.
According to certain embodiments of the first and/or second aspect of the invention, the antibiotic drug is CFZ and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 3134050.
According to certain embodiments of the first and/or second aspect of the invention, the antibiotic drug is A/S and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 2241215, 2309519, 2314492, 4595667, 87512, 4502174, 2543074, 1159191, 2325584, 4707429, 3134050, 1655745, 949314, 1427818, 2636621, 1429397, 1690889, 4856482, 1723070, 4857482, 4705326, 862961, 2876708, 980321, 4429947, 4431976, 3376607, 4565385, 750613, 1781189, 4821412, 2533864, 1511130, 1513422, 1901825, 1990704, 2506615, 2539010, 2803779, 2877526, 3239464, 545532, 1257335, 906989, 3237959, 60014, 921872, 983367, 1182426, 1200860.
According to certain embodiments of the first and/or second aspect of the invention, the antibiotic drug is CRM and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 2241215, 2309519, 2314492, 4595667, 87512, 4502174, 2543074, 1159191, 2325584, 4707429, 3134050, 1655745, 949314, 1427818, 2636621, 1429397, 1690889, 4856482, 1723070, 4857482, 4705326, 862961, 2876708, 980321, 4429947, 4431976, 3376607, 4565385, 750613, 1781189, 4821412, 2533864, 1511130, 1513422, 1901825, 1990704, 2506615, 2539010, 2803779, 2877526, 3239464, 545532, 1257335, 906989, 3237959, 60014, 921872, 983367, 1182426, 1200860.
According to certain embodiments of the first and/or second aspect of the invention, the antibiotic drug is P/T and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 2241215, 4595667, 4502174, 2543074, 1159191, 2325584, 4707429, 3134050, 1655745, 949314, 1427818, 1429397, 1690889, 4856482, 1723070, 4857482, 4705326, 862961, 2876708, 980321, 4429947, 4431976, 3376607, 4565385, 750613, 1781189, 4821412, 2533864, 1511130, 1513422, 1901825, 1990704, 2506615, 2539010, 2803779, 2877526, 3239464, 545532, 1257335, 906989, 3237959, 60014, 921872, 983367, 1182426, 1200860.
According to certain embodiments of the first and/or second aspect of the invention, the antibiotic drug is AM and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 2241215, 2309519, 2314492, 4595667, 87512, 4502174, 2543074, 1159191, 2325584, 4707429, 3134050, 1655745, 949314, 1427818, 2636621, 1429397, 1690889, 4856482, 1723070, 4857482, 4705326, 862961, 2876708, 980321, 4429947, 4431976, 3376607, 4565385, 750613, 1781189, 4821412, 2533864, 1511130, 1513422, 1901825, 1990704, 2506615, 2539010, 2803779, 2877526, 3239464, 545532, 1257335, 906989, 3237959, 60014, 921872, 983367, 1182426, 1200860.
According to certain embodiments of the first and/or second aspect of the invention, the antibiotic drug is AUG and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 2309519, 2314492, 4595667, 87512, 1159191, 2325584, 4707429, 3134050, 1427818, 2636621, 1429397, 2877526, 1200860.
According to certain embodiments of the first and/or second aspect of the invention, the antibiotic drug is TE and a mutation in at least one of the following nucleotide positions is detected with regard to reference genome NC_016822: 4595667, 4707429, 1655745, 4705326, 2876708, 1990704, 2539010, 1257335.
According to certain embodiments of the first and/or second aspect of the invention, the resistance of a bacterial microorganism belonging to the species Shigella against 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 or 16, 17, 18, 19, 20 or 21 antibiotic drugs is determined.
According to certain embodiments of the first and/or second aspect of the invention, a detected mutation is a mutation leading to an altered amino acid sequence in a polypeptide derived from a respective gene in which the detected mutation is located. According to this aspect, the detected mutation thus leads to a truncated version of the polypeptide (wherein a new stop codon is created by the mutation) or a mutated version of the polypeptide having an amino acid exchange at the respective position.
According to certain embodiments of the first and/or second aspect of the invention, determining the nucleic acid sequence information or the presence of a mutation comprises determining a partial sequence or an entire sequence of the at least two genes.
According to certain embodiments of the first and/or second aspect of the invention, determining the nucleic acid sequence information or the presence of a mutation comprises determining a partial or entire sequence of the genome of the Shigella species, wherein said partial or entire sequence of the genome comprises at least a partial sequence of said at least two genes.
According to certain embodiments of the first and/or second aspect of the invention, determining the nucleic acid sequence information or the presence of a mutation comprises using a next generation sequencing or high throughput sequencing method. According to preferred embodiments of the first and/or second aspect of the invention, a partial or entire genome sequence of the bacterial organism of Shigella species is determined by using a next generation sequencing or high throughput sequencing method.
In a further, third aspect, the present invention relates to a method of determining an antimicrobial drug, e.g. antibiotic, resistance profile for bacterial microorganisms of Shigella species, comprising:
obtaining or providing a first data set of gene sequences of a plurality of clinical isolates of Shigella species;
providing a second data set of antimicrobial drug, e.g. antibiotic, resistance of the plurality of clinical isolates of Shigella species;
aligning the gene sequences of the first data set to at least one, preferably one, reference genome of Shigella, and/or assembling the gene sequence of the first data set, at least in part;
analyzing the gene sequences of the first data set for genetic variants to obtain a third data set of genetic variants;
correlating the third data set with the second data set and statistically analyzing the correlation; and
determining the genetic sites in the genome of Shigella associated with antimicrobial drug, e.g. antibiotic, resistance.
The different steps can be carried out as described with regard to the method of the first aspect of the present invention.
When referring to the second data set, wherein the second data set e.g. comprises, respectively is, a set of antimicrobial drug, e.g. antibiotic, resistances of a plurality of clinical isolates, this can, within the scope of the invention, also refer to a self-learning data base that, whenever a new sample is analyzed, can take this sample into the second data set and thus expand its data base. The second data set thus does not have to be static and can be expanded, either by external input or by incorporating new data due to self-learning. This is, however, not restricted to the third aspect of the invention, but applies to other aspects of the invention that refer to a second data set, which does not necessarily have to refer to antimicrobial drug resistance. The same applies, where applicable, to the first data set, e.g. in the third aspect.
According to certain embodiments, statistical analysis in the present methods is carried out using Fisher's test with p<10−6, preferably p<10−9, particularly p<10−10.
The method of the third aspect of the present invention, as well as related methods, e.g. according to the 7th, 10th, 11th, 14th and 16th aspect, can, according to certain embodiments, comprise correlating different genetic sites to each other, e.g. in at least two, three, four, five, six, seven, eight, nine or ten genes. This way even higher statistical significance can be achieved.
According to certain embodiments of the method of the third aspect and related methods—as above, the second data set is provided by culturing the clinical isolates of Shigella species on agar plates provided with antimicrobial drugs, e.g. antibiotics, at different concentrations and the second data is obtained by taking the minimal concentration of the plates that inhibits growth of the respective Shigella species.
According to certain embodiments of the method of the third aspect and related methods, the antibiotic is at least one selected from the group of β-lactams, β-lactam inhibitors, quinolines and derivatives thereof, aminoglycosides, tetracyclines, and folate synthesis inhibitors, preferably Amoxicillin/K Clavulanate, Ampicillin, Aztreonam, Cefazolin, Cefepime, Cefotaxime, Ceftazidime, Ceftriaxone, Cefuroxime, Cephalothin, Ciprofloxacin, Ertapenem, Gentamicin, Imipenem, Levofloxacin, Meropenem, Piperacillin/Tazobactam, Ampicillin/Sulbactam, Tetracycline, Tobramycin, and Trimethoprim/Sulfamethoxazole.
According to certain embodiments of the method of the third aspect and related methods, the gene sequences in the third data set are comprised in at least one gene from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, and pdxA, or from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and SSON53_06085, or from the genes listed in Table 5.
According to certain embodiments of the method of the third aspect and related methods, the genetic sites in the genome of Shigella associated with antimicrobial drug, e.g. antibiotic, resistance are at least comprised in one gene from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and SSON53_06085.
According to certain embodiments of the method of the third aspect and related methods, the genetic variant has a point mutation, an insertion and or deletion of up to four bases, and/or a frameshift mutation, particularly a frameshift mutation in YP_005456953.1.
A fourth aspect of the present invention relates to a method of determining an antimicrobial drug, e.g. antibiotic, resistance profile for a bacterial microorganism belonging to the species Shigella comprising the steps of
a) obtaining or providing a sample containing or suspected of containing the bacterial microorganism;
b) determining the presence of a mutation in at least one gene of the bacterial microorganism as determined by the method of the third aspect of the invention;
wherein the presence of a mutation is indicative of a resistance to an antimicrobial drug, e.g. antibiotic, drug.
Steps a) and b) can herein be carried out as described with regard to the first aspect, as well as for the following aspects of the invention.
With this method, any mutations in the genome of Shigella species correlated with antimicrobial drug, e.g. antibiotic, resistance can be determined and a thorough antimicrobial drug, e.g. antibiotic, resistance profile can be established.
A fifth aspect of the present invention relates to a diagnostic method of determining an infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment, which also can be described as method of determining an antimicrobial drug, e.g. antibiotic, resistant Shigella infection in a patient, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing a bacterial microorganism belonging to the species Shigella from the patient;
b) determining the presence of at least one mutation in at least one gene of the bacterial microorganism belonging to the species Shigella as determined by the method of the third aspect of the present invention, wherein the presence of said at least one mutation is indicative of an antimicrobial drug, e.g. antibiotic, resistant Shigella infection in said patient.
Again, steps a) and b) can herein be carried out as described with regard to the first aspect of the present invention.
According to this aspect, a Shigella infection in a patient can be determined using sequencing methods as well as a resistance to antimicrobial drugs, e.g. antibiotics, of the Shigella species be determined in a short amount of time compared to the conventional methods.
In a sixth aspect the present invention relates to a method of selecting a treatment of a patient suffering from an infection with a potentially resistant Shigella strain, e.g. an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing a bacterial microorganism belonging to the species Shigella from the patient;
b) determining the presence of at least one mutation in at least one gene of the bacterial microorganism belonging to the species Shigella as determined by the method of the third aspect of the invention, wherein the presence of said at least one mutation is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs; and
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection.
This method can be carried out similarly to the second aspect of the invention and enables a fast was to select a suitable treatment with antibiotics for any infection with an unknown Shigella species.
A seventh aspect of the present invention relates to a method of acquiring, respectively determining, an antimicrobial drug, e.g. antibiotic, resistance profile for a bacterial microorganisms of Shigella species, comprising:
obtaining or providing a first data set of gene sequences of a clinical isolate of Shigella species;
providing a second data set of antimicrobial drug, e.g. antibiotic, resistance of a plurality of clinical isolates of Shigella species;
aligning the gene sequences of the first data set to at least one, preferably one, reference genome of Shigella, and/or assembling the gene sequence of the first data set, at least in part;
analyzing the gene sequences of the first data set for genetic variants to obtain a third data set of genetic variants of the first data set;
correlating the third data set with the second data set and statistically analyzing the correlation; and
determining the genetic sites in the genome of Shigella of the first data set associated with antimicrobial drug, e.g. antibiotic, resistance.
With this method, antimicrobial drug, e.g. antibiotic, resistances in an unknown isolate of Shigella can be determined.
According to certain embodiments, the reference genome of Shigella is NC_016822 as annotated at the NCBI. According to certain embodiments, statistical analysis in the present methods is carried out using Fisher's test with p<10−6, preferably p<10−9, particularly p<10−10. Also, according to certain embodiments, the method further comprises correlating different genetic sites to each other, e.g. in at least two, three, four, five, six, seven, eight, nine or ten genes.
An eighth aspect of the present invention relates to a computer program product comprising computer executable instructions which, when executed, perform a method according to the third, fourth, fifth, sixth or seventh aspect of the present invention.
In certain embodiments the computer program product is one on which program commands or program codes of a computer program for executing said method are stored. According to certain embodiments the computer program product is a storage medium. The same applies to the computer program products of the aspects mentioned afterwards, i.e. the twelfth and fifteenth aspect of the present invention. As noted above, the computer program products of the present invention can be self-learning, e.g. with respect to the first and second data sets.
In order to obtain the best possible information from the highly complex genetic data and develop an optimum model for diagnostic and therapeutical uses as well as the methods of the present invention—which can be applied stably in clinical routine—a thorough in-silico analysis can be necessary. The proposed principle is based on a combination of different approaches, e.g. alignment with at least one, preferably more reference genomes and/or assembly of the genome and correlation of mutations found in every sample, e.g. from each patient, with all references and drugs, e.g. antibiotics, and search for mutations which occur in several drug and several strains.
Using the above steps a list of mutations as well of genes is generated. These can be stored in databases and statistical models can be derived from the databases. The statistical models can be based on at least one or more mutations at least one or more genes. Statistical models that can be trained can be combined from mutations and genes. Examples of algorithms that can produce such models are association Rules, Support Vector Machines, Decision Trees, Decision Forests, Discriminant-Analysis, Cluster-Methods, and many more.
The goal of the training is to allow a reproducible, standardized application during routine procedures.
For this, for example, a genome or parts of the genome of a microorganism can be sequenced from a patient to be diagnosed. Afterwards, core characteristics can be derived from the sequence data, which can be used to predict resistance. These are the points in the database used for the final model, i.e. at least one mutation or at least one gene, but also combinations of mutations, etc.
The corresponding characteristics can be used as input for the statistical model and thus enable a prognosis for new patients. Not only the information regarding all resistances of all microorganisms, e.g. of Shigella species, against all drugs, e.g. antibiotics, can be integrated in a computer decision support tool, but also corresponding directives (e.g. EUCAST) so that only treatment proposals are made that are in line with the directives.
A ninth aspect of the present invention relates to the use of the computer program product according to the eighth aspect for acquiring an antimicrobial drug, e.g. antibiotic, resistance profile for bacterial microorganisms of Shigella species or in a method of the third aspect of the invention.
According to a tenth aspect, the present invention relates to a method of acquiring, respectively determining, an antimicrobial drug, e.g. antibiotic, resistance profile for a microorganism, particularly a bacterial microorganism, comprising:
obtaining or providing a first data set of gene sequences of a clinical isolate of the microorganism;
providing a second data set of antimicrobial drug, e.g. antibiotic resistance, of a plurality of clinical isolates of the microorganism;
aligning the gene sequences of the first data set to at least one, preferably one, reference genome of the microorganism, and/or assembling the gene sequences of the first data set, at least in part;
analyzing the gene sequences of the first data set for genetic variants to obtain a third data set of genetic variants of the first data set;
correlating the third data set with the second data set of antimicrobial drug, e.g. antibiotic, resistance of a plurality of clinical isolates of the microorganism and statistically analyzing the correlation; and
determining the genetic sites in the genome of the clinical isolate of the microorganism of the first data set associated with antimicrobial drug, e.g. antibiotic, resistance.
Again, the steps can be carried out as described for corresponding steps in other aspects before. This method describes the general concept of acquiring an antimicrobial drug, e.g. antibiotic, resistance profile for any microorganism.
According to an eleventh aspect, a method of acquiring, respectively determining, an antimicrobial drug, e.g. antibiotic, resistance profile for at least one microorganism, particularly at least one bacterial microorganism, obtained from a patient, comprising:
acquiring a first data set of gene sequences of the at least one microorganism;
providing at least one second data set of antimicrobial drug, e.g. antibiotic, resistance, of a plurality of clinical isolates of the microorganism;
aligning the gene sequences of the first data set to at least one, preferably one, reference genome of the microorganism, and/or assembling the gene sequence of the first data set, at least in part;
analyzing the gene sequences of the first data set for genetic variants to obtain a third data set of genetic variants of the first data set;
correlating the third data set with the at least one second data set of antimicrobial drug, e.g. antibiotic, resistance of a plurality of clinical isolates of the microorganism and statistically analyzing the correlation; and
determining the genetic sites in the genome of the at least one microorganism of the first data set associated with antimicrobial drug, e.g. antibiotic, resistance is described.
Again, the steps defined can be carried out as similar steps described in other methods before. This method discloses the general identification of antimicrobial drug, e.g. antibiotic, resistances in any sample, particularly by a general computer program product.
A simple read out concept for a diagnostic test as described in this aspect is shown schematically in
According to
According to certain embodiments, statistical analysis in the present method is carried out using Fisher's test with p<10−6, preferably p<10−9, particularly p<10−10. Also, according to certain embodiments, the method further comprises correlating different genetic sites to each other.
A twelfth aspect of the present invention relates to a computer program product comprising computer executable instructions which, when executed, perform a method according to the 10th or 11th aspect of the present invention. These computer program products thus relate to general computer program products referring to any microorganisms, applying the principles which also apply for microorganisms of Shigella species as laid out before.
According to a thirteenth aspect the use of the computer program products of the thirteenth aspect in a method of the first, second, fourth, fifth or sixth aspect is disclosed.
In a fourteenth aspect a method of selecting a treatment of a patient having an infection with a microorganism, preferably a bacterial microorganism, particularly preferably a bacterial microorganism of Shigella species, comprising:
obtaining or providing a first data set comprising a gene sequence of at least one clinical isolate of the microorganism from the patient;
providing a second data set of antimicrobial drug, e.g. antibiotic, resistance of a plurality of clinical isolates of the microorganism;
aligning the gene sequences of the first data set to at least one, preferably one, reference genome of the microorganism, and/or assembling the gene sequence of the first data set, at least in part;
analyzing the gene sequences of the first data set for genetic variants to obtain a third data set of genetic variants of the first data set;
correlating the third data set with the second data set of antimicrobial drug, e.g. antibiotic, resistance of a plurality of clinical isolates of the microorganism and statistically analyzing the correlation;
determining the genetic sites in the genome of the clinical isolate of the microorganism of the first data set associated with antimicrobial drug, e.g. antibiotic, resistance; and
selecting a treatment of the patient with one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in the determination of the genetic sites associated with antimicrobial drug, e.g. antibiotic, resistance is disclosed.
Again, the steps can be carried out as similar steps before. In this method, as well as similar ones, no aligning is necessary, as the unknown sample can be directly correlated, after the genome or genome sequences are produced, with the second data set and thus mutations and antimicrobial drug, e.g. antibiotic, resistances can be determined. The first data set can be assembled, for example, using known techniques.
According to certain embodiments, statistical analysis in the present method is carried out using Fisher's test with p<10−6, preferably p<10−9, particularly p<10−10. Also, according to certain embodiments, the method further comprises correlating different genetic sites to each other.
A fifteenth aspect of the present invention is directed to a computer program product comprising computer executable instructions which, when executed, perform a method according to the fourteenth aspect.
In a sixteenth aspect the present invention relates to a method of finding a pharmaceutical target for treating an infection with at least one microorganism, particularly at least one bacterial microorganism, comprising:
obtaining or providing a first data set of gene sequences of a plurality of clinical isolates of the at least one microorganism;
providing a second data set of sequences of the at least one microorganism representing genes relevant to a pharmaceutical compound;
aligning the gene sequences of the first data set to at least one, preferably one, reference genome each of the at least one microorganism, and/or assembling the gene sequence of the first data set, at least in part;
analyzing the gene sequences of the first data set for genetic variants to obtain a third data set of genetic variants of the first data set;
correlating the third data set to the second data set of sequences of the at least one microorganism representing genes relevant to a pharmaceutical compound; and
determining the genetic sites in the genome of the clinical isolate of the microorganism of the first data set associated relevant to a pharmaceutical compound.
According to certain embodiments the sequences of the microorganism representing genes relevant to a pharmaceutical compound represent binding sites of the microorganism. These are likely to be targets for drugs, e.g. antibiotics.
Again, the principle of the method can be conducted using a computer program product.
Algorithms that can be used in such a computer program product showed for cases of single base exchanges (e.g. gyrase) as well as deletion of whole genes or open reading frames (e.g. beta lactamase) that known targets for drugs applied successfully on the marked can be found. It is to be expected that the data is perfectly suitable for finding further drug targets/pharmaceutical targets (many different organisms that are not monoclonal but very diverse in temporal and geographical aspects). For this purpose it can be tested for all highly significant mutations if the mutation is in a gene relevant for drugs, e.g. the binding site of a medicine). The results are of diagnostic interest. The findings can be prioritized. For example, for single base exchanges preferably results will be considered which are in binding sites. The respective binding sites most likely have a function in the metabolism of the drug, so that a mutation will cause an increased or decreased or lost effectiveness. Corresponding genes or even binding partners of the genes present presumed new targets for drugs, which can enhance effects of existing therapies or work as single therapies.
In order to limit the number of candidates, all proteins with unknown function or hypothetical proteins can be excluded, concentrating on genes that are found in several organisms, e.g. more than ten, independently of one another.
According to a seventeenth aspect of the present invention, a diagnostic method of determining an infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment, which can also be described as a method of determining an antimicrobial drug, e.g. antibiotic, resistant Shigella infection of a patient is disclosed, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least two genes from the group of genes listed in Table 5, wherein the presence of said at least two mutations is indicative of an antimicrobial drug, e.g. antibiotic, resistant Shigella infection in said patient.
An eighteenth aspect of the invention discloses a method of selecting a treatment of a patient suffering from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least two genes from the group of genes listed in Table 5, wherein the presence of said at least two mutations is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs; and
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection.
Again, the steps can be carried out as in similar methods before, e.g. as in the first and second aspect of the invention. In the seventeenth and eighteenth aspect of the invention, as well as also in the twenty-third aspect of the invention, all classes of antibiotics considered in the present method are covered.
Herein, the genes in Table 5 are the following: SSON53_16770, SSON53_19730, SSON53_14680, SSON53_09620, SSON53_25600, metH, SSON53_24790, SSON53_13530, SSON53_12100, SSON53_12655, SSON53_04070, SSON53_21285, SSON53_01970, held, SSON53_08660, SSON53_14430, SSON53_17515, SSON53_11020, SSON53_00015, SSON53_01080, SSON53_12830, SSON53_13420, SSON53_13780, SSON53_16000, SSON53_16880, SSON53_18540, SSON53_22785, SSON53_24145, SSON53_25200, SSON53_25460, SSON53_03770, SSON53_18890, SSON53_16535, SSON53_19285, SSON53_03650, SSON53_15285, SSON53_18775, SSON53_17120, SSON53_24410, SSON53_26285, tig, SSON53_07220, SSON53_08070, SSON53_17780, gcvT, SSON53_20935, pyrB, SSON53_12595, SSON53_04605, SSON53_07910, SSON53_13390, queF, SSON53_18505, pdxA, SSON53_00645, SSON53_01515, SSON53_02505, SSON53_03060, dpiB, SSON53_03230, SSON53_03255, SSON53_04495, SSON53_05165, SSON53_05595, SSON53_06535, SSON53_06585, SSON53_07385, SSON53_07390, SSON53_08640, SSON53_08945, SSON53_10840, SSON53_12340, pbpG, SSON53_13270, SSON53_14190, SSON53_14265, SSON53_14460, murP, SSON53_14705, SSON53_14835, SSON53_15125, SSON53_16590, SSON53_16955, rumA, recD, SSON53_17670, SSON53_17690, SSON53_17705, SSON53_17955, SSON53_18080, SSON53_18610, SSON53_18735, SSON53_18885, SSON53_00360, thiP, SSON53_00595, SSON53_02980, SSON53_03520, SSON53_03560, bioD, SSON53_04325, SSON53_04710, putA, SSON53_05585, rne, SSON53_07205, SSON53_08060, SSON53_08955, SSON53_08960, SSON53_09205, SSON53_10160, SSON53_11100, dacD, SSON53_13065, fadJ, SSON53_14085, SSON53_14455, murQ, SSON53_16065, SSON53_16755, SSON53_18520, SSON53_20320, SSON53_21970, fadB, SSON53_24280, SSON53_25150, SSON53_25295, SSON53_25400, SSON53_00025, SSON53_04285, SSON53_04300, SSON53_06415, SSON53_07990, SSON53_08770, SSON53_09680, SSON53_10655, SSON53_11810, SSON53_12330, SSON53_12740, SSON53_12795, SSON53_13290, SSON53_13415, SSON53_13555, SSON53_15280, SSON53_16180, SSON53_19150, SSON53_22470, SSON53_24775, SSON53_24780, SSON53_02695, SSON53_04540, wcaD, SSON53_14020, SSON53_14025, SSON53_21860, and SSON53_22935.
According to certain embodiments, mutations in at least two, three, four, five, six, seven, eight, nine or ten genes are determined in any of the methods of the present invention, e.g. in at least two genes or in at least three genes. Instead of testing only single genes or mutants, a combination of several variant positions can improve the prediction accuracy and further reduce false positive findings that are influenced by other factors. Therefore, it is in particular preferred to determine the presence of a mutation in 2, 3, 4, 5, 6, 7, 8 or 9 (or more) genes selected from Table 5.
Further, according to certain embodiments, the reference genome of Shigella is again NC_016822 as annotated at the NCBI. According to certain embodiments, statistical analysis in the present methods is carried out using Fisher's test with p<10−6, preferably p<10−9, particularly p<10−10. Also, according to certain embodiments, the method further comprises correlating different genetic sites to each other. Also the other aspects of the embodiments of the first and second aspect of the invention apply.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antimicrobial drug is an antibiotic. According to certain embodiments, the antibiotic is a lactam antibiotic and a mutation in at least one of the genes listed in Table 6 or 7 is detected, or a mutation in at least one of the positions (denoted POS in the tables) listed in Table 6 or 7.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is at least one of CF, CFZ, CRM, P/T, AM, A/S, and AUG and a mutation in at least one of the genes of SSON53_16770, SSON53_19730, SSON53_14680, SSON53_09620, SSON53_25600 is detected, or a mutation in at least one of the positions of 3134050, 3722816, 2748494, 1816305, 4860933.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is at least one of CF, A/S, CRM, AM, and AUG and a mutation in at least one of the genes of metH, SSON53_24790, SSON53_13530, SSON53_12100, SSON53_12655, SSON53_04070, SSON53_21285, SSON53_01970, held, SSON53_08660, SSON53_14430, SSON53_17515, SSON53_11020, SSON53_00015, SSON53_01080 is detected, or a mutation in at least one of the positions of 4595667, 4707429, 2533894, 2239494, 2352355, 802768, 4005875, 399928, 1045235, 1633619, 2697830, 3283400, 2053683, 3915, 226692.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is P/T_and a mutation in at least one of the genes of metH, SSON53_24790, SSON53_13530, SSON53_12100, SSON53_12655, SSON53_04070, SSON53_21285 is detected, or a mutation in at least one of the positions of 4595667, 4707429, 2533894, 2239494, 2352355, 802768, 4005875.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is CFZ and a mutation in at least one of the genes of SSON53_01970, held, SSON53_08660, SSON53_14430, SSON53_17515, SSON53_11020, SSON53_00015, SSON53_01080 is detected, or a mutation in at least one of the positions of 399928, 1045235, 1633619, 2697830, 3283400, 2053683, 3915, 226692.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is a quinolone antibiotic and a mutation in at least one of the genes listed in Table 8 or 9 is detected, or a mutation in at least one of the positions (denoted POS in the tables) listed in Table 8 or 9.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is at least one of CP and LVX and a mutation in at least one of the genes of SSON53_04605, SSON53_07910, SSON53_13390, queF, SSON53_18505 is detected, or a mutation in at least one of the positions of 905340, 1506051, 2503366, 3228225, 3476068.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is CP and a mutation in at least one of the genes of SSON53_13390, SSON53_19730, SSON53_16955, SSON53_05595, SSON53_07390, SSON53_15125, dpiB, SSON53_07385, SSON53_17670, SSON53_14460, SSON53_21860 is detected, or a mutation in at least one of the positions of 2503749, 3723920, 3723919, 3723923, 3723925, 3164176, 1110022, 1422899, 2846241, 619830, 1421710, 3315886, 2704526, 3722696, 4105284.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is an aminoglycoside antibiotic and a mutation in at least one of the genes listed in Table 10 is detected, or a mutation in at least one of the positions (denoted POS in the tables) listed in Table 10.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is TO and a mutation in at least one of the genes of SSON53_04325, SSON53_04710, SSON53_14455, SSON53_16755, SSON53_03520, SSON53_16065, SSON53_00360, SSON53_00595, SSON53_25150, fadB, SSON53_10160, SSON53_02980, SSON53_02980, SSON53_08960, SSON53_03560 is detected, or a mutation in at least one of the positions of 847461, 927380, 2702371, 3130278, 707965, 3013091, 74735, 124368, 4777150, 4366586, 1917725, 600143, 600209, 1693067, 714500.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is an polyketide antibiotic and a mutation in at least one of the genes listed in Table 11 is detected, or a mutation in at least one of the positions (denoted POS in the tables) listed in Table 11.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is TE and a mutation in at least one of the genes of metH, SSON53_24790, SSON53_08770, SSON53_24780, SSON53_15280, SSON53_10655, SSON53_13555, SSON53_06415, SSON53_04285, SSON53_04300, SSON53_13530, SSON53_19150, SSON53_24775, SSON53_12100, SSON53_12795 is detected, or a mutation in at least one of the positions of 4595667, 4707429, 1655745, 1655894, 4705326, 2876708, 1990704, 2539010, 1257335, 840523, 842535, 2533894, 2539007, 3604284, 4703410, 2239494, 2378416.
According to certain embodiments of the method of the seventeenth and/or eighteenth aspect of the present invention, as well as also of the twenty-third aspect of the present invention, the antibiotic is T/S and a mutation in at least one of the genes listed in Table 12 is detected, or a mutation in at least one of the positions (denoted POS in the tables) listed in Table 12.
A nineteenth aspect of the present invention is directed to a diagnostic method of determining an infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment, which can also be described as method of determining an antimicrobial drug, e.g. antibiotic, resistant Shigella infection of a patient, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least one gene from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, and pdxA, or from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12475, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and SSON53_06085, wherein the presence of said at least one mutation is indicative of an antimicrobial drug, e.g. antibiotic, resistant Shigella infection in said patient.
A twentieth aspect of the present invention is directed to a method of selecting a treatment of a patient suffering from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least one gene from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, and pdxA, or from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12475, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and SSON53_06085, wherein the presence of said at least one mutation is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs; and
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection.
Again, in the nineteenth and the twentieth aspect the steps correspond to those in the first or second aspect, although only a mutation in at least one gene is determined.
A twenty-first aspect of the present invention is directed to a method of treating a patient suffering from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least one gene from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, and pdxA, or from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12475, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and SSON53_06085, wherein the presence of said at least one mutation is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs;
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection; and
e) treating the patient with said one or more antimicrobial, e.g. antibiotic, drugs.
A twenty-second aspect of the present invention is directed to method of treating a patient suffering from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least two genes from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12070, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_06660, SSON53_24790, SSON53_16770, SSON53_24480, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, rimO, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, and pdxA, or from the group of genes consisting of SSON53_12105, SSON53_12455, SSON53_12475, metH, SSON53_00425, SSON53_23830, SSON53_13575, flgE, SSON53_12500, SSON53_24790, SSON53_16770, SSON53_08770, SSON53_04795, SSON53_07430, SSON53_14085, SSON53_07440, SSON53_08950, SSON53_25555, SSON53_09145, SSON53_25565, SSON53_24780, SSON53_04400, SSON53_15280, SSON53_04930, SSON53_23390, rhaB, SSON53_17960, thiH, SSON53_03725, SSON53_09500, SSON53_25405, SSON53_13530, astD, SSON53_07945, SSON53_10080, SSON53_10655, SSON53_13400, SSON53_13555, SSON53_14945, SSON53_15285, fucI, SSON53_02755, SSON53_06415, SSON53_04615, SSON53_17330, pdxA, SSON53_04690, SSON53_04945, SSON53_06005, and SSON53_06085, wherein the presence of said at least two mutations is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs;
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection; and
e) treating the patient with said one or more antimicrobial, e.g. antibiotic, drugs.
A twenty-third aspect of the present invention is directed to a method of treating a patient suffering from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least two genes from the group of genes listed in Table 5, wherein the presence of said at least two mutations is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs;
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection; and
e) treating the patient with said one or more antimicrobial, e.g. antibiotic, drugs.
A twenty-fourth aspect of the present invention is directed to a method of treating a patient suffering from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least one gene from the group of genes listed in Table 13, preferably from the group of genes listed in Table 14, wherein the presence of said at least one mutation is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs;
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection; and
e) treating the patient with said one or more antimicrobial, e.g. antibiotic, drugs.
Also in the twenty-first to twenty-fourth aspect of the invention, steps a) to d) are analogous to the steps in the method of the second aspect of the present invention. Step e) can be sufficiently carried out without being restricted and can be done e.g. non-invasively.
A twenty-fifth aspect of the present invention is directed to a diagnostic method of determining an infection of a patient with Shigella species potentially resistant to antimicrobial drug treatment, which can also be described as method of determining an antimicrobial drug, e.g. antibiotic, resistant Shigella infection of a patient, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least one gene from the group of genes listed in Table 13, preferably from the group of genes listed in Table 14, wherein the presence of said at least one mutation is indicative of an antimicrobial drug, e.g. antibiotic, resistant Shigella infection in said patient.
A twenty-sixth aspect of the present invention is directed to a method of selecting a treatment of a patient suffering from an antimicrobial drug, e.g. antibiotic, resistant Shigella infection, comprising the steps of:
a) obtaining or providing a sample containing or suspected of containing at least one Shigella species from the patient;
b) determining the presence of at least one mutation in at least one gene from the group of genes listed in Table 13, preferably from the group of genes listed in Table 14, wherein the presence of said at least one mutation is indicative of a resistance to one or more antimicrobial, e.g. antibiotic, drugs;
c) identifying said at least one or more antimicrobial, e.g. antibiotic, drugs; and
d) selecting one or more antimicrobial, e.g. antibiotic, drugs different from the ones identified in step c) and being suitable for the treatment of a Shigella infection.
Again, in the twenty-fifth and the twenty-sixth aspect the steps correspond to those in the first or second aspect, although only a mutation in at least one gene is determined.
The present invention will now be described in detail with reference to several examples thereof. However, these examples are illustrative and do not limit the scope of the invention.
Whole genome sequencing was carried out in addition to classical antimicrobial susceptibility testing of the same isolates for a cohort of 470 specimens. This allowed performing genome wide correlation studies to find genetic variants (e.g. point mutations, small insertions and deletion, larger structural variants, plasmid copy number gains, gene dosage effects) in the genome and plasmids that are significantly correlated to the resistance against one or several drugs. The approach also allows for comparing the relevant sites in the genome to each other.
In the approach the different sources of genetic resistance as well as the different ways of how bacteria can become resistant were covered. By measuring clinical isolates collected in a broad geographical area and across a broad time span of three decades a complete picture going far beyond the rather artificial step of laboratory generated resistance mechanisms was tried to be generated.
To this end, a set of 21 clinically relevant antimicrobial agents with 5 different modes of action was put together, and the minimally inhibitory concentration (MIC) of the 21 drugs for the Shigella isolates was measured.
The detailed procedure is given in the following:
Bacterial Strains
The inventors selected 470 Shigella strains from the microbiology strain collection at Siemens Healthcare Diagnostics (West Sacramento, Calif.) for susceptibility testing and whole genome sequencing.
Antimicrobial Susceptibility Testing (AST) Panels Frozen reference AST panels were prepared following Clinical Laboratory Standards Institute (CLSI) recommendations. The following antimicrobial agents (with μg/ml concentrations shown in parentheses) were included in the panels: Amoxicillin/K Clavulanate (0.5/0.25-64/32), Ampicillin (0.25-128), Ampicillin/Sulbactam (0.5/0.25-64/32), Aztreonam (0.25-64), Cefazolin (0.5-32), Cefepime (0.25-64), Cefotaxime (0.25-128), Ceftazidime (0.25-64), Ceftriaxone (0.25-128), Cefuroxime (1-64), Cephalothin (1-64), Ciprofloxacin (0.015-8), Ertepenem (0.12-32), Gentamicin (0.12-32), Imipenem (0.25-32), Levofloxacin (0.25-16), Meropenem (0.12-32), Piperacillin/Tazobactam (0.25/4-256/4), Tetracycline (0.5-64), Tobramycin (0.12-32), and Trimethoprim/Sulfamethoxazole (0.25/4.7-32/608). Prior to use with clinical isolates, AST panels were tested with QC strains. AST panels were considered acceptable for testing with clinical isolates when the QC results met QC ranges described by CLSI16.
Inoculum Preparation
Isolates were cultured on trypticase soy agar with 5% sheep blood (BBL, Cockeysville, Md.) and incubated in ambient air at 35±1° C. for 18-24 h. Isolated colonies (4-5 large colonies or 5-10 small colonies) were transferred to a 3 ml Sterile Inoculum Water (Siemens) and emulsified to a final turbidity of a 0.5 McFarland standard. 2 ml of this suspension was added to 25 ml Inoculum Water with Pluronic-F (Siemens). Using the Inoculator (Siemens) specific for frozen AST panels, 5 μl of the cell suspension was transferred to each well of the AST panel. The inoculated AST panels were incubated in ambient air at 35±1° C. for 16-20 h. Panel results were read visually, and minimal inhibitory concentrations (MIC) were determined.
DNA Extraction
Four streaks of each Gram-negative bacterial isolate cultured on trypticase soy agar containing 5% sheep blood and cell suspensions were made in sterile 1.5 ml collection tubes containing 50 μl Nuclease-Free Water (AM9930, Life Technologies). Bacterial isolate samples were stored at −20° C. until nucleic acid extraction. The Tissue Preparation System (TPS) (096D0382-02_01_B, Siemens) and the VERSANT® Tissue Preparation Reagents (TPR) kit (10632404B, Siemens) were used to extract DNA from these bacterial isolates. Prior to extraction, the bacterial isolates were thawed at room temperature and were pelleted at 2000 G for 5 seconds. The DNA extraction protocol DNAext was used for complete total nucleic acid extraction of 48 isolate samples and eluates, 50 μl each, in 4 hours. The total nucleic acid eluates were then transferred into 96-Well qPCR Detection Plates (401341, Agilent Technologies) for RNase A digestion, DNA quantitation, and plate DNA concentration standardization processes. RNase A (AM2271, Life Technologies) which was diluted in nuclease-free water following manufacturer's instructions was added to 50 μl of the total nucleic acid eluate for a final working concentration of 20 μg/ml. Digestion enzyme and eluate mixture were incubated at 37° C. for 30 minutes using Siemens VERSANT® Amplification and Detection instrument. DNA from the RNase digested eluate was quantitated using the Quant-iT™ PicoGreen dsDNA Assay (P11496, Life Technologies) following the assay kit instruction, and fluorescence was determined on the Siemens VERSANT® Amplification and Detection instrument. Data analysis was performed using Microsoft® Excel 2007. 25 μl of the quantitated DNA eluates were transferred into a new 96-Well PCR plate for plate DNA concentration standardization prior to library preparation. Elution buffer from the TPR kit was used to adjust DNA concentration. The standardized DNA eluate plate was then stored at −80° C. until library preparation.
Next Generation Sequencing
Prior to library preparation, quality control of isolated bacterial DNA was conducted using a Qubit 2.0 Fluorometer (Qubit dsDNA BR Assay Kit, Life Technologies) and an Agilent 2200 TapeStation (Genomic DNA ScreenTape, Agilent Technologies). NGS libraries were prepared in 96 well format using NexteraXT DNA Sample Preparation Kit and NexteraXT Index Kit for 96 Indexes (Illumina) according to the manufacturer's protocol. The resulting sequencing libraries were quantified in a qPCR-based approach using the KAPA SYBR FAST qPCR MasterMix Kit (Peqlab) on a ViiA 7 real time PCR system (Life Technologies). 96 samples were pooled per lane for paired-end sequencing (2×100 bp) on Illumina Hiseq2000 or Hiseq2500 sequencers using TruSeq PE Cluster v3 and TruSeq SBS v3 sequencing chemistry (Illumina). Basic sequencing quality parameters were determined using the FastQC quality control tool for high throughput sequence data (Babraham Bioinformatics Institute).
Data Analysis
Raw paired-end sequencing data for the 470 Shigella samples were mapped against the Shigella reference (NC_016822) with BWA 0.6.1.20. The resulting SAM files were sorted, converted to BAM files, and PCR duplicates were marked using the Picard tools package 1.104 (http://picard.sourceforge.net/). The Genome Analysis Toolkit 3.1.1 (GATK) 21 was used to call SNPs and indels for blocks of 200 Shigella samples (parameters: -ploidy 1 -glm BOTH -stand_call_conf 30 -stand_emit_conf 10). VCF files were combined into a single file and quality filtering for SNPs was carried out (QD<2.0∥FS>60.0∥MQ<40.0) and indels (QD<2.0∥FS>200.0). Detected variants were annotated with SnpEff22 to predict coding effects. For each annotated position, genotypes of all Shigella samples were considered. Shigella samples were split into two groups, low resistance group (having lower MIC concentration for the considered drug), and high resistance group (having higher MIC concentrations) with respect to a certain MIC concentration (breakpoint). To find the best breakpoint all thresholds were evaluated and p-values were computed with Fisher's exact test relying on a 2×2 contingency table (number of Shigella samples having the reference or variant genotype vs. number of samples belonging to the low and high resistance group). The best computed breakpoint was the threshold yielding the lowest p-value for a certain genomic position and drug. For further analyses positions with non-synonymous alterations and p-value <10−9 were considered.
Since a potential reason for drug resistance is gene duplication, gene dose dependency was evaluated. For each sample the genomic coverage for each position was determined using BED Tools. 23 Gene ranges were extracted from the reference assembly NC_016822.gff and the normalized median coverage per gene was calculated. To compare low- and high-resistance isolates the best area under the curve (AUC) value was computed. Groups of at least 20% of all samples having a median coverage larger than zero for that gene and containing more than 15 samples per group were considered in order to exclude artifacts and cases with AUC >0.75 were further evaluated.
To include data on the different ways how resistance mechanisms are acquired Shigella isolates collected over more than three decades were analyzed such that also horizontal gene transfer could potentially be discovered.
In detail, the following steps were carried out: Shigella strains to be tested were seeded on agar plates and incubated under growth conditions for 24 hours. Then, colonies were picked and incubated in growth medium in the presence of a given antibiotic drug in dilution series under growth conditions for 16-20 hours. Bacterial growth was determined by observing turbidity.
Next mutations were searched that are highly correlated with the results of the phenotypic resistance test.
For sequencing, samples were prepared using a Nextera library preparation, followed by multiplexed sequencing using the Illuminat HiSeq 2500 system, paired end sequencing. Data were mapped with BWA (Li H. and Durbin R. (2010) Fast and accurate long-read alignment with Burrows-Wheeler Transform. Bioinformatics, Epub. [PMID: 20080505]) and SNP were called using samtools (Li H.*, Handsaker B.*, Wysoker A., Fennell T., Ruan J., Homer N., Marth G., Abecasis G., Durbin R. and 1000 Genome Project Data Processing Subgroup (2009) The Sequence alignment/map (SAM) format and SAMtools. Bioinformatics, 25, 2078-9. [PMID: 19505943]).
As reference genome, NC_016822 as annotated at the NCBI was determined as best suited.
The mutations were matched to the genes and the amino acid changes were calculated. Using different algorithms (SVM, homology modeling) mutations leading to amino acid changes with likely pathogenicity/resistance were calculated.
In total, whole genomes and plasmids of 470 different clinical isolates of Shigella species, particularly Shigella boydii, Shigella dysenteriae, Shigella flexneri, Shigella sonnei and other Shigella species, were sequenced, and classical antimicrobial susceptibility testing (AST) against 21 therapy forms as described above was performed for all organisms. From the classical AST a table with 470 rows (isolates) and 21 columns (MIC values for 21 drugs) was obtained. Each table entry contained the MIC for the respective isolate and the respective drug. The genetic data were mapped to different reference genomes of Shigella that have been annotated at the NCBI (http://www.ncbi.nlm.nih.gov/), and the best reference was chosen as template for the alignment—NC_016822 as annotated at the NCBI, found as SEQ ID NO 1 in the sequence listing. Additionally, assemblies were carried out and it was verified that the sequenced genomes fulfill all quality criteria to become reference genomes.
Next, genetic variants were evaluated. This approach resulted in a table with the genetic sites in columns and the same isolates in 470 rows. Each table entry contained the genetic determinant at the respective site (A, C, T, G, small insertions and deletions, . . . ) for the respective isolate.
In a next step different statistical tests were carried out
From the data, first the 50 genes with the best p-value were chosen for the list of mutations as well as the list of correlated antibiotic resistance, representing Tables 1 and 2. As can be seen from Tables 1 and 2, slight differences between the tables can be observed, showing the necessity to carry out both steps for determining statistical significant data for antimicrobial drug, e.g. antibiotic, resistance profiles.
A full list of all genetic sites, drugs, drug classes, affected genes etc. is provided in Tables 3 and 4a, 4b and 4c, wherein Table 3 corresponds to Table 1 and represents the genes having the lowest p-values after determining mutations in the genes, and Table 4, respectively Tables 4a, 4b and 4c correspond to Table 2 and represent the genes having the lowest p-values after correlating the mutations with antibiotic resistance for the respective antibiotics.
In addition, the data with the best p-values for each antibiotic class with the most antibiotic drugs as well as each antibiotic, respectively, were evaluated, being disclosed in Tables 5-12.
In Tables 3-12 the columns are designated as follows:
Gene name: affected gene;
POS: genomic position of the SNP/variant in the Shigella reference genome (see above);
P-value: significance value calculated using Fishers exact test;
genbank protein accession number: (NCBI) Accession number of the corresponding protein of the genes
Also the antibiotic/drug classes, the number of significant antibiotics correlated to the mutations (over all antibiotics or over certain classes), as well as the correlated antibiotics are denoted in the Tables.
The p-value was calculated using the Fisher exact test based on contingency table with 4 fields: #samples Resistant/wild type; #samples Resistant/mutant; #samples not Resistant/wild type; #samples not Resistant/mutant
The test is based on the distribution of the samples in the 4 fields. Even distribution indicates no significance, while clustering into two fields indicates significance.
The following results were obtained
Amoxicillin/Clavulanate, Ampicillin, Ampicillin/Sulbactam, Aztreonam, Cefazolin, Cefepime, Ceftazidime, Cefuroxime, Cephalothin, Imipenem, Piperacillin/Tazobactam, Ciprofloxacin, Levofloxacin, Gentamycin, Tobramycin, Tetracycline, Trimethoprim/Sulfamethoxazol
While in the tables only the best mutations in each gene are represented, a manifold of different SNPs has been found for each gene. Examples for multiple SNPs for two of the genes given in Table 3 are shown in the following Tables 15 and 16.
Similar results were obtained for other genes but are omitted for the sake of brevity.
Further, a synergistic effect of individual SNPs was demonstrated by exhaustively comparing significance levels for association of single SNPs with antibiotic susceptibility/resistance and significance levels for association of combinations of SNPs with antibiotic susceptibility/resistance. For a representative example of 2 SNPs
the significance level for synergistic association of two SNPs was at least 158% of the association of either SNP alone, corresponding to an increase of 58% to the original diagnostic information content. Again, similar results were obtained for other SNPs in respective genes.
Interestingly, it was also observed that the synergistic effect is enhanced for a combination of SNPs in different genes compared to SNPs from the same gene. Specifically, the above shown example reaching the increase to 158% of the original performance was obtained by combining two mutations from two different genes.
A genetic test for the combined pathogen identification and antimicrobial susceptibility testing direct from the patient sample can reduce the time-to actionable result significantly from several days to hours, thereby enabling targeted treatment. Furthermore, this approach will not be restricted to central labs, but point of care devices can be developed that allow for respective tests. Such technology along with the present methods and computer program products could revolutionize the care, e.g. in intense care units or for admissions to hospitals in general. Furthermore, even applications like real time outbreak monitoring can be achieved using the present methods.
Instead of using only single variants, a combination of several variant positions can improve the prediction accuracy and further reduce false positive findings that are influenced by other factors.
Compared to approaches using MALDI-TOF MS, the present approach has the advantage that it covers almost the complete genome and thus enables us to identify the potential genomic sites that might be related to resistance. While MALDI-TOF MS can also be used to identify point mutations in bacterial proteins, this technology only detects a subset of proteins and of these not all are equally well covered. In addition, the identification and differentiation of certain related strains is not always feasible.
The present method allows computing a best breakpoint for the separation of isolates into resistant and susceptible groups. The inventors designed a flexible software tool that allows to consider—besides the best breakpoints—also values defined by different guidelines (e.g. European and US guidelines), preparing for an application of the GAST in different countries.
The inventors demonstrate that the present approach is capable of identifying mutations in genes that are already known as drug targets, as well as detecting potential new target sites.
The current approach enables
| Number | Date | Country | Kind |
|---|---|---|---|
| PCT/EP2015/062202 | Jun 2015 | EP | regional |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2016/062379 | 6/1/2016 | WO | 00 |
