Patent Application
20230349002
This application incorporates by reference the Sequence Listing contained in the following ASCII text file being submitted concurrently herewith:
File name: 4906-00100 1811067PPUS Sequence Listing; created on Feb. 15, 2023; and having a files size of 7 KB.
The information in the Sequence Listing is incorporated herein in its entirety for all purposes.
The invention to which this application relates is a new diagnostic methodology and primers and/or drug susceptibility testing (DST) assay. In particular, the present invention relates to novel primers, and their use in a method of identifying and/or detecting the presence of drug resistance mutations in a sample from subjects with suspected or confirmed Tuberculosis.
Tuberculosis (TB), caused primarily by Mycobacterium tuberculosis1,2, is a disease of global health importance3-5. Mycobacterium tuberculosis and related bacteria in the Mycobacterium tuberculosis complex (MTBc) emerged at least 11,000 years ago and have been coevolving with their hosts since6,7. This history has resulted in a highly transmissible taxon of bacteria with longevity within their host and advanced methods of immune system evasion7.
Due to this coevolution, modern M. tuberculosis and members of the MTBc share numerous characteristics and are found in every known environment (except in the polar regions) along with members of the Non-Tuberculous mycobacterium (NTM) group7,8. The MTBc is made up of 10 mycobacterium capable of causing TB or TB-like disease within their hosts, with the three specialized human TB species being Mycobacterium tuberculosis sensu stricto, Mycobacterium canettis and Mycobacterium africanum1,7,9. Additionally, zoonotic TB transfer is well documented from cattle (Mycobacterium bovis), goats and sheep (Mycobacterium caprae), seals and sea lions (Mycobacterium pinnipedii), and rodents (Mycobacterium microti) into humans and vice versa4,6,7. Recently, three more species have been added; Mycobacterium orygis in cattle and antelope7,11, Mycobacterium suricattae in meerkats7,11, and Mycobacterium mungi in mongeese7,12.
Current research demonstrates MTBc members are highly genetically homogenous with up to 99.7% nucleotide identity and having identical 16S sequences7. MTBc members are primarily clonal with little horizontal gene transfer making differentiation between species difficult at the genetic level and impossible using microscopic methods2,4,6,13.
Mycobacteria are gram-positive acid-fast bacilli approximately 2 µm long, which are primarily transmitted via aerosols; they are strictly intracellular, and do not have a known environmental reservoir outside of their endemic hosts1,7,14. Lipid-rich cellular walls and layers of peptidoglycan, lipoglycan, mycolic acids, and waxes create an extremely hardy microbe7,14. A defining characteristic of many mycobacteria, and all members of the MTBc, is fastidiousness and slow rate of growth in culture and in vivo2,6,13,16.
Tuberculosis most commonly presents as a pulmonary disease (around 80% of cases), although extrapulmonary and disseminated disease presentations do also occur1,2,17. Mycobacterial diseases cause a high burden of disease in low- and middle-income and developing countries (LMICs) around the world3,6,18. It is estimated that one-third of the human population harbour latent TB (LTBI) and there are between nine and eleven million incident TB cases annually, according to the World Health Organization (WHO)19. The number of annual fatalities attributed to TB has been estimated at 1.5-2 million deaths globally, making TB the greatest single threat for infection associated mortality6,20,21.
The WHO defines drug resistance as a microorganism’s resistance to an antimicrobial drug that was once able to treat an infection by that microorganism. The emergence of drug resistant (DR) strains of TB is largely a result of inconsistent practice of treatment protocols, delayed treatment and/or patients defaulting on lengthy treatment courses, leading to positive selection for drug-resistance and a higher incidence of resistant strain transfer between hosts3,22,23.
There are currently several types of drug-resistant TB: multidrug-resistant (MDR) which is resistant to at least rifampicin and isoniazid; extensively drug-resistant (XDR) which has added resistances to any fluoroquinolone and at least one second-line injectable medication beyond what is found in MDR; extremely drug-resistant (XXDR) which is resistant to all first- and second-line medications; and totally drug-resistant (TDR) which has resistance to all current TB medications16,24. Additionally, some species within the MTBc have lineage specific inherent resistances, e.g. M. bovis and M. canettii, which if misdiagnosed can complicate resistance-control methods2,22,24.
Drug-resistant TB (DR-TB) is a growing issue globally as it increases in incidence21,22,25. Concerns are that drug-resistant strains will reverse the progress made towards TB eradication6,22,23. The incidence of drug resistant-TB worldwide has increased at least 10-fold in the past decade, with only 4.9% of patients demonstrating drug resistance in 2009 compared to 51% in 201819. In 2018 nearly 500,000 of approximately 10.5 million TB cases in the world were MDR and of those 31,000 (6.2%) were XDR19.
MDR-TB is the most common type of resistance16,24. MDR is defined as a TB strain which is resistant to isoniazid and rifampicin25. MDR-TB strains are typically treated with traditional WHO endorsed drug regimens which require a 6-month course of first- and second-line antibiotics. XDR-TB is an MDR strain with additional resistance to the second-line medications of any fluoroquinolones and amikacin, capreomycin, or kanamycin25,26. The specific regimen chosen to treat XDR-TB can be guided by culture or molecular (e.g. GenoType MTBDRsl - Bruker) drug susceptibility testing (DST) assays6,26,27where available. Due to difficulties in diagnosing and treating MDR and XDR strains of TB, the mortality rates in these cases are high with approximately 50% mortality MDR and over 70% in XDR-TB infections 25.
The first line treatment for TB is a combination of antibiotics; rifampicin, isoniazid, ethambutol, and pyrazinamide over 6 months. Resistance to these antibiotic therapies leads to the use of second-line antibiotics (fluoroquinolones, amikacin, capreomycin, and kanamycin), which are less effective and more toxic24,25. These therapeutics often require injections which necessitate more advanced medical infrastructure and oversight for treatment24.
Drug resistance in Mycobacteria is mutational, rather than transferrable, and numerous single nucleotide polymorphisms (SNPs) have been reported to be associated with drug-resistance over the past decades - however, not all have sufficient evidence in the literature to support this association. The World Health Organisation (WHO) and others have graded reported drug-resistance SNPs into high, moderate and low confidence brackets 28,29
The WHO has announced a goal to effectively eradicate TB by 2035 and released guidelines on how to achieve that goal in 2015 22,23,25,30. Central to the WHO defined eradication strategy was a call for new diagnostic technologies and more rapid drug-susceptibility testing (DST) capabilities23,30-32. Further was the requirement that these technologies should be effective for use in high-incidence, low-resource countries where the TB burden is high and medical infrastructure is generally lacking6,21,30.
The non-molecular ‘gold-standard’ for detection of MTb and investigation of antibiotic resistance is culturing of a sample from a patient. However, culturing requires trained lab technicians and is typically extremely slow. The current ‘gold-standard’ molecular assay for detection of MTb and investigation of rifampicin (RIF) resistance (a surrogate marker for MDR-TB) is the Xpert MTB/RIF assay, a cartridge-based nucleic acid amplification test which can give rapid results. This test is easy to use, however, it can only identify RIF resistance so cannot diagnose XDR-TB 33.
The FIND (Foundation for Innovative New Diagnostics) Seq&Treat programme (https://www.finddx.org/tb/seq-treat/) specifically called for the development of targeted next generation sequencing (tNGS) based tests for DR-TB that that could be evaluated by FIND and potentially endorsed by the WHO. Sequencing-based tests have the potential to detect all resistance associated SNPs, thereby determine which drugs will work best against the MTB strain infecting the patient (Kayomo et al. Sci Rep 10, 10786 (2020). https://doi.org/10.1038/s41598-020-67479-4).
tNGS allows sequencing of specific areas of the genome using next generation sequencing to detect variants within the regions of interest. There are different approaches to targeted sequencing, the most common being amplicon sequencing, which uses PCR primers to amplify the sequence/s of interest.
When multiple genes are to be targeted, multiplex polymerase chain reactions (multiplex PCRs) may be used to amplify several different DNA target sequences simultaneously. This process amplifies DNA in samples using multiple primers and a temperature-mediated DNA polymerase in a thermal cycler.
As drug-resistant SNPs are present at multiple sites across the genome, multiple regions need to be targeted by PCR. Multiplex PCR offers substantial advantages over amplification of single regions in separate reactions including higher throughput, cost savings (fewer deoxyribonucleotide triphosphates, enzymes, and other consumables required), turnaround time and production of more data from limited starting material.
Primer design for multiplexed PCR is, however, complex. The primers must have similar annealing temperatures, each pair needs to be specific for its target, and primer pairs should amplify similar sized PCR product to ensure similar amplification efficiency between the multiple targets in the reaction. In addition, as interaction between primers in multiplex reactions can reduce efficiency of amplification and the more primers in a reaction, the more likely this will occur. Designing efficient, sensitive and specific multiplex PCRs is challenging, and success is not assured.
Deeplex® Myc-TB, developed by Genoscreen, is an example of a targeted DR-TB test for prediction of resistance to 15 anti-tuberculous drugs, based on Illumina short read sequencing 34,35(other tests have been developed but all have similar sensitivity and turnaround time). This test takes approximately 2 days to perform and has a limit of detection of ~1000 MTB cells. There remains a need for a more rapid and sensitive test.
It is an aim of the present invention to address the abovementioned problems and meet the abovementioned needs. Accordingly, it is an aim of the present invention to provide a method for rapidly and accurately detecting and/or identifying the presence of drug resistant mutations in a sample from subjects with suspected or confirmed TB using tNGS. It is a further aim to develop primers for achieving this objective, with a focus on the development of primers for use in multiplex PCRs. It is a further aim of the present invention to provide an assay or kit that addresses the abovementioned problems.
Single nucleotide polymorphisms (SNPs) known to confer resistance to first and second-line anti-TB drugs were selected, and primers developed for the selected targets and optimized for use in multiplex PCR. The gene targets were: eis, embB, rrs, rv0678, fabG1, gyrA, rpoB, ethA, rplC, katG, gidB, inha, rrl, pncA, rpsL, tlyA.
Accordingly, in a first aspect there is provided one or more oligonucleotide primer sets for amplifying a portion of one or more genes from M. tuberculosis and/or related bacteria in the M. tuberculosis complex selected from the group comprising or consisting of one or more of eis, embB, ethA, fabG1, gidB, gyrA, inha, katG, pncA, rrl, rplC, rpoB, rpsL, rrs, rv0678 and thyA, wherein each set comprises a pair of forward and reverse primers specific for said portion, wherein each primer has a sequence as set out in SEQ ID Nos. 1-33. Preferably, the one or more sets of primers are selected from SEQ ID Nos. 1-32.
In some embodiments, the oligonucleotide primer sets comprise or consist of one or more of SEQ ID Nos. 1 and 2; 3 and 4; 5 and 6; 7 and 8; 9 and 10; 11 and 12; 13 and 14; 15 and 16; 17 and 18; 19 and 20; 21 and 22; 23 and 24; 25 and 26; 27 and 28; 29 and 30; 31 and 32; and 19 and 33.
In some embodiments, the portion of the one or more genes contains one or more mutations, preferably one or more mutations that confer antibiotic resistance, preferably wherein the one or mutations are one or more single nucleotide polymorphisms that confer antibiotic resistance. In some embodiments, the antibiotic resistance is to one or more of ethambutol, isoniazid, pyrazinamide, rifampicin, streptomycin, amikacin, bedaquiline, capreomycin, ciprofloxacin, clofazimine, ethionamide, kanamycin, linezolid, moxifloxacin, ofloxacin and quinolones.
In some embodiments, the one or more genes are from the MTBc.
In some embodiments, the sets of oligonucleotide primers can be used for multiplex PCR. Sets of primers can thus be grouped into multiplex groups. In some embodiments, one or more multiplex groups can be formed. In some embodiments, multiplex groups can be formed each comprising one or more oligonucleotide primer sets as set out in SEQ ID Nos. 1-33, preferably SEQ ID Nos. 1-32. In some embodiments, one or more multiplex groups can be formed, each comprising oligonucleotide primer sets comprising or consisting of one or more of SEQ ID Nos. 1 and 2; 3 and 4; 5 and 6; 7 and 8; 9 and 10; 11 and 12; 13 and 14; 15 and 16; 17 and 18; 19 and 20; 21 and 22; 23 and 24; 25 and 26; 27 and 28; 29 and 30; 31 and 32; and 19 and 33.
In some embodiments, a multiplex group can comprise oligonucleotide primer sets for amplifying a portion of eis, embB, rrs, nv0678, and fabG1 (Group 1). In a further embodiment, a multiplex group can comprise oligonucleotide primer sets for amplifying a portion of gyrA, rpoB, ethA, rplC, and katG (Group 2). In a further embodiment, a multiplex group can comprise oligonucleotide primer sets for amplifying a portion of gidB, inhA, rrl, pncA, rpsL, and tlyA (Group 3). Accordingly, in some embodiments, groups of oligonucleotide primer sets comprise or consist of one or more of SEQ ID Nos. 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 (Group 1 in Table 7); one or more of SEQ ID Nos. 11, 12, 13, 14, 15, 16, 17, 18, 19 and 20 (Group 2 in Table 7); and/or one or more of SEQ ID Nos. 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 and 32 (Group 3 in Table 7).
Accordingly, in one embodiment there is provided one or more multiplex groups of oligonucleotide primer sets for amplifying a portion of genes from M. tuberculosis and/or related bacteria in the MTBc selected from the group comprising or consisting of one or more of eis, embB, ethA,fabG1, gidB, gyrA, inhA, katG, pncA,rrl, rplC, rpoB, rpsL, rrs, rv0678, tlyA, wherein each oligonucleotide primer set comprises or consists of a pair of forward and reverse primers specific for said portion, wherein the multiplex groups of oligonucleotide primer sets comprise or consist of one or more of SEQ ID Nos. 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 (Group 1 in Table 7); one or more of SEQ ID Nos. 11, 12, 13, 14, 15, 16, 17, 18, 19 and 20 (Group 2 in Table 7); and/or one or more of SEQ ID Nos. 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 and 32 (Group 3 in Table 7).
In some such embodiments, the multiplex groups of oligonucleotide primer sets comprise or consist of one or more of SEQ ID Nos. 1 and 2; 3 and 4; 5 and 6; 7 and 8; and 9 and 10 (Group 1 in Table 7); one or more of SEQ ID Nos. 11 and 12; 13 and 14; 15 and 16; 17 and 18; and 19 and 20 (Group 2 in Table 7); and/or one or more of SEQ ID Nos. 21 and 22; 23 and 24; 25 and 26; 27 and 28; 29 and 30; and 31 and 32 (Group 3 in Table 7). In some embodiments, a multiplex group of oligonucleotide primer sets comprises or consists of one or more of SEQ ID Nos. 1 and 2; 3 and 4; 5 and 6; 7 and 8; and 9 and 10 (Group 1 in Table 7). In some embodiments, a multiplex group of oligonucleotide primer sets comprises or consists of one or more of SEQ ID Nos. 11 and 12; 13 and 14; 15 and 16; 17 and 18; and 19 and 20 (Group 2 in Table 7). In some embodiments, a multiplex group of oligonucleotide primer sets comprises or consists of one or more of SEQ ID Nos. 21 and 22; 23 and 24; 25 and 26; 27 and 28; 29 and 30; and 31 and 32 (Group 3 in Table 7).
In a second aspect there is provided a multiplex PCR reaction mixture comprising one or more groups of oligonucleotide primer sets for amplifying a portion of one or more genes from M. tuberculosis and/or related bacteria in the M. tuberculosis complex selected from the group comprising or consisting of one or more of eis, embB,ethA, fabG1, gidB, gyrA, inha, katG, pncA, rrl, rplC, rpoB,rpsL, rrs, rv0678, tlyA, wherein each set comprises a pair of forward and reverse primers specific for said portion, wherein the groups of oligonucleotide primer sets comprise or consist of one or more of SEQ ID Nos. 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 (Group 1 in Table 7); one or more of SEQ ID Nos. 11, 12, 13, 14, 15, 16, 17, 18, 19 and 20 (Group 2 in Table 7); and/or one or more of SEQ ID Nos. 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 and 32 (Group 3 in Table 7).
In some embodiments, a multiplex PCR reaction mixture comprises a group of oligonucleotide primer sets comprising or consisting of one or more of SEQ ID Nos. 1 and 2; 3 and 4; 5 and 6; 7 and 8; and 9 and 10 (Group 1 in Table 7); one or more of SEQ ID Nos. 11 and 12; 13 and 14; 15 and 16; 17 and 18; and 19 and 20 (Group 2 in Table 7); and/or one or more of SEQ ID Nos. 21 and 22; 23 and 24; 25 and 26; 27 and 28; 29 and 30; and 31 and 32 (Group 3 in Table 7). In one embodiment, a multiplex PCR reaction mixture comprises a group of oligonucleotide primer sets comprising or consisting of SEQ ID Nos. 1 and 2; 3 and 4; 5 and 6; 7 and 8; and 9 and 10 (Group 1 in Table 7). In a further embodiment, a multiplex PCR reaction mixture comprises a group of oligonucleotide primer sets comprising or consisting of one or more of SEQ ID Nos. 11 and 12; 13 and 14; 15 and 16; 17 and 18; and 19 and 20 (Group 2 in Table 7). In a further embodiment, a multiplex PCR reaction mixture comprises a group of oligonucleotide primer sets comprising or consisting of SEQ ID Nos. 21 and 22; 23 and 24; 25 and 26; 27 and 28; 29 and 30; and 31 and 32 (Group 3 in Table 7).
The multiplex PCR reaction mixture may comprise further ingredients and reagents required to perform multiplex PCR, such as buffers, deoxynucleotide triphosphates (dNTPs), DMSO, water and DNA polymerase.
In some multiplex embodiments, said primers may be mixed to a working concentration of 0.2 µM. Further typically with the exception of tlyA which requires a working concentration of 0.3 µM, for consistent target amplification.
In some embodiments, the portion of the one or more genes from M. tuberculosis and/or related bacteria in the M. tuberculosis complex is obtained from a sample from a subject suspected or confirmed to have TB. The sample may be one or more tissues and/or bodily fluids obtained from the subject, including one or more of sputum; urine; blood; plasma; serum; synovial fluid; pus; cerebrospinal fluid; pleural fluid; pericardial fluid; ascitic fluid; sweat; saliva; tears; vaginal fluid; semen; interstitial fluid; bronchoalveolar lavage; bronchial wash; gastric lavage; gastric wash; a transtracheal or transbronchial fine needle aspiration; bone marrow; pleural tissue; tissue from a lymph node, mediastinoscopy, thoracoscopy or transbronchial biopsy; or combinations thereof; or a culture specimen of one or more tissues and/or bodily fluids obtained from a subject suspected of having or confirmed to have TB. Typically, the sample includes cells and/or DNA from M. tuberculosis and/or related bacteria in the M. tuberculosis complex.
In a third aspect there is provided a method of detecting and/or identifying the presence of one or more mutations that confer antibiotic resistance in a sample comprising DNA from Mycobacterium tuberculosis and/or related bacteria in the M. tuberculosis complex, said method including the steps of;
In some embodiments, the mutations are within one or more genes selected from the group consisting of one or more of eis, embB, ethA, fabG1, gidB, gyrA, inha, katG, pncA, rrl, rplC, rpoB, rpsL, rrs, rv0678 and tlyA.
In some embodiments the mutations are one or more single nucleotide polymorphisms.
In some embodiments, the antibiotic resistance is to one or more of ethambutol, isoniazid, pyrazinamide, rifampicin, streptomycin, amikacin, bedaquiline, capreomycin, ciprofloxacin, clofazimine, ethionamide, kanamycin, linezolid, moxifloxacin, ofloxacin and quinolones.
The amplification step uses one or more groups of oligonucleotide primer sets. In some embodiments, the groups of oligonucleotide primer sets comprise or consist of one or more forward and reverse primer pairs selected from SEQ ID Nos. 1 and 2; 3 and 4; 5 and 6; 7 and 8; 9 and 10; 11 and 12; 13 and 14; 15 and 16; 17 and 18; 19 and 20; 21 and 22; 23 and 24; 25 and 26; 27 and 28; 29 and 30; 31 and 32 and 19 and 33.
In some embodiments, the one or more groups of oligonucleotide primer sets comprise or consist of one or more of SEQ ID Nos. 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 (Group 1 in Table 7); one or more of SEQ ID Nos. 11, 12, 13, 14, 15, 16, 17, 18, 19 and 20 (Group 2 in Table 7) and/or one or more of SEQ ID Nos. 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 and 32 (Group 3 in Table 7). In some embodiments, the amplification step uses a group of oligonucleotide primer sets consisting of SEQ ID Nos. 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 (Group 1). In some embodiments, the amplification step uses a group of oligonucleotide primer sets consisting of SEQ ID Nos. 11, 12, 13, 14, 15, 16, 17, 18, 19 and 20 (Group 2 in Table 7). In some embodiments the amplification step uses a group of oligonucleotide primer sets consisting of SEQ ID Nos. 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 and 32 (Group 3 in Table 7).
Detection of a mutation is indicative of antibiotic resistance. Identification of the mutation informs or allows identification of the nature of the antibiotic resistance (i.e. the antibiotic to which the bacteria is resistant).
Accordingly, in a fourth aspect there is provided a method of predicting whether a patient suffering from tuberculosis will respond to treatment with one or more of ethambutol, isoniazid, pyrazinamide, rifampicin, streptomycin, amikacin, bedaquiline, capreomycin, ciprofloxacin, clofazimine, ethionamide, kanamycin, linezolid and moxifloxacin, said method comprising a step of determining the presence of one or more drug resistant mutations in one or more genes selected from the group comprising one or more of eis, embB, ethA, f fabG1, gidB, gyrA, inha, katG, pncA, rrl, rplC, rpoB, rpsL, rrs, rv0678 and tlyA in DNA obtained from a sample from the patient, the method comprising:
In some embodiments, the method is for predicting whether a patient suffering from tuberculosis will respond to treatment with one or more of ethambutol, isoniazid, streptomycin, amikacin, bedaquiline, capreomycin, clofazimine, ethionamide, kanamycin, wherein the one or more genes are eis, embB, rrs, rv0678, and fabG1; and the group of oligonucleotide primer sets consists of SEQ ID Nos. 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 (Group 1 in Table 7).
In some embodiments, the method is for predicting whether a patient suffering from tuberculosis will respond to treatment with one or more of isoniazid, rifampicin, ciprofloxacin, ethionamide, linezolid, moxifloxacin, ofloxacin and quinolones, wherein the one or more genes are gyrA, rpoB, ethA, rplC, and katG; and the group of oligonucleotide primer sets consists of SEQ ID Nos. 11, 12, 13, 14, 15, 16, 17, 18, 19 and 20 (Group 2 in Table 7).
In some embodiments, the method is for predicting whether a patient suffering from tuberculosis will respond to treatment with one or more of pyrazinamide, streptomycin, capreomycin and ethionamide, wherein the one or more genes are gidB, inha, rrl, pncA, rpsL, and tlyA; and the group of oligonucleotide primer sets consists of SEQ ID Nos. 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 and 32 (Group 3 in Table 7).
In some embodiments according to the third or fourth aspect, the DNA is from M. tuberculosis.
In some embodiments according to the third or fourth aspect, the sample is a clinical sample. The sample may be one or more tissues and/or bodily fluids obtained from a subjected suspected of having or confirmed to have TB, including one or more of sputum; urine; blood; plasma; serum; synovial fluid; pus; cerebrospinal fluid; pleural fluid; pericardial fluid; ascitic fluid; sweat; saliva; tears; vaginal fluid; semen; interstitial fluid; bronchoalveolar lavage; bronchial wash; gastric lavage; gastric wash; a transtracheal or transbronchial fine needle aspiration; bone marrow; pleural tissue; tissue from a lymph node, mediastinoscopy, thoracoscopy or transbronchial biopsy; or combinations thereof; or a culture specimen of one or more tissues and/or bodily fluids obtained from a subject suspected of having or confirmed to have TB. Typically, the sample includes cells and/or DNA from M. tuberculosis and/or related bacteria in the M. tuberculosis complex. In some embodiments, the sample is a sputum sample from a subject suspected or confirmed to have TB.
In some embodiments, the samples undergo mechanical disruption in order to disrupt the cells in the sample and achieve cell lysis. Any suitable means may be used, for example bead beating.
The step of isolating or extracting DNA from the sample may be carried out by any suitable means, including by the use of an appropriate kit, using given or standard protocols. For example, a Maxwell RSC PureFood Pathogen Kit from Promega AS1660, with instructions for use. In some embodiments, a Maxwell RSC PureFood Pathogen Kit from Promega AS1660 may be used. In some such embodiments, the following modifications were made from the kit instructions: The kit teaches use of a 800 µl sample; in some embodiments, a 400 µl sample after bead beating was used. The kit teaches adding 200 µl lysis buffer A and incubating at 56° C. for 4 min with shaking; in some embodiments, 200 µl lysis buffer A was added together with 40 µl Proteinase k, with incubation at 65° C. for 10 min. The kit teaches addition of 300 µl of lysis buffer and then placing the sample on the robot; in some embodiments, 300 µl lysis buffer was added together with 400 µl PBS and the sample was then placed on the robot.
In embodiments according to the third or fourth aspect wherein more than one group of primer sets are used for the amplification step, each group may be run as a separate multiplex group template.
Labelled nucleotides or labelled primers may be used in the amplification of the DNA for the purpose of, for example, quality control. For example, a fluorescent DNA-binding dye may be added to enable DNA quantitation. Any suitable dyes or probes with dyes may be used, such as probes with fluorescent dyes, such as use of a sybr green assay such as Roche Lightcycler® 480 SYBR Green I master.
In embodiments wherein more than one group of primer sets are used for the amplification step and each group is run as a separate multiplex group template, one or more multiplex group templates may be pooled to make a single template for DNA quantitation and/or sequencing.
Samples may then undergo barcode ligation and adaptor ligation to create a library for sequencing. Barcoding can be used when the amount of data required per sample is less than the total amount of data that can be generated: it allows pooling of multiple samples and sequencing of them together. Any suitable means may be used, including the use of barcoding kits, using given or standard protocols. For example, Oxford Nanopore Technologies provides amplicon barcoding with native barcoding expansion 96 (EXP-NBD196 and SQK-LSK109), including instructions for use. In some embodiments, the Oxford Nanopore Technologies amplicon barcoding with native barcoding expansion 96 (EXP-NBD196 and SQK-LSK109) may be used following the instructions for use provided.
The DNA sequencing step may be carried out by any suitable means. In preferred embodiments, the DNA sequencing is tNGS or third-generation sequencing (also known as long-read sequencing). Third-generation sequencing may be carried out using Oxford Nanopore Technologies’ MinION, or PacBio’s sequencing platform of single molecule real time sequencing (SMRT). Oxford Nanopore’s sequencing technology is based on detecting the changes in electrical current passing through a nanopore as a piece of DNA moves through the pore. The current measurably changes as the bases G, A, T and C pass through the pore in different combinations. SMRT is based on the properties of zero-mode waveguides. Signals in the form of fluorescent light emission from each nucleotide are incorporated by a DNA polymerase bound to the bottom of the zL well. In preferred embodiments the sequencing is long-read nanopore sequencing.
The step of detecting of one or more mutations may be carried out by any suitable method, such as suitable bioinformatics tools and programmes. In some embodiments, the Oxford nanopore technologies workflow for TB may be used in desktop program EPI2ME with the FASTQ TB RESISTANCE PROFILE v2020.03.11.
The oligonucleotide primer sets of the first aspect, the PCR reaction mixture of the second aspect and/or the method of the third aspect can be used to identify both the presence and identity of drug resistance mutations in the genes of TB bacteria from a particular subject. Such information informs decisions regarding drug administration and allows a tailored treatment regime to be determined for the patient depending upon the identified mutations.
As such, in a fifth aspect, there is provided a method for determining an appropriate antibiotic treatment regime for a patient with tuberculosis, comprising detecting and/or identifying the presence of one or more mutations that confer antibiotic resistance in a sample from the patient according to the third aspect, and determining an appropriate antibiotic regime on the basis of the mutations detected/identified. The disclosure herein also provides a method of assigning a patient with tuberculosis to one of a certain number of treatment pathways comprising detecting and/or identifying the presence of one or more mutations that confer antibiotic resistance in a sample from the patient using a method according to the third aspect, and assigning the patient to a treatment regime on the basis of the mutations detected/identified.
In a sixth aspect there is provided a kit comprising one or more oligonucleotide primer sets or groups of oligonucleotide primer sets according to the first aspect. The kit may be used to carry out a method according to one or more of steps (a) (b) or (c) of the third aspect. The kit may further comprise ingredients and reagents required to carry out the method according to one or more of steps (a) (b) or (c) of the third aspect, including buffers, DNA polymerase and nucleotides. In some embodiments, the kit further comprises reagents required for the amplification of the gene regions between the primers. The kit may further comprise a sample collection container for receiving the sample. Samples may be processed according to the method of the third aspect immediately, alternatively they may be stored at low temperatures, for example in a fridge or freezer before the method is carried out. The sample may be processed before the method is carried out. For instance, a sedimentation assay may be carried out, and/or a preservative and/or dilutant may be added. Thus, the sample collection container may contain suitable processing solutions, such as buffers, preservative and dilutants.
Gene targets and their corresponding primer pairs according to the disclosure herein are as shown in Table 1.
For a better understanding of the invention, and to show how embodiments of the same may be carried into effect, reference will now be made, by way of example, to the accompanying Figures in which:
Selected target single nucleotide polymorphisms (SNPs) that confer resistance to first and second-line anti-TB drugs were chosen primarily from WHO/FIND evidence published in the WHO next-generation sequencing technical guide 36. The targets for rpsL were selected from prior literature by Karimi, et al. and Meier, et al37,38. Targets for gidB were selected on evidence from Villellas, et al39. Targets for ethA were selected on evidence from Morlock, et al40. Targets for embB were selected on evidence from Zhao, et al41. Finally, targets for tlyA were selected from prior literature by Maus, et al42.
Base positions and genes as listed are based on the H37Rv M. tuberculosis reference genome available through the NCBI database (NC_000962.3)43. Targeted mutations were identified either as their codon location or their nucleotide location. Mutations were identified by the codon which they effect when the SNP occurs within an annotated gene region and the prior literature explicitly states the altered amino acid. Targets were listed by nucleotide mutation in the event they occur within a gene promoter region or the supporting literature does not explicitly identify the amino acid mutation. These promoter region SNPs are further identified by a “-” prior to its position indicating it occurs before the annotated gene. The effect of the mutated base is also included; e.g. Asparagine to Histidine or nucleotide A to nucleotide C (Table A, appended).
Primers were developed for the chosen gene/promotor targets (n=16; Table 2) that amplified ~1000 bp regions containing the targeted SNPs of interest. As discussed above, interaction between primers in multiplex reactions can reduce efficiency of amplification and the more primers in a reaction, the more likely this will occur. Therefore designing efficient, sensitive and specific multiplex PCRs is complex.
The following genes were targeted in the DR-TB sequencing assay: eis, embB, ethA, fabG1, gidB, gyrA, inha, katG, pncA, rrl, rplC, rpoB, rpsL, rrs, rv0678, tlyA. Initially, gene target primer pairs were grouped into 5 sets of three (Table 3). DNA was extracted from M. bovis BCG and used to test the specificity and sensitivity of the triplex assays.
The multiplex PCRs were performed as follows:
Nested qPCR was performed on the amplified products from the multiplex PCR to evaluate the amplification of all the targets. Nested PCR on all amplified products resulted in very similar Ct values, indicating the same amplification efficiency across all primers (
While the triplex assays worked well, the requirement for 5 PCR reactions was considered too laborious and expensive for the tNGS assay. Hence, the primer pairs were combined in a new format to make three groups (two 5-plex and one 6-plex reaction), in order to simplify the assay. Multiplex efficiency was again measured by nested qPCR (
Formulations 1-6 had multiple late Cts and/or total dropouts indicative of inhibition and competition within the multiplex groups. Version 7 showed multiplex groups 2 and 3 had Ct ranges <1.5 while group 1 had a range of approximately 15Cts (
Concurrently to optimising the group formulations, various primers were redesigned to overcome primer interactions. In total there were 48 multiplex primer combinations with >300 primer designs (Table 5) before the optimal sequences were determined.
After testing ~400 samples provided by FIND in a lab validation study (described below), a redesign was required for the katG reverse primer to avoid a common non-resistance conferring SNP in the primer binding site. To overcome this, five new reverse primers were tested where each primer was shifted towards the 3′ 1 bp at a time (up to 5bp shift) (Table 6). Option 5 was selected for the final assay as the mutation site was avoided and the performance of the assay wasn’t negatively affected.
The final optimal iteration of primers consisted of two 5-plex groups and one 6-plex group (Table 7).
Visualized target regions are shown as either the parent or complement strand depending on gene orientation. Target regions were designed to be 900-1100 bp long as this is a good size for PCR and nanopore sequencing. Keeping the PCR products a uniform size reduces bias toward certain targets in multiplex PCR and sequencing reactions.
The target region for identified eis mutations encompasses the promoter region, denoted in bold text, of the 1,209 base pair eis gene. The eis gene is on the complement strand. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-TGTCGGGTACCTTTCGAGC-3′ [Sequence ID No. 1]
Reverse Primer: 5′-TCCATGTACAGCGCCATCC-3′ [Sequence ID No. 2]
TGTCGGGTACCTTTCGAGCCGCCGAGCTGACCGCGGCGGAACTAGGGTCC
The embB target region on the parent strand is a subsection of the overall 3,297 base pair embB gene. The region chosen contains all the high confidence SNPS and the majority of known embB SNPs. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-CGCCGTGGTGATATTCGGC-3′ [Sequence ID No. 3]
Reverse Primer: 5′-GCACACCGTAGCTGGAGAC-3′ [Sequence ID No. 4]
CGCCGTGGTGATATTCGGCTTCCTGCTCTGGCATGTCATCGGCGCGAATT
The rrs primers target includes a subset of the 1,537 base pair rrs gene on the parent strand and some sequence outside the gene at the 3′ end as some of the target SNPs are at the 3′ end of the gene. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-CTCTGGGCAGTAACTGACGC-3′ [Sequence ID No. 5]
Reverse Primer: 5′-GAGTGTTGCCTCAGGACCC-3′ [Sequence ID No. 6]
CTCTGGGCAGTAACTGACGCTGAGGAGCGAAAGCGTGGGGAGCGAACAGG
The rv0678 target region contains the entire 498 base pair rv0678 gene on the parent strand along with intergenic regions on either side. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-GCTCGTCCTTCACTTCGCC-3′ [Sequence ID No. 7]
Reverse Primer: 5′-ATCAGTCGTCCTCTCCGGT-3′ [Sequence ID No. 8]
GCTCGTCCTTCACTTCGCCATCGACGGTGATTCGGCAGGTGATGGAAGTG
The fabG1 target region covers the 744 bp fabG1 gene on the parent strand along the gene promoter region (denoted in bold), targeting the high confidence SNPs located therein, and some intergenic sequence at the 3′ end. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-CTTTTGCACGCAATTGCGC-3′ [Sequence ID No. 9]
Reverse Primer: 5′-AGCAGTCCTGTCATGTGCG-3′ [Sequence ID No. 10]AAGTGTGC
CTTTTHCACGAATTGCGCGGYCAGTTCCACACCCTGCGGCACGTACACGT
The gyrA target region is a subset of the overall 2,517 bp gyrA gene on the parent strand. This target region was designed to encompass all the high confidence gyrA resistance-conferring SNPs.. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-TGACAGACACGACGTTGCC-3′ [Sequence ID No. 11]
Reverse Primer: 5′-CGATCGCTAGCATGTTGGC-3′ [Sequence ID No. 12]
TGACAGACACGACGTTGCCGCCTGACGACTCGCTCGACCGGATCGAACCG
The rpoB target region is a subset of the 3,519 bp rpoB gene on the parent strand. This target region was designed to encompass all the high confidence rpoB resistance-conferring SNPs. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-TCATCATCAACGGGACCGAG-3′ [Sequence ID No. 13]
Reverse Primer: 5′-ACACGATCTCGTCGCTAACC-3′ [Sequence ID No. 14]
TCATCATCAACGGGACCGAGCGTGTGGTGGTCAGCCAGCTGGTGCGGTCG
The ethA target region covers a subset of the 1470 base pair ethA gene on the complement strand. This section was chosen to cover the high confidence SNPs located at the 5′ end of the gene. Sequence outside the annotated gene is underlined. Forward and reverse primer locations are written italics.
Forward Primer: 5′- TGGATCCATGACCGAGCAC -3′ [Sequence ID No. 15]
Reverse Primer: 5′- GTCCAGGAGGCATTGGTGT -3′ [Sequence ID No. 16]
TGGATCCATGACCGAGCACCTCGACGTTGTCATCGTGGGCGCTGGAATCT
The rplC target region contains the entire 654 bp rplC gene on the parent strand along with intergenic regions on the 5′ and 3′ ends. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-AGTACAAGGACTCGCGGGA-3′ [Sequence ID No. 17]
Reverse Primer: 5′-TCGAGTGGGTACCCTGGC-3′ [Sequence ID No. 18]
AGTACAAGGAGTCGCGGGAGCACTTCGAGATGCGCACACACAAGCGGTTG
The katG target region is a subset of the 2,223 base pair katG gene, which is on the complement strand. The region was chosen to cover all high confidence SNPs. Forward and reverse primer locations are highlighted in italics.
Forward Primer: 5′- CTGTGGCCGGTCAAGAAGA -3′ [Sequence ID No.19]
Reverse Primer: 5′- TGCCCGGATCTGGCTCTTA -3′ [Sequence ID No.33]
CTGTGGCCGGTCAAGAAGAAGTACGGCAAGAAGCTCTCATGGGCGGACCT
The katG target region is a subset of the 2,223 bp katG gene, which is on the complement strand. The region was chosen to cover all the high confidence SNPs. Forward and reverse primer locations are written in italics.
Forward Primer: 5′- CTGTGGCCGGTCAAGAAGA -3′ [Sequence ID No. 19]
Reverse Primer: 5′- GGATCTGGCTCTTAAGGCTGG -3′ [Sequence ID No. 20]
CTGTGGCCGGTCAAGAAGAAGTACGGCAAGAAGCTCTCATGGGCGGACCT
The gidB target region contains the entire 675 bp gidB gene on the parent strand along with intergenic sequence on the 5′ and 3′ ends. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-TGACACAGACCTCACGAGC-3′ [Sequence ID No. 21]
Reverse Primer: 5′-GCCCTTCTGATTCGCGATG-3′ [Sequence ID No. 22]
TGACACAGACCTCAGGAGCCGGCGGAGTGCGTAATGTCTCCGATCGAGCC
The inha target region contains a subset of the inhA 810 bp gene on the parent strand along with the promoter region, denoted in bold, to cover all the high confidence SNPs in the gene and promotor. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are highlighted in italics.
Forward Primer: 5′-GGGCGCTGCAATTTATCCC-3′ [Sequence ID No. 23]
Reverse Primer: 5′-GGCGTAGATGATGTCACCC-3′ [Sequence ID No. 24]
GGGCGCTGCAATTTATCCCAGCGAAGCGGGTCGGCACCCCCGCCGAGGTC
The rrl target region is a subsection of the overall 3,138 bp rrl gene on the parent strand, targeting all the high confidence SNPs. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-GGTCCGTGCGAAGTCGC-3′ [Sequence ID No. 25]
Reverse Primer: 5′-TGAACCCGTGTTCTGCGG-3′ [Sequence ID No. 26]
GGTCCGTGCGAAGTCGCAAGACGATGTATACGGACTGACGCCTGCCCGGT
The pncA target region contains the entire 561 base pair pncA gene on the complement strand along with intergenic regions at the 5′ and 3′ ends. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-TCACCGGACGGATTTGTCG-3′ [Sequence ID No. 27]
Reverse Primer: 5′-TCCAGATCGCGATGGAACG-3′ [Sequence ID No. 28]
TCACCGGACGGATTTGTCGCTCACTACATCACCGGCGTGATCTATCCCGC
The rpsL target region contains the entire 375 bp rpsL gene on the parent strand along with intergenic regions at the 5′ and 3′ ends. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-GCGGCGGGTATTGTGGTT-3′ [Sequence ID No. 29]
Reverse Primer: 5′-TAACCGGCGCTTCTCACC-3′ [Sequence ID No. 30]
GCGGCGGGTATTGTGGTTGCTCGTGCCTGGCGGCTTACGCTTGATGTAGG
The tlyA target region contains the entire 807 base pair tfyA gene on the parent strand along with intergenic regions at the 5′ and 3′ ends. Sequence outside the annotated gene is highlighted in grey. Forward and reverse primer locations are written in italics.
Forward Primer: 5′-CGTTGATGCGCAGCGATC-3′ [Sequence ID No. 31]
Reverse Primer: 5′-GGTCTCGGTGGCTTCGTC-3′ [Sequence ID No. 32]
CGTTGATGCGCAGCGATCATCCGGTGACTAGCGTAGGAACGCAATGACCA
The present disclosure provides a means of accurately and rapidly identifying the presence of multiple drug resistance mutations in a sample from a patient with suspected or confirmed Tuberculosis. Such information informs decisions regarding drug administration, and allows a tailored regimen to be determined for the patient depending upon the identified mutations. Furthermore, the disclosed methods can be successfully carried out on samples taken directly from patients, such as sputum, thereby adding to their potential for use in lower and middle income and developing countries. The development of optimised primers for this purpose means the advantages of using a multiplex assay can be realised. The disclosed methods are highly sensitive (<100 MTB cells), rapid (taking approximately 8 hours) and can detect a broad range of mutations, and thus represent a major improvement over current culture, molecular (e.g. GenoType MTBDRsl line probe assay) and tNGS based tests. This allows the correct treatment pathway to be determined and for patients to commence treatment promptly and not be lost to follow-up (a major problem in developing countries). This reduces the spread of disease and helps prevent the development of drug-resistant bacterial strains.
Wherever the term ‘comprising’ is used herein we also contemplate options wherein the terms ‘consisting of or ‘consisting essentially of are used instead. In addition, any and all liquid compositions described herein can be aqueous solutions. Note too that whenever the phrase “one or more” is used for a range, for example in relation to a number of sequences W, X, Y and Z (“one or more of SEQ ID Nos. W, X, Y and Z”) this is a disclosure of each value alone (SEQ ID No. W; SEQ ID No. X; SEQ ID No. Y; SEQ ID No. Z), or in combination, e.g. SEQ ID Nos. W and X and SEQ ID No. Y and Z). Similarly, whenever the phrase “one or more” is used in relation to a range of pairs, for example in relation to a number of pairs of sequences (“one or more of SEQ ID Nos. W and X; and Y and Z”) this is a disclosure of each pair alone (SEQ ID No. W and X) or in combination (e.g. SEQ ID Nos. W and X and SEQ ID Nos. Y and Z).
The following Examples are provided to illustrate embodiments of the present invention and should not be construed as limiting thereof.
A study was conducted using sputum spiked with well characterized M. tuberculosis isolates (whole-genome sequence and culture confirmed resistance profiles) to evaluate the developed primers and method. DNA was extracted on the MagNA Pure Compact, PCR amplified in 3 multiplex reactions per sample, pooled, washed, barcoded, and sequenced on the MinION in batches of 80 as described below.
DNA Extraction:
Multiplex PCR:
1. Prepare 3 multiplex 10x primer mixes as follows:
2. In MSC-II mix PCR Master Mix (Qiagen Multiplex PCR kit) for each multiplex primer group in the following ratio per sample:
3. In MSC-II add 45 µL mastermix to 0.2 mL thin-walled PCR tubes.
a. Each sample requires three tubes, one for each Multiplex Primer Group.
4. In MSC-II carefully add 5 µL extracted DNA to PCR tubes.
5. In MSC-II seal PCR tubes tightly and vortex.
6. In MSC-II briefly spin down PCR tubes and remove bubbles.
7. Load PCR tubes into a thermocycler and run an amplification protocol with the following parameters:
8. Carefully remove PCR tubes and return to MSC-II.
9. In MSC-II transfer PCR product to clean PCR tubes.
10. Submerge clean PCR tubes in a 1:16 dilution of Bioguard for minimum 30 seconds for removal from CL3.
The three multiplex reactions for each sample are then pooled as follows:
The pooled samples were then prepared for nanopore sequencing as follows:
1. Transfer 45 µL of pooled DNA to a thin-walled PCR plate
2. Add following reagents to the DNA
3. Mix by pipette
4. Spin down tube and incubate for 5 minutes at 20° C. followed by 5 minutes at 65° C.
5. Transfer samples to a clean 96-well plate
6. Perform a 1x bead wash by adding 60 µL AMPure XP Beads
7. Incubate sample for 5 minutes on a hula mixer
8. Briefly spin down plate
9. Place plate on magnet-rack and let incubate for 5 minutes
10. Remove supernatant
11. Wash bead pellet with 180 µL 70% ethanol
12. Remove supernatant
13. Wash bead pellet with 180 µL 70% ethanol
14. Remove supernatant
15. Briefly spin down plate and return to magnet-rack
16. Remove residual supernatant
17. Air dry pellet for approximately 30 seconds
18. Resuspend pellet in 31 µL nuclease free H2O
19. Incubate samples for 2 minutes at room temperature
20. Return plate to magnet-rack and pellet beads for 2 minutes
21. Carefully remove eluted supernatant and transfer 30 µL to a clean 96-well plate
1. In a thin-walled PCR plate combine the following:
2. Briefly vortex
3. Spin down samples
4. PCR amplify using the following cycling conditions
5. Perform 0.8x bead wash (40 µL) using AMPure XP beads as described above
6. Resuspend pellet in 45 µL nuclease free H2O
7. Incubate samples for 2 minutes at room temperature
8. Return plate to magnet-rack and pellet beads for 2 minutes
9. Carefully remove eluted supernatant and transfer to a clean 96-well plate.
10. Quantify as described above
11. Pool each barcoded sample equimolar into a fresh 1.5 mL Eppendorf
12. Perform 0.8x bead wash using AMPure XP beads on pooled samples as described above and resuspend in 45 µL nuclease free H2O
1. In a 0.2 mL thin walled PCR tube combine the following:
2. Vortex and briefly spin down
3. Incubate for 5 minutes at 20° C. followed by 5 minutes at 65° C.
4. Transfer sample to a clean 1.5 mL Eppendorf
5. Perform a 0.8x bead wash (48 µL) using AMPure XP beads as described above
6. Resuspend pellet in 61 µL nuclease free H2O
7. Incubate samples for 2 minutes at room temperature
8. Return plate to magnet-rack and pellet beads for 2 minutes
9. Carefully remove eluted supernatant and transfer to a clean 1.5 mL Eppendorf.
1. Thaw and spin down Adapter Mix (AMX), T4 Ligase, Ligation Buffer (LNB), and Elution Buffer (EB) (Oxford Nanopore Technologies Ligation Sequencing Kit SQK-LSK109).
2. Place thawed and vortexed reagents on ice
3. Thaw one tube of Short Fragment Buffer (SFB) at room temperature
a. Vortex and spin down before placing on ice
4. Mix the following in a 1.5 mL Eppendorf in order:
5. Gently mix tube by flicking and spin down
6. Incubate for 10 minutes at room temperature
7. Perform a 0.6x bead wash (60 µL) using AMPure XP beads
8. Incubate samples for 5 minutes on a hula mixer
9. Briefly spin down samples
10. Place tube on magnet-rack and let incubate for 5 minutes
11. Remove supernatant
12. Resuspend pellet in 125 µL SFB
13. Place tube on magnet-rack and let incubate for 10 minutes
14. Carefully remove supernatant
15. Resuspend pellet in 125 µL SFB
16. Place tube on magnet-rack and let incubate for 10 minutes
17. Carefully remove supernatant
18. Briefly spin down tube and return to magnet-rack
19. Remove residual supernatant
20. Air dry pellet for approximately 30 seconds
21. Resuspend pellet in 15 µL EB
22. Incubate at room temperature for 10 minutes
23. Place tube on magnet-rack until elute is clear and colourless
24. Carefully remove and retain 15 µL eluted supernatant in clean 1.5 mL Eppendorf
25. Perform Qubit HS Assay on 1 µL elute
Sequencing library loading on MinION
Resistance to first- and second-line anti-TB drugs was identified using the ONT Epi2Me FastQ TB Resistance Profile pipeline. Wild-type and mutant nucleotides were reported for all drug resistance associated SNP sites detected within the PCR product fastQ sequences. The presence of SNPs in specific target genes indicated resistance to specific anti-TB drugs (Table 8).
This method also allowed for identification of heteroresistance by comparison of the relative number of reads for wild-type compared to the number of reads for mutants (Table 9). Heteroresistance was called when >15% and <80% mutant bases were detected.
Table 9: Example heteroresistance detection results from two sequenced samples. Boxes with vertical stripes signify >80% of reads at that site are resistant associated mutants (resistant, no heteroresistance). Boxes with diagonal stripes signify 51%-79% of reads at that site are resistance associated mutants (heteroresistant, majority resistant bases). Black boxes signify 20%-50% of reads at that site are resistance associated mutants (heteroresistant, majority wild-type bases).
Raw read numbers could also be visualised, providing a more detailed analysis if required (Table 10). These results display the codon or nucleotide location within the annotated gene as well as the number of wild-type or mutant bases recorded at that location.
Following on from Example 1, a set of samples were processed with an altered DNA extraction and simplified library preparation method. Here, DNA was extracted instead using the Promega Maxwell RSC 48 with the PureFood Pathogen kit and within the library preparation alterations were made to the end-prep and barcode/adapter ligation reactions. The resistance profile was compared between methods to ensure the same profile was identified. Details of the method alterations are below:
DNA Extraction:
End Prep
Barcode Ligation
Adapter Ligation:
Resistance to ‘first- and second-line anti-TB drugs was identified using the ONT Epi2Me FastQ TB Resistance Profile pipeline. Wild-type and mutant nucleotides were reported for all drug resistance associated SNP sites detected within the PCR product fastQ sequences. The presence of SNPs (>15% mutant bases) in specific target genes indicated resistance to specific anti-TB drugs (Table 11).
Raw read numbers could also be visualised, providing a more detailed analysis if required (Table 12) e.g. for identifying heteroresistance. These results display the codon or nucleotide location within the annotated gene as well as the number of wild-type or mutant bases recorded at that location.
As can be seen from both results tables the alterations in methodology did not change the resistance profile of this sample. Therefore the optimised method (using the Promega Maxwell and simplified library preparation) would be the method of choice for this assay.
1. Coscolla M, Gagneux S. Seminars in Immunology Consequences of genomic diversity in Mycobacterium tuberculosis. Semin. Immunol. 2014;26(6):431-444. Available at http://dx.doi.org/10.1016/j.smim.2014.09.012.
2. Doughty EL, Sergeant MJ, Adetifa I, Antonio M, Pallen MJ. Culture-independent detection and characterisation of Mycobacterium tuberculosis and M. africanum in sputum samples using shotgun metagenomics on a benchtop sequencer. PeerJ. 2014;2:1-18.
3. Chatterjee A, Nilgiriwala K, Saranath D, Rodrigues C, Mistry N. Whole genome sequencing of clinical strains of Mycobacterium tuberculosis from Mumbai, India : A potential tool for determining drug-resistance and strain lineage. Tuberculosis. 2017;107:63-72. Available at https://doi.org/10.1016/j.tube.2017.08.002.
4. Costa P, Botelho A, Couto I, Viveiros M, Inácio J. Standing of nucleic acid testing strategies in veterinary diagnosis laboratories to uncover Mycobacterium tuberculosis complex members. Front. Mol. Biosci. 2014;1 (October):1-11.
5. Gupta S, Kakkar V. Biosensors and Bioelectronics Recent technological advancements in tuberculosis diagnostics - A review. Biosens. Bioelectron. 2018;115(May):14-29. Available at https://doi.org/10.1016/j.bios.2018.05.017.
6. Wlodarska M, Johnston JC, Gardy JL. A Microbiological Revolution Meets an Ancient Disease : Improving the Management of Tuberculosis with Genomics. 2015;28(2):523-539.
7. Jagielski T, Minias A, Ingen J Van, Rastogi N, Brzostek A. Methodological and Clinical Aspects of the Molecular Epidemiology of Mycobacterium tuberculosis and Other Mycobacteria. Clin. Microbiol. Rev. 2016;29(2):239-290.
8. N′Dira Sanoussi C, Affolabi D, Rigouts L, Anagonou S, Jong B de. Genotypic characterization directly applied to sputum improves the detection of Mycobacterium africanum West African 1, under-represented in positive cultures. PL0S Negl. Trop. Dis. 2017:1-13.
9. Rue-albrecht K, Magee DA, Killick KE, et al. Comparative functional genomics and the bovine macrophage response to strains of the Mycobacterium genus. Front. Immunol. 2014;5(November):1-14.
10. IngenJ Van, Rahim Z, Mulder A, et al. Characterization of Mycobacterium orygis as M tuberculosis Complex Subspecies. Emerg. Infect. Dis. 2012;18(4):653--655.
11. Dippenaar A, David S, Parsons C, et al. Whole genome sequence analysis of Mycobacterium suricattae. Tuberculosis. 2015;95(6):682-688. Available at http://dx.doi.org/10.1016/j.tube.2015.10.001.
12. Alexander KA, Laver PN, Williams MC, et al. Pathology of the Emerging Mycobacterium tuberculosis Complex Pathogen, Mycobacterium mungi, in the Banded Mongoose ( Mungos mungo ). 2018;55(2):303-309.
13. Guthrie JL, Gardy JL. A brief primer on genomic epidemiology: lessons learned from Mycobacterium tuberculosis. Ann. N. Y. Acad. Sci. 2016:59-78.
14. Mcnerney R, Clark TG, Campino S, et al. International Journal of Infectious Diseases Removing the bottleneck in whole genome sequencing of Mycobacterium tuberculosis for rapid drug resistance analysis: a call to action. Int. J. Infect. Dis. 2017;56:130-135. Available at http://dx.doi.org/10.1016/j.ijid.2016.11.422.
15. Pankhurst LJ, Elias O, Votintseva AA, et al. Rapid, comprehensive, and aff ordable mycobacterial diagnosis with whole-genome sequencing: a prospective study. Lancet Respir. 4(1):49-58. Available at http://dx.doi.org/10.1016/S2213-2600(15)00466-X.
16. Brown AC, Bryant JM, Einer-jensen K, et al. Rapid Whole-Genome Sequencing of Mycobacterium tuberculosis Isolates Directly from Clinical Samples. J. Clin. Microbiol. 2015;53(7):2230-2237.
17. Kulchavenya E. Extrapulmonary tuberculosis: are statistical reports accurate? Ther. Adv. Infect. Dis. 2014;2(2):61-70.
18. Fisher M, Dolby T, Surtie S, et al. Improved method for collection of sputum for tuberculosis testing to ensure adequate sample volumes for molecular diagnostic testing. J. Microbiol. Methods. 2017;135:35-40. Available at: http://dx.doi.org/10.1016/j.mimet.2017.01.011.
19. World Health Organization. Global Tuberculosis Report. 2019.
20. Quan TP, Bawa Z, Foster D, et al. Evaluation of Whole-Genome Sequencing for Mycobacterial Species Identification and Drug Susceptibility Testing in a Clinical Setting: a Large-Scale Prospective Assessment of Performance against Line Probe Assays and Phenotyping. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2018;56(2):1-14.
21. Zumla A, Al-Tawfiq JA, Enne VI, et al. Rapid point of care diagnostic tests for viral and bacterial respiratory tract infections-needs, advances, and future prospects. Lancet Infect. Dis. 2014;14(11):1123-1135.
22. Walker TM, Kohl TA, Omar S V, et al. Whole-genome sequencing for prediction of Mycobacterium tuberculosis drug susceptibility and resistance: a retrospective cohort study. Lancet Infect. Dis. 2015;15:1193-1202.
23. Gardy JL. Towards genomic prediction of drug resistance in tuberculosis. Lancet Infect. Dis. 2015;15(10):1124-1125. Available at http://dx.doi.org/10.1016/S1473-3099(15)00088-2.
24. Bradley P, Gordon NC, Walker TM, et al. Rapid antibiotic-resistance predictions from genome sequence data for Staphylococcus aureus and Mycobacterium tuberculosis. Nat. Commun. 2015;6:1-14. Available at http://dx.doi.org/10.1038/ncomms10063.
25. Papaventsis D, Casali N, Kontsevaya I, et al. Whole genome sequencing of Mycobacterium tuberculosis for detection of drug resistance : a systematic review. Clin. Microbiol. Infect. 2017;23(2):61-68. Available at: http://dx.doi.org/10.1016/j.cmi.2016.09.008.
26. Nimmo C, Doyle R, Burgess C, et al. International Journal of Infectious Diseases Rapid identi fi cation of a Mycobacterium tuberculosis full genetic drug resistance pro fi le through whole genome sequencing directly from sputum. Int. J. Infect Dis. 2017;62:44-46. Available at http://dx.doi.org/10.1016/j.ijid.2017.07.007.
27. Linger Y, Knickerbocker C, Sipes D, et al. Genotyping Multidrug-Resistant Mycobacterium tuberculosis from Primary Sputum and Decontaminated Sediment with an Integrated Microfluidic Amplification Microarray Test. J. Clin. Microbiol. 2018;56(3):1-11.
28. Miotto P, Tessema B, Tagliani E, et al. A standardised method for interpreting the association between mutations and phenotypic drug resistance in Mycobacterium tuberculosis. Eur. Respir. J. 2017;50. Available at: http://dx.doi.org/10.1183/13993003.01354-2017.
29. World Health Organization. The use of next-generation sequencing technologies for the detection of mutations associated with drug resistance in Mycobacterium tuberculosis complex: technical guide. 2018.
30. Votintseva AA, Bradley P, Pankhurst LJ, et al. Same-Day Diagnostic and Surveillance Data for Tuberculosis via Whole-Genome Sequencing of Direct Respiratory Samples. J. Clin. Microbiol. 2017;55(5):1285-1298.
31. Haas CT, Roe JK, Pollara G, Mehta M, Noursadeghi M. Diagnostic ‘omics’ for active tuberculosis. BMC Med. 2016. Available at: http://dx.doi.org/10.1186/s12916-016-0583-9.
32. Lee RS, Pai M. Real-Time Sequencing of Mycobacterium tuberculosis: Are We There Yet? J. Clin. Microbiol. 2017;55(5):1249-1254.
33. Allahyartorkaman M, Mirsaeidi M, Hamzehloo G, et al. Low diagnostic accuracy of Xpert MTB /RIF assay for extrapulmonary tuberculosis : A multicenter surveillance. Sci. Rep. 2019;9:1-6. Available at: http://dx.doi.org/10.1038/s41598-019-55112-y.
34. Jouet A, Gaudin C, Badalato N, et al. free prediction of susceptibility or resistance to 13 anti-tuberculous drugs. Eur. Respir. J. 2020;(June 2020). Available at http://dx.doi.org/10.1183/13993003.02338-2020.
35. Feuerriegel S, Kohl TA, Utpatel C, et al. Early View Rapid genomic first- and second-line drug resistance prediction from clinical Mycobacterium tuberculosis specimens using Deeplex ® -MycTB. Eur. Respir. J. 2020.
36. World Health Organization. The Use of Next-Generation Sequencing Technologies for the Detection of Mutations Associated with Drug Resistance in Mycobacterium tuberculosis Complex: Technical Guide. 2018.
37. Meier A, Kirschner P, Bange FC, Vogel U, Bottger EC. Genetic alterations in streptomycin-resistant Mycobacterium tuberculosis: Mapping of mutations conferring resistance. Antimicrob. Agents Chemother. 1994;38(2):228-233.
38. Karimi, S., Mirhendi, H., Zaniani F., Manesh, S., Salehi, M., Esfahani B. Rapid detection of streptomycin-resistant Mycobacterium tuberculosis by rpsL-restriction fragment length polymorphism. Adv. Biomed. Res. 2017;6(126).
39. Villellas C, Aristimuño L, Vitoria MA, et al. Analysis of mutations in streptomycin-resistant strains reveals a simple and reliable genetic marker for identification of the Mycobacterium tuberculosis Beijing genotype. J. Clin. Microbiol 2013;51(7):2124-2130.
40. Morlock GP, Metchock B, Sikes D, Crawford JT, Cooksey RC. ethA, inhA and katG Loci of ethionamide- resistant Clinical MTB isolates. Antimicrob. Agents Chemother. 2003;47(12):3799-3805.
41. Zhao L, Sun Q, Liu H, et al. Analysis of embCAB Mutations Associated with Ethambutol Resistance in Multidrug-Resistant Mycobacterium tuberculosis Isolates from China. Antimicrob. Agents Chemother. 2015;59(4):2045-2050.
42. Maus CE, Plikaytis BB, Shinnick TM. Mutation of tlyA confers capreomycin resistance in Mycobacterium tuberculosis. Antimicrob. Agents Chemother. 2005;49(2):571-7. Available at http://www.ncbi.nlm.nih.gov/pubmed/15673735%0Ahttp://www.pubmedcentral.nih.gov/articlere nder.fcgi?artid=PMC547314.
43. (NCBI) NC for BI. Mycobacterium tuberculosis. Available at https://www.ncbi.nlm.nih.gov/genome/?term=h37rv [Accessed Jul. 17, 2020].
Design changes occurred only in Multiplex Group 1. Groups 2 and 3 remained unchanged during this period.
Design changes occurred only in Multiplex Group 2. Groups 1 and 3 remained unchanged during this period.
Design changes occurred only in Multiplex Group 3. Groups 1 and 2 remained unchanged ed during this period.
| Number | Date | Country | Kind |
|---|---|---|---|
| 2013928.3 | Sep 2020 | GB | national |
The present application is a filing under 35 U.S.C. 371 as the National Stage of International Application No. PCT/GB2021/052121, filed Aug. 16, 2021, entitled “METHOD AND COMPOSITIONS FOR DRUG RESISTANCE SCREENING,” which claims priority to United Kingdom Application No. 2013928.3filed with the Intellectual Property Office of the United Kingdom on Sep. 4, 2020, both of which are incorporated herein by reference in their entirety for all purposes.
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/GB2021/052121 | 8/16/2021 | WO |
